Skip to content

Changed DAE solver from IDA to ODASSL as assimulo does not support ID… #15

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.

Sign up for GitHub

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

Open
wredsen wants to merge 1 commit into
base: main
Choose a base branch
from
Open

Changed DAE solver from IDA to ODASSL as assimulo does not support ID… #15

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 4 additions & 4 deletions docs/demo_notebooks/modeltools/four-bar_linkage_model_and_simulation.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -37,8 +37,8 @@
"from symbtools.modeltools import Rz # Rotationsmatrix in the xy plane (c, -s; s, c)\n",
"import scipy.integrate as sc_integrate\n",
"import matplotlib.pyplot as pl\n",
"from assimulo.solvers import IDA # Imports the solver IDA from Assimulo\n",
"from assimulo.problem import Implicit_Problem # Imports the problem formulation from Assimulo\n",
"from assimulo.solvers import ODASSL # Imports the solver ODASSL from Assimulo\n",
"from assimulo.problem import Overdetermined_Problem # Imports the problem formulation from Assimulo\n",
"import ipydex\n",
"plt = pl\n",
"\n",
Expand Down Expand Up @@ -867,13 +867,13 @@
"metadata": {},
"outputs": [],
"source": [
"model = Implicit_Problem(dae.model_func, yy0, yyd0, t0)\n",
"model = Overdetermined_Problem(dae.model_func, yy0, yyd0, t0)\n",
"model.name = 'fourbar linkage'\n",
"\n",
"# indicate which components of y are differential (1) and which are purely algebraic (0)\n",
"# model.algvar = dae.diff_alg_vars ##:\n",
"\n",
"sim = IDA(model)\n",
"sim = ODASSL(model)\n",
"# sim.suppress_alg = True #Necessary to set to True for >1 index problem\n",
"\n",
"# values <= 30 produce lots of output\n",
Expand Down