diff options
Diffstat (limited to 'src/All_Coeffs_mod.F90')
-rw-r--r-- | src/All_Coeffs_mod.F90 | 106 |
1 files changed, 106 insertions, 0 deletions
diff --git a/src/All_Coeffs_mod.F90 b/src/All_Coeffs_mod.F90 index 1b773ad..539e46e 100644 --- a/src/All_Coeffs_mod.F90 +++ b/src/All_Coeffs_mod.F90 @@ -76,6 +76,18 @@ module All_Coeffs_mod end subroutine coeffs_2_2_1 + subroutine coeffs_2_2 ( a, q ) + + CCTK_REAL, dimension(2), intent(OUT) :: a + CCTK_REAL, dimension(1,1), intent(OUT) :: q + + a(1) = -2.0_wp + a(2) = 1.0_wp + + q(1,1) = 0.0_wp + + end subroutine coeffs_2_2 + subroutine coeffs_1_4_2 ( a, q ) CCTK_REAL, dimension(2), intent(OUT) :: a @@ -267,6 +279,25 @@ module All_Coeffs_mod end subroutine coeffs_2_4_2_opt + subroutine coeffs_2_4 ( a, q ) + + CCTK_REAL, dimension(3), intent(OUT) :: a + CCTK_REAL, dimension(3,2), intent(OUT) :: q + + a(1) = -2.5_wp + a(2) = 1.333333333333333333333333333333333333333_wp + a(3) = -0.08333333333333333333333333333333333333333_wp + + q(1,1) = zero + q(2,1) = zero + q(3,1) = zero + + q(1,2) = 1.0_wp + q(2,2) = -2.0_wp + q(3,2) = 1.0_wp + + end subroutine coeffs_2_4 + subroutine coeffs_1_6_3 ( a, q ) CCTK_REAL, dimension(3), intent(OUT) :: a @@ -622,6 +653,36 @@ module All_Coeffs_mod end subroutine coeffs_2_6_3 + subroutine coeffs_2_6 ( a, q ) + + CCTK_REAL, dimension(4), intent(OUT) :: a + CCTK_REAL, dimension(5,3), intent(OUT) :: q + + a(1) = -2.722222222222222222222222222222222222222_wp + a(2) = 1.500000000000000000000000000000000000000_wp + a(3) = -0.1500000000000000000000000000000000000000_wp + a(4) = 0.01111111111111111111111111111111111111111_wp + + q(1,1) = zero + q(2,1) = zero + q(3,1) = zero + q(4,1) = zero + q(5,1) = zero + + q(1,2) = 1.0_wp + q(2,2) = -2.0_wp + q(3,2) = 1.0_wp + q(4,2) = zero + q(5,2) = zero + + q(1,3) = -0.08333333333333333333333333333333333333333_wp + q(2,3) = 1.333333333333333333333333333333333333333_wp + q(3,3) = -2.5_wp + q(4,3) = 1.333333333333333333333333333333333333333_wp + q(5,3) = -0.08333333333333333333333333333333333333333_wp + + end subroutine coeffs_2_6 + subroutine coeffs_1_8_4 ( a, q ) CCTK_REAL, dimension(4), intent(OUT) :: a @@ -1250,6 +1311,51 @@ module All_Coeffs_mod end subroutine coeffs_2_8_4 + subroutine coeffs_2_8 ( a, q ) + + CCTK_REAL, dimension(5), intent(OUT) :: a + CCTK_REAL, dimension(7,4), intent(OUT) :: q + + a(1) = -2.847222222222222222222222222222222222222_wp + a(2) = 1.600000000000000000000000000000000000000_wp + a(3) = -0.2000000000000000000000000000000000000000_wp + a(4) = 0.02539682539682539682539682539682539682540_wp + a(5) = -0.001785714285714285714285714285714285714286_wp + + q(1,1) = zero + q(2,1) = zero + q(3,1) = zero + q(4,1) = zero + q(5,1) = zero + q(6,1) = zero + q(7,1) = zero + + q(1,2) = 1.0_wp + q(2,2) = -2.0_wp + q(3,2) = 1.0_wp + q(4,2) = zero + q(5,2) = zero + q(6,2) = zero + q(7,2) = zero + + q(1,3) = -0.08333333333333333333333333333333333333333_wp + q(2,3) = 1.333333333333333333333333333333333333333_wp + q(3,3) = -2.5_wp + q(4,3) = 1.333333333333333333333333333333333333333_wp + q(5,3) = -0.08333333333333333333333333333333333333333_wp + q(6,3) = zero + q(7,3) = zero + + q(1,4) = 0.01111111111111111111111111111111111111111_wp + q(2,4) = -0.15_wp + q(3,4) = 1.5_wp + q(4,4) = -2.722222222222222222222222222222222222222_wp + q(5,4) = 1.5_wp + q(6,4) = -0.15_wp + q(7,4) = 0.01111111111111111111111111111111111111111_wp + + end subroutine coeffs_2_8 + subroutine coeffs_1_4_3 ( a, q ) CCTK_REAL, dimension(2), intent(OUT) :: a |