Skip to content

Commit

Permalink
Revert "modernize tableaus"
Browse files Browse the repository at this point in the history
This reverts commit b13d4d5.
  • Loading branch information
ChrisRackauckas committed Jul 30, 2018
1 parent 3e13840 commit 696b119
Show file tree
Hide file tree
Showing 8 changed files with 1,572 additions and 1,572 deletions.
1,280 changes: 640 additions & 640 deletions src/tableaus/feagin_tableaus.jl

Large diffs are not rendered by default.

308 changes: 154 additions & 154 deletions src/tableaus/high_order_rk_tableaus.jl

Large diffs are not rendered by default.

126 changes: 63 additions & 63 deletions src/tableaus/low_order_rk_tableaus.jl
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,7 @@ constructBogakiShampine3()
Constructs the tableau object for the Bogakai-Shampine Order 2/3 method.
"""
function BS3ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function BS3ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
a21 = T(0.5)
a32 = T(0.75)
a41 = T(0.2222222222222222)
Expand Down Expand Up @@ -80,7 +80,7 @@ struct OwrenZen3ConstantCache{T,T2} <: OrdinaryDiffEqConstantCache
r32::T
end

function OwrenZen3ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function OwrenZen3ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
a21 = T(0.5217391304347826)
a31 = T(-0.18133333333333335)
a32 = T(0.9813333333333333)
Expand All @@ -104,7 +104,7 @@ function OwrenZen3ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T
r13,r12,r23,r22,r33,r32)
end

function OwrenZen3ConstantCache(::Type{T},::Type{T2}) where {T,T2}
Base.@pure function OwrenZen3ConstantCache(::Type{T},::Type{T2}) where {T,T2}
a21 = T(12//23)
a31 = T(-68//375)
a32 = T(368//375)
Expand Down Expand Up @@ -168,7 +168,7 @@ struct OwrenZen4ConstantCache{T,T2} <: OrdinaryDiffEqConstantCache
r62::T
end

function OwrenZen4ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function OwrenZen4ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
a21 = T(0.16666666666666666)
a31 = T(0.03214024835646457)
a32 = T(0.26515704894083275)
Expand Down Expand Up @@ -215,7 +215,7 @@ function OwrenZen4ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T
r54,r53,r52,r64,r63,r62)
end

function OwrenZen4ConstantCache(::Type{T},::Type{T2}) where {T,T2}
Base.@pure function OwrenZen4ConstantCache(::Type{T},::Type{T2}) where {T,T2}
a21 = T(1//6)
a31 = T(44//1369)
a32 = T(363//1369)
Expand Down Expand Up @@ -331,7 +331,7 @@ struct OwrenZen5ConstantCache{T,T2} <: OrdinaryDiffEqConstantCache
r82::T
end

function OwrenZen5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function OwrenZen5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
a21 = T(0.16666666666666666)
a31 = T(0.0625)
a32 = T(0.1875)
Expand Down Expand Up @@ -412,7 +412,7 @@ function OwrenZen5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T
r62,r75,r74,r73,r72,r85,r84,r83,r82)
end

function OwrenZen5ConstantCache(::Type{T},::Type{T2}) where {T,T2}
Base.@pure function OwrenZen5ConstantCache(::Type{T},::Type{T2}) where {T,T2}
a21 = T(1//6)
a31 = T(1//16)
a32 = T(3//16)
Expand Down Expand Up @@ -553,7 +553,7 @@ struct Tsit5ConstantCache{T,T2} <: OrdinaryDiffEqConstantCache
r74::T
end

function Tsit5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function Tsit5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
c1 = T2(0.161)
c2 = T2(0.327)
c3 = T2(0.9)
Expand Down Expand Up @@ -604,43 +604,43 @@ function Tsit5ConstantCache(T::Type,T2::Type)
c1 = T2(161//1000)
c2 = T2(327//1000)
c3 = T2(9//10)
c4 = T2(big".9800255409045096857298102862870245954942137979563024768854764293221195950761080302604")
c4 = T2(parse(BigFloat,".9800255409045096857298102862870245954942137979563024768854764293221195950761080302604"))
c5 = T2(1)
c6 = T2(1)
a21= T(161//1000)
a31= T(big"-.8480655492356988544426874250230774675121177393430391537369234245294192976164141156943e-2")
a32= T(big".3354806554923569885444268742502307746751211773934303915373692342452941929761641411569")
a41= T(big"2.897153057105493432130432594192938764924887287701866490314866693455023795137503079289")
a42= T(big"-6.359448489975074843148159912383825625952700647415626703305928850207288721235210244366")
a43= T(big"4.362295432869581411017727318190886861027813359713760212991062156752264926097707165077")
a51= T(big"5.325864828439256604428877920840511317836476253097040101202360397727981648835607691791")
a52= T(big"-11.74888356406282787774717033978577296188744178259862899288666928009020615663593781589")
a53= T(big"7.495539342889836208304604784564358155658679161518186721010132816213648793440552049753")
a54= T(big"-.9249506636175524925650207933207191611349983406029535244034750452930469056411389539635e-1")
a61= T(big"5.861455442946420028659251486982647890394337666164814434818157239052507339770711679748")
a62= T(big"-12.92096931784710929170611868178335939541780751955743459166312250439928519268343184452")
a63= T(big"8.159367898576158643180400794539253485181918321135053305748355423955009222648673734986")
a64= T(big"-.7158497328140099722453054252582973869127213147363544882721139659546372402303777878835e-1")
a65= T(big"-.2826905039406838290900305721271224146717633626879770007617876201276764571291579142206e-1")
a71= T(big".9646076681806522951816731316512876333711995238157997181903319145764851595234062815396e-1")
a31= T(parse(BigFloat,"-.8480655492356988544426874250230774675121177393430391537369234245294192976164141156943e-2"))
a32= T(parse(BigFloat,".3354806554923569885444268742502307746751211773934303915373692342452941929761641411569"))
a41= T(parse(BigFloat,"2.897153057105493432130432594192938764924887287701866490314866693455023795137503079289"))
a42= T(parse(BigFloat,"-6.359448489975074843148159912383825625952700647415626703305928850207288721235210244366"))
a43= T(parse(BigFloat,"4.362295432869581411017727318190886861027813359713760212991062156752264926097707165077"))
a51= T(parse(BigFloat,"5.325864828439256604428877920840511317836476253097040101202360397727981648835607691791"))
a52= T(parse(BigFloat,"-11.74888356406282787774717033978577296188744178259862899288666928009020615663593781589"))
a53= T(parse(BigFloat,"7.495539342889836208304604784564358155658679161518186721010132816213648793440552049753"))
a54= T(parse(BigFloat,"-.9249506636175524925650207933207191611349983406029535244034750452930469056411389539635e-1"))
a61= T(parse(BigFloat,"5.861455442946420028659251486982647890394337666164814434818157239052507339770711679748"))
a62= T(parse(BigFloat,"-12.92096931784710929170611868178335939541780751955743459166312250439928519268343184452"))
a63= T(parse(BigFloat,"8.159367898576158643180400794539253485181918321135053305748355423955009222648673734986"))
a64= T(parse(BigFloat,"-.7158497328140099722453054252582973869127213147363544882721139659546372402303777878835e-1"))
a65= T(parse(BigFloat,"-.2826905039406838290900305721271224146717633626879770007617876201276764571291579142206e-1"))
a71= T(parse(BigFloat,".9646076681806522951816731316512876333711995238157997181903319145764851595234062815396e-1"))
a72= T(1//100)
a73= T(big".4798896504144995747752495322905965199130404621990332488332634944254542060153074523509")
a74= T(big"1.379008574103741893192274821856872770756462643091360525934940067397245698027561293331")
a75= T(big"-3.290069515436080679901047585711363850115683290894936158531296799594813811049925401677")
a76= T(big"2.324710524099773982415355918398765796109060233222962411944060046314465391054716027841")
# b1 = T(big".9468075576583945807478876255758922856117527357724631226139574065785592789071067303271e-1")
# b2 = T(big".9183565540343253096776363936645313759813746240984095238905939532922955247253608687270e-2")
# b3 = T(big".4877705284247615707855642599631228241516691959761363774365216240304071651579571959813")
# b4 = T(big"1.234297566930478985655109673884237654035539930748192848315425833500484878378061439761")
# b5 = T(big"-2.707712349983525454881109975059321670689605166938197378763992255714444407154902012702")
# b6 = T(big"1.866628418170587035753719399566211498666255505244122593996591602841258328965767580089")
a73= T(parse(BigFloat,".4798896504144995747752495322905965199130404621990332488332634944254542060153074523509"))
a74= T(parse(BigFloat,"1.379008574103741893192274821856872770756462643091360525934940067397245698027561293331"))
a75= T(parse(BigFloat,"-3.290069515436080679901047585711363850115683290894936158531296799594813811049925401677"))
a76= T(parse(BigFloat,"2.324710524099773982415355918398765796109060233222962411944060046314465391054716027841"))
# b1 = T(parse(BigFloat,".9468075576583945807478876255758922856117527357724631226139574065785592789071067303271e-1"))
# b2 = T(parse(BigFloat,".9183565540343253096776363936645313759813746240984095238905939532922955247253608687270e-2"))
# b3 = T(parse(BigFloat,".4877705284247615707855642599631228241516691959761363774365216240304071651579571959813"))
# b4 = T(parse(BigFloat,"1.234297566930478985655109673884237654035539930748192848315425833500484878378061439761"))
# b5 = T(parse(BigFloat,"-2.707712349983525454881109975059321670689605166938197378763992255714444407154902012702"))
# b6 = T(parse(BigFloat,"1.866628418170587035753719399566211498666255505244122593996591602841258328965767580089"))
# b7 = T(1//66)
btilde1 = T(big"-1.780011052225771443378550607539534775944678804333659557637450799792588061629796e-03")
btilde2 = T(big"-8.164344596567469032236360633546862401862537590159047610940604670770447527463931e-04")
btilde3 = T(big"7.880878010261996010314727672526304238628733777103128603258129604952959142646516e-03")
btilde4 = T(big"-1.44711007173262907537165147972635116720922712343167677619514233896760819649515e-01")
btilde5 = T(big"5.823571654525552250199376106520421794260781239567387797673045438803694038950012e-01")
btilde6 = T(big"-4.580821059291869466616365188325542974428047279788398179474684434732070620889539e-01")
btilde1 = T(parse(BigFloat,"-1.780011052225771443378550607539534775944678804333659557637450799792588061629796e-03"))
btilde2 = T(parse(BigFloat,"-8.164344596567469032236360633546862401862537590159047610940604670770447527463931e-04"))
btilde3 = T(parse(BigFloat,"7.880878010261996010314727672526304238628733777103128603258129604952959142646516e-03"))
btilde4 = T(parse(BigFloat,"-1.44711007173262907537165147972635116720922712343167677619514233896760819649515e-01"))
btilde5 = T(parse(BigFloat,"5.823571654525552250199376106520421794260781239567387797673045438803694038950012e-01"))
btilde6 = T(parse(BigFloat,"-4.580821059291869466616365188325542974428047279788398179474684434732070620889539e-01"))
btilde7 = T(1//66)

r11,r12,r13,r14,r22,r23,r24,r32,r33,r34,r42,r43,r44,r52,r53,r54,r62,r63,r64,r72,r73,r74 = Tsit5Interp(T)
Expand Down Expand Up @@ -703,30 +703,30 @@ simplifying assumption
Ch. Tsitouras
"""
function Tsit5Interp(T::Type)
r11 = T(big"0.999999999999999974283372471559910888475488471328")
r12 = T(big"-2.763706197274825911336735930481400260916070804192")
r13 = T(big"2.91325546182191274375068099306808")
r11 = T(parse(BigFloat,"0.999999999999999974283372471559910888475488471328"))
r12 = T(parse(BigFloat,"-2.763706197274825911336735930481400260916070804192"))
r13 = T(parse(BigFloat,"2.91325546182191274375068099306808"))
r14 = T(-1.0530884977290216)

