Opened 12 years ago

Closed 11 years ago

#1810 closed defect (invalid)

omc simulation runtime c crashes on the HumMod model

Reported by: Pavol Privitzer Owned by: Willi Braun
Priority: blocker Milestone: 1.9.0
Component: Backend Version: trunk
Keywords: Cc: Jens Frenkel

Description

Omc crashes on the HumMod model translation.
The test is in the testsuite/3rdParty/HumMod/translateHumMod.mos .
Note that the .mo files are taken from HumMod_OMC/ directory which is created by the "failingtest" of the testing Makefile.
(there are some other issues with the model coming from Dymola/OMC incompatibility, which will be reported after this one is resolved; they are "preprocessed" by the testing Makefile).

See the translateHumMod.mos for the output from omc after it crashes (on Windows).

Change History (8)

comment:1 by Pavol Privitzer, 12 years ago

Two more notes:

  • it is preferable to translate the model with +d=nogen option (as is done in the test Makefile), otherwise omc is trying to read initialization data from file (it is not there for this test) and on every read (on Windows) it is creating a .dll for external functions
  • after the backend will come into generation phase, it is expected another crash beacause of the "madness" function extractValuesAndDims() in SimCodeUtil(as we call it; search for the madness word in the SimCodeUtil) ... HumMod model contains some large strong components with more discrete unknowns

comment:2 by Jens Frenkel, 12 years ago

  1. test

IndexReduction.differentiateEqns failed for eqn 2572:
gases.carbonDioxyd.veinsCO2.cdCO2 * 10.0 (gases.acidBase.Veins_pH.k - gases.carbonDioxyd.veinsCO2.pK) = water.busConnector.CO2Veins_cHCO3
IndexReduction.differentiateEqns failed for eqn 2565:
der(gases.carbonDioxyd.veinsCO2.tCO2_P) = der(water.busConnector.CO2Veins_cHCO3) + der(gases.carbonDioxyd.veinsCO2.cdCO2)
IndexReduction.differentiateEqns failed for eqn 2565:
der(gases.carbonDioxyd.veinsCO2.tCO2_P) = der(water.busConnector.CO2Veins_cHCO3) + der(gases.carbonDioxyd.veinsCO2.cdCO2)

and some more

in reply to:  2 comment:3 by Jens Frenkel, 12 years ago

Replying to jfrenkel:

  1. test

IndexReduction.differentiateEqns failed for eqn 2572:
gases.carbonDioxyd.veinsCO2.cdCO2 * 10.0 (gases.acidBase.Veins_pH.k - gases.carbonDioxyd.veinsCO2.pK) = water.busConnector.CO2Veins_cHCO3

461: gases.carbonDioxyd.veinsCO2.cdCO2:STATE(stateSelect=StateSelect.never ) .HumMod.Main.test.Main_test, .HumMod.Gases.Gases, .HumMod.Gases.CO2.CO2, .HumMod.Gases.CO2.MeassureBloodCO2, .Real type: Real
1281: gases.acidBase.Veins_pH.k:PARAM(fixed = true ) = Physiolibrary.Utilities.readRealParameter("setup/default.txt", gases.acidBase.Veins_pH.varName) .HumMod.Main.test.Main_test, .HumMod.Gases.Gases, .HumMod.Gases.AcidBase.AcidBase_variables, .HumMod.Gases.AcidBase.AcidBase_variables.Variable, .Real type: Real (fixed = true )
463: gases.carbonDioxyd.veinsCO2.pK:STATE(stateSelect=StateSelect.never ) .HumMod.Main.test.Main_test, .HumMod.Gases.Gases, .HumMod.Gases.CO2.CO2, .HumMod.Gases.CO2.MeassureBloodCO2, .Real type: Real
75: water.busConnector.CO2Veins_cHCO3:STATE(stateSelect=StateSelect.never ) .HumMod.Main.test.Main_test type: Real

IndexReduction.differentiateEqns failed for eqn 2565:
der(gases.carbonDioxyd.veinsCO2.tCO2_P) = der(water.busConnector.CO2Veins_cHCO3) + der(gases.carbonDioxyd.veinsCO2.cdCO2)
IndexReduction.differentiateEqns failed for eqn 2565:
der(gases.carbonDioxyd.veinsCO2.tCO2_P) = der(water.busConnector.CO2Veins_cHCO3) + der(gases.carbonDioxyd.veinsCO2.cdCO2)

and some more

comment:4 by Jens Frenkel, 12 years ago

Owner: changed from somebody to Jens Frenkel
Status: newassigned

comment:5 by Jens Frenkel, 12 years ago

I added rules for derive rx
and log10(x). This solves the causalization problems.

comment:6 by Jens Frenkel, 12 years ago

