Comments (4)
One sanity check here is to solve the system with a different method.
from worlddynamics.jl.
The initialization inaccuracy persists after trying different methods.
In the following, we provide the code snippet for Fig. 4.69d.
fig_4_69d_variables = [
(dlm.fr, 0, 4, "fr"),
(dlm.falm, 0, 0.1, "falm"),
(lfd.lfert, 0, 600, "lfert"),
(lfd.lfd, 0, 80, "lfd"),
(lfr.lfr, 0, 80, "lfr"),
(lfr.lfrt, 0, 20, "lfrt"),
]
system = agriculture_historicalrun()
sol_4_69d = WorldDynamics.solve(system, (1900, 2100))
plotvariables(sol_4_69d, (t, 1900, 1970), fig_4_69d_variables, name="Fig. 4.69d", showlegend=true, showaxis=true, colored=true)
After solving the system, we get the following initial values of the flm
variable :
sol_4_69d[dlm.falm][1:5] = (0.08, 0.04002382698512611, 0.040040848269282996, 0.040175000234290166)
The origin of the update is coming from the interpolation function :
falm ~ interpolate(pfr, tables[:falm], ranges[:falm])
where
tables[:falm] = (0.0, 0.04, 0.07, 0.09, 0.1)
ranges[:falm] = (0, 4)
The initial values of pfr
are reported below :
sol_4_69d[dlm.pfr][1:4] = (1.0, 1.0110931817510707, 1.0217794925728327, 1.068709296379689)
from worlddynamics.jl.
As a sanity check, test if changing the solver as suggested in #60 makes any difference.
from worlddynamics.jl.
Unfortunately, it doesn't.
from worlddynamics.jl.
Related Issues (20)
- Improve E4A variable names HOT 4
- Implement figures of E4A (GL scenario) HOT 3
- Implement scenario and figures of E4A (TLTL scenario) HOT 1
- Fix documentation version
- Fix bug in interpolate function
- Missing World1, World3-91, and World3-03 figures documentation in webpage
- Improve some docstrings HOT 7
- Small improvement of function interpolate with single value HOT 1
- JuliaFormatter
- Add documentation link HOT 2
- Compatibility issues HOT 8
- DELAYN function HOT 1
- Inconsistency in variable initializations for Pop15
- Reorganization of scenarios
- Check missing and unnecessary `using` statements
- Improve documentation for `interpolate` function
- Remove unnecessary type restrictions
- Update default solver for each model to match Vensim's
- JOSS Review HOT 7
- Compatibility with ModelingToolkit 9 HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from worlddynamics.jl.