diff options
Diffstat (limited to 'src/Derivatives_6_3_min_err_coeff.F90')
| -rw-r--r-- | src/Derivatives_6_3_min_err_coeff.F90 | 193 |
1 files changed, 6 insertions, 187 deletions
diff --git a/src/Derivatives_6_3_min_err_coeff.F90 b/src/Derivatives_6_3_min_err_coeff.F90 index b2d00e7..bf8b8f2 100644 --- a/src/Derivatives_6_3_min_err_coeff.F90 +++ b/src/Derivatives_6_3_min_err_coeff.F90 @@ -4,13 +4,13 @@ subroutine deriv_gf_6_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,63 +28,7 @@ subroutine deriv_gf_6_3_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) = -1.582533518939116418785258993332844897062_wp - q(2,1) = 2.033426786468126253898161347360808173712_wp - q(3,1) = -0.1417052898146741610733887894481170575600_wp - q(4,1) = -0.4501096599735708523162117824920488989702_wp - q(5,1) = 0.1042956382142412661862395105494407610836_wp - q(6,1) = 0.03662604404499391209045870736276191879693_wp - q(7,1) = 0.0_wp - q(8,1) = 0.0_wp - q(9,1) = 0.0_wp - q(1,2) = -0.4620701275035953590186631853846278325646_wp - q(2,2) = 0.0_wp - q(3,2) = 0.2873679417026202568532985205129449923126_wp - q(4,2) = 0.2585974499280928196267362923074433487080_wp - q(5,2) = -0.06894808744606961472005221923058251153103_wp - q(6,2) = -0.01494717668104810274131940820517799692506_wp - q(7,2) = 0.0_wp - q(8,2) = 0.0_wp - q(9,2) = 0.0_wp - q(1,3) = 0.07134398748360337973038301686379010397038_wp - q(2,3) = -0.6366933020423417826592908754928085932593_wp - q(3,3) = 0.0_wp - q(4,3) = 0.6067199374180168986519150843189505198519_wp - q(5,3) = -0.02338660408468356531858175098561718651857_wp - q(6,3) = -0.01798401877459493040442547470431484404443_wp - q(7,3) = 0.0_wp - q(8,3) = 0.0_wp - q(9,3) = 0.0_wp - q(1,4) = 0.1146397975178068401430112823144985150596_wp - q(2,4) = -0.2898424301162697370942324201800071793273_wp - q(3,4) = -0.3069262456316931913128086944558079603132_wp - q(4,4) = 0.0_wp - q(5,4) = 0.5203848121857539166740071338174418292578_wp - q(6,4) = -0.05169127637022742348368508279860701098408_wp - q(7,4) = 0.01343534241462959507370778130248180630715_wp - q(8,4) = 0.0_wp - q(9,4) = 0.0_wp - q(1,5) = -0.03614399304268576976452921364705641609825_wp - q(2,5) = 0.1051508663818248421520867474440761344449_wp - q(3,5) = 0.01609777419666805778308369351834662756172_wp - q(4,5) = -0.7080721616106272031118456849378369336023_wp - q(5,5) = 0.0_wp - q(6,5) = 0.7692160858661111736140494493705980473867_wp - q(7,5) = -0.1645296432652024882569506157166433921544_wp - q(8,5) = 0.01828107147391138758410562396851593246160_wp - q(9,5) = 0.0_wp - q(1,6) = -0.01141318406360863692889821914555232596651_wp - q(2,6) = 0.02049729840293952857599941220163960606616_wp - q(3,6) = 0.01113095018331244864875173213474522093204_wp - q(4,6) = 0.06324365883611076515355091406993789453750_wp - q(5,6) = -0.6916640154753724474963890679085181638850_wp - q(6,6) = 0.0_wp - q(7,6) = 0.7397091390607520376247117645715851236273_wp - q(8,6) = -0.1479418278121504075249423529143170247255_wp - q(9,6) = 0.01643798086801671194721581699047966941394_wp - + call coeffs_1_6_3_opt ( a, q ) first = .false. end if @@ -295,13 +239,13 @@ end subroutine deriv_gf_6_3_opt subroutine up_deriv_gf_6_3_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, 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 @@ -320,132 +264,7 @@ subroutine up_deriv_gf_6_3_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva logical, save :: first = .true. if ( first ) then - a1(-3) = -0.03333333333333333333333333333333333333333_wp - a1(-2) = 0.2500000000000000000000000000000000000000_wp - a1(-1) = -1.000000000000000000000000000000000000000_wp - a1(0) = 0.3333333333333333333333333333333333333333_wp - a1(1) = 0.5000000000000000000000000000000000000000_wp - a1(2) = -0.05000000000000000000000000000000000000000_wp - a1(3) = 0.0_wp - - q1(1,1) = -1.529782401641145871492417026888416733827_wp - q1(2,1) = 1.875173434574214612019635448027523684006_wp - q1(3,1) = 0.01654806207923748080513710988516743214616_wp - q1(4,1) = -0.5028607772715413996090537489364770622056_wp - q1(5,1) = 0.1042956382142412661862395105494407610836_wp - q1(6,1) = 0.03662604404499391209045870736276191879693_wp - q1(7,1) = 0.0_wp - q1(8,1) = 0.0_wp - q1(9,1) = 0.0_wp - q1(1,2) = -0.4980311697078740570957463452947252270539_wp - q1(2,2) = 0.1198701406809289935902771997003246482977_wp - q1(3,2) = 0.1435237728855054645449658808725554143554_wp - q1(4,2) = 0.3305195343366502157809026121276381376866_wp - q1(5,2) = -0.08093510151416251407907993920061497636080_wp - q1(6,2) = -0.01494717668104810274131940820517799692506_wp - q1(7,2) = 0.0_wp - q1(8,2) = 0.0_wp - q1(9,2) = 0.0_wp - q1(1,3) = 0.1510193840162481602541749755506215314879_wp - q1(2,3) = -0.9553948881729209047544587102401343033294_wp - q1(3,3) = 0.5046108447067502766506824050165990409443_wp - q1(4,3) = 0.2083429547547929960329552908847933822643_wp - q1(5,3) = 0.1359641889806059957290021663880456685165_wp - q1(6,3) = -0.04454248428547652391235612759992531988360_wp - q1(7,3) = 0.0_wp - q1(8,3) = 0.0_wp - q1(9,3) = 0.0_wp - q1(1,4) = 0.1012044551031772450693035010120167087525_wp - q1(2,4) = -0.2092303756284921666519857323651163414844_wp - q1(3,4) = -0.5084563818511371174184254139930350549204_wp - q1(4,4) = 0.2687068482925919014741556260496361261429_wp - q1(5,4) = 0.3188546759663099905683904142802147346506_wp - q1(6,4) = 0.02892077811755014695856160501628382685880_wp - q1(7,4) = 0.0_wp - q1(8,4) = 0.0_wp - q1(9,4) = 0.0_wp - q1(1,5) = -0.03614399304268576976452921364705641609825_wp - q1(2,5) = 0.08686979490791345456798112347556020198327_wp - q1(3,5) = 0.1257842030401363832877174373294422223313_wp - q1(4,5) = -0.9822882337192980168734300444655759205262_wp - q1(5,5) = 0.3656214294782277516821124793703186492319_wp - q1(6,5) = 0.4950000137574403598524650898428590604628_wp - q1(7,5) = -0.05484321442173416275231687190554779738479_wp - q1(8,5) = 0.0_wp - q1(9,5) = 0.0_wp - q1(1,6) = -0.01141318406360863692889821914555232596651_wp - q1(2,6) = 0.02049729840293952857599941220163960606616_wp - q1(3,6) = -0.005307030684704263298464084855734448481896_wp - q1(4,6) = 0.1618715440442110368368458160128159110211_wp - q1(5,6) = -0.9382337284956231267046263227657132050941_wp - q1(6,6) = 0.3287596173603342389443163398095933882788_wp - q1(7,6) = 0.4931394260405013584164745097143900824182_wp - q1(8,6) = -0.04931394260405013584164745097143900824182_wp - q1(9,6) = 0.0_wp - - a2(-3) = 0.0_wp - a2(-2) = 0.05000000000000000000000000000000000000000_wp - a2(-1) = -0.5000000000000000000000000000000000000000_wp - a2(0) = -0.3333333333333333333333333333333333333333_wp - a2(1) = 1.000000000000000000000000000000000000000_wp - a2(2) = -0.2500000000000000000000000000000000000000_wp - a2(3) = 0.03333333333333333333333333333333333333333_wp - - q2(1,1) = -1.635284636237086966078100959777273060297_wp - q2(2,1) = 2.191680138362037895776687246694092663418_wp - q2(3,1) = -0.2999586417085858029519146887814015472662_wp - q2(4,1) = -0.3973585426756003050233698160476207357348_wp - q2(5,1) = 0.1042956382142412661862395105494407610836_wp - q2(6,1) = 0.03662604404499391209045870736276191879693_wp - q2(7,1) = 0.0_wp - q2(8,1) = 0.0_wp - q2(9,1) = 0.0_wp - q2(1,2) = -0.4261090852993166609415800254745304380753_wp - q2(2,2) = -0.1198701406809289935902771997003246482977_wp - q2(3,2) = 0.4312121105197350491616311601533345702699_wp - q2(4,2) = 0.1866753655195354234725699724872485597294_wp - q2(5,2) = -0.05696107337797671536102449926055004670126_wp - q2(6,2) = -0.01494717668104810274131940820517799692506_wp - q2(7,2) = 0.0_wp - q2(8,2) = 0.0_wp - q2(9,2) = 0.0_wp - q2(1,3) = -0.008331409049041400793408941823041323547141_wp - q2(2,3) = -0.3179917159117626605641230407454828831892_wp - q2(3,3) = -0.5046108447067502766506824050165990409443_wp - q2(4,3) = 1.005096920081240801270874877753107657440_wp - q2(5,3) = -0.1827373971499731263661656683592800415536_wp - q2(6,3) = 0.008574446736286663103505178191295631794744_wp - q2(7,3) = 0.0_wp - q2(8,3) = 0.0_wp - q2(9,3) = 0.0_wp - q2(1,4) = 0.1280751399324364352167190636169803213668_wp - q2(2,4) = -0.3704544846040473075364791079948980171701_wp - q2(3,4) = -0.1053961094122492652071919749185808657060_wp - q2(4,4) = -0.2687068482925919014741556260496361261429_wp - q2(5,4) = 0.7219149484051978427796238533546689238650_wp - q2(6,4) = -0.1323033308580049939259317706134978488270_wp - q2(7,4) = 0.02687068482925919014741556260496361261429_wp - q2(8,4) = 0.0_wp - q2(9,4) = 0.0_wp - q2(1,5) = -0.03614399304268576976452921364705641609825_wp - q2(2,5) = 0.1234319378557362297361923714125920669065_wp - q2(3,5) = -0.09358865464680026772155005029274896720786_wp - q2(4,5) = -0.4338560895019563893502613254100979466783_wp - q2(5,5) = -0.3656214294782277516821124793703186492319_wp - q2(6,5) = 1.043432157974781987375633808898337034311_wp - q2(7,5) = -0.2742160721086708137615843595277389869240_wp - q2(8,5) = 0.03656214294782277516821124793703186492319_wp - q2(9,5) = 0.0_wp - q2(1,6) = -0.01141318406360863692889821914555232596651_wp - q2(2,6) = 0.02049729840293952857599941220163960606616_wp - q2(3,6) = 0.02756893105132916059596754912522489034598_wp - q2(4,6) = -0.03538422637198950652974398787294012194614_wp - q2(5,6) = -0.4450943024551217682881518130513231226759_wp - q2(6,6) = -0.3287596173603342389443163398095933882788_wp - q2(7,6) = 0.9862788520810027168329490194287801648364_wp - q2(8,6) = -0.2465697130202506792082372548571950412091_wp - q2(9,6) = 0.03287596173603342389443163398095933882788_wp - + call coeffs_up_6_3_opt ( a1, q1, a2, q2 ) first = .false. end if |