-
Notifications
You must be signed in to change notification settings - Fork 6
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Fix Pop15
test for Julia1.9
#187
Conversation
@@ -2,7 +2,7 @@ _inits = Dict{Symbol, Float64}( | |||
:p1 => 5.3e7, | |||
:p2 => 2.1e8, | |||
:p3 => 2.2e8, | |||
:p4 => 1.4e8, | |||
:p4 => 1.8e8, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Checked: page 173, Appendix C of DGFW.
function historicalrunsolution() | ||
isdefined(@__MODULE__, :_solution_historicalrun) && return _solution_historicalrun | ||
global _solution_historicalrun = solve(historicalrun(), (1900, 2100)) | ||
global _solution_historicalrun = WorldDynamics.solve(historicalrun(), (1900, 2100); solver = Rodas5()) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This PR removes the AutoVern9
method on top of Rodas5
. Keeping it raises a mass matrix issue, but removing it results in oscillating plots. We're keeping AutoVern9
in all other plots. We keep investigating the issue in #60.
@variables fpu(t) = inits[:fpu] | ||
@variables lmf(t) = inits[:lmf] | ||
@variables cmi(t) = inits[:cmi] | ||
@variables cdr(t) le(t) hsapc(t) lmhs(t) lmhs1(t) lmhs2(t) lmc(t) lmp(t) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Doing this here makes the plots for pop15
work, but the other plots were working already with
@variables cdr(t) le(t) lmf(t) hsapc(t) lmhs(t) lmhs1(t) lmhs2(t) fpu(t) cmi(t) lmc(t) lmp(t)
It is unclear why this is the case and we opened issue #188 to investigate it.
Issue #182