reorganize common to Lagrange and Hermite interpolators

(i.e. code depending only on molecule size/shape)
--> now there are separate *functions* rather than just code fragments
--> ../template.c should hopefully compile in finite time.....


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@109 df1f8a13-aa1d-4dd4-9681-27ded5b42416

1 files changed, 76 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/common/store.h b/src/GeneralizedPolynomial-Uniform/common/store.h
new file mode 100644
index 0000000..e0911af
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/store.h
@@ -0,0 +1,76 @@
+/* store.h -- store molecule coefficients in Jacobian */
+/* $Header$ */
+
+/*
+ * prerequisite headers:
+ *	"cctk.h"
+ *	"../InterpLocalUniform.h"
+ *	"structs.h"
+ */
+
+void LocalInterp_store_1d_cube2
+      (fp factor, const struct coeffs_struct_1d_cube_size2 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride);
+void LocalInterp_store_1d_cube3
+      (fp factor, const struct coeffs_struct_1d_cube_size3 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride);
+void LocalInterp_store_1d_cube4
+      (fp factor, const struct coeffs_struct_1d_cube_size4 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride);
+void LocalInterp_store_1d_cube5
+      (fp factor, const struct coeffs_struct_1d_cube_size5 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride);
+void LocalInterp_store_1d_cube6
+      (fp factor, const struct coeffs_struct_1d_cube_size6 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride);
+void LocalInterp_store_1d_cube7
+      (fp factor, const struct coeffs_struct_1d_cube_size7 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride);
+
+void LocalInterp_store_2d_cube2
+      (fp factor, const struct coeffs_struct_2d_cube_size2 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride);
+void LocalInterp_store_2d_cube3
+      (fp factor, const struct coeffs_struct_2d_cube_size3 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride);
+void LocalInterp_store_2d_cube4
+      (fp factor, const struct coeffs_struct_2d_cube_size4 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride);
+void LocalInterp_store_2d_cube5
+      (fp factor, const struct coeffs_struct_2d_cube_size5 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride);
+void LocalInterp_store_2d_cube6
+      (fp factor, const struct coeffs_struct_2d_cube_size6 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride);
+
+void LocalInterp_store_3d_cube2
+      (fp factor, const struct coeffs_struct_3d_cube_size2 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride);
+void LocalInterp_store_3d_cube3
+      (fp factor, const struct coeffs_struct_3d_cube_size3 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride);
+void LocalInterp_store_3d_cube4
+      (fp factor, const struct coeffs_struct_3d_cube_size4 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride);
+void LocalInterp_store_3d_cube5
+      (fp factor, const struct coeffs_struct_3d_cube_size5 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride);
+void LocalInterp_store_3d_cube6
+      (fp factor, const struct coeffs_struct_3d_cube_size6 *coeffs,
+       fp Jacobian_ptr[],
+       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride);