aboutsummaryrefslogtreecommitdiff
path: root/src/mtx_spherical_harmonics/chebyshev12.h
diff options
context:
space:
mode:
Diffstat (limited to 'src/mtx_spherical_harmonics/chebyshev12.h')
-rw-r--r--src/mtx_spherical_harmonics/chebyshev12.h36
1 files changed, 36 insertions, 0 deletions
diff --git a/src/mtx_spherical_harmonics/chebyshev12.h b/src/mtx_spherical_harmonics/chebyshev12.h
new file mode 100644
index 0000000..0da747a
--- /dev/null
+++ b/src/mtx_spherical_harmonics/chebyshev12.h
@@ -0,0 +1,36 @@
+/*
+ Evaluates all circular harmonics
+ at the angles phi up to the order nmax.
+ using the recurrence for the Chebyshev
+ polynomials of the first and second kind
+ T has the dimensions length(phi) x 2nmax+1
+
+ Implementation by Franz Zotter, Institute of Electronic Music and Acoustics
+ (IEM), University of Music and Dramatic Arts (KUG), Graz, Austria
+ http://iem.at/Members/zotter, 2008.
+
+ This code is published under the Gnu General Public License, see
+ "LICENSE.txt"
+*/
+
+#ifndef _chebyshev12_h__
+#define _chebyshev12_h__
+#include <math.h>
+#include <stdlib.h>
+
+typedef struct _Cheby12WorkSpace_
+{
+ size_t nmax;
+ size_t l;
+ double *t;
+} Cheby12WorkSpace;
+
+Cheby12WorkSpace *chebyshev12_alloc(const size_t nmax, const size_t l);
+
+void chebyshev12_free(Cheby12WorkSpace *wc);
+
+void chebyshev12(double *phi, Cheby12WorkSpace *wc);
+
+#endif
+
+