No. of Equations: 27713
No. of Variables: 27713
##########################################################
Statistics
##########################################################
Number of independent Subsystems: 2
selected States: cardioVascularSystem.heart.RightAtrium.stateValue, cardioVascularSystem.heart.rightVentricle.ventricle.stateValue, cardioVascularSystem.heart.LeftAtrium.stateValue, cardioVascularSystem.heart.leftVentricle.ventricle.stateValue, cardioVascularSystem.heart.atriopeptin.ANPPool.stateValue, cardioVascularSystem.heart.atriopeptin.integrator.stateValue, cardioVascularSystem.heart.atriopeptin.integrator1.stateValue, cardioVascularSystem.pulmonaryCirculation.pulmCapys.stateValue, cardioVascularSystem.pulmonaryCirculation.pulmArty.stateValue, cardioVascularSystem.pulmonaryCirculation.pulmVeins.stateValue, cardioVascularSystem.systemicCirculation.systemicArtys.stateValue, cardioVascularSystem.systemicCirculation.peripheral.bone.Vasculature_Bone.integrator.stateValue, cardioVascularSystem.systemicCirculation.peripheral.brain.Vasculature_Brain.integrator.stateValue, cardioVascularSystem.systemicCirculation.peripheral.fat.Vasculature_Fat.integrator.stateValue, cardioVascularSystem.systemicCirculation.peripheral.kidney.AfferentArtery.stateValue, cardioVascularSystem.systemicCirculation.peripheral.kidney.AfferentArtery.kidneyMyogenic.stateValue, cardioVascularSystem.systemicCirculation.peripheral.kidney.TGFEffect.integrator.stateValue, cardioVascularSystem.systemicCirculation.peripheral.skin.Vasculature_Skin.integrator.stateValue, cardioVascularSystem.systemicCirculation.peripheral.otherTissue.Vasculature_OtherTissue.integrator.stateValue, cardioVascularSystem.systemicCirculation.peripheral.skeletalMuscle.Vasculature_skeletalMuscle.integrator.stateValue, cardioVascularSystem.systemicCirculation.peripheral.skeletalMuscle.metabolicVasolidation.effect.stateValue, cardioVascularSystem.systemicCirculation.peripheral.respiratoryMuscle.Vasculature.integrator.stateValue, cardioVascularSystem.systemicCirculation.leftCororaryCirculation.Vasculature.integrator.stateValue, cardioVascularSystem.systemicCirculation.rightCororaryCirculation.Vasculature.integrator.stateValue, cardioVascularSystem.systemicCirculation.splanchnicCirculation.portalVein.stateValue, cardioVascularSystem.systemicCirculation.splanchnicCirculation.GITract.Vasculature.integrator.stateValue, cardioVascularSystem.systemicCirculation.sequesteredBlood.pressureControledCompartment.stateValue, cardioVascularSystem.systemicCirculation.sequesteredBlood1.pressureControledCompartment.stateValue, cardioVascularSystem.redCells.RBCVolume.stateValue, cardioVascularSystem.redCells.EPOEffect.integrator.stateValue, nutrientsAndMetabolism.tissuesMetabolism.respiratoryMuscle.glycogen.stateValue, nutrientsAndMetabolism.tissuesMetabolism.respiratoryMuscle.InsulinEffect.integrator.stateValue, nutrientsAndMetabolism.tissuesMetabolism.skeletalMuscle.glycogen.stateValue, nutrientsAndMetabolism.tissuesMetabolism.skeletalMuscle.InsulinEffect.integrator.stateValue, nutrientsAndMetabolism.liverMetabolism.Glycogen.stateValue, nutrientsAndMetabolism.Glucose.GlucosePool.stateValue, nutrientsAndMetabolism.Glucose.bladder.Bladder.stateValue, nutrientsAndMetabolism.KetoAcids.bladder.Bladder.stateValue, nutrientsAndMetabolism.CellProteins.CellProteins.stateValue, nutrientsAndMetabolism.Urea.Medulla.stateValue, nutrientsAndMetabolism.Urea.Urea.stateValue, nutrientsAndMetabolism.Urea.UreaCell.stateValue, nutrientsAndMetabolism.Urea.bladder.Bladder.stateValue, nutrientsAndMetabolism.Lipids.LipidDeposits.stateValue, nutrientsAndMetabolism.Lipids.FAPool.stateValue, nutrientsAndMetabolism.lactate.Lactate1.stateValue, nutrientsAndMetabolism.glucagon.GlucagonPool.stateValue, nutrientsAndMetabolism.insulin.InsulinPool.stateValue, nutrientsAndMetabolism.insulin.InsulinStorage.stateValue, nutrientsAndMetabolism.insulin.integrator.stateValue, nutrientsAndMetabolism.insulin.integrator1.stateValue, nutrientsAndMetabolism.thyroxine.ThyroidPool.stateValue, nutrientsAndMetabolism.leptin.LeptinPool.stateValue, electrolytes.sodium.NaPool.stateValue, electrolytes.sodium.gILumen.GILumen.stateValue, electrolytes.sodium.bladder.Bladder.stateValue, electrolytes.sodium.kidneyNa.Medulla.stateValue, electrolytes.potassium.KPool.stateValue, electrolytes.potassium.gILumenK.GILumen.stateValue, electrolytes.potassium.bladder.Bladder.stateValue, electrolytes.chloride.gILumenCl.GILumen.stateValue, electrolytes.chloride.bladder.Bladder.stateValue, electrolytes.phosphate.PO4Pool.stateValue, electrolytes.phosphate.bladder.Bladder.stateValue, electrolytes.sulphate.SO4Pool.stateValue, electrolytes.sulphate.bladder.Bladder.stateValue, hormones.aldosterone.AldoPool.stateValue, hormones.aldosterone.AldoDelayed.stateValue, hormones.renin.ReninPool.stateValue, hormones.renin.ReninFree.stateValue, hormones.renin.ReninGranules.stateValue, hormones.renin.integrator.stateValue, hormones.catechols.epinephrine.EpiPool.stateValue, hormones.catechols.norepinephrine.NEPool.stateValue, hormones.erythropoietin.EPOPool.stateValue, nerves.lowPressureReceptors.adaptivePressure.stateValue, nerves.Baroreflex.adaptivePressure.stateValue, water.plasma.waterColloidOsmoticCompartment.stateValue, water.UpperTorso.extravascularH2O.stateValue, water.MiddleTorso.extravascularH2O.stateValue, water.LowerTorso.extravascularH2O.stateValue, water.GILumen.GILumen.stateValue, water.Kidney.Nephron.H2OChannels.stateValue, water.Bladder.volumeCompartement.stateValue, water.antidiureticHormone.ADHPool.stateValue, water.antidiureticHormone.SlowMass.stateValue, water.antidiureticHormone.FastMass.stateValue, water.antidiureticHormone.integrator1.stateValue, proteins.plasma.stateValue, proteins.UpperTorso.stateValue, proteins.MiddleTorso.stateValue, proteins.LowerTorso.stateValue, proteins.peritoneum.stateValue, proteins.synthesis.synthetized, proteins.degradation.degraded, proteins.Bladder.stateValue, status.skeletalMuscle.integrator.y, status.bone.integrator.y, status.fat.integrator.y, status.brain.integrator.y, status.rightHeart.integrator.y, status.respiratoryMuscle.integrator.y, status.otherTissue.integrator.y, status.liver.integrator.y, status.leftHeart.integrator.y, status.kidney.integrator.y, status.GITract.integrator.y, status.skin.integrator.y, gases.oxygen.artys.stateValue, gases.oxygen.veins.stateValue, gases.carbonDioxyd.veins.stateValue, gases.carbonDioxyd.artys.stateValue, gases.ventilation.respiratoryNeuralDrive2_1.peripheralChemoreceptors.integrator.stateValue, heat.bladderHeat.Bladder.stateValue, heat.gILumenHeat.GILumen.stateValue, heat.skinHeat.skin.stateValue, heat.muscleHeat.skeletalMuscle.stateValue, heat.hypothalamus.HypothalamusShiveringAcclimation.integrator.stateValue, heat.hypothalamus.HypothalamusSweatingAcclimation1.integrator.stateValue, heat.sweatGland.AcclimationEffect.integrator.stateValue, heat.sweatGland.SweatFuel.stateValue, heat.core1.core.stateValue
Number of States: 122
Toplevel Inputs: 0

