Split the basis sets definitions into a separate file.

1 files changed, 27 insertions, 0 deletions

diff --git a/src/maximal_slicing_axi.h b/src/maximal_slicing_axi.h
new file mode 100644
index 0000000..b15420c
--- /dev/null
+++ b/src/maximal_slicing_axi.h
@@ -0,0 +1,27 @@
+
+#define SQR(x) ((x) * (x))
+#define SGN(x) ((x) >= 0.0 ? 1.0 : -1.0)
+
+#define SCALE_FACTOR scale_factor
+
+/* a set of basis functions */
+typedef struct BasisSet {
+    /* evaluate the idx-th basis function at the specified point*/
+    double (*eval)      (double coord, int idx);
+    /* evaluate the first derivative of the idx-th basis function at the specified point*/
+    double (*eval_diff1)(double coord, int idx);
+    /* evaluate the second derivative of the idx-th basis function at the specified point*/
+    double (*eval_diff2)(double coord, int idx);
+    /**
+     * Get the idx-th collocation point for the specified order.
+     * idx runs from 0 to order - 1 (inclusive)
+     */
+    double (*colloc_point)(int order, int idx);
+} BasisSet;
+
+extern const BasisSet msa_cheb_basis;
+extern const BasisSet msa_full_basis;
+extern const BasisSet msa_tb_even_basis;
+extern const BasisSet msa_sb_even_basis;
+
+extern double scale_factor;