diff options
Diffstat (limited to 'src/All_Coeffs_mod.F90')
| -rw-r--r-- | src/All_Coeffs_mod.F90 | 148 |
1 files changed, 129 insertions, 19 deletions
diff --git a/src/All_Coeffs_mod.F90 b/src/All_Coeffs_mod.F90 index 69e58a9..1b773ad 100644 --- a/src/All_Coeffs_mod.F90 +++ b/src/All_Coeffs_mod.F90 @@ -202,25 +202,25 @@ module All_Coeffs_mod q(2,1) = -5.000000000000000000000000000000000000000_wp q(3,1) = 4.000000000000000000000000000000000000000_wp q(4,1) = -1.000000000000000000000000000000000000000_wp - q(5,1) = 0_wp - q(6,1) = 0_wp + q(5,1) = zero + q(6,1) = zero q(1,2) = 1.000000000000000000000000000000000000000_wp q(2,2) = -2.000000000000000000000000000000000000000_wp q(3,2) = 1.000000000000000000000000000000000000000_wp - q(4,2) = 0_wp - q(5,2) = 0_wp - q(6,2) = 0_wp + q(4,2) = zero + q(5,2) = zero + q(6,2) = zero q(1,3) = -0.09302325581395348837209302325581395348837_wp q(2,3) = 1.372093023255813953488372093023255813953_wp q(3,3) = -2.558139534883720930232558139534883720930_wp q(4,3) = 1.372093023255813953488372093023255813953_wp q(5,3) = -0.09302325581395348837209302325581395348837_wp - q(6,3) = 0_wp + q(6,3) = zero q(1,4) = -0.02040816326530612244897959183673469387755_wp - q(2,4) = 0_wp + q(2,4) = zero q(3,4) = 1.204081632653061224489795918367346938776_wp q(4,4) = -2.408163265306122448979591836734693877551_wp q(5,4) = 1.306122448979591836734693877551020408163_wp @@ -566,9 +566,9 @@ module All_Coeffs_mod 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(7,1) = zero + q(8,1) = zero + q(9,1) = zero q(1,2) = 1.053412969283276450511945392491467576792_wp q(2,2) = -2.350398179749715585893060295790671217292_wp @@ -576,9 +576,9 @@ module All_Coeffs_mod 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(7,2) = zero + q(8,2) = zero + q(9,2) = zero q(1,3) = -0.6441780400836099840157383499323742776343_wp q(2,3) = 4.137556867084716586745358416328538054838_wp @@ -586,9 +586,9 @@ module All_Coeffs_mod 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(7,3) = zero + q(8,3) = zero + q(9,3) = zero q(1,4) = 0.2133575915904708589911053057162405921503_wp q(2,4) = -1.159078186228774025004665049449524164956_wp @@ -597,8 +597,8 @@ module All_Coeffs_mod 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(8,4) = zero + q(9,4) = zero q(1,5) = -0.1790783293131903008759680081249206550717_wp q(2,5) = 0.9156510515847827006897719097795268926410_wp @@ -608,7 +608,7 @@ module All_Coeffs_mod q(6,5) = 1.735582497566755532986331513689644957894_wp q(7,5) = -0.1645296432652024882569506157166433921544_wp q(8,5) = 0.01218738098260759172273708264567728830773_wp - q(9,5) = 0_wp + q(9,5) = zero q(1,6) = 0.04002001476374207590389869333272451922711_wp q(2,6) = -0.1875223548929628699496967344733377472356_wp @@ -1140,6 +1140,116 @@ module All_Coeffs_mod end subroutine coeffs_up_8_4_opt + subroutine coeffs_2_8_4 ( a, q ) + + CCTK_REAL, dimension(5), intent(OUT) :: a + CCTK_REAL, dimension(12,8), intent(OUT) :: q + + a(1) = -2.847222222222222222222222222222222222222_wp + a(2) = 1.600000000000000000000000000000000000000_wp + a(3) = -0.2000000000000000000000000000000000000000_wp + a(4) = 0.02539682539682539682539682539682539682540_wp + a(5) = -0.001785714285714285714285714285714285714286_wp + + q(1,1) = 3.700125961665716499242827594045243158419_wp + q(2,1) = -12.53408910332763232879029889760479228385_wp + q(3,1) = 17.08522275831908082197574724401198070962_wp + q(4,1) = -12.00251923331432998485655188090486316838_wp + q(5,1) = 4.335222758319080821975747244011980709621_wp + q(6,1) = -0.5340891033276323287902988976047922838482_wp + q(7,1) = -0.04987403833428350075717240595475684158085_wp + q(8,1) = zero + q(9,1) = zero + q(10,1) = zero + q(11,1) = zero + q(12,1) = zero + q(1,2) = 0.8545561806165568782871327076593746414458_wp + q(2,2) = -1.385651335219520350522366987732872958967_wp + q(3,2) = 0.03489488503609432577108173221703695010854_wp + q(4,2) = 0.6174959481995095555971180534964638500300_wp + q(5,2) = -0.01537226034806520970159454957687817246322_wp + q(6,2) = -0.1687175231686941483830240392722911697288_wp + q(7,2) = 0.07110835640429802975122382498579196986742_wp + q(8,2) = -0.008314251520179080799570741776625110292494_wp + q(9,2) = zero + q(10,2) = zero + q(11,2) = zero + q(12,2) = zero + q(1,3) = 0.1486597947505276786174351637623847280273_wp + q(2,3) = 0.2067945264072191961776578930368936098873_wp + q(3,3) = -0.6126228482043964280320857972314409478969_wp + q(4,3) = 0.6747701953118921125519965255870017768863_wp + q(5,3) = -1.911835643616456845050516722454935307278_wp + q(6,3) = 2.589340655152612946529422151593721302602_wp + q(7,3) = -1.360526641378450595342844756571493806946_wp + q(8,3) = 0.2654199615770519345489355422778686447177_wp + q(9,3) = zero + q(10,3) = zero + q(11,3) = zero + q(12,3) = zero + q(1,4) = -0.1146188247411133220554202699779770824365_wp + q(2,4) = 0.5239525746096007180691923905405359669894_wp + q(3,4) = 0.09661320523910878808200465962763626256375_wp + q(4,4) = -1.414029112733921772841533040350356265715_wp + q(5,4) = 1.222591772291551121101947203543540982713_wp + q(6,4) = -0.6020259924428416149507501533753687531597_wp + q(7,4) = 0.3679434182813619635245551059859820833417_wp + q(8,4) = -0.08042704050374588092999589599399319429637_wp + q(9,4) = zero + q(10,4) = zero + q(11,4) = zero + q(12,4) = zero + q(1,5) = 0.06885441184905793064871518875563556079618_wp + q(2,5) = -0.05682897197659991639833322507928914093129_wp + q(3,5) = -1.192628972259900927264800832893135534458_wp + q(4,5) = 5.326668514854093024959085487041627095051_wp + q(5,5) = -8.658036133958342620625148610263805114859_wp + q(6,5) = 6.351205788124211735819368535961188726091_wp + q(7,5) = -2.217166245530019638125083881812697165497_wp + q(8,5) = 0.3822584308534509596846452244576658437999_wp + q(9,5) = -0.004326821955950548698447886167190269992355_wp + q(10,5) = zero + q(11,5) = zero + q(12,5) = zero + q(1,6) = 0.03324412543986457361212616687247594270027_wp + q(2,6) = -0.2013444667830157913459721318051963127270_wp + q(3,6) = 0.5214244967518298047750888317723818148279_wp + q(4,6) = -0.8467142789066131514432621699356515897162_wp + q(5,6) = 2.050234909410534381962583240886249492721_wp + q(6,6) = -2.978298213712682019535273900331397359519_wp + q(7,6) = 1.578055555463864252208238842945777112098_wp + q(8,6) = -0.1750701723207555800161607707603235528915_wp + q(9,6) = 0.01986478753019001522837715937418159597272_wp + q(10,6) = -0.001396742873216485445745269018497143466832_wp + q(11,6) = zero + q(12,6) = zero + q(1,7) = -0.01592923134237101475574629012747169731725_wp + q(2,7) = 0.1175046087844397827007651199218350175991_wp + q(3,7) = -0.3793709248668048299249134817817685155210_wp + q(4,7) = 0.7165680376452736733332948640656278517079_wp + q(5,7) = -0.9910611454636313458152579274571837912909_wp + q(6,7) = 2.185128799543527842172503734448943570849_wp + q(7,7) = -3.127418850946242350478043794675628822387_wp + q(8,7) = 1.685621415770389829447415862701279666019_wp + q(9,7) = -0.2166153552278720352907291696202456084323_wp + q(10,7) = 0.02750671177496787749723545011050737884854_wp + q(11,7) = -0.001934065671677428886524367585895050075288_wp + q(12,7) = zero + q(1,8) = zero + q(2,8) = -0.01256791782038527350173477151453027289777_wp + q(3,8) = 0.06770121895926724008958693099893390610910_wp + q(4,8) = -0.1432796115427160189070360250265966875352_wp + q(5,8) = 0.1563022312812120974195915804261720302791_wp + q(6,8) = -0.2217546957854982665994400260882731886207_wp + q(7,8) = 1.541934599402382617061776818198730931660_wp + q(8,8) = -2.798782027001839321767451883871457955840_wp + q(9,8) = 1.585203927110954750232022339670564355947_wp + q(10,8) = -0.1981504908888693437790027924588205444934_wp + q(11,8) = 0.02516196709699928174971464031223118025313_wp + q(12,8) = -0.001769200811507761998026810646953754861548_wp + + end subroutine coeffs_2_8_4 + subroutine coeffs_1_4_3 ( a, q ) CCTK_REAL, dimension(2), intent(OUT) :: a |