QuasinormalModes.jl: A Julia package for computing discrete eigenvalues of second order ODEs

TL;DR

QuasinormalModes.jl is a versatile Julia package that efficiently computes discrete eigenvalues of second order ODEs, including those in General Relativity and quantum mechanics, with flexible precision and parallel processing options.

Contribution

The package offers a user-friendly, adaptable API for computing eigenvalues of second order ODEs, supporting arbitrary precision and parallel computation, tailored for research integration.

Findings

Supports a wide range of second order ODE eigenvalue problems.

Allows flexible precision and parallel processing options.

Designed for easy integration into existing research workflows.

Abstract

QuasinormalModes.jl is a Julia package for computing the quasinormal modes of any General Relativity model whose perturbation equation can be expressed as second order homogeneous ODE. Not only that, the package can be used to compute the discrete eigenvalues of any second order homogeneous ODE (such as the energy eigenstates of the time independent Schr\"odinger equation) provided that these eigenvalues actually exist. The package features a flexible and user friendly API where the user simply needs to provide the coefficients of the problem ODE after incorporating boundary and asymptotic conditions on it. The user can also choose to use machine or arbitrary precision arithmetic for the underlying floating point operations involved and whether or not to do computations sequentially or in parallel using threads. The API also tries not to force any particular workflow on the users so that they can incorporate and adapt the existing functionality on their research pipelines without unwanted intrusions.