sxs-collaboration / spectre

SpECTRE is a code for multi-scale, multi-physics problems in astrophysics and gravitational physics.
https://spectre-code.org
Other
159 stars 189 forks source link

FindHorizon.py: also compute horizon quantities #6024

Closed nilsvu closed 4 months ago

nilsvu commented 4 months ago

Proposed changes

This allows to get the horizon area, mass, spin, etc from BBH initial data just by calling a Python function.

Upgrade instructions

The FindHorizons3D executable was removed. You can find horizons in 3D volume data in Python like this:

from spectre.Pipelines.Bbh.FindHorizon import find_horizon
from spectre.SphericalHarmonics import Strahlkorper

horizon, quantities = find_horizon(
    h5_files, subfile_name, obs_id, obs_time,
    initial_guess=Strahlkorper(l_max, m_max, radius, center)
)

mass = quantities["ChristodoulouMass"]
spin = quantities["DimensionlessSpinMagnitude"]

Code review checklist

Further comments