diff options
Diffstat (limited to 'src/Derivatives_6_5_min_err_coeff.F90')
| -rw-r--r-- | src/Derivatives_6_5_min_err_coeff.F90 | 78 |
1 files changed, 3 insertions, 75 deletions
diff --git a/src/Derivatives_6_5_min_err_coeff.F90 b/src/Derivatives_6_5_min_err_coeff.F90 index f3735ab..7f9e209 100644 --- a/src/Derivatives_6_5_min_err_coeff.F90 +++ b/src/Derivatives_6_5_min_err_coeff.F90 @@ -4,13 +4,13 @@ subroutine deriv_gf_6_5_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,79 +28,7 @@ subroutine deriv_gf_6_5_opt ( var, ni, nj, nk, dir, bb, gsize, delta, dvar ) 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) = -2.465354921110524023660777656111276003457_wp - q(2,1) = 6.092129526663144141964665936667656020742_wp - q(3,1) = -7.730323816657860354911664841669140051855_wp - q(4,1) = 6.973765088877147139882219788892186735807_wp - q(5,1) = -3.980323816657860354911664841669140051855_wp - q(6,1) = 1.292129526663144141964665936667656020742_wp - q(7,1) = -0.1820215877771906903274443227779426701237_wp - q(8,1) = 0.0_wp - q(9,1) = 0.0_wp - q(10,1) = 0.0_wp - q(1,2) = -0.2234725650784319828746535134412736890421_wp - q(2,2) = -0.9329308121107134563129925525068570679651_wp - q(3,2) = 1.586820596545839371759081303802027231274_wp - q(4,2) = -0.3647002340377160216914505558624668821400_wp - q(5,2) = -0.2666957784872806143914117440166232718819_wp - q(6,2) = 0.3112949048634705032101261273629794071371_wp - q(7,2) = -0.1404504214762266650000768489896480092493_wp - q(8,2) = 0.03488568514730479833596013512958238764128_wp - q(9,2) = -0.004964021886392518344179263072091597647654_wp - q(10,2) = 0.0002126465201465853095969115943714918742904_wp - q(1,3) = 0.1582216737061633151406179477554921935333_wp - q(2,3) = -1.137049298003377811733609086574457439398_wp - q(3,3) = 1.212364522932578587741649981040340946798_wp - q(4,3) = -0.9562288729513894906148167047868730813830_wp - q(5,3) = 1.066548057336766350478498057851678826640_wp - q(6,3) = -0.3478788551267041838265477441805600110467_wp - q(7,3) = -0.03133923293520187620333693909408071632123_wp - q(8,3) = 0.04098845955755862691072597869183962277781_wp - q(9,3) = -0.005963188634687155197078928402509551508436_wp - q(10,3) = 0.0003367341182936373038974376991292099082999_wp - q(1,4) = 0.02915734641890708196910927068736798144670_wp - q(2,4) = -0.1169665089768926152768236581512624861308_wp - q(3,4) = -0.1112219092451476301503253995474190870412_wp - q(4,4) = -0.7924486261248032107393766820001361351677_wp - q(5,4) = 1.266650704820613624987450232358951199911_wp - q(6,4) = -0.2899273290506621673153239836530375587273_wp - q(7,4) = 0.002515684257201926199329020583484434062150_wp - q(8,4) = 0.01329713961871764653006682056620518602804_wp - q(9,4) = -0.001124464399630667352932212208930962568134_wp - q(10,4) = 0.00006796268169601114882659136477742818715059_wp - q(1,5) = -0.04582150000326981674750984653096293434777_wp - q(2,5) = 0.2240986548857151482718685516611524323427_wp - q(3,5) = -0.3246718493011818141660859125588209338018_wp - q(4,5) = -0.3929792921782506986152017485694441380503_wp - q(5,5) = 0.1166355818729375628072830916953646214341_wp - q(6,5) = 0.3449626905957060254933930895775644438105_wp - q(7,5) = 0.1430419813354607083034935179267283951745_wp - q(8,5) = -0.07764802499372607792980458731991885121073_wp - q(9,5) = 0.01332439335504217034559288889042994978834_wp - q(10,5) = -0.0009426355684332077630290447720929851395193_wp - q(1,6) = 0.003172814452954821196677290327889903944225_wp - q(2,6) = 0.00001061446045061551877105554145609103530766_wp - q(3,6) = -0.08747763580209736614983637747947172321794_wp - q(4,6) = 0.3975827322299876034907453299884380895682_wp - q(5,6) = -1.148835072393422871630425744497391344782_wp - q(6,6) = 0.3583006649535242306065761818925080902380_wp - q(7,6) = 0.5647665154270147564019144982190032455071_wp - q(8,6) = -0.09698196887272109736153117076061707705561_wp - q(9,6) = 0.008843905091972988427261446924164441884143_wp - q(10,6) = 0.0006174304523363194998474898440202828786385_wp - q(1,7) = -0.008639107540858839028043929986084287776394_wp - q(2,7) = 0.04722773954485212324714352753530343274219_wp - q(3,7) = -0.1008747537650261142294540111407681552350_wp - q(4,7) = 0.08043834953845218736895768965086958762389_wp - q(5,7) = 0.1295138674713300902982857323205417604553_wp - q(6,7) = -0.7909424166489541737614153656634872155367_wp - q(7,7) = 0.03807866847647628589685997987877954466259_wp - q(8,7) = 0.7367055699548196242687865288427927434250_wp - q(9,7) = -0.1480235854665196220062411065981933720158_wp - q(10,7) = 0.01651566843542843794512095516024596165494_wp - + call coeffs_1_6_5_opt ( a, q ) first = .false. end if |