Ticket #3926: calculateParametersBug.mos

// name:      calculateParametersBug
// keywords:  parameter
// status:    correct

loadString("
model foo
  final parameter Real p1 = 1;
  final parameter Real p2 = 2;
  final parameter Real p3 = 3;
  final parameter Real p4 = 4;
  parameter Real q1 (fixed=false);
  parameter Real q2 (fixed=false);
initial equation
  -p1 = 3 * (p3 * q1 + p4 * q2);
  -p2 = 3 * (q1 * p4 - p3 * q2);
end foo;
");getErrorString();

setCommandLineOptions("--initOptModules-=calculateStrongComponentJacobians");
simulate(foo);getErrorString(); 

val(p1, 0.0);getErrorString(); 
val(p2, 0.0);getErrorString(); 
val(p3, 0.0);getErrorString();
val(p4, 0.0);getErrorString();
val(q1, 0.0);getErrorString();
val(q2, 0.0);getErrorString();

// Correct results:
// 1.0
// ""
// 2.0
// ""
// 3.0
// ""
// 4.0
// ""
// -0.1466666666666667
// ""
// 0.02666666666666667
// ""


// Result:
// true
// ""
// true
// record SimulationResult
//     resultFile = "foo_res.mat",
//     simulationOptions = "startTime = 0.0, stopTime = 1.0, numberOfIntervals = 500, tolerance = 1e-06, method = 'dassl', fileNamePrefix = 'foo', options = '', outputFormat = 'mat', variableFilter = '.*', cflags = '', simflags = ''",
//     messages = "stdout            | warning | The default linear solver fails, the fallback solver with total pivoting is started at time 0.000000. That might raise performance issues, for more information use -lv LOG_LS.
// "
// end SimulationResult;
// ""
// 1.0
// ""
// 2.0
// ""
// 3.0
// ""
// 4.0
// ""
// 0.0
// ""
// 0.0
// ""
// endResult