#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_ILU_2d_cube_o1_s0
#include "template.h"

#define N_DIMS		2
#define MOLECULE_MIN_M	0
#define MOLECULE_MAX_M	1
#define MOLECULE_SIZE	2

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

#define DATA_VAR_DCL_FILE_NAME		"2d.coeffs/2d.cube.size2/data-var.dcl.c"
#define DATA_VAR_ASSIGN_FILE_NAME	"2d.coeffs/2d.cube.size2/data-var.assign.c"

#define COEFF_I_DCL_FILE_NAME		"2d.coeffs/2d.cube.size2/coeff-I.dcl.c"
#define COEFF_DX_DCL_FILE_NAME		"2d.coeffs/2d.cube.size2/coeff-dx.dcl.c"
#define COEFF_DY_DCL_FILE_NAME		"2d.coeffs/2d.cube.size2/coeff-dy.dcl.c"

#define INTERP_I_COMPUTE_FILE_NAME	"2d.coeffs/2d.cube.size2/interp-I.compute.c"
#define INTERP_DX_COMPUTE_FILE_NAME	"2d.coeffs/2d.cube.size2/interp-dx.compute.c"
#define INTERP_DY_COMPUTE_FILE_NAME	"2d.coeffs/2d.cube.size2/interp-dy.compute.c"

#define COEFF_I_COMPUTE_FILE_NAME	"2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c"
#define COEFF_DX_COMPUTE_FILE_NAME	"2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c"
#define COEFF_DY_COMPUTE_FILE_NAME	"2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c"

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