/* $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_3
#include "../template.h"

#define N_DIMS		3
#define MOLECULE_MIN_M	-2
#define MOLECULE_MAX_M	3
#define MOLECULE_SIZE	6

/* 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.size6/data-var.dcl.c"
#define DATA_VAR_ASSIGN_FILE_NAME	"../common/3d.cube.size6/data-var.assign.c"

#define COEFF_I_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-I.store.c"
#define COEFF_DX_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-dx.store.c"
#define COEFF_DY_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-dy.store.c"
#define COEFF_DZ_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-dz.store.c"
#define COEFF_DXX_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-dxx.store.c"
#define COEFF_DXY_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-dxy.store.c"
#define COEFF_DXZ_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-dxz.store.c"
#define COEFF_DYY_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-dyy.store.c"
#define COEFF_DYZ_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-dyz.store.c"
#define COEFF_DZZ_VAR_STORE_FILE_NAME	"../common/3d.cube.size6/coeff-dzz.store.c"

#define COEFF_I_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-I.dcl.c"
#define COEFF_DX_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-dx.dcl.c"
#define COEFF_DY_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-dy.dcl.c"
#define COEFF_DZ_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-dz.dcl.c"
#define COEFF_DXX_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-dxx.dcl.c"
#define COEFF_DXY_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-dxy.dcl.c"
#define COEFF_DXZ_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-dxz.dcl.c"
#define COEFF_DYY_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-dyy.dcl.c"
#define COEFF_DYZ_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-dyz.dcl.c"
#define COEFF_DZZ_DCL_FILE_NAME		"../common/3d.cube.size6/coeff-dzz.dcl.c"

#define INTERP_I_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-I.compute.c"
#define INTERP_DX_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-dx.compute.c"
#define INTERP_DY_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-dy.compute.c"
#define INTERP_DZ_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-dz.compute.c"
#define INTERP_DXX_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-dxx.compute.c"
#define INTERP_DXY_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-dxy.compute.c"
#define INTERP_DXZ_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-dxz.compute.c"
#define INTERP_DYY_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-dyy.compute.c"
#define INTERP_DYZ_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-dyz.compute.c"
#define INTERP_DZZ_COMPUTE_FILE_NAME	"../common/3d.cube.size6/interp-dzz.compute.c"

#define COEFF_I_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-I.compute.c"
#define COEFF_DX_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-dx.compute.c"
#define COEFF_DY_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-dy.compute.c"
#define COEFF_DZ_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-dz.compute.c"
#define COEFF_DXX_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-dxx.compute.c"
#define COEFF_DXY_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-dxy.compute.c"
#define COEFF_DXZ_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-dxz.compute.c"
#define COEFF_DYY_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-dyy.compute.c"
#define COEFF_DYZ_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-dyz.compute.c"
#define COEFF_DZZ_COMPUTE_FILE_NAME	"3d.coeffs/3d.cube.order3/coeff-dzz.compute.c"

/* actual code */
#include "../template.c"