Initialization of non-fixed parameters

[Martin Sjölund]

loadString("
model M
  parameter Real r[30](each fixed=false);
  Real x = time;
initial equation
  for i in 1:30 loop
    der(r[i]) = i-1;
  end for;
end M;
");
simulate(M,numberOfIntervals=10);

generates C-code

  ...
  localData->initialResiduals[i++] = 2.0;
  ...

Which is clearly wrong

[Willi Braun]

The initial equations make no sense, since if "r" is a parameter the result of "der(r)=0".
Then you get wrong initial equations like 0=1,0=2,0=3...
May be need somewhere a check for constant equations like that, which cannot be satisfied.

[Martin Sjölund]

Oh yea... It's true der(param)=0. But there are still problems with initialization if der is not used (albeit different ones); convergence is then too slow.