id	summary	reporter	owner	description	type	status	priority	milestone	component	version	resolution	keywords	cc
2352	New Bug in model calling random function with two outputs	Peter Fritzson	Mahder Alemseged Gebremedhin	"This model has worked fine in OpenModelica for several years. It worked fine in March 2013.  It is used in several models
in DrModelica and in Peter Fritzsons book.

Now it cannot simulate. OpenModelica r17270,
see the model below.

{{{
simulate(TestRandom,stopTime=10)

record SimulationResult
    resultFile = """",
    messages = ""Simulation failed for model: TestRandom
Warning: Assuming fixed start value for the following 6 variables:
         randomSeed[1]:DISCRETE(start = 23 fixed = false ) .TestRandom, .Integer type: Integer [3]
         randomSeed[2]:DISCRETE(start = 87 fixed = false ) .TestRandom, .Integer type: Integer [3]
         randomSeed[2]:DISCRETE(start = 87 fixed = false ) .TestRandom, .Integer type: Integer [3]
         randomSeed[3]:DISCRETE(start = 187 fixed = false ) .TestRandom, .Integer type: Integer [3]
         randomSeed[1]:DISCRETE(start = 23 fixed = false ) .TestRandom, .Integer type: Integer [3]
         servTime:DISCRETE(fixed = false ) .TestRandom, .Real type: Real 
Error: Fel vid bygge av simulator. Loggfil: gcc   -falign-functions -msse2 -mfpmath=sse   -I\""C:/OpenModelica1.9.0//include/omc\"" -I.   -DOPENMODELICA_XML_FROM_FILE_AT_RUNTIME  -c -o TestRandom.o TestRandom.c
TestRandom.c: In function 'initializeStaticNLSData7':
TestRandom.c:140: error: 'INTEGER_ATTRIBUTE' has no member named 'nominal'
TestRandom.c: In function 'residualFunc7':
TestRandom.c:157: error: lvalue required as unary '&' operand
TestRandom.c:158: error: incompatible types when assigning to type 'double' from type 'integer_array'
mingw32-make: *** [TestRandom.o] Error 1
}}}




{{{#!mo
package Random

import Modelica.Math; // import might not be needed
constant Real NV_MAGICCONST=4*exp(-0.5)/sqrt(2.0);
type Seed = Real[3];

function random ""input random number generator with
external storage of the seed""
  input Seed si ""input random seed"";
  output Real x ""uniform random variate between 0 and 1"";
  output Seed so ""output random seed"";
algorithm 
  so[1] := abs(rem((171 * si[1]),30269));
  so[2] := abs(rem((172 * si[2]),30307));
  so[3] := abs(rem((170 * si[3]),30323));
  // zero is a poor Seed, therfore substitute 1;
  if so[1] == 0 then
    so[1] := 1;
  end if;
  if so[2] == 0 then
    so[2] := 1;
   end if;
  if so[3] == 0 then
    so[3] := 1;
  end if;
  x := rem((so[1]/30269.0 +so[2]/30307.0 +
    so[3]/30323.0),1.0);
end random;

function normalvariate ""normally distributed random
variable""
  input Real mu ""mean value"";
  input Real sigma ""standard deviation"";
  input Seed si ""input random seed"";
  output Real x ""gaussian random variate"";
  output Seed so ""output random seed"";
  protected 
  Seed s1, s2;
  Real z, zz, u1, u2;
  Boolean breakcond=false;
algorithm 
  s1 := si;
  u2 := 1;
  while not breakcond loop
    (u1,s2) := random(s1);
    (u2,s1) := random(s2);
    z := NV_MAGICCONST*(u1-0.5)/u2;
    zz := z*z/4.0;
    breakcond := zz <= (- Math.log(u2));
   end while;
   x := mu + z*sigma;
  so := s1;
end normalvariate;

  model TestRandom
    Real   servTime;
    parameter Real  mean =  2.0;
    parameter Real  stdev = 0.5;
     Random.Seed  randomSeed(start={23,87,187});
  equation 
    when sample(0,1) then
    (servTime,randomSeed) = Random.normalvariate(mean,stdev,pre(randomSeed));
    end when;
  end TestRandom;
end Random;
}}}"	defect	closed	high	1.9.2	Backend	trunk	fixed		Mahder Alemseged Gebremedhin