This PR adds calculations for k_mu source terms (radiation/geodesics.hpp) for arbitrary metric/geometry as it was previously specialized to Minkowski. This was causing unphysical positions when running the torus problem. It now runs beyond the problem point and a longer run in in the queue. The leptoneq problem seems reasonable.
Also removes a print statement from plot_snap2d that I left from before.
TODO:
[ ] There are redundant calculations (in Minowski all the sources are 0, in FMKS only 2 are non zero) and we are doing more integration than necessary. Performant way to specialize? Template on geometry?
[ ] Write a simple test for particles around a Kerr black hole to say for sure that the Push is doing what we want.
This PR adds calculations for
k_mu
source terms (radiation/geodesics.hpp
) for arbitrary metric/geometry as it was previously specialized to Minkowski. This was causing unphysical positions when running the torus problem. It now runs beyond the problem point and a longer run in in the queue. Theleptoneq
problem seems reasonable.Also removes a print statement from
plot_snap2d
that I left from before.TODO: