[Users] Interest in Null / Timelike geodesic integrator

[Yosef Zlochower]

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

yosef at astro.rit.edu

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.