r22 = T(big"0.13169999999999999727")
r23 = T(big"-0.22339999999999999818")
r22 = T(parse(BigFloat,"0.13169999999999999727"))
r23 = T(parse(BigFloat,"-0.22339999999999999818"))
r24 = T(0.1017)

r32 = T(big"3.93029623689475152850687446709813398")
r33 = T(big"-5.94103387213150473470249202589458001")
r34 = T(big"2.490627285651252793")
r32 = T(parse(BigFloat,"3.93029623689475152850687446709813398"))
r33 = T(parse(BigFloat,"-5.94103387213150473470249202589458001"))
r34 = T(parse(BigFloat,"2.490627285651252793"))

r42 = T(big"-12.411077166933676983734381540685453484102414134010752")
r43 = T(big"30.3381886302823215981729903691836576")
r44 = T(big"-16.54810288924490272")
r42 = T(parse(BigFloat,"-12.411077166933676983734381540685453484102414134010752"))
r43 = T(parse(BigFloat,"30.3381886302823215981729903691836576"))
r44 = T(parse(BigFloat,"-16.54810288924490272"))

r52 = T(big"37.50931341651103919496903965334519631242339792120440212")
r53 = T(big"-88.1789048947664011014276693541209817")
r54 = T(big"47.37952196281928122")
r52 = T(parse(BigFloat,"37.50931341651103919496903965334519631242339792120440212"))
r53 = T(parse(BigFloat,"-88.1789048947664011014276693541209817"))
r54 = T(parse(BigFloat,"47.37952196281928122"))

