diff options
Diffstat (limited to 'src/Derivatives_4_3_min_err_coeff.F90')
-rw-r--r-- | src/Derivatives_4_3_min_err_coeff.F90 | 43 |
1 files changed, 3 insertions, 40 deletions
diff --git a/src/Derivatives_4_3_min_err_coeff.F90 b/src/Derivatives_4_3_min_err_coeff.F90 index b42d857..7bc0a82 100644 --- a/src/Derivatives_4_3_min_err_coeff.F90 +++ b/src/Derivatives_4_3_min_err_coeff.F90 @@ -4,13 +4,13 @@ subroutine deriv_gf_4_3_opt ( var, ni, nj, nk, dir, bb, gsize, delta, dvar ) + use All_Coeffs_mod + implicit none DECLARE_CCTK_FUNCTIONS DECLARE_CCTK_PARAMETERS - CCTK_REAL, parameter :: zero = 0.0 - integer, parameter :: wp = kind(zero) CCTK_INT, intent(IN) :: ni, nj, nk CCTK_REAL, dimension(ni,nj,nk), intent(IN) :: var CCTK_INT, intent(IN) :: dir @@ -28,44 +28,7 @@ subroutine deriv_gf_4_3_opt ( var, ni, nj, nk, dir, bb, gsize, delta, dvar ) 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) = -2.09329763466349871588733_wp - q(2,1) = 4.0398572053206615302160_wp - q(3,1) = -3.0597858079809922953240_wp - q(4,1) = 1.37319053865399486354933_wp - q(5,1) = -0.25996430133016538255400_wp - q(6,1) = 0.0_wp - q(7,1) = 0.0_wp - q(1,2) = -0.31641585285940445272297_wp - q(2,2) = -0.53930788973980422327388_wp - q(3,2) = 0.98517732028644343383297_wp - q(4,2) = -0.05264665989297578146709_wp - q(5,2) = -0.113807251750624235013258_wp - q(6,2) = 0.039879767889849911803103_wp - q(7,2) = -0.0028794339334846531588787_wp - q(1,3) = 0.13026916185021164524452_wp - q(2,3) = -0.87966858995059249256890_wp - q(3,3) = 0.38609640961100070000134_wp - q(4,3) = 0.31358369072435588745988_wp - q(5,3) = 0.085318941913678384633511_wp - q(6,3) = -0.039046615792734640274641_wp - q(7,3) = 0.0034470016440805155042908_wp - q(1,4) = -0.01724512193824647912172_wp - q(2,4) = 0.16272288227127504381134_wp - q(3,4) = -0.81349810248648813029217_wp - q(4,4) = 0.13833269266479833215645_wp - q(5,4) = 0.59743854328548053399616_wp - q(6,4) = -0.066026434346299887619324_wp - q(7,4) = -0.0017244594505194129307249_wp - q(1,5) = -0.00883569468552192965061_wp - q(2,5) = 0.03056074759203203857284_wp - q(3,5) = 0.05021168274530854232278_wp - q(4,5) = -0.66307364652444929534068_wp - q(5,5) = 0.014878787464005191116088_wp - q(6,5) = 0.65882706381707471953820_wp - q(7,5) = -0.082568940408449266558615_wp - + call coeffs_1_4_3_opt ( a, q ) first = .false. end if |