@@ -553,7 +553,7 @@ function prepare_ida(instance::Instance, first_F_args, initial_bindings::Abstrac
553553 end
554554
555555 der_x0 = zeros (size (x0)) # todo: allow other initial guesses for der_x?
556- (F, eliminated_f, x0, der_x0, diffstates, params, states, state_sizes, state_offsets, eliminated_Ts)
556+ (F, eliminated_f, x0, der_x0, x_nominal, diffstates, params, states, state_sizes, state_offsets, eliminated_Ts)
557557end
558558
559559# Temporary until we can store it in SimulationModel
@@ -602,7 +602,7 @@ function simulate_ida(instance::Instance, t::Vector{Float64},
602602 initial_bindings[simulationModel_symbol] = initial_m
603603
604604 prep = prepare_ida (instance, [simulationModel_symbol], initial_bindings, store_eliminated= storeEliminated, need_eliminated_f= storeEliminated)
605- F, eliminated_f, x0, der_x0, diffstates, params, states, state_sizes, state_offsets, eliminated = prep
605+ F, eliminated_f, x0, der_x0, x_nominal, diffstates, params, states, state_sizes, state_offsets, eliminated = prep
606606
607607 eliminated_results = Vector[Vector {T} () for T in values (eliminated)]
608608 # temporary until we can store it in simulationModel
@@ -645,9 +645,9 @@ function simulate_ida(instance::Instance, t::Vector{Float64},
645645 if false
646646 m = ModiaSimulationModel (model_name_of (instance), F, x0, der_x0, jac;
647647 xw_states= diffstates, maxSparsity= maxSparsity, nc= 1 , hev= hev, nz= initial_m. nz_preInitial,
648- xNames= xNames)
648+ xNames= xNames, x_nominal = x_nominal )
649649 else
650- m = ModiaSimulationModel (string (model_name_of (instance)), F, x0, getVariableName;
650+ m = ModiaSimulationModel (string (model_name_of (instance)), F, x0, getVariableName; x_nominal = x_nominal,
651651 maxSparsity= maxSparsity, nc= 1 , nz= initial_m. nz_preInitial, hev= hev, jac= jac, x_fixed= diffstates)
652652 end
653653
0 commit comments