|
454 | 454 | function alg_cache(alg::MPE, u, rate_prototype, ::Type{uEltypeNoUnits}, |
455 | 455 | ::Type{uBottomEltypeNoUnits}, ::Type{tTypeNoUnits}, |
456 | 456 | uprev, uprev2, f, t, dt, reltol, p, calck, |
457 | | - ::Val{false}) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
| 457 | + ::Val{false}, |
| 458 | + verbose) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
458 | 459 | if !(f isa PDSFunction || f isa ConservativePDSFunction) |
459 | 460 | throw(ArgumentError("MPE can only be applied to production-destruction systems")) |
460 | 461 | end |
@@ -503,7 +504,8 @@ get_tmp_cache(integrator, ::MPE, cache::OrdinaryDiffEqMutableCache) = (cache.σ, |
503 | 504 | function alg_cache(alg::MPE, u, rate_prototype, ::Type{uEltypeNoUnits}, |
504 | 505 | ::Type{uBottomEltypeNoUnits}, ::Type{tTypeNoUnits}, |
505 | 506 | uprev, uprev2, f, t, dt, reltol, p, calck, |
506 | | - ::Val{true}) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
| 507 | + ::Val{true}, |
| 508 | + verbose) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
507 | 509 | P = p_prototype(u, f) |
508 | 510 | σ = zero(u) |
509 | 511 | tab = MPEConstantCache(alg.small_constant_function(uEltypeNoUnits)) |
|
679 | 681 | function alg_cache(alg::MPRK22, u, rate_prototype, ::Type{uEltypeNoUnits}, |
680 | 682 | ::Type{uBottomEltypeNoUnits}, ::Type{tTypeNoUnits}, |
681 | 683 | uprev, uprev2, f, t, dt, reltol, p, calck, |
682 | | - ::Val{false}) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
| 684 | + ::Val{false}, |
| 685 | + verbose) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
683 | 686 | if !(f isa PDSFunction || f isa ConservativePDSFunction) |
684 | 687 | throw(ArgumentError("MPRK22 can only be applied to production-destruction systems")) |
685 | 688 | end |
@@ -761,7 +764,8 @@ get_tmp_cache(integrator, ::MPRK22, cache::OrdinaryDiffEqMutableCache) = (cache. |
761 | 764 | function alg_cache(alg::MPRK22, u, rate_prototype, ::Type{uEltypeNoUnits}, |
762 | 765 | ::Type{uBottomEltypeNoUnits}, ::Type{tTypeNoUnits}, |
763 | 766 | uprev, uprev2, f, t, dt, reltol, p, calck, |
764 | | - ::Val{true}) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
| 767 | + ::Val{true}, |
| 768 | + verbose) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
765 | 769 | a21, b1, b2 = get_constant_parameters(alg) |
766 | 770 | tab = MPRK22ConstantCache(a21, b1, b2, alg.small_constant_function(uEltypeNoUnits)) |
767 | 771 | tmp = zero(u) |
@@ -1121,7 +1125,8 @@ end |
1121 | 1125 | function alg_cache(alg::Union{MPRK43I, MPRK43II}, u, rate_prototype, ::Type{uEltypeNoUnits}, |
1122 | 1126 | ::Type{uBottomEltypeNoUnits}, ::Type{tTypeNoUnits}, |
1123 | 1127 | uprev, uprev2, f, t, dt, reltol, p, calck, |
1124 | | - ::Val{false}) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
| 1128 | + ::Val{false}, |
| 1129 | + verbose) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
1125 | 1130 | if !(f isa PDSFunction || f isa ConservativePDSFunction) |
1126 | 1131 | throw(ArgumentError("MPRK43 can only be applied to production-destruction systems")) |
1127 | 1132 | end |
@@ -1230,7 +1235,8 @@ end |
1230 | 1235 | function alg_cache(alg::Union{MPRK43I, MPRK43II}, u, rate_prototype, ::Type{uEltypeNoUnits}, |
1231 | 1236 | ::Type{uBottomEltypeNoUnits}, ::Type{tTypeNoUnits}, |
1232 | 1237 | uprev, uprev2, f, t, dt, reltol, p, calck, |
1233 | | - ::Val{true}) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
| 1238 | + ::Val{true}, |
| 1239 | + verbose) where {uEltypeNoUnits, uBottomEltypeNoUnits, tTypeNoUnits} |
1234 | 1240 | a21, a31, a32, b1, b2, b3, c2, c3, beta1, beta2, q1, q2 = get_constant_parameters(alg) |
1235 | 1241 | tab = MPRK43ConstantCache(a21, a31, a32, b1, b2, b3, c2, c3, |
1236 | 1242 | beta1, beta2, q1, q2, |
|
0 commit comments