Extending Source Terms to SR

[hwhitehead]

Does anyone have a good understanding on how Newtonian source terms might be extended to special relativity? For instance, if I wished to study a flow moving at relativistic speeds about a Newtonian point source of gravity, would the implementation in pointmass.cpp be sufficient, or is extra caution required? To clarify, this would be an externally imposed gravity term, not from gas self-gravity.

My concern is that quantities such as the energy flux across the boundary e.g. in FLUX[X1DIR] might require additional nuance, but it is possible that this is already accounted for by the routine that populates the FLUX arrays.

If these routines are not SR compliant, it would be worth protecting against their inclusion if the RELATIVISTIC_DYNAMICS flag is not switched on.