aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat
-rw-r--r--interface.ccl18
-rw-r--r--src/Coefficients_2_1.F9062
-rw-r--r--src/Coefficients_4_2.F9086
-rw-r--r--src/Coefficients_4_3.F90133
-rw-r--r--src/Coefficients_6_3.F90112
-rw-r--r--src/Coefficients_6_5.F90162
-rw-r--r--src/Coefficients_8_4.F90231
-rw-r--r--src/get_coeffs.c135
8 files changed, 930 insertions, 9 deletions
diff --git a/interface.ccl b/interface.ccl
index 8b4304d..e3c2d98 100644
--- a/interface.ccl
+++ b/interface.ccl
@@ -4,15 +4,6 @@
implements: SummationByParts
inherits: grid
-#SUBROUTINE Diff_gf ( CCTK_REAL IN ARRAY var, \
-# CCTK_INT IN ni, \
-# CCTK_INT IN nj, \
-# CCTK_INT IN nk, \
-# CCTK_INT IN dir, \
-# CCTK_REAL IN delta, \
-# CCTK_REAL OUT ARRAY dvar )
-#PROVIDES FUNCTION Diff_gf WITH deriv_gf LANGUAGE Fortran
-
SUBROUTINE Diff_gf ( CCTK_POINTER_TO_CONST IN cctkGH, \
CCTK_INT IN dir, \
CCTK_STRING IN var_name, \
@@ -25,6 +16,15 @@ SUBROUTINE Diff_gv ( CCTK_POINTER_TO_CONST IN cctkGH, \
CCTK_REAL OUT ARRAY dvar )
PROVIDES FUNCTION Diff_gv WITH DiffGv LANGUAGE C
+SUBROUTINE Diff_coeff ( CCTK_POINTER_TO_CONST IN cctkGH, \
+ CCTK_INT IN dir, \
+ CCTK_INT IN nsize, \
+ CCTK_INT OUT ARRAY imin, \
+ CCTK_INT OUT ARRAY imax, \
+ CCTK_REAL OUT ARRAY q, \
+ CCTK_INT IN table_handle )
+PROVIDES FUNCTION Diff_coeff WITH DiffCoeff LANGUAGE C
+
CCTK_REAL FUNCTION GetScalProdCoeff ()
PROVIDES FUNCTION GetScalProdCoeff WITH GetCoeff LANGUAGE Fortran
diff --git a/src/Coefficients_2_1.F90 b/src/Coefficients_2_1.F90
new file mode 100644
index 0000000..f8fe61e
--- /dev/null
+++ b/src/Coefficients_2_1.F90
@@ -0,0 +1,62 @@
+#include "cctk.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+subroutine set_coeff_2_1 ( nsize, loc_order, bb, gsize, imin, imax, dd )
+
+ implicit none
+
+ DECLARE_CCTK_FUNCTIONS
+ DECLARE_CCTK_PARAMETERS
+
+ CCTK_REAL, parameter :: zero = 0.0
+ integer, parameter :: wp = kind(zero)
+ CCTK_INT, intent(IN) :: nsize, loc_order
+ CCTK_INT, dimension(2), intent(IN) :: bb
+ CCTK_INT, intent(IN) :: gsize
+ CCTK_INT, dimension(nsize), intent(OUT) :: imin, imax
+ CCTK_REAL, dimension(nsize,nsize), intent(OUT) :: dd
+
+ CCTK_REAL, dimension(1), save :: a
+ CCTK_REAL, dimension(3,2), save :: q
+
+ CCTK_INT :: i, il, ir
+
+ logical, save :: first = .true.
+
+ if ( first ) then
+ a(1) = 1.0_wp/2.0_wp
+
+ q(1,1) = -1.0_wp; q(2,1) = 1.0_wp; q(3,1) = zero
+ q(1,2) = -1.0_wp/2.0_wp; q(2,2) = zero; q(3,2) = 1.0_wp/2.0_wp
+ first = .false.
+ end if
+
+ dd = zero
+ imin = 0
+ imax = -1
+
+ if ( bb(1) == 0 ) then
+ il = 1 + gsize
+ else
+ dd(1,1) = q(1,1); dd(2,1) = q(2,1)
+ imin(1) = 1; imax(1) = 2
+ dd(1,2) = q(1,2); dd(3,2) = q(3,2)
+ imin(2) = 1; imax(2) = 3
+ il = 3
+ end if
+ if ( bb(2) == 0 ) then
+ ir = nsize - gsize
+ else
+ dd(nsize-2,nsize-1) = -q(3,2); dd(nsize,nsize-1) = -q(1,2)
+ imin(nsize-1) = nsize-2; imax(nsize-1) = nsize
+ dd(nsize-1,nsize) = -q(2,1); dd(nsize,nsize) = -q(1,1)
+ imin(nsize) = nsize-1; imax(nsize) = nsize
+ ir = nsize - 2
+ end if
+ do i = il, ir
+ dd(i-1,i) = -a(1); dd(i+1,i) = a(1)
+ imin(i) = i-1; imax(i) = i+1
+ end do
+
+end subroutine set_coeff_2_1
diff --git a/src/Coefficients_4_2.F90 b/src/Coefficients_4_2.F90
new file mode 100644
index 0000000..778da65
--- /dev/null
+++ b/src/Coefficients_4_2.F90
@@ -0,0 +1,86 @@
+#include "cctk.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+subroutine set_coeff_4_2 ( nsize, loc_order, bb, gsize, imin, imax, dd )
+
+ implicit none
+
+ DECLARE_CCTK_FUNCTIONS
+ DECLARE_CCTK_PARAMETERS
+
+ CCTK_REAL, parameter :: zero = 0.0
+ integer, parameter :: wp = kind(zero)
+ CCTK_INT, intent(IN) :: nsize, loc_order
+ CCTK_INT, dimension(2), intent(IN) :: bb
+ CCTK_INT, intent(IN) :: gsize
+ CCTK_INT, dimension(nsize), intent(OUT) :: imin, imax
+ CCTK_REAL, dimension(nsize,nsize), intent(OUT) :: dd
+
+ CCTK_REAL, dimension(2), save :: a
+ CCTK_REAL, dimension(6,4), save :: q
+
+ CCTK_INT :: i, il, ir
+
+ logical, save :: first = .true.
+
+ if ( first ) then
+ a(1) = 2.0_wp/3.0_wp; a(2) = -1.0_wp/12.0_wp
+
+ q(1,1) = -24.0_wp/17.0_wp; q(2,1) = 59.0_wp/34.0_wp
+ q(3,1) = -4.0_wp/17.0_wp; q(4,1) = -3.0_wp/34.0_wp
+ q(5,1) = zero; q(6,1) = zero
+ q(1,2) = -1.0_wp/2.0_wp; q(2,2) = zero
+ q(3,2) = 1.0_wp/2.0_wp; q(4,2) = zero
+ q(5,2) = zero; q(6,2) = zero
+ q(1,3) = 4.0_wp/43.0_wp; q(2,3) = -59.0_wp/86.0_wp
+ q(3,3) = zero; q(4,3) = 59.0_wp/86.0_wp
+ q(5,3) = -4.0_wp/43.0_wp; q(6,3) = zero
+ q(1,4) = 3.0_wp/98.0_wp; q(2,4) = zero
+ q(3,4) = -59.0_wp/98.0_wp; q(4,4) = zero
+ q(5,4) = 32.0_wp/49.0_wp; q(6,4) = -4.0_wp/49.0_wp
+ first = .false.
+ end if
+
+ dd = zero
+ imin = 0
+ imax = -1
+
+ if ( bb(1) == 0 ) then
+ il = 1 + gsize
+ else
+ dd(1:4,1) = q(1:4,1)
+ imin(1) = 1; imax(1) = 4
+
+ dd(1:3,2) = q(1:3,2)
+ imin(2) = 1; imax(2) = 3
+
+ dd(1:5,3) = q(1:5,3)
+ imin(3) = 1; imax(3) = 5
+
+ dd(1:6,4) = q(1:6,4)
+ imin(4) = 1; imax(4) = 6
+
+ il = 5
+ end if
+ if ( bb(2) == 0 ) then
+ ir = nsize - gsize
+ else
+ dd(nsize-5:nsize,nsize-3) = -q(6:1:-1,4)
+ imin(nsize-3) = nsize-5; imax(nsize-3) = nsize
+
+ dd(nsize-4:nsize,nsize-2) = -q(5:1:-1,3)
+ imin(nsize-2) = nsize-4; imax(nsize-2) = nsize
+
+ dd(nsize-2:nsize,nsize-1) = -q(3:1:-1,2)
+ imin(nsize-1) = nsize-2; imax(nsize-1) = nsize
+
+ dd(nsize-3:nsize,nsize) = -q(4:1:-1,1)
+ imin(nsize) = nsize-3; imax(nsize) = nsize
+ ir = nsize - 4
+ end if
+ do i = il, ir
+ dd(i-2:i-1,i) = -a(2:1:-1); dd(i+1:i+2,i) = a(1:2)
+ imin(i) = i-2; imax(i) = i+2
+ end do
+end subroutine set_coeff_4_2
diff --git a/src/Coefficients_4_3.F90 b/src/Coefficients_4_3.F90
new file mode 100644
index 0000000..89c824d
--- /dev/null
+++ b/src/Coefficients_4_3.F90
@@ -0,0 +1,133 @@
+#include "cctk.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+subroutine set_coeff_4_3 ( nsize, loc_order, bb, gsize, imin, imax, dd )
+
+ implicit none
+
+ DECLARE_CCTK_FUNCTIONS
+ DECLARE_CCTK_PARAMETERS
+
+ CCTK_REAL, parameter :: zero = 0.0
+ integer, parameter :: wp = kind(zero)
+ CCTK_INT, intent(IN) :: nsize, loc_order
+ CCTK_INT, dimension(2), intent(IN) :: bb
+ CCTK_INT, intent(IN) :: gsize
+ CCTK_INT, dimension(nsize), intent(OUT) :: imin, imax
+ CCTK_REAL, dimension(nsize,nsize), intent(OUT) :: dd
+
+ CCTK_REAL, dimension(2), save :: a
+ CCTK_REAL, dimension(7,5), save :: q
+ CCTK_REAL :: f0, f1, f2, f3
+
+ CCTK_INT :: i, il, ir
+
+ logical, save :: first = .true.
+
+ if ( first ) then
+ a(1) = 2.0_wp/3.0_wp; a(2) = -1.0_wp/12.0_wp
+
+ f0 = sqrt(26116897.0_wp)
+ f1 = -56764003702447356523.0_wp + 8154993476273221.0_wp * f0
+ f2 = -55804550303.0_wp + 9650225.0_wp * f0
+ f3 = 3262210757.0_wp + 271861.0_wp * f0
+
+ q(1,1) = -11._wp/6.0_wp; q(2,1) = 3.0_wp
+ q(3,1) = -3.0_wp/2.0_wp; q(4,1) = 1.0_wp/3.0_wp
+ q(5,1) = zero; q(6,1) = zero; q(7,1) = zero
+
+ q(1,2) = -24.0_wp * ( -779042810827742869.0_wp + &
+ 104535124033147.0_wp * f0 ) / f1
+ q(2,2) = -( -176530817412806109689.0_wp + &
+ 29768274816875927.0_wp * f0 ) / ( 6.0_wp * f1 )
+ q(3,2) = 343.0_wp * ( -171079116122226871.0_wp + &
+ 27975630462649.0_wp * f0 ) / f1
+ q(4,2) = -3.0_wp * ( -7475554291248533227.0_wp + &
+ 1648464218793925.0_wp * f0 ) / ( 2.0_wp * f1 )
+ q(5,2) = ( -2383792768180030915.0_wp + &
+ 1179620587812973.0_wp * f0 ) / ( 3.0_wp * f1 )
+ q(6,2) = -1232.0_wp * ( -115724529581315.0_wp + 37280576429.0_wp * f0 ) / f1
+ q(7,2) = zero
+
+ q(1,3) = -12.0_wp * ( -380966843.0_wp + 86315.0_wp * f0 ) / f2
+ q(2,3) = ( 5024933015.0_wp + 2010631.0_wp * f0 ) / ( 3.0_wp * f2 )
+ q(3,3) = -231.0_wp * ( -431968921.0_wp + 86711.0_wp * f0 ) / ( 2.0_wp * f2 )
+ q(4,3) = ( -65931742559.0_wp + 12256337.0_wp * f0 ) / f2
+ q(5,3) = -( -50597298167.0_wp + 9716873.0_wp * f0 ) / ( 6.0_wp * f2 )
+ q(6,3) = -88.0_wp * ( -15453061.0_wp + 2911.0_wp * f0 ) / f2
+ q(7,3) = zero
+
+ q(1,4) = 48.0_wp * ( -56020909845192541.0_wp + &
+ 9790180507043.0_wp * f0 ) / f1
+ q(2,4) = ( -9918249049237586011.0_wp + &
+ 1463702013196501.0_wp * f0 ) / ( 6.0_wp * f1)
+ q(3,4) = -13.0_wp * ( -4130451756851441723.0_wp + &
+ 664278707201077.0_wp * f0 ) / f1
+ q(4,4) = 3.0_wp * ( -26937108467782666617.0_wp + &
+ 5169063172799767.0_wp * f0 ) / ( 2.0_wp * f1 )
+ q(5,4) = -( 6548308508012371315.0_wp + &
+ 3968886380989379.0_wp * f0 ) / ( 3.0_wp * f1 )
+ q(6,4) = 88.0_wp * ( -91337851897923397.0_wp + &
+ 19696768305507.0_wp * f0 ) / f1
+ q(7,4) = 242.0_wp * ( -120683.0_wp + 15.0_wp * f0 ) / f3
+
+ q(1,5) = 264.0_wp * ( -120683.0_wp + 15.0_wp * f0 ) / f3
+ q(2,5) = ( -43118111.0_wp + 23357.0_wp * f0 ) / ( 3.0_wp * f3 )
+ q(3,5) = -47.0_wp * ( -28770085.0_wp + 2259.0_wp * f0 ) / ( 2.0_wp * f3 )
+ q(4,5) = -3.0_wp * ( 1003619433.0_wp + 11777.0_wp * f0 ) / f3
+ q(5,5) = -11.0_wp * ( -384168269.0_wp + 65747.0_wp * f0 ) / ( 6.0_wp * f3 )
+ q(6,5) = 22.0_wp * ( 87290207.0_wp + 10221.0_wp * f0 ) / f3
+ q(7,5) = -66.0_wp * ( 3692405.0_wp + 419.0_wp * f0 ) / f3
+ first = .false.
+ end if
+
+ dd = zero
+ imin = 0
+ imax = -1
+
+ if ( bb(1) == 0 ) then
+ il = 1 + gsize
+ else
+ dd(1:4,1) = q(1:4,1)
+ imin(1) = 1; imax(1) = 4
+
+ dd(1:6,2) = q(1:6,2)
+ imin(2) = 1; imax(2) = 6
+
+ dd(1:6,3) = q(1:6,3)
+ imin(3) = 1; imax(3) = 6
+
+ dd(1:7,4) = q(1:7,4)
+ imin(4) = 1; imax(4) = 7
+
+ dd(1:7,5) = q(1:7,5)
+ imin(5) = 1; imax(5) = 7
+
+ il = 6
+ end if
+ if ( bb(2) == 0 ) then
+ ir = nsize - gsize
+ else
+ dd(nsize-6:nsize,nsize-4) = -q(7:1:-1,5)
+ imin(nsize-4) = nsize-6; imax(nsize-4) = nsize
+
+ dd(nsize-6:nsize,nsize-3) = -q(7:1:-1,4)
+ imin(nsize-3) = nsize-6; imax(nsize-3) = nsize
+
+ dd(nsize-5:nsize,nsize-2) = -q(6:1:-1,3)
+ imin(nsize-2) = nsize-5; imax(nsize-2) = nsize
+
+ dd(nsize-5:nsize,nsize-1) = -q(6:1:-1,2)
+ imin(nsize-1) = nsize-5; imax(nsize-1) = nsize
+
+ dd(nsize-3:nsize,nsize) = -q(4:1:-1,1)
+ imin(nsize) = nsize-3; imax(nsize) = nsize
+
+ ir = nsize - 5
+ end if
+ do i = il, ir
+ dd(i-2:i-1,i) = -a(2:1:-1); dd(i+1:i+2,i) = a(1:2)
+ imin(i) = i-2; imax(i) = i+2
+ end do
+end subroutine set_coeff_4_3
diff --git a/src/Coefficients_6_3.F90 b/src/Coefficients_6_3.F90
new file mode 100644
index 0000000..958d170
--- /dev/null
+++ b/src/Coefficients_6_3.F90
@@ -0,0 +1,112 @@
+#include "cctk.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+subroutine set_coeff_6_3 ( nsize, loc_order, bb, gsize, imin, imax, dd )
+
+ implicit none
+
+ DECLARE_CCTK_FUNCTIONS
+ DECLARE_CCTK_PARAMETERS
+
+ CCTK_REAL, parameter :: zero = 0.0
+ integer, parameter :: wp = kind(zero)
+ CCTK_INT, intent(IN) :: nsize, loc_order
+ CCTK_INT, dimension(2), intent(IN) :: bb
+ CCTK_INT, intent(IN) :: gsize
+ CCTK_INT, dimension(nsize), intent(OUT) :: imin, imax
+ CCTK_REAL, dimension(nsize,nsize), intent(OUT) :: dd
+
+ CCTK_REAL, dimension(3), save :: a
+ CCTK_REAL, dimension(9,6), save :: q
+
+ CCTK_INT :: i, il, ir
+
+ logical, save :: first = .true.
+
+ if ( first ) then
+ a(1) = 3.0_wp/4.0_wp; a(2) = -3.0_wp/20.0_wp; a(3) = 1.0_wp/60.0_wp
+
+ q(1,1) = -21600.0_wp/13649.0_wp; q(2,1) = 81763.0_wp/40947.0_wp
+ q(3,1) = 131.0_wp/27298.0_wp; q(4,1) = -9143.0_wp/13649.0_wp
+ q(5,1) = 20539.0_wp/81894.0_wp; q(6,1) = zero;
+ q(7,1) = zero; q(8,1) = zero; q(9,1) = zero
+ q(1,2) = -81763.0_wp/180195.0_wp; q(2,2) = zero
+ q(3,2) = 7357.0_wp/36039.0_wp; q(4,2) = 30637.0_wp/72078.0_wp
+ q(5,2) = -2328.0_wp/12013.0_wp; q(6,2) = 6611.0_wp/360390.0_wp
+ q(7,2) = zero; q(8,2) = zero; q(9,2) = zero
+ q(1,3) = -131.0_wp/54220.0_wp; q(2,3) = -7357.0_wp/16266.0_wp
+ q(3,3) = zero; q(4,3) = 645.0_wp/2711.0_wp
+ q(5,3) = 11237.0_wp/32532.0_wp; q(6,3) = -3487.0_wp/27110.0_wp
+ q(7,3) = zero; q(8,3) = zero; q(9,3) = zero
+ q(1,4) = 9143.0_wp/53590.0_wp; q(2,4) = -30637.0_wp/64308.0_wp
+ q(3,4) = -645.0_wp/5359.0_wp; q(4,4) = zero
+ q(5,4) = 13733.0_wp/32154.0_wp; q(6,4) = -67.0_wp/4660.0_wp
+ q(7,4) = 72.0_wp/5359.0_wp; q(8,4) = zero; q(9,4) = zero;
+ q(1,5) = -20539.0_wp/236310.0_wp; q(2,5) = 2328.0_wp/7877.0_wp
+ q(3,5) = -11237.0_wp/47262.0_wp; q(4,5) = -13733.0_wp/23631.0_wp
+ q(5,5) = zero; q(6,5) = 89387.0_wp/118155.0_wp
+ q(7,5) = -1296.0_wp/7877.0_wp; q(8,5) = 144.0_wp/7877.0_wp; q(9,5) = zero
+ q(1,6) = zero; q(2,6) = -6611.0_wp/262806.0_wp
+ q(3,6) = 3487.0_wp/43801.0_wp; q(4,6) = 1541.0_wp/87602.0_wp
+ q(5,6) = -89387.0_wp/131403.0_wp; q(6,6) = zero
+ q(7,6) = 32400.0_wp/43801.0_wp; q(8,6) = -6480.0_wp/43801.0_wp
+ q(9,6) = 720.0_wp/43801.0_wp
+ first = .false.
+ end if
+
+ dd = zero
+ imin = 0
+ imax = -1
+
+ if ( bb(1) == 0 ) then
+ il = 1 + gsize
+ else
+ dd(1:5,1) = q(1:5,1)
+ imin(1) = 1; imax(1) = 5
+
+ dd(1:6,2) = q(1:6,2)
+ imin(2) = 1; imax(2) = 6
+
+ dd(1:6,3) = q(1:6,3)
+ imin(3) = 1; imax(3) = 6
+
+ dd(1:7,4) = q(1:7,4)
+ imin(4) = 1; imax(4) = 7
+
+ dd(1:8,5) = q(1:8,5)
+ imin(5) = 1; imax(5) = 8
+
+ dd(2:9,6) = q(2:9,6)
+ imin(6) = 2; imax(6) = 9
+
+ il = 7
+ end if
+ if ( bb(2) == 0 ) then
+ ir = nsize - gsize
+ else
+ dd(nsize-8:nsize-1,nsize-5) = -q(9:2:-1,6)
+ imin(nsize-5) = nsize-8; imax(nsize-5) = nsize-1
+
+ dd(nsize-7:nsize,nsize-4) = -q(8:1:-1,5)
+ imin(nsize-4) = nsize-7; imax(nsize-4) = nsize
+
+ dd(nsize-6:nsize,nsize-3) = -q(7:1:-1,4)
+ imin(nsize-3) = nsize-6; imax(nsize-3) = nsize
+
+ dd(nsize-5:nsize,nsize-2) = -q(6:1:-1,3)
+ imin(nsize-2) = nsize-5; imax(nsize-2) = nsize
+
+ dd(nsize-5:nsize,nsize-1) = -q(6:1:-1,2)
+ imin(nsize-1) = nsize-5; imax(nsize-1) = nsize
+
+ dd(nsize-4:nsize,nsize) = -q(5:1:-1,1)
+ imin(nsize) = nsize-4; imax(nsize) = nsize
+
+ ir = nsize - 6
+ end if
+ do i = il, ir
+ dd(i-3:i-1,i) = -a(3:1:-1); dd(i+1:i+3,i) = a(1:3)
+ imin(i) = i-3; imax(i) = i+3
+ end do
+end subroutine set_coeff_6_3
diff --git a/src/Coefficients_6_5.F90 b/src/Coefficients_6_5.F90
new file mode 100644
index 0000000..361ab5b
--- /dev/null
+++ b/src/Coefficients_6_5.F90
@@ -0,0 +1,162 @@
+#include "cctk.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+subroutine set_coeff_6_5 ( nsize, loc_order, bb, gsize, imin, imax, dd )
+
+ implicit none
+
+ DECLARE_CCTK_FUNCTIONS
+ DECLARE_CCTK_PARAMETERS
+
+ CCTK_REAL, parameter :: zero = 0.0
+ integer, parameter :: wp = kind(zero)
+ CCTK_INT, intent(IN) :: nsize, loc_order
+ CCTK_INT, dimension(2), intent(IN) :: bb
+ CCTK_INT, intent(IN) :: gsize
+ CCTK_INT, dimension(nsize), intent(OUT) :: imin, imax
+ CCTK_REAL, dimension(nsize,nsize), intent(OUT) :: dd
+
+ CCTK_REAL, dimension(3), save :: a
+ CCTK_REAL, dimension(10,7), save :: q
+
+ CCTK_INT :: i, il, ir
+
+ logical, save :: first = .true.
+
+ if ( first ) then
+ a(1) = 3.0_wp/4.0_wp; a(2) = -3.0_wp/20.0_wp; a(3) = 1.0_wp/60.0_wp
+
+ q(1,1) = -137.0_wp/60.0_wp; q(2,1) = 5.0_wp; q(3,1) = -5.0_wp;
+ q(4,1) = 10.0_wp/3.0_wp; q(5,1) = -5.0_wp/4.0_wp; q(6,1) = 1.0_wp/5.0_wp;
+ q(7,1) = zero; q(8,1) = zero; q(9,1) = zero; q(10,1) = zero;
+
+ q(1,2) = -0.3209098489859668118347_wp
+ q(2,2) = -0.359441670715467075894_wp
+ q(3,2) = 0.195756852998105503889_wp
+ q(4,2) = 1.394685510250317033169_wp
+ q(5,2) = -1.448965775497476572821_wp
+ q(6,2) = 0.6519476244467816674835_wp
+ q(7,2) = -0.1115052611983591304249_wp
+ q(8,2) = -0.00156743129793461356831_wp
+ q(9,2) = zero; q(10,2) = zero
+
+ q(1,3) = 0.0980796259621024870281_wp
+ q(2,3) = -0.784752101846956387607_wp
+ q(3,3) = 0.364206221878666751300_wp
+ q(4,3) = 0.102097753636344358379_wp
+ q(5,3) = 0.377167650934576412896_wp
+ q(6,3) = -0.173241400242683302020_wp
+ q(7,3) = 0.0062120424243610783650_wp
+ q(8,3) = 0.01153111791917461507898_wp
+ q(9,3) = -0.00130091066558601341956_wp
+ q(10,3) = zero
+
+ q(1,4) = 0.0771947815148234674634_wp
+ q(2,4) = -0.4124719432622107102831_wp
+ q(3,4) = 0.654798717867622763906_wp
+ q(4,4) = -1.873474657556531216053_wp
+ q(5,4) = 2.161961221935152129759_wp
+ q(6,4) = -0.729147814207596325277_wp
+ q(7,4) = 0.1292509053479867634674_wp
+ q(8,4) = -0.01092898521375837831427_wp
+ q(9,4) = 0.00316983426828354261009_wp
+ q(10,4) = -0.000352060693772037278516_wp
+
+ q(1,5) = -0.01135836398682271847110_wp
+ q(2,5) = 0.0114974985688122284189_wp
+ q(3,5) = 0.228765565489867047381_wp
+ q(4,5) = -1.179122428698534075609_wp
+ q(5,5) = 0.774618306056840861215_wp
+ q(6,5) = 0.016612403904600041087_wp
+ q(7,5) = 0.2399304273448163041220_wp
+ q(8,5) = -0.0959180672407791626759_wp
+ q(9,5) = 0.01612460763754670343830_wp
+ q(10,5) = -0.001149949076347228906108_wp
+
+ q(1,6) = -0.01912046962023373131259_wp
+ q(2,6) = 0.1569245474161812198493_wp
+ q(3,6) = -0.567250839387340928718_wp
+ q(4,6) = 1.230411310316952671353_wp
+ q(5,6) = -2.048795717924254691221_wp
+ q(6,6) = 0.982715833107347313532_wp
+ q(7,6) = 0.2876584784911885215053_wp
+ q(8,6) = -0.0207657038198929292487_wp
+ q(9,6) = -0.00335290957957120811713_wp
+ q(10,6) = 0.001575470999623762377395_wp
+
+ q(1,7) = -0.0076072322555736921424_wp
+ q(2,7) = 0.034503414345314189270_wp
+ q(3,7) = -0.044377233426986476366_wp
+ q(4,7) = -0.049408595803628143797_wp
+ q(5,7) = 0.304693537653175477283_wp
+ q(6,7) = -0.935863321236636595622_wp
+ q(7,7) = 0.111375967229143365899_wp
+ q(8,7) = 0.7149322477536284813995_wp
+ q(9,7) = -0.1444768161656941519630_wp
+ q(10,7) = 0.01622803190725754603783_wp
+
+ first = .false.
+ end if
+
+ dd = zero
+ imin = 0
+ imax = -1
+
+ if ( bb(1) == 0 ) then
+ il = 1 + gsize
+ else
+ dd(1:6,1) = q(1:6,1)
+ imin(1) = 1; imax(1) = 6
+
+ dd(1:8,2) = q(1:8,2)
+ imin(2) = 1; imax(2) = 8
+
+ dd(1:9,3) = q(1:9,3)
+ imin(3) = 1; imax(3) = 9
+
+ dd(1:10,4) = q(1:10,4)
+ imin(4) = 1; imax(4) = 10
+
+ dd(1:10,5) = q(1:10,5)
+ imin(5) = 1; imax(5) = 10
+
+ dd(1:10,6) = q(1:10,6)
+ imin(6) = 1; imax(6) = 10
+
+ dd(1:10,7) = q(1:10,7)
+ imin(7) = 1; imax(7) = 10
+
+ il = 8
+ end if
+ if ( bb(2) == 0 ) then
+ ir = nsize - gsize
+ else
+ dd(nsize-9:nsize,nsize-6) = -q(10:1:-1,7)
+ imin(nsize-6) = nsize-9; imax(nsize-6) = nsize
+
+ dd(nsize-9:nsize,nsize-5) = -q(10:1:-1,6)
+ imin(nsize-5) = nsize-9; imax(nsize-5) = nsize
+
+ dd(nsize-9:nsize,nsize-4) = -q(10:1:-1,5)
+ imin(nsize-4) = nsize-9; imax(nsize-4) = nsize
+
+ dd(nsize-9:nsize,nsize-3) = -q(10:1:-1,4)
+ imin(nsize-3) = nsize-9; imax(nsize-3) = nsize
+
+ dd(nsize-8:nsize,nsize-2) = -q(9:1:-1,3)
+ imin(nsize-2) = nsize-8; imax(nsize-2) = nsize
+
+ dd(nsize-7:nsize,nsize-1) = -q(8:1:-1,2)
+ imin(nsize-1) = nsize-7; imax(nsize-1) = nsize
+
+ dd(nsize-5:nsize,nsize) = -q(6:1:-1,1)
+ imin(nsize) = nsize-5; imax(nsize) = nsize
+
+ ir = nsize - 7
+ end if
+ do i = il, ir
+ dd(i-3:i-1,i) = -a(3:1:-1); dd(i+1:i+3,i) = a(1:3)
+ imin(i) = i-3; imax(i) = i+3
+ end do
+end subroutine set_coeff_6_5
diff --git a/src/Coefficients_8_4.F90 b/src/Coefficients_8_4.F90
new file mode 100644
index 0000000..eb029ba
--- /dev/null
+++ b/src/Coefficients_8_4.F90
@@ -0,0 +1,231 @@
+#include "cctk.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+subroutine set_coeff_8_4 ( nsize, loc_order, bb, gsize, imin, imax, dd )
+
+ implicit none
+
+ DECLARE_CCTK_FUNCTIONS
+ DECLARE_CCTK_PARAMETERS
+
+ CCTK_REAL, parameter :: zero = 0.0
+ integer, parameter :: wp = kind(zero)
+ CCTK_INT, intent(IN) :: nsize, loc_order
+ CCTK_INT, dimension(2), intent(IN) :: bb
+ CCTK_INT, intent(IN) :: gsize
+ CCTK_INT, dimension(nsize), intent(OUT) :: imin, imax
+ CCTK_REAL, dimension(nsize,nsize), intent(OUT) :: dd
+
+ CCTK_REAL, dimension(4), save :: a
+ CCTK_REAL, dimension(12,8), save :: q
+ CCTK_REAL, parameter :: x1 = 0.541_wp, x2 = -0.0675_wp, x3 = 0.748_wp
+
+ CCTK_INT :: i, il, ir
+
+ logical, save :: first = .true.
+
+ if ( first ) then
+ a(1) = 4.0_wp/5.0_wp; a(2) = -1.0_wp/5.0_wp
+ a(3) = 4.0_wp/105.0_wp; a(4) = -1.0_wp/280.0_wp
+
+ q(1,1) = -2540160.0_wp/1498139.0_wp
+ q(2,1) = 9.0_wp * ( 2257920.0_wp*x1 + 11289600.0_wp*x2 + &
+ 22579200.0_wp*x3 - 15849163.0_wp ) / 5992556.0_wp
+ q(3,1) = 3.0_wp * ( -33868800.0_wp*x1 - 162570240.0_wp*x2 - &
+ 304819200.0_wp*x3 + 235236677.0_wp ) / 5992556.0_wp
+ q(4,1) = ( 609638400.0_wp*x1 + 2743372800.0_wp*x2 + &
+ 4572288000.0_wp * x3 - 3577778591.0_wp ) / 17977668.0_wp
+ q(5,1) = 3.0_wp * ( -16934400*x1 - 67737600.0_wp*x2 - &
+ 84672000.0_wp*x3 + 67906303.0_wp ) / 1498139.0_wp
+ q(6,1) = 105.0_wp * ( 967680.0_wp*x1 + 2903040.0_wp*x2 - &
+ 305821.0_wp ) / 5992556.0_wp
+ q(7,1) = 49.0_wp * ( -1244160.0_wp*x1 + 18662400.0_wp * x3 - &
+ 13322233.0_wp ) / 17977668.0_wp
+ q(8,1) = 3.0_wp * ( -6773760.0_wp*x2 - 33868800.0_wp*x3 + &
+ 24839327.0_wp ) / 5992556.0_wp
+ q(9,1) = zero; q(10,1) = zero; q(11,1) = zero; q(12,1) = zero
+
+ q(1,2) = 9.0_wp * ( -2257920.0_wp*x1 - 11289600.0_wp*x2 - &
+ 22579200.0_wp*x3 + 15849163.0_wp ) / 31004596.0_wp
+ q(2,2) = zero
+ q(3,2) = 3.0_wp * ( 7257600.0_wp*x1 + 33868800.0_wp*x2 + &
+ 60963840.0_wp*x3 - 47167457.0_wp ) / 2214614.0_wp
+ q(4,2) = 3.0_wp * ( -9676800.0_wp*x1 - 42336000.0_wp*x2 - &
+ 67737600.0_wp*x3 + 53224573.0_wp ) / 1107307.0_wp
+ q(5,2) = 7.0_wp * ( 55987200.0_wp*x1 + 217728000.0_wp*x2 + &
+ 261273600.0_wp*x3 - 211102099.0_wp ) / 13287684.0_wp
+ q(6,2) = 3.0_wp * ( -11612160.0_wp*x1 - 33868800.0_wp*x2 + &
+ 3884117.0_wp ) / 2214614.0_wp
+ q(7,2) = 150.0_wp * ( 24192.0_wp*x1 - 338688.0_wp*x3 + &
+ 240463.0_wp ) / 1107307.0_wp
+ q(8,2) = ( 152409600.0_wp*x2 + 731566080.0_wp*x3 - &
+ 536324953.0_wp ) / 46506894.0_wp
+ q(9,2) = zero; q(10,2) = zero; q(11,2) = zero; q(12,2) = zero
+
+ q(1,3) = ( 33868800.0_wp*x1 + 162570240.0_wp*x2 + &
+ 304819200.0_wp*x3 - 235236677.0_wp ) / 1743924.0_wp
+ q(2,3) = ( -7257600.0_wp*x1 - 33868800.0_wp*x2 - &
+ 60963840.0_wp*x3 + 47167457.0_wp ) / 124566.0_wp
+ q(3,3) = zero
+ q(4,3) = ( 24192000.0_wp*x1 + 101606400.0_wp*x2 + &
+ 152409600.0_wp*x3 - 120219461.0_wp ) / 124566.0_wp
+ q(5,3) = ( -72576000.0_wp*x1 - 270950400.0_wp*x2 - &
+ 304819200.0_wp*x3 + 249289259.0_wp ) / 249132.0_wp
+ q(6,3) = 9.0_wp * ( 806400.0_wp*x1 + 2257920.0_wp*x2 - &
+ 290167.0_wp ) / 41522.0_wp
+ q(7,3) = 6.0_wp * ( -134400.0_wp*x1 + 1693440.0_wp*x3 - &
+ 1191611.0_wp ) / 20761.0_wp
+ q(8,3) = 5.0_wp * ( -2257920.0_wp*x2 - 10160640.0_wp*x3 + &
+ 7439833.0_wp ) / 290654.0_wp
+ q(9,3) = zero; q(10,3) = zero; q(11,3) = zero; q(12,3) = zero
+
+ q(1,4) = ( -609638400.0_wp*x1 - 2743372800.0_wp*x2 - &
+ 4572288000.0_wp*x3 + 3577778591.0_wp ) / 109619916.0_wp
+ q(2,4) = 3.0_wp * ( 9676800.0_wp*x1 + 42336000.0_wp*x2 + &
+ 67737600.0_wp*x3 - 53224573.0_wp ) / 1304999.0_wp
+ q(3,4) = 3.0_wp * ( -24192000.0_wp*x1 - 101606400.0_wp*x2 - &
+ 152409600.0_wp*x3 + 120219461.0_wp ) / 2609998.0_wp
+ q(4,4) = zero
+ q(5,4) = 9.0_wp * ( 16128000.0_wp*x1 + 56448000.0_wp*x2 + &
+ 56448000.0_wp*x3 - 47206049.0_wp ) / 5219996.0_wp
+ q(6,4) = 3.0_wp * ( -19353600.0_wp*x1 - 50803200.0_wp*x2 + &
+ 7628371.0_wp ) / 2609998.0_wp
+ q(7,4) = 2.0_wp * ( 10886400.0_wp*x1 - 114307200.0_wp*x3 + &
+ 79048289.0_wp ) / 3914997.0_wp
+ q(8,4) = 75.0_wp * ( 1354752.0_wp*x2 + 5419008.0_wp*x3 - &
+ 3952831.0_wp ) / 18269986.0_wp
+ q(9,4) = zero; q(10,4) = zero; q(11,4) = zero; q(12,4) = zero
+
+ q(1,5) = 3.0_wp * ( 16934400.0_wp*x1 + 67737600.0_wp*x2 + &
+ 84672000.0_wp*x3 - 67906303.0_wp ) / 2096689.0_wp
+ q(2,5) = 7.0_wp * ( -55987200.0_wp*x1 - 217728000.0_wp*x2 - &
+ 261273600.0_wp*x3 + 211102099.0_wp ) / 3594324.0_wp
+ q(3,5) = 3.0_wp * ( 72576000.0_wp*x1 + 270950400.0_wp*x2 + &
+ 304819200.0_wp*x3 - 249289259.0_wp ) / 1198108.0_wp
+ q(4,5) = 9.0_wp * ( -16128000.0_wp*x1 - 56448000.0_wp*x2 - &
+ 56448000.0_wp*x3 + 47206049.0_wp ) / 1198108.0_wp
+ q(5,5) = zero
+ q(6,5) = 105.0_wp * ( 414720.0_wp*x1 + 967680.0_wp*x2 - &
+ 165527.0_wp ) / 1198108.0_wp
+ q(7,5) = 15.0_wp * ( -967680.0_wp*x1 + 6773760.0_wp*x3 - &
+ 4472029.0_wp ) / 1198108.0_wp
+ q(8,5) = ( -304819200.0_wp*x2 - 914457600.0_wp*x3 + &
+ 657798011.0_wp ) / 25160268.0_wp
+ q(9,5) = -2592.0_wp/299527.0_wp;
+ q(10,5) = zero; q(11,5) = zero; q(12,5) = zero
+
+ q(1,6) = 5.0_wp * ( -967680.0_wp*x1 - 2903040.0_wp*x2 + &
+ 305821.0_wp ) / 1237164.0_wp
+ q(2,6) = ( 11612160.0_wp*x1 + 33868800.0_wp*x2 - &
+ 3884117.0_wp ) / 618582.0_wp
+ q(3,6) = 9.0_wp * ( -806400.0_wp*x1 - 2257920.0_wp*x2 + &
+ 290167.0_wp ) / 206194.0_wp
+ q(4,6) = ( 19353600.0_wp*x1 + 50803200.0_wp*x2 - &
+ 7628371.0_wp ) / 618582.0_wp
+ q(5,6) = 35.0_wp * ( -414720.0_wp*x1 - 967680.0_wp*x2 + &
+ 165527.0_wp ) / 1237164.0_wp
+ q(6,6) = zero; q(7,6) = 80640.0_wp*x1 / 103097.0_wp
+ q(8,6) = 80640.0_wp*x2 / 103097.0_wp; q(9,6) = 3072.0_wp / 103097.0_wp
+ q(10,6) = -288.0_wp/103097.0_wp; q(11,6) = zero; q(12,6) = zero
+
+ q(1,7) = 7.0_wp * ( 1244160.0_wp*x1 - 18662400.0_wp*x3 + &
+ 13322233.0_wp ) / 8041092.0_wp
+ q(2,7) = 150.0_wp * ( -24192.0_wp*x1 + 338688.0_wp*x3 - &
+ 240463.0_wp ) / 670091.0_wp
+ q(3,7) = 54.0_wp * ( 134400.0_wp*x1 - 1693440.0_wp*x3 + &
+ 1191611.0_wp ) / 670091.0_wp
+ q(4,7) = 2.0_wp * ( -10886400.0_wp*x1 + 114307200.0_wp*x3 - &
+ 79048289.0_wp ) / 2010273.0_wp
+ q(5,7) = 15.0_wp * ( 967680.0_wp*x1 - 6773760.0_wp*x3 + &
+ 4472029.0_wp ) / 2680364.0_wp
+ q(6,7) = -725760.0_wp*x1 / 670091.0_wp; q(7,7) = zero
+ q(8,7) = 725760.0_wp*x3 / 670091.0_wp; q(9,7) = -145152.0_wp/670091.0_wp
+ q(10,7) = 27648.0_wp/670091.0_wp; q(11,7) = -2592.0_wp/670091.0_wp
+ q(12,7) = zero
+
+ q(1,8) = 3.0_wp * ( 6773760.0_wp*x2 + 33868800.0_wp*x3 - &
+ 24839327.0_wp ) / 20510956.0_wp
+ q(2,8) = ( -152409600.0_wp*x2 - 731566080.0_wp*x3 + &
+ 536324953.0_wp ) / 30766434.0_wp
+ q(3,8) = 45.0_wp * ( 2257920.0_wp*x2 + 10160640.0_wp*x3 - &
+ 7439833.0_wp ) / 10255478.0_wp
+ q(4,8) = 75.0_wp * ( -1354752.0_wp*x2 - 5419008.0_wp*x3 + &
+ 3952831.0_wp ) / 10255478.0_wp
+ q(5,8) = ( 304819200.0_wp*x2 + 914457600.0_wp*x3 - &
+ 657798011.0_wp ) / 61532868.0_wp
+ q(6,8) = -5080320.0_wp*x2 / 5127739.0_wp
+ q(7,8) = -5080320.0_wp*x3 / 5127739.0_wp; q(8,8) = zero
+ q(9,8) = 4064256.0_wp/5127739.0_wp; q(10,8) = -1016064.0_wp/5127739.0_wp
+ q(11,8) = 193536.0_wp/5127739.0_wp; q(12,8) = -18144.0_wp/5127739.0_wp
+
+ first = .false.
+ end if
+
+ dd = zero
+ imin = 0
+ imax = -1
+
+ if ( bb(1) == 0 ) then
+ il = 1 + gsize
+ else
+ dd(1:8,1) = q(1:8,1)
+ imin(1) = 1; imax(1) = 8
+
+ dd(1:8,2) = q(1:8,2)
+ imin(2) = 1; imax(2) = 8
+
+ dd(1:8,3) = q(1:8,3)
+ imin(3) = 1; imax(3) = 8
+
+ dd(1:8,4) = q(1:8,4)
+ imin(4) = 1; imax(4) = 8
+
+ dd(1:9,5) = q(1:9,5)
+ imin(5) = 1; imax(5) = 9
+
+ dd(1:10,6) = q(1:10,6)
+ imin(6) = 1; imax(6) = 10
+
+ dd(1:11,7) = q(1:11,7)
+ imin(7) = 1; imax(7) = 11
+
+ dd(1:12,8) = q(1:12,8)
+ imin(8) = 1; imax(8) = 12
+
+ il = 9
+ end if
+ if ( bb(2) == 0 ) then
+ ir = nsize - gsize
+ else
+ dd(nsize-11:nsize,nsize-7) = -q(12:1:-1,8)
+ imin(nsize-7) = nsize-11; imax(nsize-7) = nsize
+
+ dd(nsize-10:nsize,nsize-6) = -q(11:1:-1,7)
+ imin(nsize-6) = nsize-10; imax(nsize-6) = nsize
+
+ dd(nsize-9:nsize,nsize-5) = -q(10:1:-1,6)
+ imin(nsize-5) = nsize-9; imax(nsize-5) = nsize
+
+ dd(nsize-8:nsize,nsize-4) = -q(9:1:-1,5)
+ imin(nsize-4) = nsize-8; imax(nsize-4) = nsize
+
+ dd(nsize-7:nsize,nsize-3) = -q(8:1:-1,4)
+ imin(nsize-3) = nsize-7; imax(nsize-3) = nsize
+
+ dd(nsize-7:nsize,nsize-2) = -q(8:1:-1,3)
+ imin(nsize-2) = nsize-7; imax(nsize-2) = nsize
+
+ dd(nsize-7:nsize,nsize-1) = -q(8:1:-1,2)
+ imin(nsize-1) = nsize-7; imax(nsize-1) = nsize
+
+ dd(nsize-7:nsize,nsize) = -q(8:1:-1,1)
+ imin(nsize) = nsize-7; imax(nsize) = nsize
+
+ ir = nsize - 8
+ end if
+ do i = il, ir
+ dd(i-4:i-1,i) = -a(4:1:-1); dd(i+1:i+4,i) = a(1:4)
+ imin(i) = i-4; imax(i) = i+4
+ end do
+end subroutine set_coeff_8_4
diff --git a/src/get_coeffs.c b/src/get_coeffs.c
new file mode 100644
index 0000000..0d8e365
--- /dev/null
+++ b/src/get_coeffs.c
@@ -0,0 +1,135 @@
+#include "cctk.h"
+#include "cctk_Parameters.h"
+#include "cctk_Arguments.h"
+#include "util_Table.h"
+#include <stdio.h>
+
+#include <assert.h>
+
+void DiffCoeff ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir,
+ const CCTK_INT nsize, CCTK_INT *imin, CCTK_INT *imax,
+ CCTK_REAL *q, const CCTK_INT table_handle )
+{
+ cGH const * restrict const cctkGH = cctkGH_;
+ DECLARE_CCTK_PARAMETERS
+ DECLARE_CCTK_ARGUMENTS
+
+ CCTK_INT ni, nj, nk, gsize, loc_order;
+ CCTK_INT ierr;
+ int nelements;
+ CCTK_INT lsh[3], bbox[6], bb[2], nghostzones[3];
+ void CCTK_FCALL CCTK_FNAME(set_coeff_2_1)(const CCTK_INT *nsize,
+ const CCTK_INT *loc_order,
+ const CCTK_INT *bb,
+ const CCTK_INT *gsize,
+ CCTK_INT *imin,
+ CCTK_INT *imax,
+ CCTK_REAL *q);
+ void CCTK_FCALL CCTK_FNAME(set_coeff_4_2)(const CCTK_INT *nsize,
+ const CCTK_INT *loc_order,
+ const CCTK_INT *bb,
+ const CCTK_INT *gsize,
+ CCTK_INT *imin,
+ CCTK_INT *imax,
+ CCTK_REAL *q);
+ void CCTK_FCALL CCTK_FNAME(set_coeff_6_3)(const CCTK_INT *nsize,
+ const CCTK_INT *loc_order,
+ const CCTK_INT *bb,
+ const CCTK_INT *gsize,
+ CCTK_INT *imin,
+ CCTK_INT *imax,
+ CCTK_REAL *q);
+ void CCTK_FCALL CCTK_FNAME(set_coeff_8_4)(const CCTK_INT *nsize,
+ const CCTK_INT *loc_order,
+ const CCTK_INT *bb,
+ const CCTK_INT *gsize,
+ CCTK_INT *imin,
+ CCTK_INT *imax,
+ CCTK_REAL *q);
+ void CCTK_FCALL CCTK_FNAME(set_coeff_4_3)(const CCTK_INT *nsize,
+ const CCTK_INT *loc_order,
+ const CCTK_INT *bb,
+ const CCTK_INT *gsize,
+ CCTK_INT *imin,
+ CCTK_INT *imax,
+ CCTK_REAL *q);
+ void CCTK_FCALL CCTK_FNAME(set_coeff_6_5)(const CCTK_INT *nsize,
+ const CCTK_INT *loc_order,
+ const CCTK_INT *bb,
+ const CCTK_INT *gsize,
+ CCTK_INT *imin,
+ CCTK_INT *imax,
+ CCTK_REAL *q);
+
+ ni = cctk_lsh[0]; nj = cctk_lsh[1]; nk = cctk_lsh[2];
+
+ if ( table_handle >=0 ) {
+ nelements = Util_TableGetInt ( table_handle, &loc_order, "order" );
+ assert(nelements>0);
+ } else {
+ loc_order = order;
+ }
+// printf("loc_order = %d\n",loc_order);
+// printf("dir = %d, nsize = %d\n",dir,nsize);
+
+
+ switch(dir) {
+ case 0: {
+ assert(nsize==ni);
+ bb[0] = cctk_bbox[0]; bb[1] = cctk_bbox[1];
+ gsize = cctk_nghostzones[0];
+ break;
+ }
+ case 1: {
+ assert(nsize==nj);
+ bb[0] = cctk_bbox[2]; bb[1] = cctk_bbox[3];
+ gsize = cctk_nghostzones[1];
+ break;
+ }
+ case 2: {
+ assert(nsize==nk);
+ bb[0] = cctk_bbox[4]; bb[1] = cctk_bbox[5];
+ gsize = cctk_nghostzones[2];
+ break;
+ }
+ default:
+ assert(0);
+ }
+// printf("bbox = [%d,%d], gsize = %d\n",bb[0],bb[1],gsize);
+
+ if ( CCTK_Equals(norm_type,"Diagonal") ) {
+ switch(loc_order) {
+ case 2: {
+ CCTK_FNAME(set_coeff_2_1)(&nsize,&loc_order,bb,&gsize,imin,imax,q);
+ break;
+ }
+ case 4: {
+ CCTK_FNAME(set_coeff_4_2)(&nsize,&loc_order,bb,&gsize,imin,imax,q);
+ break;
+ }
+ case 6: {
+ CCTK_FNAME(set_coeff_6_3)(&nsize,&loc_order,bb,&gsize,imin,imax,q);
+ break;
+ }
+ case 8: {
+ CCTK_FNAME(set_coeff_8_4)(&nsize,&loc_order,bb,&gsize,imin,imax,q);
+ break;
+ }
+ default:
+ assert(0);
+ }
+ } else {
+ switch(loc_order) {
+ case 4: {
+ CCTK_FNAME(set_coeff_4_3)(&nsize,&loc_order,bb,&gsize,imin,imax,q);
+ break;
+ }
+ case 6: {
+ CCTK_FNAME(set_coeff_6_5)(&nsize,&loc_order,bb,&gsize,imin,imax,q);
+ break;
+ }
+ default:
+ assert(0);
+ }
+ }
+}
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-12 22:01:52 +0000