diff options
Diffstat (limited to 'src/All_Coeffs_mod.F90')
| -rw-r--r-- | src/All_Coeffs_mod.F90 | 41 |
1 files changed, 40 insertions, 1 deletions
diff --git a/src/All_Coeffs_mod.F90 b/src/All_Coeffs_mod.F90 index ef965d8..76dfcc9 100644 --- a/src/All_Coeffs_mod.F90 +++ b/src/All_Coeffs_mod.F90 @@ -194,6 +194,45 @@ module All_Coeffs_mod CCTK_REAL, dimension(3), intent(OUT) :: a CCTK_REAL, dimension(6,4), intent(OUT) :: q + a(1) = -2.500000000000000000000000000000000000000_wp + a(2) = 1.333333333333333333333333333333333333333_wp + a(3) = -0.08333333333333333333333333333333333333333_wp + + q(1,1) = 2.000000000000000000000000000000000000000_wp + q(2,1) = -5.000000000000000000000000000000000000000_wp + q(3,1) = 4.000000000000000000000000000000000000000_wp + q(4,1) = -1.000000000000000000000000000000000000000_wp + q(5,1) = 0_wp + q(6,1) = 0_wp + + q(1,2) = 1.000000000000000000000000000000000000000_wp + q(2,2) = -2.000000000000000000000000000000000000000_wp + q(3,2) = 1.000000000000000000000000000000000000000_wp + q(4,2) = 0_wp + q(5,2) = 0_wp + q(6,2) = 0_wp + + q(1,3) = -0.09302325581395348837209302325581395348837_wp + q(2,3) = 1.372093023255813953488372093023255813953_wp + q(3,3) = -2.558139534883720930232558139534883720930_wp + q(4,3) = 1.372093023255813953488372093023255813953_wp + q(5,3) = -0.09302325581395348837209302325581395348837_wp + q(6,3) = 0_wp + + q(1,4) = -0.02040816326530612244897959183673469387755_wp + q(2,4) = 0_wp + q(3,4) = 1.204081632653061224489795918367346938776_wp + q(4,4) = -2.408163265306122448979591836734693877551_wp + q(5,4) = 1.306122448979591836734693877551020408163_wp + q(6,4) = -0.08163265306122448979591836734693877551020_wp + + end subroutine coeffs_2_4_2 + + subroutine coeffs_2_4_2_opt ( a, q ) + + CCTK_REAL, dimension(3), intent(OUT) :: a + CCTK_REAL, dimension(6,4), intent(OUT) :: q + a(1) = -2.5_wp a(2) = 1.333333333333333333333333333333333333333_wp a(3) = -0.08333333333333333333333333333333333333333_wp @@ -226,7 +265,7 @@ module All_Coeffs_mod q(5,4) = -0.9258553965216912759079296546576971585129_wp q(6,4) = 0.4150007042443253806547133319378374518503_wp - end subroutine coeffs_2_4_2 + end subroutine coeffs_2_4_2_opt subroutine coeffs_1_6_3 ( a, q ) |