diff options
author | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2008-04-18 21:40:46 +0000 |
---|---|---|
committer | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2008-04-18 21:40:46 +0000 |
commit | 6fdbf297784da5711d8801877103cd7b4f43f383 (patch) | |
tree | 368fb8943e537bb957a6f9a7af5ea3bb2f7fd4f6 /src/All_Coeffs_mod.F90 | |
parent | 0616ca3051e71b2a90c641755ea0e20e20e0f51e (diff) |
Add stencils for non-SBP operators.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@104 f69c4107-0314-4c4f-9ad4-17e986b73f4a
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 |