|
| 1 | +module DiffEqDevTools |
| 2 | + |
| 3 | +using DiffEqBase: AbstractODEAlgorithm |
| 4 | +using DiffEqBase, RecipesBase, RecursiveArrayTools, DiffEqNoiseProcess, StructArrays |
| 5 | +using NLsolve, LinearAlgebra, RootedTrees |
| 6 | + |
| 7 | +using LinearAlgebra, Distributed |
| 8 | + |
| 9 | +using Statistics |
| 10 | + |
| 11 | +import Base: length |
| 12 | + |
| 13 | +import DiffEqBase: AbstractODEProblem, AbstractDDEProblem, AbstractDDEAlgorithm, |
| 14 | + AbstractODESolution, AbstractRODEProblem, AbstractSDEProblem, |
| 15 | + AbstractSDDEProblem, AbstractEnsembleProblem, |
| 16 | + AbstractDAEProblem, AbstractBVProblem, @def, ConvergenceSetup, |
| 17 | + DEAlgorithm, |
| 18 | + ODERKTableau, AbstractTimeseriesSolution, ExplicitRKTableau, |
| 19 | + ImplicitRKTableau |
| 20 | + |
| 21 | +import LinearAlgebra: norm, I |
| 22 | + |
| 23 | +const TIMESERIES_ERRORS = Set([:l2, :l∞, :L2, :L∞]) |
| 24 | +const DENSE_ERRORS = Set([:L2, :L∞]) |
| 25 | +const WEAK_TIMESERIES_ERRORS = Set([:weak_l2, :weak_l∞]) |
| 26 | +const WEAK_DENSE_ERRORS = Set([:weak_L2, :weak_L∞]) |
| 27 | +const WEAK_ERRORS = union(Set([:weak_final]), |
| 28 | + WEAK_TIMESERIES_ERRORS, WEAK_DENSE_ERRORS) |
| 29 | +const ALL_ERRORS = union([:final], |
| 30 | + TIMESERIES_ERRORS, DENSE_ERRORS, WEAK_TIMESERIES_ERRORS, WEAK_DENSE_ERRORS, WEAK_ERRORS) |
| 31 | + |
| 32 | +include("benchmark.jl") |
| 33 | +include("convergence.jl") |
| 34 | +include("plotrecipes.jl") |
| 35 | +include("test_solution.jl") |
| 36 | +include("ode_tableaus.jl") |
| 37 | +include("tableau_info.jl") |
| 38 | + |
| 39 | +export ConvergenceSimulation, Shootout, ShootoutSet, TestSolution |
| 40 | + |
| 41 | +#Benchmark Functions |
| 42 | +export Shootout, ShootoutSet, WorkPrecision, WorkPrecisionSet |
| 43 | + |
| 44 | +export test_convergence, analyticless_test_convergence, appxtrue!, appxtrue |
| 45 | + |
| 46 | +export get_sample_errors |
| 47 | + |
| 48 | +#Tab Functions |
| 49 | +export stability_region, residual_order_condition, check_tableau, imaginary_stability_interval |
| 50 | + |
| 51 | +#Tableaus |
| 52 | +export deduce_Butcher_tableau |
| 53 | +export constructEuler, constructKutta3, constructRK4, constructRK438Rule, |
| 54 | + constructImplicitEuler, constructMidpointRule, constructTrapezoidalRule, |
| 55 | + constructLobattoIIIA4, constructLobattoIIIB2, constructLobattoIIIB4, |
| 56 | + constructLobattoIIIC2, constructLobattoIIIC4, constructLobattoIIICStar2, |
| 57 | + constructLobattoIIICStar4, constructLobattoIIID2, constructLobattoIIID4, |
| 58 | + constructRadauIA3, constructRadauIA5, |
| 59 | + constructRadauIIA3, constructRadauIIA5, |
| 60 | + constructRalston, constructRalston4, constructHeun, constructRKF5, |
| 61 | + constructBogakiShampine3, |
| 62 | + constructCashKarp, constructRKF8, constructDormandPrince8, |
| 63 | + constructFeagin10, constructFeagin12, constructFeagin14, |
| 64 | + constructDormandPrince8_64bit, constructRKF5, constructRungeFirst5, |
| 65 | + constructCassity5, constructLawson5, |
| 66 | + constructLutherKonen5, constructLutherKonen52, |
| 67 | + constructLutherKonen53, constructPapakostasPapaGeorgiou5, |
| 68 | + constructPapakostasPapaGeorgiou52, |
| 69 | + constructTsitouras5, constructBogakiShampine5, constructSharpSmart5, |
| 70 | + constructButcher6, constructButcher7, |
| 71 | + constructDverk, constructClassicVerner6, |
| 72 | + constructClassicVerner7, constructClassicVerner8, |
| 73 | + constructVernerRobust7, constructEnrightVerner7, |
| 74 | + constructTanakaYamashitaStable7, |
| 75 | + constructTanakaYamashitaEfficient7, constructSharpSmart7, |
| 76 | + constructSharpVerner7, |
| 77 | + constructVernerEfficient7, constructCooperVerner8, |
| 78 | + constructCooperVerner82, |
| 79 | + constructTsitourasPapakostas8, constructdverk78, constructEnrightVerner8, |
| 80 | + constructCurtis8, constructVernerRobust9, constructVernerEfficient9, |
| 81 | + constructSharp9, constructTsitouras9, |
| 82 | + constructTsitouras92, |
| 83 | + constructOno12, constructCurtis10, |
| 84 | + constructOno10, |
| 85 | + constructCurtis10, constructBaker10, |
| 86 | + constructHairer10, constructButcher63, |
| 87 | + constructButcher6, constructButcher62, |
| 88 | + constructVerner6, constructDormandPrince6, |
| 89 | + constructSharpVerner6, constructVerner9162, |
| 90 | + constructVerner916, constructVernerRobust6, |
| 91 | + constructVernerEfficient6, constructPapakostas6, constructLawson6, |
| 92 | + constructTsitourasPapakostas6, constructDormandLockyerMcCorriganPrince6, |
| 93 | + constructTanakaKasugaYamashitaYazaki6D, |
| 94 | + constructTanakaKasugaYamashitaYazaki6C, |
| 95 | + constructTanakaKasugaYamashitaYazaki6B, |
| 96 | + constructTanakaKasugaYamashitaYazaki6A, |
| 97 | + constructMikkawyEisa, constructChummund6, constructChummund62, |
| 98 | + constructHuta62, constructHuta6, constructRKF4, |
| 99 | + constructVerner7, constructVerner8, |
| 100 | + constructVerner6, constructSSPRK22, constructSSPRK33, |
| 101 | + constructSSPRK43, constructSSPRK104, |
| 102 | + constructRKO65 |
| 103 | + |
| 104 | +end # module |
0 commit comments