PTAfast calculates the overlap reduction function in Pulsar Timing Array produced by the stochastic gravitational wave background for arbitrary polarizations, propagation velocities, and pulsar distances.
https://ui.adsabs.harvard.edu/abs/2022arXiv220812538B and https://ui.adsabs.harvard.edu/abs/2022JCAP...11..046B