%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  Computation of QNMs for slowly rotating Kerr-de Sitter black holes
%  by means of a quantization condition
%  (c) Semyon Dyatlov, UC Berkeley, 2010 <dyatlov@berkeley.edu>
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% info.txt - description of the project

The front end consists of the following functions:
getQNM - compute the QNM using WKB approximation
outputQNM - output QNMs into a file
loadQNM - load QNMs from a file

Some information on the inside workings:

Jet - a matrix containing the values of a function and several of its
derivatives on the grid. If A is a jet, then A(i, j) is the value of the
(j - 1)'st derivative of the function at the grid point i, divided by (j - 1)!

Polynomials are always assumed to be row vectors, with A(i) being the
    coefficient next to x^(i - 1).

KDSData - a structure containing information about the current constants of
the problem and the QNMs. The fields are:
a, Lambda, M, alpha - the corresponding constants of the black hole
deltarc - the coefficient vector of the degree 4 polynomial Delta_r
omegaR, omegaI, k - the spectral parameters currently being used.
    Should only exist inside the computation functions
    omega = omegaR + i * omegaI
rgrid - the current grid of values of r
    Should only exist inside the radial computation functions