/* $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"
#include "../common/structs.h"
#include "../common/fetch.h"
#include "../common/evaluate.h"
#include "../common/store.h"

/* function prototype */
#define FUNCTION_NAME	LocalInterp_U_LagrTP_1d_cube60
#include "../template.h"

#define N_DIMS		1
#define MOLECULE_MIN_M	-3
#define MOLECULE_MAX_M	3
#define MOLECULE_SIZE	7

/* which derivative ops do we support? */
#define HAVE_OP_I
#define HAVE_OP_DX
#define HAVE_OP_DXX

#define XYZ				x
#define FP_XYZ				fp x
#define STRIDE_IJK			stride_i
#define JACOBIAN_MIJK_STRIDE		Jacobian_mi_stride

#define DATA_STRUCT			data_struct_1d_cube_size7
#define COEFFS_STRUCT			coeffs_struct_1d_cube_size7

#define FETCH_DATA_REAL			LocalInterp_fetch_1d_cube7_r
#define FETCH_DATA_REAL4		LocalInterp_fetch_1d_cube7_r4
#define FETCH_DATA_REAL8		LocalInterp_fetch_1d_cube7_r8
#define FETCH_DATA_REAL16		LocalInterp_fetch_1d_cube7_r16
#define FETCH_DATA_COMPLEX		LocalInterp_fetch_1d_cube7_c
#define FETCH_DATA_COMPLEX8		LocalInterp_fetch_1d_cube7_c8
#define FETCH_DATA_COMPLEX16		LocalInterp_fetch_1d_cube7_c16
#define FETCH_DATA_COMPLEX32		LocalInterp_fetch_1d_cube7_c32

#define EVALUATE_MOLECULE		LocalInterp_eval_1d_cube7

#define STORE_COEFFS			LocalInterp_store_1d_cube7

/* note pathnames are all relative to "../template.c" */
#define COEFFS_I_COMPUTE_FILE_NAME	"Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-I.compute.c"
#define COEFFS_DX_COMPUTE_FILE_NAME	"Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-dx.compute.c"
#define COEFFS_DXX_COMPUTE_FILE_NAME	"Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-dxx.compute.c"

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