Opened 9 years ago

Last modified 8 years ago

#3926 accepted defect

Parameter calculation with numerical Jacobian leads to wrong results — at Version 1

Reported by: Patrick Täuber Owned by: Willi Braun
Priority: normal Milestone: Future
Component: Backend Version:
Keywords: Cc: Francesco Casella, Lennart Ochel

Description (last modified by Martin Sjölund)

The non-fixed parameters of the following model are calculated to 0 if calculateStrongComponentJacobians is skipped for initialization, which is obviously wrong:

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;

A testing script is attached.

Change History (2)

by Patrick Täuber, 9 years ago

Attachment: calculateParametersBug.mos​ added

comment:1 by Martin Sjölund, 9 years ago

Description: modified (diff)
Note: See TracTickets for help on using tickets.