diff options
Diffstat (limited to 'src/All_Coeffs_mod.F90')
| -rw-r--r-- | src/All_Coeffs_mod.F90 | 77 |
1 files changed, 77 insertions, 0 deletions
diff --git a/src/All_Coeffs_mod.F90 b/src/All_Coeffs_mod.F90 index 76dfcc9..69e58a9 100644 --- a/src/All_Coeffs_mod.F90 +++ b/src/All_Coeffs_mod.F90 @@ -497,6 +497,7 @@ module All_Coeffs_mod q2(7,1) = zero q2(8,1) = zero q2(9,1) = zero + q2(1,2) = -0.4261090852993166609415800254745304380753_wp q2(2,2) = -0.1198701406809289935902771997003246482977_wp q2(3,2) = 0.4312121105197350491616311601533345702699_wp @@ -506,6 +507,7 @@ module All_Coeffs_mod q2(7,2) = zero q2(8,2) = zero q2(9,2) = zero + q2(1,3) = -0.008331409049041400793408941823041323547141_wp q2(2,3) = -0.3179917159117626605641230407454828831892_wp q2(3,3) = -0.5046108447067502766506824050165990409443_wp @@ -515,6 +517,7 @@ module All_Coeffs_mod q2(7,3) = zero q2(8,3) = zero q2(9,3) = zero + q2(1,4) = 0.1280751399324364352167190636169803213668_wp q2(2,4) = -0.3704544846040473075364791079948980171701_wp q2(3,4) = -0.1053961094122492652071919749185808657060_wp @@ -524,6 +527,7 @@ module All_Coeffs_mod q2(7,4) = 0.02687068482925919014741556260496361261429_wp q2(8,4) = zero q2(9,4) = zero + q2(1,5) = -0.03614399304268576976452921364705641609825_wp q2(2,5) = 0.1234319378557362297361923714125920669065_wp q2(3,5) = -0.09358865464680026772155005029274896720786_wp @@ -533,6 +537,7 @@ module All_Coeffs_mod q2(7,5) = -0.2742160721086708137615843595277389869240_wp q2(8,5) = 0.03656214294782277516821124793703186492319_wp q2(9,5) = zero + q2(1,6) = -0.01141318406360863692889821914555232596651_wp q2(2,6) = 0.02049729840293952857599941220163960606616_wp q2(3,6) = 0.02756893105132916059596754912522489034598_wp @@ -545,6 +550,78 @@ module All_Coeffs_mod end subroutine coeffs_up_6_3_opt + subroutine coeffs_2_6_3 ( a, q ) + + CCTK_REAL, dimension(4), intent(OUT) :: a + CCTK_REAL, dimension(9,6), 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) = 2.788238454587637677973966346740908979901_wp + q(2,1) = -8.024525606271521723203165067037878232838_wp + q(3,1) = 8.215717879209710113072996800742423132342_wp + q(4,1) = -3.382384545876376779739663467409089799008_wp + q(5,1) = 0.2745256062715217232031650670378782328376_wp + q(6,1) = 0.1284282120790289886927003199257576867658_wp + q(7,1) = 0_wp + q(8,1) = 0_wp + q(9,1) = 0_wp + + q(1,2) = 1.053412969283276450511945392491467576792_wp + q(2,2) = -2.350398179749715585893060295790671217292_wp + q(3,2) = 1.867463026166097838452787258248009101251_wp + q(4,2) = -1.034129692832764505119453924914675767918_wp + q(5,2) = 0.6003981797497155858930602957906712172924_wp + q(6,2) = -0.1367463026166097838452787258248009101251_wp + q(7,2) = 0_wp + q(8,2) = 0_wp + q(9,2) = 0_wp + + q(1,3) = -0.6441780400836099840157383499323742776343_wp + q(2,3) = 4.137556867084716586745358416328538054838_wp + q(3,3) = -8.108447067502766506824050165990409443010_wp + q(4,3) = 6.941780400836099840157383499323742776343_wp + q(5,3) = -2.887556867084716586745358416328538054838_wp + q(6,3) = 0.5608447067502766506824050165990409443010_wp + q(7,3) = 0_wp + q(8,3) = 0_wp + q(9,3) = 0_wp + + q(1,4) = 0.2133575915904708589911053057162405921503_wp + q(2,4) = -1.159078186228774025004665049449524164956_wp + q(3,4) = 3.511693723953473906823412328170678609193_wp + q(4,4) = -4.723144865335572557069104932512284630217_wp + q(5,4) = 2.489690240716551595446911737264415002799_wp + q(6,4) = -0.3414753996392361759034645767245132798408_wp + q(7,4) = 0.008956894943086396715805187534987870871431_wp + q(8,4) = 0_wp + q(9,4) = 0_wp + + q(1,5) = -0.1790783293131903008759680081249206550717_wp + q(2,5) = 0.9156510515847827006897719097795268926410_wp + q(3,5) = -1.987601032541999915365409843002835258770_wp + q(4,5) = 3.387647581566586263806017519360162498413_wp + q(5,5) = -3.719859506580339384706529558630612331260_wp + q(6,5) = 1.735582497566755532986331513689644957894_wp + q(7,5) = -0.1645296432652024882569506157166433921544_wp + q(8,5) = 0.01218738098260759172273708264567728830773_wp + q(9,5) = 0_wp + + q(1,6) = 0.04002001476374207590389869333272451922711_wp + q(2,6) = -0.1875223548929628699496967344733377472356_wp + q(3,6) = 0.3471267779274445788908929019885390744504_wp + q(4,6) = -0.4177910702190969764769449708149737829426_wp + q(5,6) = 1.560601736642238000654475164189554271972_wp + q(6,6) = -2.684870208442729618045250108445012670944_wp + q(7,6) = 1.479418278121504075249423529143170247255_wp + q(8,6) = -0.1479418278121504075249423529143170247255_wp + q(9,6) = 0.01095865391201114129814387799365311294263_wp + + end subroutine coeffs_2_6_3 + subroutine coeffs_1_8_4 ( a, q ) CCTK_REAL, dimension(4), intent(OUT) :: a |