Single Equations: 2133
Array Equations: 0
Algorithms: 0
Complex Equations: 0

Equationsystems with constant Jacobian: 0 {}
Equationsystems with time varying Jacobian: 10 {2, 2, 2, 2, 2, 2, 2, 2, 2, 2}
Equationsystems with nonlinear Jacobian: 1 {2}
Equationsystems without analytic Jacobian: 0 {}

mixed Equationsystems with Single Equation: 0 {}
mixed Equationsystems with Array Equation: 0 {}
mixed Equationsystems with Algorithm: 0 {}
mixed Equationsystems with Complex Equation: 0 {}
mixed Equationsystems with constant Jacobian: 0 {}
mixed Equationsystems with time varying Jacobian: 0 {}
mixed Equationsystems with nonlinear Jacobian: 0 {}
mixed Equationsystems without analytic Jacobian: 0 {}
torn linear Equationsystems: 4 {1 2, 1 3, 1 2, 2 6}
torn nonlinear Equationsystems: 16 {23 145, 9 32, 9 32, 5 31, 9 33, 9 34, 9 35, 9 33, 7 28, 11 66, 1 4, 7 28, 9 34, 7 28, 6 15, 6 15}
##########################################################
"SimCode: The model has been translated"

comment:7 by Jens Frenkel, 12 years ago

Cc: Jens Frenkel added
Owner: changed from Jens Frenkel to Willi Braun
Summary: omc backend crashes on the HumMod modelomc simulation runtime c crashes on the HumMod model

comment:8 by Martin Sjölund, 11 years ago

Resolution: invalid
Status: assignedclosed

Contains expressions like:
blood_in.m_flow * ((if (blood_in.m_flow>0) then inStream(blood_in.C_outflow) else blood_in.C_outflow)[i]) which are not legal Modelica.

Version 0, edited 11 years ago by Martin Sjölund (next)
Note: See TracTickets for help on using tickets.