diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/3d.cube.order1.smooth0.c')
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange/3d.cube.order1.smooth0.c | 57 |
1 files changed, 33 insertions, 24 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.cube.order1.smooth0.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.cube.order1.smooth0.c index f97eb83..7c06d13 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.cube.order1.smooth0.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.cube.order1.smooth0.c @@ -9,9 +9,13 @@ #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_ULagrange_3dcube_10 +#define FUNCTION_NAME LocalInterp_ULagrange_3d_cube10 #include "../template.h" #define N_DIMS 3 @@ -24,29 +28,34 @@ #define HAVE_OP_DX #define HAVE_OP_DY #define HAVE_OP_DZ - -#define DATA_VAR_DCL_FILE_NAME "../common/3d.cube.size2/data-var.dcl.c" -#define DATA_VAR_ASSIGN_FILE_NAME "../common/3d.cube.size2/data-var.assign.c" - -#define COEFF_I_VAR_STORE_FILE_NAME "../common/3d.cube.size2/coeff-I.store.c" -#define COEFF_DX_VAR_STORE_FILE_NAME "../common/3d.cube.size2/coeff-dx.store.c" -#define COEFF_DY_VAR_STORE_FILE_NAME "../common/3d.cube.size2/coeff-dy.store.c" -#define COEFF_DZ_VAR_STORE_FILE_NAME "../common/3d.cube.size2/coeff-dz.store.c" - -#define COEFF_I_DCL_FILE_NAME "../common/3d.cube.size2/coeff-I.dcl.c" -#define COEFF_DX_DCL_FILE_NAME "../common/3d.cube.size2/coeff-dx.dcl.c" -#define COEFF_DY_DCL_FILE_NAME "../common/3d.cube.size2/coeff-dy.dcl.c" -#define COEFF_DZ_DCL_FILE_NAME "../common/3d.cube.size2/coeff-dz.dcl.c" - -#define INTERP_I_COMPUTE_FILE_NAME "../common/3d.cube.size2/interp-I.compute.c" -#define INTERP_DX_COMPUTE_FILE_NAME "../common/3d.cube.size2/interp-dx.compute.c" -#define INTERP_DY_COMPUTE_FILE_NAME "../common/3d.cube.size2/interp-dy.compute.c" -#define INTERP_DZ_COMPUTE_FILE_NAME "../common/3d.cube.size2/interp-dz.compute.c" - -#define COEFF_I_COMPUTE_FILE_NAME "3d.coeffs/3d.cube.order1.smooth0/coeff-I.compute.c" -#define COEFF_DX_COMPUTE_FILE_NAME "3d.coeffs/3d.cube.order1.smooth0/coeff-dx.compute.c" -#define COEFF_DY_COMPUTE_FILE_NAME "3d.coeffs/3d.cube.order1.smooth0/coeff-dy.compute.c" -#define COEFF_DZ_COMPUTE_FILE_NAME "3d.coeffs/3d.cube.order1.smooth0/coeff-dz.compute.c" +/* n.b. no 2nd derivatives for linear interpolation! */ + +#define XYZ x, y, z +#define FP_XYZ fp x, fp y, fp z +#define STRIDE_IJK stride_i, stride_j, stride_k +#define JACOBIAN_MIJK_STRIDE Jacobian_mi_stride, Jacobian_mj_stride, Jacobian_mk_stride + +#define DATA_STRUCT data_struct_3d_cube_size2 +#define COEFFS_STRUCT coeffs_struct_3d_cube_size2 + +#define FETCH_DATA_REAL LocalInterp_fetch_3d_cube2_r +#define FETCH_DATA_REAL4 LocalInterp_fetch_3d_cube2_r4 +#define FETCH_DATA_REAL8 LocalInterp_fetch_3d_cube2_r8 +#define FETCH_DATA_REAL16 LocalInterp_fetch_3d_cube2_r16 +#define FETCH_DATA_COMPLEX LocalInterp_fetch_3d_cube2_c +#define FETCH_DATA_COMPLEX8 LocalInterp_fetch_3d_cube2_c8 +#define FETCH_DATA_COMPLEX16 LocalInterp_fetch_3d_cube2_c16 +#define FETCH_DATA_COMPLEX32 LocalInterp_fetch_3d_cube2_c32 + +#define EVALUATE_MOLECULE LocalInterp_eval_3d_cube2 + +#define STORE_COEFFS LocalInterp_store_3d_cube2 + +/* note pathnames are all relative to "../template.c" */ +#define COEFFS_I_COMPUTE_FILE_NAME "Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeffs-I.compute.c" +#define COEFFS_DX_COMPUTE_FILE_NAME "Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeffs-dx.compute.c" +#define COEFFS_DY_COMPUTE_FILE_NAME "Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeffs-dy.compute.c" +#define COEFFS_DZ_COMPUTE_FILE_NAME "Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeffs-dz.compute.c" /* actual code */ #include "../template.c" |