diff options
author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-09-01 18:10:41 +0000 |
---|---|---|
committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-09-01 18:10:41 +0000 |
commit | f66e9bbb167e3e3e98b0647e30fa22139d6f28f0 (patch) | |
tree | 01b5d8d23679ef302f5df71c1a62fe1bb1b6628e /src/GeneralizedPolynomial-Uniform/common/3d.cube.size5/store-coeffs.c | |
parent | 5ae03060c304e68a184efaf4d77a6d46d3bd24ce (diff) |
reorganize common to Lagrange and Hermite interpolators
(i.e. code depending only on molecule size/shape)
--> now there are separate *functions* rather than just code fragments
--> ../template.c should hopefully compile in finite time.....
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@109 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/common/3d.cube.size5/store-coeffs.c')
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/common/3d.cube.size5/store-coeffs.c | 125 |
1 files changed, 125 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/common/3d.cube.size5/store-coeffs.c b/src/GeneralizedPolynomial-Uniform/common/3d.cube.size5/store-coeffs.c new file mode 100644 index 0000000..40c057c --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/common/3d.cube.size5/store-coeffs.c @@ -0,0 +1,125 @@ +COEFF(-2,-2,-2) = factor * coeffs->coeff_m2_m2_m2; +COEFF(-1,-2,-2) = factor * coeffs->coeff_m1_m2_m2; +COEFF(0,-2,-2) = factor * coeffs->coeff_0_m2_m2; +COEFF(1,-2,-2) = factor * coeffs->coeff_p1_m2_m2; +COEFF(2,-2,-2) = factor * coeffs->coeff_p2_m2_m2; +COEFF(-2,-1,-2) = factor * coeffs->coeff_m2_m1_m2; +COEFF(-1,-1,-2) = factor * coeffs->coeff_m1_m1_m2; +COEFF(0,-1,-2) = factor * coeffs->coeff_0_m1_m2; +COEFF(1,-1,-2) = factor * coeffs->coeff_p1_m1_m2; +COEFF(2,-1,-2) = factor * coeffs->coeff_p2_m1_m2; +COEFF(-2,0,-2) = factor * coeffs->coeff_m2_0_m2; +COEFF(-1,0,-2) = factor * coeffs->coeff_m1_0_m2; +COEFF(0,0,-2) = factor * coeffs->coeff_0_0_m2; +COEFF(1,0,-2) = factor * coeffs->coeff_p1_0_m2; +COEFF(2,0,-2) = factor * coeffs->coeff_p2_0_m2; +COEFF(-2,1,-2) = factor * coeffs->coeff_m2_p1_m2; +COEFF(-1,1,-2) = factor * coeffs->coeff_m1_p1_m2; +COEFF(0,1,-2) = factor * coeffs->coeff_0_p1_m2; +COEFF(1,1,-2) = factor * coeffs->coeff_p1_p1_m2; +COEFF(2,1,-2) = factor * coeffs->coeff_p2_p1_m2; +COEFF(-2,2,-2) = factor * coeffs->coeff_m2_p2_m2; +COEFF(-1,2,-2) = factor * coeffs->coeff_m1_p2_m2; +COEFF(0,2,-2) = factor * coeffs->coeff_0_p2_m2; +COEFF(1,2,-2) = factor * coeffs->coeff_p1_p2_m2; +COEFF(2,2,-2) = factor * coeffs->coeff_p2_p2_m2; +COEFF(-2,-2,-1) = factor * coeffs->coeff_m2_m2_m1; +COEFF(-1,-2,-1) = factor * coeffs->coeff_m1_m2_m1; +COEFF(0,-2,-1) = factor * coeffs->coeff_0_m2_m1; +COEFF(1,-2,-1) = factor * coeffs->coeff_p1_m2_m1; +COEFF(2,-2,-1) = factor * coeffs->coeff_p2_m2_m1; +COEFF(-2,-1,-1) = factor * coeffs->coeff_m2_m1_m1; +COEFF(-1,-1,-1) = factor * coeffs->coeff_m1_m1_m1; +COEFF(0,-1,-1) = factor * coeffs->coeff_0_m1_m1; +COEFF(1,-1,-1) = factor * coeffs->coeff_p1_m1_m1; +COEFF(2,-1,-1) = factor * coeffs->coeff_p2_m1_m1; +COEFF(-2,0,-1) = factor * coeffs->coeff_m2_0_m1; +COEFF(-1,0,-1) = factor * coeffs->coeff_m1_0_m1; +COEFF(0,0,-1) = factor * coeffs->coeff_0_0_m1; +COEFF(1,0,-1) = factor * coeffs->coeff_p1_0_m1; +COEFF(2,0,-1) = factor * coeffs->coeff_p2_0_m1; +COEFF(-2,1,-1) = factor * coeffs->coeff_m2_p1_m1; +COEFF(-1,1,-1) = factor * coeffs->coeff_m1_p1_m1; +COEFF(0,1,-1) = factor * coeffs->coeff_0_p1_m1; +COEFF(1,1,-1) = factor * coeffs->coeff_p1_p1_m1; +COEFF(2,1,-1) = factor * coeffs->coeff_p2_p1_m1; +COEFF(-2,2,-1) = factor * coeffs->coeff_m2_p2_m1; +COEFF(-1,2,-1) = factor * coeffs->coeff_m1_p2_m1; +COEFF(0,2,-1) = factor * coeffs->coeff_0_p2_m1; +COEFF(1,2,-1) = factor * coeffs->coeff_p1_p2_m1; +COEFF(2,2,-1) = factor * coeffs->coeff_p2_p2_m1; +COEFF(-2,-2,0) = factor * coeffs->coeff_m2_m2_0; +COEFF(-1,-2,0) = factor * coeffs->coeff_m1_m2_0; +COEFF(0,-2,0) = factor * coeffs->coeff_0_m2_0; +COEFF(1,-2,0) = factor * coeffs->coeff_p1_m2_0; +COEFF(2,-2,0) = factor * coeffs->coeff_p2_m2_0; +COEFF(-2,-1,0) = factor * coeffs->coeff_m2_m1_0; +COEFF(-1,-1,0) = factor * coeffs->coeff_m1_m1_0; +COEFF(0,-1,0) = factor * coeffs->coeff_0_m1_0; +COEFF(1,-1,0) = factor * coeffs->coeff_p1_m1_0; +COEFF(2,-1,0) = factor * coeffs->coeff_p2_m1_0; +COEFF(-2,0,0) = factor * coeffs->coeff_m2_0_0; +COEFF(-1,0,0) = factor * coeffs->coeff_m1_0_0; +COEFF(0,0,0) = factor * coeffs->coeff_0_0_0; +COEFF(1,0,0) = factor * coeffs->coeff_p1_0_0; +COEFF(2,0,0) = factor * coeffs->coeff_p2_0_0; +COEFF(-2,1,0) = factor * coeffs->coeff_m2_p1_0; +COEFF(-1,1,0) = factor * coeffs->coeff_m1_p1_0; +COEFF(0,1,0) = factor * coeffs->coeff_0_p1_0; +COEFF(1,1,0) = factor * coeffs->coeff_p1_p1_0; +COEFF(2,1,0) = factor * coeffs->coeff_p2_p1_0; +COEFF(-2,2,0) = factor * coeffs->coeff_m2_p2_0; +COEFF(-1,2,0) = factor * coeffs->coeff_m1_p2_0; +COEFF(0,2,0) = factor * coeffs->coeff_0_p2_0; +COEFF(1,2,0) = factor * coeffs->coeff_p1_p2_0; +COEFF(2,2,0) = factor * coeffs->coeff_p2_p2_0; +COEFF(-2,-2,1) = factor * coeffs->coeff_m2_m2_p1; +COEFF(-1,-2,1) = factor * coeffs->coeff_m1_m2_p1; +COEFF(0,-2,1) = factor * coeffs->coeff_0_m2_p1; +COEFF(1,-2,1) = factor * coeffs->coeff_p1_m2_p1; +COEFF(2,-2,1) = factor * coeffs->coeff_p2_m2_p1; +COEFF(-2,-1,1) = factor * coeffs->coeff_m2_m1_p1; +COEFF(-1,-1,1) = factor * coeffs->coeff_m1_m1_p1; +COEFF(0,-1,1) = factor * coeffs->coeff_0_m1_p1; +COEFF(1,-1,1) = factor * coeffs->coeff_p1_m1_p1; +COEFF(2,-1,1) = factor * coeffs->coeff_p2_m1_p1; +COEFF(-2,0,1) = factor * coeffs->coeff_m2_0_p1; +COEFF(-1,0,1) = factor * coeffs->coeff_m1_0_p1; +COEFF(0,0,1) = factor * coeffs->coeff_0_0_p1; +COEFF(1,0,1) = factor * coeffs->coeff_p1_0_p1; +COEFF(2,0,1) = factor * coeffs->coeff_p2_0_p1; +COEFF(-2,1,1) = factor * coeffs->coeff_m2_p1_p1; +COEFF(-1,1,1) = factor * coeffs->coeff_m1_p1_p1; +COEFF(0,1,1) = factor * coeffs->coeff_0_p1_p1; +COEFF(1,1,1) = factor * coeffs->coeff_p1_p1_p1; +COEFF(2,1,1) = factor * coeffs->coeff_p2_p1_p1; +COEFF(-2,2,1) = factor * coeffs->coeff_m2_p2_p1; +COEFF(-1,2,1) = factor * coeffs->coeff_m1_p2_p1; +COEFF(0,2,1) = factor * coeffs->coeff_0_p2_p1; +COEFF(1,2,1) = factor * coeffs->coeff_p1_p2_p1; +COEFF(2,2,1) = factor * coeffs->coeff_p2_p2_p1; +COEFF(-2,-2,2) = factor * coeffs->coeff_m2_m2_p2; +COEFF(-1,-2,2) = factor * coeffs->coeff_m1_m2_p2; +COEFF(0,-2,2) = factor * coeffs->coeff_0_m2_p2; +COEFF(1,-2,2) = factor * coeffs->coeff_p1_m2_p2; +COEFF(2,-2,2) = factor * coeffs->coeff_p2_m2_p2; +COEFF(-2,-1,2) = factor * coeffs->coeff_m2_m1_p2; +COEFF(-1,-1,2) = factor * coeffs->coeff_m1_m1_p2; +COEFF(0,-1,2) = factor * coeffs->coeff_0_m1_p2; +COEFF(1,-1,2) = factor * coeffs->coeff_p1_m1_p2; +COEFF(2,-1,2) = factor * coeffs->coeff_p2_m1_p2; +COEFF(-2,0,2) = factor * coeffs->coeff_m2_0_p2; +COEFF(-1,0,2) = factor * coeffs->coeff_m1_0_p2; +COEFF(0,0,2) = factor * coeffs->coeff_0_0_p2; +COEFF(1,0,2) = factor * coeffs->coeff_p1_0_p2; +COEFF(2,0,2) = factor * coeffs->coeff_p2_0_p2; +COEFF(-2,1,2) = factor * coeffs->coeff_m2_p1_p2; +COEFF(-1,1,2) = factor * coeffs->coeff_m1_p1_p2; +COEFF(0,1,2) = factor * coeffs->coeff_0_p1_p2; +COEFF(1,1,2) = factor * coeffs->coeff_p1_p1_p2; +COEFF(2,1,2) = factor * coeffs->coeff_p2_p1_p2; +COEFF(-2,2,2) = factor * coeffs->coeff_m2_p2_p2; +COEFF(-1,2,2) = factor * coeffs->coeff_m1_p2_p2; +COEFF(0,2,2) = factor * coeffs->coeff_0_p2_p2; +COEFF(1,2,2) = factor * coeffs->coeff_p1_p2_p2; +COEFF(2,2,2) = factor * coeffs->coeff_p2_p2_p2; |