The current sink code is incorrect:
The loops are not treated in a symmetric way. By convention
iact(a,b) has to yield the same result as
I have fixed the code in
src/runner_doiact_functions_hydro.h to make sure the symmetric calls are made.
But there is still a mistake in the GEAR implementation as the symmetric function is empty.