-
-
Notifications
You must be signed in to change notification settings - Fork 195
/
OrdinaryDiffEq.jl
235 lines (179 loc) · 8.29 KB
/
OrdinaryDiffEq.jl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
module OrdinaryDiffEq
using Reexport
@reexport using DiffEqBase
using Logging
using MuladdMacro
using LinearAlgebra
# Interfaces
import DiffEqBase: solve!, step!, initialize!, isadaptive
# Internal utils
import DiffEqBase: ODE_DEFAULT_NORM, ODE_DEFAULT_ISOUTOFDOMAIN, ODE_DEFAULT_PROG_MESSAGE, ODE_DEFAULT_UNSTABLE_CHECK
using DiffEqOperators: DiffEqArrayOperator, DEFAULT_UPDATE_FUNC
import RecursiveArrayTools: chain, recursivecopy!
using Parameters, GenericSVD, ForwardDiff, RecursiveArrayTools,
NLsolve, DataStructures, DiffEqDiffTools
import ForwardDiff.Dual
using ExponentialUtilities
# Required by temporary fix in not in-place methods with 12+ broadcasts
# `MVector` is used by Nordsieck forms
import StaticArrays: SArray, MVector, SVector, @SVector, StaticArray
# Integrator Interface
import DiffEqBase: resize!,deleteat!,addat!,full_cache,user_cache,u_cache,du_cache,
resize_non_user_cache!,deleteat_non_user_cache!,addat_non_user_cache!,
terminate!,get_du, get_dt,get_proposed_dt,set_proposed_dt!,
u_modified!,savevalues!,add_tstop!,add_saveat!,set_reltol!,
set_abstol!, postamble!, last_step_failed,
isautodifferentiable
using DiffEqBase: check_error!, @def, @..
macro tight_loop_macros(ex)
:($(esc(ex)))
end
const CompiledFloats = Union{Float32,Float64,
ForwardDiff.Dual{ForwardDiff.Tag{T,W},K,3} where {T,W<:Union{Float64,Float32},
K<:Union{Float64,Float32}}}
include("misc_utils.jl")
include("algorithms.jl")
include("nlsolve/type.jl")
include("nlsolve/newton.jl")
include("nlsolve/functional.jl")
include("nlsolve/utils.jl")
include("caches/basic_caches.jl")
include("caches/low_order_rk_caches.jl")
include("caches/high_order_rk_caches.jl")
include("caches/low_storage_rk_caches.jl")
include("caches/ssprk_caches.jl")
include("caches/feagin_caches.jl")
include("caches/verner_caches.jl")
include("caches/sdirk_caches.jl")
include("caches/firk_caches.jl")
include("caches/kencarp_kvaerno_caches.jl")
include("caches/generic_implicit_caches.jl")
include("caches/linear_caches.jl")
include("caches/linear_nonlinear_caches.jl")
include("caches/symplectic_caches.jl")
include("caches/rosenbrock_caches.jl")
include("caches/rkn_caches.jl")
include("caches/adams_bashforth_moulton_caches.jl")
include("caches/nordsieck_caches.jl")
include("caches/bdf_caches.jl")
include("caches/rkc_caches.jl")
include("caches/extrapolation_caches.jl")
include("caches/prk_caches.jl")
include("cache_utils.jl")
include("alg_utils.jl")
include("tableaus/low_order_rk_tableaus.jl")
include("tableaus/high_order_rk_tableaus.jl")
include("tableaus/symplectic_tableaus.jl")
include("tableaus/verner_tableaus.jl")
include("tableaus/feagin_tableaus.jl")
include("tableaus/rosenbrock_tableaus.jl")
include("tableaus/sdirk_tableaus.jl")
include("tableaus/firk_tableaus.jl")
include("tableaus/rkn_tableaus.jl")
include("tableaus/rkc_tableaus.jl")
include("integrators/type.jl")
include("integrators/controllers.jl")
include("integrators/integrator_utils.jl")
include("integrators/integrator_interface.jl")
include("perform_step/fixed_timestep_perform_step.jl")
include("perform_step/symplectic_perform_step.jl")
include("perform_step/rkn_perform_step.jl")
include("perform_step/split_perform_step.jl")
include("perform_step/linear_perform_step.jl")
include("perform_step/iif_perform_step.jl")
include("perform_step/exponential_rk_perform_step.jl")
include("perform_step/explicit_rk_perform_step.jl")
include("perform_step/low_order_rk_perform_step.jl")
include("perform_step/high_order_rk_perform_step.jl")
include("perform_step/verner_rk_perform_step.jl")
include("perform_step/feagin_rk_perform_step.jl")
include("perform_step/low_storage_rk_perform_step.jl")
include("perform_step/ssprk_perform_step.jl")
include("perform_step/sdirk_perform_step.jl")
include("perform_step/kencarp_kvaerno_perform_step.jl")
include("perform_step/firk_perform_step.jl")
include("perform_step/generic_implicit_perform_step.jl")
include("perform_step/rosenbrock_perform_step.jl")
include("perform_step/threaded_rk_perform_step.jl")
include("perform_step/composite_perform_step.jl")
include("perform_step/adams_bashforth_moulton_perform_step.jl")
include("perform_step/nordsieck_perform_step.jl")
include("perform_step/bdf_perform_step.jl")
include("perform_step/rkc_perform_step.jl")
include("perform_step/extrapolation_perform_step.jl")
include("perform_step/prk_perform_step.jl")
include("dense/generic_dense.jl")
include("dense/interpolants.jl")
include("dense/rosenbrock_interpolants.jl")
include("dense/stiff_addsteps.jl")
include("dense/low_order_rk_addsteps.jl")
include("dense/verner_addsteps.jl")
include("dense/high_order_rk_addsteps.jl")
include("derivative_utils.jl")
include("nordsieck_utils.jl")
include("adams_utils.jl")
include("bdf_utils.jl")
include("rkc_utils.jl")
include("derivative_wrappers.jl")
include("iterator_interface.jl")
include("constants.jl")
include("composite_solution.jl")
include("solve.jl")
include("initdt.jl")
include("interp_func.jl")
include("composite_algs.jl")
#General Functions
export solve, solve!, init, step!
export OrdinaryDiffEqAlgorithm
#Callback Necessary
export ode_addsteps!, addsteps!, ode_interpolant,
terminate!, savevalues!, copyat_or_push!, isfsal
export constructDormandPrince
# Reexport the Alg Types
export FunctionMap, Euler, Heun, Ralston, Midpoint, RK4, ExplicitRK, OwrenZen3, OwrenZen4, OwrenZen5,
BS3, BS5, RK46NL, DP5, DP5Threaded, Tsit5, DP8, Vern6, Vern7, Vern8, TanYam7, TsitPap8,
Vern9, Feagin10, Feagin12, Feagin14, CompositeAlgorithm, Anas5
export SSPRK22, SSPRK33, SSPRK53, SSPRK53_2N1, SSPRK53_2N2, SSPRK63, SSPRK73, SSPRK83, SSPRK432,
SSPRKMSVS32, SSPRKMSVS43, SSPRK932, SSPRK54, SSPRK104
export ORK256, CarpenterKennedy2N54, HSLDDRK64, DGLDDRK73_C, DGLDDRK84_C, DGLDDRK84_F, NDBLSRK124, NDBLSRK134, NDBLSRK144,
CFRLDDRK64, TSLDDRK74,
CKLLSRK43_2,CKLLSRK54_3C,CKLLSRK95_4S,CKLLSRK95_4C,CKLLSRK95_4M,
CKLLSRK54_3C_3R, CKLLSRK54_3M_3R, CKLLSRK54_3N_3R, CKLLSRK85_4C_3R, CKLLSRK85_4M_3R, CKLLSRK85_4P_3R,
CKLLSRK54_3N_4R, CKLLSRK54_3M_4R, CKLLSRK65_4M_4R, CKLLSRK85_4FM_4R, CKLLSRK75_4M_5R,
ParsaniKetchesonDeconinck3S32, ParsaniKetchesonDeconinck3S82,
ParsaniKetchesonDeconinck3S53, ParsaniKetchesonDeconinck3S173,
ParsaniKetchesonDeconinck3S94, ParsaniKetchesonDeconinck3S184,
ParsaniKetchesonDeconinck3S105, ParsaniKetchesonDeconinck3S205,
KYK2014DGSSPRK_3S2
export RadauIIA5
export ImplicitEuler, ImplicitMidpoint, Trapezoid, TRBDF2, SDIRK2, Kvaerno3,
KenCarp3, Cash4, Hairer4, Hairer42, SSPSDIRK2, Kvaerno4, Kvaerno5,
KenCarp4, KenCarp5, ESDIRK54I8L2SA
export GenericImplicitEuler, GenericTrapezoid
export MidpointSplitting, LinearExponential
export Rosenbrock23, Rosenbrock32, RosShamp4, Veldd4, Velds4, GRK4T, GRK4A,
Ros4LStab, ROS3P, Rodas3, Rodas4, Rodas42, Rodas4P, Rodas5, RosenbrockW6S4OS
export GenericIIF1, GenericIIF2
export LawsonEuler, NorsettEuler, ETD1, ETDRK2, ETDRK3, ETDRK4, HochOst4, Exp4, EPIRK4s3A, EPIRK4s3B,
EPIRK5s3, EXPRB53s3, EPIRK5P1, EPIRK5P2, ETD2, Exprb32, Exprb43
export SymplecticEuler, VelocityVerlet, VerletLeapfrog, PseudoVerletLeapfrog,
McAte2, Ruth3, McAte3, CandyRoz4, McAte4, McAte42, McAte5,
CalvoSanz4, Yoshida6, KahanLi6, McAte8, KahanLi8, SofSpa10
export SplitEuler
export Nystrom4, Nystrom4VelocityIndependent, Nystrom5VelocityIndependent,
IRKN3, IRKN4, DPRKN6, DPRKN8, DPRKN12, ERKN4, ERKN5
export ROCK2, ROCK4, RKC, IRKC, ESERK5
export AB3, AB4, AB5, ABM32, ABM43, ABM54
export VCAB3, VCAB4, VCAB5, VCABM3, VCABM4, VCABM5
export VCABM
export IMEXEuler, CNAB2, CNLF2
export AN5, JVODE, JVODE_Adams, JVODE_BDF
export ABDF2, QNDF1, QBDF1, QNDF2, QBDF2, QNDF, QBDF
export SBDF2, SBDF3, SBDF4
export MEBDF2
export AutoSwitch, AutoTsit5, AutoDP5,
AutoVern6, AutoVern7, AutoVern8, AutoVern9
export AitkenNeville
export KuttaPRK2p5
export NLNewton, NLFunctional, NLAnderson
end # module