Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Add parse(BigFloat, ...) back #453

Merged
merged 2 commits into from
Jul 24, 2018
Merged

Add parse(BigFloat, ...) back #453

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
4bf1435
Add `parse(BigFloat, ...)` back
YingboMa Jul 23, 2018
33650e7
Fix Nordsieck test
YingboMa Jul 23, 2018
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
1,294 changes: 647 additions & 647 deletions src/tableaus/feagin_tableaus.jl
View file
Open in desktop

Large diffs are not rendered by default.

298 changes: 149 additions & 149 deletions src/tableaus/high_order_rk_tableaus.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,83 +768,83 @@ function DP8ConstantCache(T::Type,T2::Type)
c4 = T2((6-sqrt(6))/10 * c6)
c3 = T2(2//3 * c4)
c2 = T2(2//3 * c3)
b1 = T( 5.42937341165687622380535766363e-2)
b6 = T( 4.45031289275240888144113950566)
b7 = T( 1.89151789931450038304281599044)
b8 = T(-5.8012039600105847814672114227)
b9 = T( 3.1116436695781989440891606237e-1)
b10 = T(-1.52160949662516078556178806805e-1)
b11 = T( 2.01365400804030348374776537501e-1)
b12 = T( 4.47106157277725905176885569043e-2)
# bhh1 = T(0.244094488188976377952755905512)
# bhh2 = T(0.733846688281611857341361741547)
# bhh3 = T(0.220588235294117647058823529412e-01)
btilde1 = T(-1.898007540724076157147023288757e-1)
btilde6 = T( 4.45031289275240888144113950566)
btilde7 = T( 1.89151789931450038304281599044)
btilde8 = T(-5.8012039600105847814672114227)
btilde9 = T(-4.22682321323791962932445679177e-1)
btilde10 = T(-1.52160949662516078556178806805e-1)
btilde11 = T( 2.01365400804030348374776537501e-1)
btilde12 = T(2.26517921983608258118062039631e-2)
er1 = T( 0.1312004499419488073250102996e-01)
er6 = T(-0.1225156446376204440720569753e+01)
er7 = T(-0.4957589496572501915214079952)
er8 = T( 0.1664377182454986536961530415e+01)
er9 = T(-0.3503288487499736816886487290)
er10 = T( 0.3341791187130174790297318841)
er11 = T( 0.8192320648511571246570742613e-01)
er12 = T(-0.2235530786388629525884427845e-01)
a0201 = T( 5.26001519587677318785587544488e-2)
a0301 = T( 1.97250569845378994544595329183e-2)
a0302 = T( 5.91751709536136983633785987549e-2)
a0401 = T( 2.95875854768068491816892993775e-2)
a0403 = T( 8.87627564304205475450678981324e-2)
a0501 = T( 2.41365134159266685502369798665e-1)
a0503 = T(-8.84549479328286085344864962717e-1)
a0504 = T( 9.24834003261792003115737966543e-1)
a0601 = T( 3.7037037037037037037037037037e-2)
a0604 = T( 1.70828608729473871279604482173e-1)
a0605 = T( 1.25467687566822425016691814123e-1)
a0701 = T( 3.7109375e-2)
a0704 = T( 1.70252211019544039314978060272e-1)
a0705 = T( 6.02165389804559606850219397283e-2)
a0706 = T(-1.7578125e-2)
a0801 = T( 3.70920001185047927108779319836e-2)
a0804 = T( 1.70383925712239993810214054705e-1)
a0805 = T( 1.07262030446373284651809199168e-1)
a0806 = T(-1.53194377486244017527936158236e-2)
a0807 = T( 8.27378916381402288758473766002e-3)
a0901 = T( 6.24110958716075717114429577812e-1)
a0904 = T(-3.36089262944694129406857109825)
a0905 = T(-8.68219346841726006818189891453e-1)
a0906 = T( 2.75920996994467083049415600797e1)
a0907 = T( 2.01540675504778934086186788979e1)
a0908 = T(-4.34898841810699588477366255144e1)
a1001 = T( 4.77662536438264365890433908527e-1)
a1004 = T(-2.48811461997166764192642586468e0)
a1005 = T(-5.90290826836842996371446475743e-1)
a1006 = T( 2.12300514481811942347288949897e1)
a1007 = T( 1.52792336328824235832596922938e1)
a1008 = T(-3.32882109689848629194453265587e1)
a1009 = T(-2.03312017085086261358222928593e-2)
a1101 = T(-9.3714243008598732571704021658e-1)
a1104 = T( 5.18637242884406370830023853209e0)
a1105 = T( 1.09143734899672957818500254654e0)
a1106 = T(-8.14978701074692612513997267357e0)
a1107 = T(-1.85200656599969598641566180701e1)
a1108 = T( 2.27394870993505042818970056734e1)
a1109 = T( 2.49360555267965238987089396762e0)
a1110 = T(-3.0467644718982195003823669022e0)
a1201 = T( 2.27331014751653820792359768449e0)
a1204 = T( -1.05344954667372501984066689879e1)
a1205 = T( -2.00087205822486249909675718444e0)
a1206 = T( -1.79589318631187989172765950534e1)
a1207 = T( 2.79488845294199600508499808837e1)
a1208 = T( -2.85899827713502369474065508674e0)
a1209 = T( -8.87285693353062954433549289258e0)
a1210 = T( 1.23605671757943030647266201528e1)
a1211 = T( 6.43392746015763530355970484046e-1)
b1 = T(parse(BigFloat," 5.42937341165687622380535766363e-2"))
b6 = T(parse(BigFloat," 4.45031289275240888144113950566"))
b7 = T(parse(BigFloat," 1.89151789931450038304281599044"))
b8 = T(parse(BigFloat,"-5.8012039600105847814672114227"))
b9 = T(parse(BigFloat," 3.1116436695781989440891606237e-1"))
b10 = T(parse(BigFloat,"-1.52160949662516078556178806805e-1"))
b11 = T(parse(BigFloat," 2.01365400804030348374776537501e-1"))
b12 = T(parse(BigFloat," 4.47106157277725905176885569043e-2"))
# bhh1 = T(parse(BigFloat,"0.244094488188976377952755905512"))
# bhh2 = T(parse(BigFloat,"0.733846688281611857341361741547"))
# bhh3 = T(parse(BigFloat,"0.220588235294117647058823529412e-01"))
btilde1 = T(parse(BigFloat,"-1.898007540724076157147023288757e-1"))
btilde6 = T(parse(BigFloat," 4.45031289275240888144113950566"))
btilde7 = T(parse(BigFloat," 1.89151789931450038304281599044"))
btilde8 = T(parse(BigFloat,"-5.8012039600105847814672114227"))
btilde9 = T(parse(BigFloat,"-4.22682321323791962932445679177e-1"))
btilde10 = T(parse(BigFloat,"-1.52160949662516078556178806805e-1"))
btilde11 = T(parse(BigFloat," 2.01365400804030348374776537501e-1"))
btilde12 = T(parse(BigFloat,"2.26517921983608258118062039631e-2"))
er1 = T(parse(BigFloat," 0.1312004499419488073250102996e-01"))
er6 = T(parse(BigFloat,"-0.1225156446376204440720569753e+01"))
er7 = T(parse(BigFloat,"-0.4957589496572501915214079952"))
er8 = T(parse(BigFloat," 0.1664377182454986536961530415e+01"))
er9 = T(parse(BigFloat,"-0.3503288487499736816886487290"))
er10 = T(parse(BigFloat," 0.3341791187130174790297318841"))
er11 = T(parse(BigFloat," 0.8192320648511571246570742613e-01"))
er12 = T(parse(BigFloat,"-0.2235530786388629525884427845e-01"))
a0201 = T(parse(BigFloat," 5.26001519587677318785587544488e-2"))
a0301 = T(parse(BigFloat," 1.97250569845378994544595329183e-2"))
a0302 = T(parse(BigFloat," 5.91751709536136983633785987549e-2"))
a0401 = T(parse(BigFloat," 2.95875854768068491816892993775e-2"))
a0403 = T(parse(BigFloat," 8.87627564304205475450678981324e-2"))
a0501 = T(parse(BigFloat," 2.41365134159266685502369798665e-1"))
a0503 = T(parse(BigFloat,"-8.84549479328286085344864962717e-1"))
a0504 = T(parse(BigFloat," 9.24834003261792003115737966543e-1"))
a0601 = T(parse(BigFloat," 3.7037037037037037037037037037e-2"))
a0604 = T(parse(BigFloat," 1.70828608729473871279604482173e-1"))
a0605 = T(parse(BigFloat," 1.25467687566822425016691814123e-1"))
a0701 = T(parse(BigFloat," 3.7109375e-2"))
a0704 = T(parse(BigFloat," 1.70252211019544039314978060272e-1"))
a0705 = T(parse(BigFloat," 6.02165389804559606850219397283e-2"))
a0706 = T(parse(BigFloat,"-1.7578125e-2"))
a0801 = T(parse(BigFloat," 3.70920001185047927108779319836e-2"))
a0804 = T(parse(BigFloat," 1.70383925712239993810214054705e-1"))
a0805 = T(parse(BigFloat," 1.07262030446373284651809199168e-1"))
a0806 = T(parse(BigFloat,"-1.53194377486244017527936158236e-2"))
a0807 = T(parse(BigFloat," 8.27378916381402288758473766002e-3"))
a0901 = T(parse(BigFloat," 6.24110958716075717114429577812e-1"))
a0904 = T(parse(BigFloat,"-3.36089262944694129406857109825"))
a0905 = T(parse(BigFloat,"-8.68219346841726006818189891453e-1"))
a0906 = T(parse(BigFloat," 2.75920996994467083049415600797e1"))
a0907 = T(parse(BigFloat," 2.01540675504778934086186788979e1"))
a0908 = T(parse(BigFloat,"-4.34898841810699588477366255144e1"))
a1001 = T(parse(BigFloat," 4.77662536438264365890433908527e-1"))
a1004 = T(parse(BigFloat,"-2.48811461997166764192642586468e0"))
a1005 = T(parse(BigFloat,"-5.90290826836842996371446475743e-1"))
a1006 = T(parse(BigFloat," 2.12300514481811942347288949897e1"))
a1007 = T(parse(BigFloat," 1.52792336328824235832596922938e1"))
a1008 = T(parse(BigFloat,"-3.32882109689848629194453265587e1"))
a1009 = T(parse(BigFloat,"-2.03312017085086261358222928593e-2"))
a1101 = T(parse(BigFloat,"-9.3714243008598732571704021658e-1"))
a1104 = T(parse(BigFloat," 5.18637242884406370830023853209e0"))
a1105 = T(parse(BigFloat," 1.09143734899672957818500254654e0"))
a1106 = T(parse(BigFloat,"-8.14978701074692612513997267357e0"))
a1107 = T(parse(BigFloat,"-1.85200656599969598641566180701e1"))
a1108 = T(parse(BigFloat," 2.27394870993505042818970056734e1"))
a1109 = T(parse(BigFloat," 2.49360555267965238987089396762e0"))
a1110 = T(parse(BigFloat,"-3.0467644718982195003823669022e0"))
a1201 = T(parse(BigFloat," 2.27331014751653820792359768449e0"))
a1204 = T(parse(BigFloat," -1.05344954667372501984066689879e1"))
a1205 = T(parse(BigFloat," -2.00087205822486249909675718444e0"))
a1206 = T(parse(BigFloat," -1.79589318631187989172765950534e1"))
a1207 = T(parse(BigFloat," 2.79488845294199600508499808837e1"))
a1208 = T(parse(BigFloat," -2.85899827713502369474065508674e0"))
a1209 = T(parse(BigFloat," -8.87285693353062954433549289258e0"))
a1210 = T(parse(BigFloat," 1.23605671757943030647266201528e1"))
a1211 = T(parse(BigFloat," 6.43392746015763530355970484046e-1"))
c14,c15,c16,a1401,a1407,a1408,a1409,a1410,a1411,a1412,a1413,a1501,a1506,a1507,a1508,a1511,a1512,a1513,a1514,a1601,a1606,a1607,a1608,a1609,a1613,a1614,a1615 = DP8Interp(T,T2)
d401,d406,d407,d408,d409,d410,d411,d412,d413,d414,d415,d416,d501,d506,d507,d508,d509,d510,d511,d512,d513,d514,d515,d516,d601,d606,d607,d608,d609,d610,d611,d612,d613,d614,d615,d616,d701,d706,d707,d708,d709,d710,d711,d712,d713,d714,d715,d716 = DP8Interp_polyweights(T)
DP8ConstantCache(c7,c8,c9,c10,c11,c6,c5,c4,c3,c2,b1,b6,b7,b8,b9,b10,b11,b12,btilde1,btilde6,btilde7,btilde8,btilde9,btilde10,btilde11,btilde12,er1,er6,er7,er8,er9,er10,er11,er12,a0201,a0301,a0302,a0401,a0403,a0501,a0503,a0504,a0601,a0604,a0605,a0701,a0704,a0705,a0706,a0801,a0804,a0805,a0806,a0807,a0901,a0904,a0905,a0906,a0907,a0908,a1001,a1004,a1005,a1006,a1007,a1008,a1009,a1101,a1104,a1105,a1106,a1107,a1108,a1109,a1110,a1201,a1204,a1205,a1206,a1207,a1208,a1209,a1210,a1211,c14,c15,c16,a1401,a1407,a1408,a1409,a1410,a1411,a1412,a1413,a1501,a1506,a1507,a1508,a1511,a1512,a1513,a1514,a1601,a1606,a1607,a1608,a1609,a1613,a1614,a1615,d401,d406,d407,d408,d409,d410,d411,d412,d413,d414,d415,d416,d501,d506,d507,d508,d509,d510,d511,d512,d513,d514,d515,d516,d601,d606,d607,d608,d609,d610,d611,d612,d613,d614,d615,d616,d701,d706,d707,d708,d709,d710,d711,d712,d713,d714,d715,d716)
Expand Down Expand Up @@ -886,30 +886,30 @@ function DP8Interp(T::Type,T2::Type)
c14 = T2(1//10)
c15 = T2(2//10)
c16 = T2(7//9)
a1401 = T( 5.61675022830479523392909219681e-2)
a1407 = T( 2.53500210216624811088794765333e-1)
a1408 = T(-2.46239037470802489917441475441e-1)
a1409 = T(-1.24191423263816360469010140626e-1)
a1410 = T( 1.5329179827876569731206322685e-1)
a1411 = T( 8.20105229563468988491666602057e-3)
a1412 = T( 7.56789766054569976138603589584e-3)
a1413 = T(-8.298e-3)
a1501 = T( 3.18346481635021405060768473261e-2)
a1506 = T( 2.83009096723667755288322961402e-2)
a1507 = T( 5.35419883074385676223797384372e-2)
a1508 = T(-5.49237485713909884646569340306e-2)
a1511 = T(-1.08347328697249322858509316994e-4)
a1512 = T( 3.82571090835658412954920192323e-4)
a1513 = T(-3.40465008687404560802977114492e-4)
a1514 = T( 1.41312443674632500278074618366e-1)
a1601 = T(-4.28896301583791923408573538692e-1)
a1606 = T(-4.69762141536116384314449447206e0)
a1607 = T( 7.68342119606259904184240953878e0)
a1608 = T( 4.06898981839711007970213554331e0)
a1609 = T( 3.56727187455281109270669543021e-1)
a1613 = T(-1.39902416515901462129418009734e-3)
a1614 = T( 2.9475147891527723389556272149e0)
a1615 = T(-9.15095847217987001081870187138e0)
a1401 = T(parse(BigFloat," 5.61675022830479523392909219681e-2"))
a1407 = T(parse(BigFloat," 2.53500210216624811088794765333e-1"))
a1408 = T(parse(BigFloat,"-2.46239037470802489917441475441e-1"))
a1409 = T(parse(BigFloat,"-1.24191423263816360469010140626e-1"))
a1410 = T(parse(BigFloat," 1.5329179827876569731206322685e-1"))
a1411 = T(parse(BigFloat," 8.20105229563468988491666602057e-3"))
a1412 = T(parse(BigFloat," 7.56789766054569976138603589584e-3"))
a1413 = T(parse(BigFloat,"-8.298e-3"))
a1501 = T(parse(BigFloat," 3.18346481635021405060768473261e-2"))
a1506 = T(parse(BigFloat," 2.83009096723667755288322961402e-2"))
a1507 = T(parse(BigFloat," 5.35419883074385676223797384372e-2"))
a1508 = T(parse(BigFloat,"-5.49237485713909884646569340306e-2"))
a1511 = T(parse(BigFloat,"-1.08347328697249322858509316994e-4"))
a1512 = T(parse(BigFloat," 3.82571090835658412954920192323e-4"))
a1513 = T(parse(BigFloat,"-3.40465008687404560802977114492e-4"))
a1514 = T(parse(BigFloat," 1.41312443674632500278074618366e-1"))
a1601 = T(parse(BigFloat,"-4.28896301583791923408573538692e-1"))
a1606 = T(parse(BigFloat,"-4.69762141536116384314449447206e0"))
a1607 = T(parse(BigFloat," 7.68342119606259904184240953878e0"))
a1608 = T(parse(BigFloat," 4.06898981839711007970213554331e0"))
a1609 = T(parse(BigFloat," 3.56727187455281109270669543021e-1"))
a1613 = T(parse(BigFloat,"-1.39902416515901462129418009734e-3"))
a1614 = T(parse(BigFloat," 2.9475147891527723389556272149e0"))
a1615 = T(parse(BigFloat,"-9.15095847217987001081870187138e0"))

return c14,c15,c16,a1401,a1407,a1408,a1409,a1410,a1411,a1412,a1413,a1501,a1506,a1507,a1508,a1511,a1512,a1513,a1514,a1601,a1606,a1607,a1608,a1609,a1613,a1614,a1615
end
Expand Down Expand Up @@ -969,54 +969,54 @@ end


function DP8Interp_polyweights(T::Type)
d401 = T(-0.84289382761090128651353491142e+01)
d406 = T( 0.56671495351937776962531783590e+00)
d407 = T(-0.30689499459498916912797304727e+01)
d408 = T( 0.23846676565120698287728149680e+01)
d409 = T( 0.21170345824450282767155149946e+01)
d410 = T(-0.87139158377797299206789907490e+00)
d411 = T( 0.22404374302607882758541771650e+01)
d412 = T( 0.63157877876946881815570249290e+00)
d413 = T(-0.88990336451333310820698117400e-01)
d414 = T( 0.18148505520854727256656404962e+02)
d415 = T(-0.91946323924783554000451984436e+01)
d416 = T(-0.44360363875948939664310572000e+01)
d501 = T( 0.10427508642579134603413151009e+02)
d506 = T( 0.24228349177525818288430175319e+03)
d507 = T( 0.16520045171727028198505394887e+03)
d508 = T(-0.37454675472269020279518312152e+03)
d509 = T(-0.22113666853125306036270938578e+02)
d510 = T( 0.77334326684722638389603898808e+01)
d511 = T(-0.30674084731089398182061213626e+02)
d512 = T(-0.93321305264302278729567221706e+01)
d513 = T( 0.15697238121770843886131091075e+02)
d514 = T(-0.31139403219565177677282850411e+02)
d515 = T(-0.93529243588444783865713862664e+01)
d516 = T( 0.35816841486394083752465898540e+02)
d601 = T( 0.19985053242002433820987653617e+02)
d606 = T(-0.38703730874935176555105901742e+03)
d607 = T(-0.18917813819516756882830838328e+03)
d608 = T( 0.52780815920542364900561016686e+03)
d609 = T(-0.11573902539959630126141871134e+02)
d610 = T( 0.68812326946963000169666922661e+01)
d611 = T(-0.10006050966910838403183860980e+01)
d612 = T( 0.77771377980534432092869265740e+00)
d613 = T(-0.27782057523535084065932004339e+01)
d614 = T(-0.60196695231264120758267380846e+02)
d615 = T( 0.84320405506677161018159903784e+02)
d616 = T( 0.11992291136182789328035130030e+02)
d701 = T(-0.25693933462703749003312586129e+02)
d706 = T(-0.15418974869023643374053993627e+03)
d707 = T(-0.23152937917604549567536039109e+03)
d708 = T( 0.35763911791061412378285349910e+03)
d709 = T( 0.93405324183624310003907691704e+02)
d710 = T(-0.37458323136451633156875139351e+02)
d711 = T( 0.10409964950896230045147246184e+03)
d712 = T( 0.29840293426660503123344363579e+02)
d713 = T(-0.43533456590011143754432175058e+02)
d714 = T( 0.96324553959188282948394950600e+02)
d715 = T(-0.39177261675615439165231486172e+02)
d716 = T(-0.14972683625798562581422125276e+03)
d401 = T(parse(BigFloat,"-0.84289382761090128651353491142e+01"))
d406 = T(parse(BigFloat," 0.56671495351937776962531783590e+00"))
d407 = T(parse(BigFloat,"-0.30689499459498916912797304727e+01"))
d408 = T(parse(BigFloat," 0.23846676565120698287728149680e+01"))
d409 = T(parse(BigFloat," 0.21170345824450282767155149946e+01"))
d410 = T(parse(BigFloat,"-0.87139158377797299206789907490e+00"))
d411 = T(parse(BigFloat," 0.22404374302607882758541771650e+01"))
d412 = T(parse(BigFloat," 0.63157877876946881815570249290e+00"))
d413 = T(parse(BigFloat,"-0.88990336451333310820698117400e-01"))
d414 = T(parse(BigFloat," 0.18148505520854727256656404962e+02"))
d415 = T(parse(BigFloat,"-0.91946323924783554000451984436e+01"))
d416 = T(parse(BigFloat,"-0.44360363875948939664310572000e+01"))
d501 = T(parse(BigFloat," 0.10427508642579134603413151009e+02"))
d506 = T(parse(BigFloat," 0.24228349177525818288430175319e+03"))
d507 = T(parse(BigFloat," 0.16520045171727028198505394887e+03"))
d508 = T(parse(BigFloat,"-0.37454675472269020279518312152e+03"))
d509 = T(parse(BigFloat,"-0.22113666853125306036270938578e+02"))
d510 = T(parse(BigFloat," 0.77334326684722638389603898808e+01"))
d511 = T(parse(BigFloat,"-0.30674084731089398182061213626e+02"))
d512 = T(parse(BigFloat,"-0.93321305264302278729567221706e+01"))
d513 = T(parse(BigFloat," 0.15697238121770843886131091075e+02"))
d514 = T(parse(BigFloat,"-0.31139403219565177677282850411e+02"))
d515 = T(parse(BigFloat,"-0.93529243588444783865713862664e+01"))
d516 = T(parse(BigFloat," 0.35816841486394083752465898540e+02"))
d601 = T(parse(BigFloat," 0.19985053242002433820987653617e+02"))
d606 = T(parse(BigFloat,"-0.38703730874935176555105901742e+03"))
d607 = T(parse(BigFloat,"-0.18917813819516756882830838328e+03"))
d608 = T(parse(BigFloat," 0.52780815920542364900561016686e+03"))
d609 = T(parse(BigFloat,"-0.11573902539959630126141871134e+02"))
d610 = T(parse(BigFloat," 0.68812326946963000169666922661e+01"))
d611 = T(parse(BigFloat,"-0.10006050966910838403183860980e+01"))
d612 = T(parse(BigFloat," 0.77771377980534432092869265740e+00"))
d613 = T(parse(BigFloat,"-0.27782057523535084065932004339e+01"))
d614 = T(parse(BigFloat,"-0.60196695231264120758267380846e+02"))
d615 = T(parse(BigFloat," 0.84320405506677161018159903784e+02"))
d616 = T(parse(BigFloat," 0.11992291136182789328035130030e+02"))
d701 = T(parse(BigFloat,"-0.25693933462703749003312586129e+02"))
d706 = T(parse(BigFloat,"-0.15418974869023643374053993627e+03"))
d707 = T(parse(BigFloat,"-0.23152937917604549567536039109e+03"))
d708 = T(parse(BigFloat," 0.35763911791061412378285349910e+03"))
d709 = T(parse(BigFloat," 0.93405324183624310003907691704e+02"))
d710 = T(parse(BigFloat,"-0.37458323136451633156875139351e+02"))
d711 = T(parse(BigFloat," 0.10409964950896230045147246184e+03"))
d712 = T(parse(BigFloat," 0.29840293426660503123344363579e+02"))
d713 = T(parse(BigFloat,"-0.43533456590011143754432175058e+02"))
d714 = T(parse(BigFloat," 0.96324553959188282948394950600e+02"))
d715 = T(parse(BigFloat,"-0.39177261675615439165231486172e+02"))
d716 = T(parse(BigFloat,"-0.14972683625798562581422125276e+03"))

return d401,d406,d407,d408,d409,d410,d411,d412,d413,d414,d415,d416,d501,d506,d507,d508,d509,d510,d511,d512,d513,d514,d515,d516,d601,d606,d607,d608,d609,d610,d611,d612,d613,d614,d615,d616,d701,d706,d707,d708,d709,d710,d711,d712,d713,d714,d715,d716
end