add Maple code and C skeletons for Hermite polynomial interpolation

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

1 files changed, 82 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/3d.cube.order2.c b/src/GeneralizedPolynomial-Uniform/Hermite/3d.cube.order2.c
new file mode 100644
index 0000000..c22c56d
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/3d.cube.order2.c
@@ -0,0 +1,82 @@
+/* $Header$ */
+
+#include <math.h>
+#include <limits.h>
+#include <stdlib.h>
+#include <string.h>
+#include <stdio.h>
+
+#include "util_ErrorCodes.h"
+#include "cctk.h"
+#include "../InterpLocalUniform.h"
+
+/* function prototype */
+#define FUNCTION_NAME	LocalInterp_UHermite_3dcube_2
+#include "../template.h"
+
+#define N_DIMS		3
+#define MOLECULE_MIN_M	-1
+#define MOLECULE_MAX_M	2
+#define MOLECULE_SIZE	4
+
+/* which derivative ops do we support? */
+#define HAVE_OP_I
+#define HAVE_OP_DX
+#define HAVE_OP_DY
+#define HAVE_OP_DZ
+#define HAVE_OP_DXX
+#define HAVE_OP_DXY
+#define HAVE_OP_DXZ
+#define HAVE_OP_DYY
+#define HAVE_OP_DYZ
+#define HAVE_OP_DZZ
+
+#define DATA_VAR_DCL_FILE_NAME		"../common/3d.cube.size4/data-var.dcl.c"
+#define DATA_VAR_ASSIGN_FILE_NAME	"../common/3d.cube.size4/data-var.assign.c"
+
+#define COEFF_I_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-I.store.c"
+#define COEFF_DX_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-dx.store.c"
+#define COEFF_DY_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-dy.store.c"
+#define COEFF_DZ_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-dz.store.c"
+#define COEFF_DXX_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-dxx.store.c"
+#define COEFF_DXY_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-dxy.store.c"
+#define COEFF_DXZ_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-dxz.store.c"
+#define COEFF_DYY_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-dyy.store.c"
+#define COEFF_DYZ_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-dyz.store.c"
+#define COEFF_DZZ_VAR_STORE_FILE_NAME	"../common/3d.cube.size4/coeff-dzz.store.c"
+
+#define COEFF_I_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-I.dcl.c"
+#define COEFF_DX_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-dx.dcl.c"
+#define COEFF_DY_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-dy.dcl.c"
+#define COEFF_DZ_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-dz.dcl.c"
+#define COEFF_DXX_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-dxx.dcl.c"
+#define COEFF_DXY_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-dxy.dcl.c"
+#define COEFF_DXZ_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-dxz.dcl.c"
+#define COEFF_DYY_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-dyy.dcl.c"
+#define COEFF_DYZ_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-dyz.dcl.c"
+#define COEFF_DZZ_DCL_FILE_NAME		"../common/3d.cube.size4/coeff-dzz.dcl.c"
+
+#define INTERP_I_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-I.compute.c"
+#define INTERP_DX_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-dx.compute.c"
+#define INTERP_DY_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-dy.compute.c"
+#define INTERP_DZ_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-dz.compute.c"
+#define INTERP_DXX_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-dxx.compute.c"
+#define INTERP_DXY_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-dxy.compute.c"
+#define INTERP_DXZ_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-dxz.compute.c"
+#define INTERP_DYY_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-dyy.compute.c"
+#define INTERP_DYZ_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-dyz.compute.c"
+#define INTERP_DZZ_COMPUTE_FILE_NAME	"../common/3d.cube.size4/interp-dzz.compute.c"
+
+#define COEFF_I_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-I.compute.c"
+#define COEFF_DX_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-dx.compute.c"
+#define COEFF_DY_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-dy.compute.c"
+#define COEFF_DZ_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-dz.compute.c"
+#define COEFF_DXX_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-dxx.compute.c"
+#define COEFF_DXY_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-dxy.compute.c"
+#define COEFF_DXZ_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-dxz.compute.c"
+#define COEFF_DYY_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-dyy.compute.c"
+#define COEFF_DYZ_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-dyz.compute.c"
+#define COEFF_DZZ_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order2/coeff-dzz.compute.c"
+
+/* actual code */
+#include "../template.c"