#1152: Add option to compute center of mass to hydro_analysis and lots of smaller changes ------------------------------------+--------------------------------------- Reporter: rhaas | Owner: Type: enhancement | Status: review Priority: minor | Milestone: Component: EinsteinToolkit thorn | Version: Resolution: | Keywords: Hydro_Analysis ------------------------------------+---------------------------------------
Comment (by rhaas):
myrho is declared by DECLARE_MY_VARIABLE(rho) in line 44. I am fine renaming the routine to "centroid" rather than "center of mass". Since this is only used to track the star and move the grids along using just rho is fine. For the actual center of mass one would either have to compute sqrt(det(g_ij)) * w_lorentz * rho or integrate dens. This would be more expensive. It is very easy to change the code so that it takes a parameter specifying what "rho" should be "HYDROBASE::rho" would be the current situation and "GRHYDRO::dens" would be the actual center of mass. Would that be acceptable?