Skip to content

Commit 4283554

Browse files
fix to new ODE problem setup
1 parent 5e346d6 commit 4283554

File tree

1 file changed

+8
-8
lines changed

1 file changed

+8
-8
lines changed

src/ode_premade_problems.jl

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -139,7 +139,7 @@ Lotka-Voltera Equations
139139
140140
with initial condition ``x=y=1``
141141
"""
142-
prb_ode_lotkavoltera = ODEProblem(lotka,[1;1],[0;1.0])
142+
prb_ode_lotkavoltera = ODEProblem(lotka,[1;1],(0.0,1.0))
143143

144144
## Fitzhugh-Nagumo
145145

@@ -157,7 +157,7 @@ Fitzhugh-Nagumo
157157
158158
with initial condition ``v=w=1``
159159
"""
160-
prob_ode_fitzhughnagumo = ODEProblem(fitz,[1;1],[0;1.0])
160+
prob_ode_fitzhughnagumo = ODEProblem(fitz,[1;1],(0.0,1.0))
161161

162162
#Van der Pol Equations
163163
van = @ode_def_noinvhes VanDerPol begin
@@ -179,7 +179,7 @@ with ``μ=1.0`` and ``u0=[0,\\sqrt{3}]``
179179
180180
Non-stiff parameters.
181181
"""
182-
prob_ode_vanderpol = ODEProblem(van,[0;sqrt(3)],[0;1.0])
182+
prob_ode_vanderpol = ODEProblem(van,[0;sqrt(3)],(0.0,1.0))
183183

184184
van_stiff = VanDerPol=1e6)
185185
"""Van der Pol Equations
@@ -195,7 +195,7 @@ with ``μ=10^6`` and ``u0=[0,\\sqrt{3}]``
195195
196196
Stiff parameters.
197197
"""
198-
prob_ode_vanderpol_stiff = ODEProblem(van_stiff,[0;sqrt(3)],[0;1.0])
198+
prob_ode_vanderpol_stiff = ODEProblem(van_stiff,[0;sqrt(3)],(0.0,1.0))
199199

200200
# ROBER
201201

@@ -222,7 +222,7 @@ Hairer Norsett Wanner Solving Ordinary Differential Euations I - Nonstiff Proble
222222
223223
Usually solved on `[0,1e11]`
224224
"""
225-
prob_ode_rober = ODEProblem(rober,[1.0;0.0;0.0],[0,1e11])
225+
prob_ode_rober = ODEProblem(rober,[1.0;0.0;0.0],(0.0,1e11))
226226

227227
# Three Body
228228
const threebody_μ = parse(BigFloat,"0.012277471"); const threebody_μ′ = 1 - threebody_μ
@@ -258,7 +258,7 @@ From Hairer Norsett Wanner Solving Ordinary Differential Euations I - Nonstiff P
258258
Usually solved on `t₀ = 0.0`; `T = parse(BigFloat,"17.0652165601579625588917206249")`
259259
Periodic with that setup.
260260
"""
261-
prob_ode_threebody = ODEProblem(threebody,[0.994, 0.0, 0.0, parse(BigFloat,"-2.00158510637908252240537862224")],[0.0;parse(BigFloat,"17.0652165601579625588917206249")])
261+
prob_ode_threebody = ODEProblem(threebody,[0.994, 0.0, 0.0, parse(BigFloat,"-2.00158510637908252240537862224")],(parse(BigFloat,"0.0"),parse(BigFloat,"17.0652165601579625588917206249")))
262262

263263
# Rigid Body Equations
264264

@@ -289,7 +289,7 @@ or Hairer Norsett Wanner Solving Ordinary Differential Euations I - Nonstiff Pro
289289
290290
Usually solved from 0 to 20.
291291
"""
292-
prob_ode_rigidbody = ODEProblem(rigid,[1.0,0.0,0.9],[0;20])
292+
prob_ode_rigidbody = ODEProblem(rigid,[1.0,0.0,0.9],(0.0,20.0))
293293

294294
# Pleiades Problem
295295

@@ -363,4 +363,4 @@ From Hairer Norsett Wanner Solving Ordinary Differential Euations I - Nonstiff P
363363
364364
Usually solved from 0 to 3.
365365
"""
366-
prob_ode_pleides = ODEProblem(pleides,[3.0,3.0,-1.0,-3.0,2.0,-2.0,2.0,3.0,-3.0,2.0,0,0,-4.0,4.0,0,0,0,0,0,1.75,-1.5,0,0,0,-1.25,1,0,0],[0;3])
366+
prob_ode_pleides = ODEProblem(pleides,[3.0,3.0,-1.0,-3.0,2.0,-2.0,2.0,3.0,-3.0,2.0,0,0,-4.0,4.0,0,0,0,0,0,1.75,-1.5,0,0,0,-1.25,1,0,0],(0.0,3.0))

0 commit comments

Comments
 (0)