CUDA版で --mutual を明示的に指定すると結果がずれる

mutualについての扱いでGPU版の細部で不整合があるので結果がずれるこの扱いはどうするべきか？GPU版だとright hand side側への書き込みが遅いが、カーネルを書くのはユーザー側なのでユーザー側で対処すべき

$ MYTH_WORKER_NUM=1 ./cuda_tapas --numBodies 10000
Rel. L2 Error (pot)  : 3.6826631e-03
Rel. L2 Error (acc)  : 2.9516280e-03

$ MYTH_WORKER_NUM=1 ./cuda_tapas --numBodies 10000 --mutual 0
Rel. L2 Error (pot)  : 9.0998369e-01
Rel. L2 Error (acc)  : 2.1257406e-01

$ MYTH_WORKER_NUM=1 ./cuda_tapas --numBodies 10000 --mutual 1
Rel. L2 Error (pot)  : 3.6826631e-03
Rel. L2 Error (acc)  : 2.9516280e-03
# これはユーザーが自分で実行時エラーにすべき

keisukefukuda / tapas

CUDA版で --mutual を明示的に指定すると結果がずれる #13