diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.cube.order4.smooth0.c')
| -rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.cube.order4.smooth0.c | 71 |
1 files changed, 71 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.cube.order4.smooth0.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.cube.order4.smooth0.c new file mode 100644 index 0000000..4842b2f --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.cube.order4.smooth0.c @@ -0,0 +1,71 @@ +/* $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_3d_cube40 +#include "../template.h" + +#define N_DIMS 3 +#define MOLECULE_MIN_M -2 +#define MOLECULE_MAX_M 2 +#define MOLECULE_SIZE 5 + +/* 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 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_size5 +#define COEFFS_STRUCT coeffs_struct_3d_cube_size5 + +#define FETCH_DATA_REAL LocalInterp_fetch_3d_cube5_r +#define FETCH_DATA_REAL4 LocalInterp_fetch_3d_cube5_r4 +#define FETCH_DATA_REAL8 LocalInterp_fetch_3d_cube5_r8 +#define FETCH_DATA_REAL16 LocalInterp_fetch_3d_cube5_r16 +#define FETCH_DATA_COMPLEX LocalInterp_fetch_3d_cube5_c +#define FETCH_DATA_COMPLEX8 LocalInterp_fetch_3d_cube5_c8 +#define FETCH_DATA_COMPLEX16 LocalInterp_fetch_3d_cube5_c16 +#define FETCH_DATA_COMPLEX32 LocalInterp_fetch_3d_cube5_c32 + +#define EVALUATE_MOLECULE LocalInterp_eval_3d_cube5 + +#define STORE_COEFFS LocalInterp_store_3d_cube5 + +#define COEFFS_I_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-I.compute.c" +#define COEFFS_DX_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-dx.compute.c" +#define COEFFS_DY_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-dy.compute.c" +#define COEFFS_DZ_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-dz.compute.c" +#define COEFFS_DXX_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-dxx.compute.c" +#define COEFFS_DXY_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-dxy.compute.c" +#define COEFFS_DXZ_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-dxz.compute.c" +#define COEFFS_DYY_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-dyy.compute.c" +#define COEFFS_DYZ_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-dyz.compute.c" +#define COEFFS_DZZ_COMPUTE_FILE_NAME "Lagrange-tensor-product/3d.coeffs/3d.cube.order4.smooth0/coeffs-dzz.compute.c" + +/* actual code */ +#include "../template.c" |