Skip to content

Commit

Permalink
More documentation, remove testing notebooks
Browse files Browse the repository at this point in the history
  • Loading branch information
duetosymmetry committed Mar 10, 2019
1 parent da6d5eb commit 99e9563
Show file tree
Hide file tree
Showing 8 changed files with 54 additions and 3,220 deletions.
31 changes: 21 additions & 10 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -5,9 +5,20 @@


# Welcome to qnm
Python implementation of Cook-Zalutskiy spectral approach to computing Kerr QNM frequencies.
Python implementation of Cook-Zalutskiy spectral approach to computing
Kerr quasinormal frequencies (QNMs).

TODO basic info
With this python package, you can compute the QNMs labeled by
different (s,l,m,n), at a desired dimensionless spin parameter 0≤a<1.
The angular sector is treated as a spectral decomposition of
spin-weighted *spheroidal* harmonics into spin-weighted spherical
harmonics. Therefore the spherical-spheroidal decomposition
coefficients come for free when solving for ω and A.

We have precomputed a large number of low-lying modes (s=-2 and s=-1,
all l<8, all n<7). These can be automatically installed with a single
function call, and interpolated for good initial guesses for
root-finding at some value of a.

## Installation

Expand Down Expand Up @@ -95,21 +106,21 @@ for ind in mode_list:

plt.figure(figsize=(16,8))

plt.subplot(1, 2, 1)
plt.subplot(1, 2, 1)
for mode, seq in modes.iteritems():
plt.plot(np.real(seq.omega),np.imag(seq.omega))


modestr = "{},{},{},n".format(s,l,m)
plt.xlabel(r'$\textrm{Re}[\omega_{' + modestr + r'}]$', fontsize=16)
plt.ylabel(r'$\textrm{Im}[\omega_{' + modestr + r'}]$', fontsize=16)
plt.gca().tick_params(labelsize=16)
plt.gca().invert_yaxis()

plt.subplot(1, 2, 2)
plt.subplot(1, 2, 2)
for mode, seq in modes.iteritems():
plt.plot(np.real(seq.A),np.imag(seq.A))

plt.xlabel(r'$\textrm{Re}[A_{' + modestr + r'}]$', fontsize=16)
plt.ylabel(r'$\textrm{Im}[A_{' + modestr + r'}]$', fontsize=16)
plt.gca().tick_params(labelsize=16)
Expand All @@ -129,21 +140,21 @@ for ind in mode_list:

plt.figure(figsize=(16,8))

plt.subplot(1, 2, 1)
plt.subplot(1, 2, 1)
for mode, seq in modes.iteritems():
plt.plot(np.real(seq.omega),np.imag(seq.omega))


modestr = "{},{},m,0".format(s,l)
plt.xlabel(r'$\textrm{Re}[\omega_{' + modestr + r'}]$', fontsize=16)
plt.ylabel(r'$\textrm{Im}[\omega_{' + modestr + r'}]$', fontsize=16)
plt.gca().tick_params(labelsize=16)
plt.gca().invert_yaxis()

plt.subplot(1, 2, 2)
plt.subplot(1, 2, 2)
for mode, seq in modes.iteritems():
plt.plot(np.real(seq.A),np.imag(seq.A))

plt.xlabel(r'$\textrm{Re}[A_{' + modestr + r'}]$', fontsize=16)
plt.ylabel(r'$\textrm{Im}[A_{' + modestr + r'}]$', fontsize=16)
plt.gca().tick_params(labelsize=16)
Expand Down
2 changes: 1 addition & 1 deletion docs/index.rst
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
qnm's documentation: index
==========================

TODO point people in the right direction here
For a quick start, take a look at the `usage section <README.html#usage>`_ of the Welcome node.

.. toctree::
:maxdepth: 2
Expand Down
1 change: 1 addition & 0 deletions docs/notebooks/example_22n.png
1 change: 1 addition & 0 deletions docs/notebooks/example_2m0.png
259 changes: 0 additions & 259 deletions notebooks/Continued fraction testing.ipynb

This file was deleted.

Loading

0 comments on commit 99e9563

Please sign in to comment.