Hi, I have been working on an integrator for timelike and null geodesics that evolve the geodesics using MoL. With an RK4 integration, and a fourth-order interpolation stencil, the code should be third-order accurate (due to the lower-order accuracy obtained when interpolating derivatives).
The code appears to work well with a typical moving boxes BBH simulation (I haven't actually tried null geodesics yet). It however, will not scale well to arbitrary numbers of geodesics as all processes need to know the locations of each geodesic. I'd be happy to share the code if there is some interest. -- Dr. Yosef Zlochower Center for Computational Relativity and Gravitation Associate Professor School of Mathematical Sciences Rochester Institute of Technology 85 Lomb Memorial Drive Rochester, NY 14623 Office:74-2067 Phone: +1 585-475-6103 [email protected] CONFIDENTIALITY NOTE: The information transmitted, including attachments, is intended only for the person(s) or entity to which it is addressed and may contain confidential and/or privileged material. Any review, retransmission, dissemination or other use of, or taking of any action in reliance upon this information by persons or entities other than the intended recipient is prohibited. If you received this in error, please contact the sender and destroy any copies of this information. _______________________________________________ Users mailing list [email protected] http://lists.einsteintoolkit.org/mailman/listinfo/users
