search

znasipak / pybhpt

A python package for solving problems in black hole perturbation theory
GNU General Public License v3.0
7 stars 1 forks source link

Use `spherical` package for Wigner 3j and spherical harmonics? #1

Open duetosymmetry opened 1 year ago

duetosymmetry commented 1 year ago

Mike Boyle (@moble) has taken great care in the package spherical to make sure that Wigner 3j symbols can be evaluated without using factorials, to avoid overflows and loss of precision when there are large cancellations.

There is also a bunch of overlap with the my qnm package in computing the decomposition of spheroidal harmonics into spherical harmonics, but this is not a big deal.

znasipak commented 1 year ago

These are great suggestions. This code was designed originally as a C/C++ library that was focused on just running generic-orbit self-force calculations. This package is an attempt to make that code more accessible through Python. But now that it is actually in Python, I should incorporate these awesome tools that you and Mike Boyle built. I expect that they are more efficient than what I have currently implemented.