r62 = T(big"-27.896526289197287805948263144598643896")
r63 = T(big"65.09189467479367152629021928716553658")
r64 = T(big"-34.87065786149660974")
r62 = T(parse(BigFloat,"-27.896526289197287805948263144598643896"))
r63 = T(parse(BigFloat,"65.09189467479367152629021928716553658"))
r64 = T(parse(BigFloat,"-34.87065786149660974"))

r72 = T(1.5)
r73 = T(-4)
Expand Down Expand Up @@ -866,7 +866,7 @@ end
An Efficient Runge-Kutta (4,5) Pair by P.Bogacki and L.F.Shampine
Computers and Mathematics with Applications, Vol. 32, No. 6, 1996, pages 15 to 28
"""
function BS5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function BS5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
c1 =T2(0.16666666666666666)
c2 =T2(0.2222222222222222)
c3 =T2(0.42857142857142855)
Expand Down Expand Up @@ -999,7 +999,7 @@ Used in the lazy construction of the dense output
k9, k10, k11 are not computed until called in the dense routine
"""
function BS5Interp(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function BS5Interp(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}

c6 = T2(0.5)
c7 = T2(0.8333333333333334)
Expand Down Expand Up @@ -1269,7 +1269,7 @@ struct DP5ConstantCache{T,T2} <: OrdinaryDiffEqConstantCache
d7::T
end

function DP5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function DP5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
a21 = T(0.2)
a31 = T(0.075)
a32 = T(0.225)
Expand Down Expand Up @@ -1312,7 +1312,7 @@ function DP5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:Com
DP5ConstantCache(a21,a31,a32,a41,a42,a43,a51,a52,a53,a54,a61,a62,a63,a64,a65,a71,a73,a74,a75,a76,btilde1,btilde3,btilde4,btilde5,btilde6,btilde7,c1,c2,c3,c4,c5,c6,d1,d3,d4,d5,d6,d7)
end

function DP5_dense_ds(::Type{T}) where T<:CompiledFloats
Base.@pure function DP5_dense_ds(::Type{T}) where T<:CompiledFloats
d1 = T(-1.1270175653862835)
d3 = T(2.675424484351598)
d4 = T(-5.685526961588504)
Expand Down Expand Up @@ -1420,7 +1420,7 @@ struct Anas5ConstantCache{T,T2} <: OrdinaryDiffEqConstantCache
b6::T
end

function Anas5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function Anas5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
a21 = T(0.1)
a31 = T(-0.2222222222222222)
a32 = T(0.5555555555555556)
Expand Down Expand Up @@ -1449,7 +1449,7 @@ function Anas5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:C
Anas5ConstantCache(a21,a31,a32,a41,a42,a43,a51,a52,a53,a54,a61,a62,a63,a64,a65,c2,c3,c4,c5,c6,b1,b3,b4,b5,b6)
end

function Anas5ConstantCache(::Type{T},::Type{T2}) where {T,T2}
Base.@pure function Anas5ConstantCache(::Type{T},::Type{T2}) where {T,T2}
a21 = T(1//10)
a31 = T(-2//9)
a32 = T(5//9)
Expand Down
2 changes: 1 addition & 1 deletion src/tableaus/rkc_tableaus.jl
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
function ROCK2ConstantCache(::Type{T}, ::Type{T2}, zprev) where {T, T2}
Base.@pure function ROCK2ConstantCache(::Type{T}, ::Type{T2}, zprev) where {T, T2}
ms = SVector{46, Int}(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,
22,24,26,28,30,33,36,39,43,47,51,56,61,66,72,78,85,93,
102,112,123,135,148,163,180,198,)
Expand Down
14 changes: 7 additions & 7 deletions src/tableaus/rkn_tableaus.jl
Original file line number Diff line number Diff line change
Expand Up @@ -9,7 +9,7 @@ struct IRKN3ConstantCache{T,T2} <: OrdinaryDiffEqConstantCache
bbar2::T
end

function IRKN3ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function IRKN3ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
bconst1 = T(1.5)
bconst2 = T(-0.5)
c1 = T2(0.5)
Expand Down Expand Up @@ -49,7 +49,7 @@ struct IRKN4ConstantCache{T,T2} <: OrdinaryDiffEqConstantCache
bbar3::T
end

function IRKN4ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function IRKN4ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
bconst1 = T(1.5)
bconst2 = T(-0.5)
c1 = T2(0.25)
Expand Down Expand Up @@ -198,7 +198,7 @@ function ERKN4ConstantCache(T::Type,T2::Type)
ERKN4ConstantCache(c1, c2, c3, a21, a31, a32, a41, a42, a43, b1, b2, b3, b4, bp1, bp2, bp3, bp4, btilde1, btilde2, btilde3, btilde4, bptilde1, bptilde2, bptilde3, bptilde4)
end

function ERKN4ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function ERKN4ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
ERKN4ConstantCache(
T2(0.25),
T2(0.7),
Expand Down Expand Up @@ -284,7 +284,7 @@ function ERKN5ConstantCache(T::Type,T2::Type)
ERKN5ConstantCache(c1, c2, c3, a21, a31, a32, a41, a42, a43, b1, b2, b3, b4, bp1, bp2, bp3, bp4, btilde1, btilde2, btilde3, btilde4)
end

function ERKN5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function ERKN5ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
ERKN5ConstantCache(
T2(0.5),
T2(0.2714285714285714),
Expand Down Expand Up @@ -398,7 +398,7 @@ struct DPRKN6ConstantCache{T,T2} <: OrdinaryDiffEqConstantCache
rp61::T
end

function DPRKN6ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function DPRKN6ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
c1 = T2(0.12929590313670442)
c2 = T2(0.25859180627340883)
c3 = T2(0.67029708261548)
Expand Down Expand Up @@ -754,7 +754,7 @@ function DPRKN8ConstantCache(T::Type,T2::Type)
DPRKN8ConstantCache(c1, c2, c3, c4, c5, c6, c7, c8, a21, a31, a32, a41, a42, a43, a51, a52, a53, a54, a61, a62, a63, a64, a65, a71, a72, a73, a74, a75, a76, a81, a82, a83, a84, a85, a86, a87, a91, a93, a94, a95, a96, a97, b1, b3, b4, b5, b6, b7, bp1, bp3, bp4, bp5, bp6, bp7, bp8, btilde1, btilde3, btilde4, btilde5, btilde6, btilde7, bptilde1, bptilde3, bptilde4, bptilde5, bptilde6, bptilde7, bptilde8, bptilde9)
end

function DPRKN8ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function DPRKN8ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
DPRKN8ConstantCache(
T2(0.05),
T2(0.1),
Expand Down Expand Up @@ -1292,7 +1292,7 @@ function DPRKN12ConstantCache(T::Type,T2::Type)
DPRKN12ConstantCache(c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11, c12, c13, c14, c15, c16, a21, a31, a32, a41, a42, a43, a51, a53, a54, a61, a63, a64, a65, a71, a73, a74, a75, a76, a81, a84, a85, a86, a87, a91, a93, a94, a95, a96, a97, a98, a101, a103, a104, a105, a106, a107, a108, a109, a111, a113, a114, a115, a116, a117, a118, a119, a1110, a121, a123, a124, a125, a126, a127, a128, a129, a1210, a1211, a131, a133, a134, a135, a136, a137, a138, a139, a1310, a1311, a1312, a141, a143, a144, a145, a146, a147, a148, a149, a1410, a1411, a1412, a1413, a151, a153, a154, a155, a156, a157, a158, a159, a1510, a1511, a1512, a1513, a1514, a161, a163, a164, a165, a166, a167, a168, a169, a1610, a1611, a1612, a1613, a1614, a1615, a171, a173, a174, a175, a176, a177, a178, a179, a1710, a1711, a1712, a1713, a1714, a1715, b1, b7, b8, b9, b10, b11, b12, b13, b14, b15, bp1, bp7, bp8, bp9, bp10, bp11, bp12, bp13, bp14, bp15, bp16, bp17, btilde1, btilde7, btilde8, btilde9, btilde10, btilde11, btilde12, btilde13, btilde14, btilde15, bptilde1, bptilde7, bptilde8, bptilde9, bptilde10, bptilde11, bptilde12, bptilde13, bptilde14, bptilde15, bptilde16, bptilde17)
end

function DPRKN12ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
Base.@pure function DPRKN12ConstantCache(::Type{T},::Type{T2}) where {T<:CompiledFloats,T2<:CompiledFloats}
DPRKN12ConstantCache(
T2(2.0e-2),
T2(4.0e-2),
Expand Down

0 comments on commit 696b119

Please sign in to comment.