root of negative value

[jan.hagemann@…]

The example
testsuite/simulation/modelica/tearing/Tearing19-cel.mos
fails with activated wrapFunctionCalls-module.

The error dump:

assert            | warning | The following assertion has been violated during initialization at time 0.000000
assert            | debug   | Model error: Argument of asin($cse6) outside the domain -1.0 <= nan <= 1.0

This is strange, since cse variables are only introduced for the simulation system. The initialization is solved propably.

[Lennart Ochel]

It seems that the problem comes from a bad value that is calculated for $Px5 in function Tearing19_eqFunction_46 in Tearing19_02nls.c.

[Lennart Ochel]

Tearing19_eqFunction_46 is one of the inner equations of a tearing set:

/*
 equation index: 46
 type: SIMPLE_ASSIGN
 x5 = (1.0 - (x4 ^ 4.0 + x0 * x1 * x2)) ^ 0.2
 */
void Tearing19_eqFunction_46(DATA *data, threadData_t *threadData)
{
  TRACE_PUSH
  const int equationIndexes[2] = {1,46};
  modelica_real tmp1054;
  tmp1054 = $Px4;
  tmp1054 *= tmp1054;
  e3 = pow(1.0 - (tmp1054 * tmp1054) + ($Px0) * (($Px1) * ($Px2)), 0.2);
  TRACE_POP
}

In this case, x5 is root of a negative number, which let function pow fail.

Probably, we should not solve explicitly for roots. Patrick, what do you think?

[Lennart Ochel]

The same issue is present in ​simulation/modelica/others/TestExpressionSolve.mos.

[Lennart Ochel]

Fixed with 63d57dd2ca84c3e72641edba77fde0494bc56b73/OMCompiler.

[Martin Sjölund]

Milestone renamed

[Martin Sjölund]

Milestone renamed