testArrayEquations.mos

// name:     testArrayEquations
// keywords: array equations, for loops
// status: correct
// teardown_command: rm -f *ArrayEquationsTest*
// cflags: -d=-newInst

setCommandLineOptions("--std=3.3"); getErrorString();
setCommandLineOptions("-d=newInst,-nfScalarize,dumpSimCode"); getErrorString();
setCommandLineOptions("--simCodeTarget=Cpp"); getErrorString();

loadString("
model ArrayEquationsTest
  function myfor
    input Real[:] u;
    input Real[size(u, 1)] previous_x;
    input Boolean isFirstTick;
    output Real[size(u, 1)] x;
  algorithm
    for i in 1:size(u,1) loop
      x[i] := if isFirstTick then previous_x[i] else previous_x[i] + u[i];
    end for;
  end myfor;

  parameter Integer n = 10;
  input Real[n] u(start = 1:n);
  Real[n] x1(each start = 1);
  Real[n] x2(each start = 2);
  Real[n] x3(each start = 3);
  Real[n] x4(each start = 4);
  output Real y1, y2, y3, y4;
equation
  when Clock(0.1) then
    // for eqation
    for i in 1:n loop
      x1[i] = if firstTick(x1[i]) then previous(x1[i]) else previous(x1[i]) + u[i];
    end for;
    // function with for loop
    x2 = myfor(u, previous(x2), firstTick());
    // array equation
    x3 = if firstTick(x3) then previous(x3) else previous(x3) + u;
  end when;
  // for equation requiring inline integration
  for i in 1:n loop
    der(x4[i]) = x3[i];
  end for;
  y1 = sum(x1);
  y2 = sum(x2);
  y3 = sum(x3);
  y4 = sum(x4);
end ArrayEquationsTest;
"); getErrorString();

simulate(ArrayEquationsTest); getErrorString();

val(y1, 1.0);
val(y2, 1.0);
val(y3, 1.0);
val(y4, 1.0);

// Result:
// true
// ""
// true
// ""
// true
// ""
// true
// ""
//
//
// *********************
// * SimCode Equations *
// *********************
//
//
// allEquations:
// ========================================
//
// ========================================
//
//
//
// odeEquations (0 systems):
// ========================================
// ========================================
//
//
//
// algebraicEquations (0 systems):
// ========================================
// ========================================
//
//
// clockPartitions (1 systems):
//
// BaseClock:
// ========================================
// Clock(RCONST 0.1
// )========================================
//
// SubPartition Vars:
// ========================================
// index:-1: $DER.x4 (no alias)  initial: 	no arrCref index:() [10]
// index:-1: x1 (no alias)  initial: 	no arrCref index:() [10] (PREVIOUS)
// index:-1: x2 (no alias)  initial: 	no arrCref index:() [10] (PREVIOUS)
// index:-1: x3 (no alias)  initial: 	no arrCref index:() [10] (PREVIOUS)
// index:-1: x4 (no alias)  initial: 	no arrCref index:() [10] (PREVIOUS)
// index:-1: y1 (no alias)  initial: 	no arrCref index:() []
// index:-1: y2 (no alias)  initial: 	no arrCref index:() []
// index:-1: y3 (no alias)  initial: 	no arrCref index:() []
// index:-1: y4 (no alias)  initial: 	no arrCref index:() []
// partition equations:
// ========================================
// 26: $CLKPRE.x1=x1 [Real[10]]
// 25: $CLKPRE.x2=x2 [Real[10]]
// 24: $CLKPRE.x3=x3 [Real[10]]
// 23: $CLKPRE.x4=x4 [Real[10]]
// 14 FOR-LOOP:  for i in (1:10) loop
// x1[i]=if firstTick() then previous(x1[i]) else previous(x1[i]) + u[i][Real]
// end for;
// 15: y1=sum(x1) [Real]
// 16: x2=ArrayEquationsTest.myfor(u, {previous(x2[1]), previous(x2[2]), previous(x2[3]), previous(x2[4]), previous(x2[5]), previous(x2[6]), previous(x2[7]), previous(x2[8]), previous(x2[9]), previous(x2[10])}, firstTick()) [Real[10]]
// 17: y2=sum(x2) [Real]
// 18: x3=if firstTick() then {previous(x3[1]), previous(x3[2]), previous(x3[3]), previous(x3[4]), previous(x3[5]), previous(x3[6]), previous(x3[7]), previous(x3[8]), previous(x3[9]), previous(x3[10])} else {previous(x3[1]), previous(x3[2]), previous(x3[3]), previous(x3[4]), previous(x3[5]), previous(x3[6]), previous(x3[7]), previous(x3[8]), previous(x3[9]), previous(x3[10])} + u [Real[10]]
// 19 FOR-LOOP:  for i in (1:10) loop
// $DER.x4[i]=x3[i][Real]
// end for;
// 20 FOR-LOOP:  for i in (1:10) loop
// x4[i]=(if firstTick() then 0.0 else $DER.x4[i]) * interval() + previous(x4[i])[Real]
// end for;
// 21: y4=sum(x4) [Real]
// 22: y3=sum(x3) [Real]removedEquations equations:
// ========================================
// SubClock:
// factor(1/1) shift(0/1) Hold Events: false
//
// ========================================
//
//
//
// initialEquations: (13)
// ========================================
// 1: $DER.x4=0.0 [Real]
// 2: x1={1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0} [Real[10]]
// 3: $CLKPRE.x1=x1 [Real[10]]
// 4: x2={2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0} [Real[10]]
// 5: $CLKPRE.x2=x2 [Real[10]]
// 6: x3={3.0, 3.0, 3.0, 3.0, 3.0, 3.0, 3.0, 3.0, 3.0, 3.0} [Real[10]]
// 7: $CLKPRE.x3=x3 [Real[10]]
// 8: x4={4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0} [Real[10]]
// 9: $CLKPRE.x4=x4 [Real[10]]
// 10: y1=0.0 [Real]
// 11: y2=0.0 [Real]
// 12: y3=0.0 [Real]
// 13: y4=0.0 [Real]
// ========================================
//
//
//
// initialEquations_lambda0: (0)
// ========================================
//
// removedInitialEquations:
// ========================================
//
// startValueEquations:
// ========================================
//
// nominalValueEquations:
// ========================================
//
// minValueEquations:
// ========================================
//
// maxValueEquations:
// ========================================
//
// parameterEquations:
// ========================================
// 27: n=10 [Integer]
//
// removedEquations:
// ========================================
//
// algorithmAndEquationAsserts:
// ========================================
//
// equationsForZeroCrossings:
// ========================================
//
// generic calls:
// ========================================
//
// jacobianEquations:
// ========================================
//
// jacobianMatrices:
// ========================================
// 	Jacobian idx: 0
//
// 	Jacobian idx: 1
//
// 	Jacobian idx: 2
//
// 	Jacobian idx: 3
//
// 	Jacobian idx: 4
//
// 	Jacobian idx: 5
//
//
// modelInfo:
// ========================================
// algVars (104)
// ----------------------
// index:0: $CLKPRE.x1 (no alias)  hideResult  initial: x1	no arrCref index:(1) [10]
// index:1: $CLKPRE.x2 (no alias)  hideResult  initial: x2	no arrCref index:(11) [10]
// index:2: $CLKPRE.x3 (no alias)  hideResult  initial: x3	no arrCref index:(21) [10]
// index:3: $CLKPRE.x4 (no alias)  hideResult  initial: x4	no arrCref index:(31) [10]
// index:4: $DER.x4 (no alias)  initial: 	no arrCref index:(41) [10]
// index:5: u (no alias)  initial: 1.0:10.0	no arrCref index:(51) [10]
// index:6: x1 (no alias)  initial: 	no arrCref index:(61) [10]
// index:7: x2 (no alias)  initial: 	no arrCref index:(71) [10]
// index:8: x3 (no alias)  initial: 	no arrCref index:(81) [10]
// index:9: x4 (no alias)  initial: 	no arrCref index:(91) [10]
// index:10: y1 (no alias)  initial: 	no arrCref index:(101) []
// index:11: y2 (no alias)  initial: 	no arrCref index:(102) []
// index:12: y3 (no alias)  initial: 	no arrCref index:(103) []
// index:13: y4 (no alias)  initial: 	no arrCref index:(104) []
// intParamVars
// ----------------------
// index:0: n (no alias)  initial: 10	no arrCref index:(105) []
// functions:
// -----------
//
// Function: ArrayEquationsTest.myfor
// 	outVars: x (Real, )  ,
// 	functionArguments: u (array(Real)[:], )  , previous_x (array(Real)[size(u, 1)], )  , isFirstTick (Boolean, )  ,
// 	variableDeclarations: x (Real, )  ,
// record SimulationResult
//     resultFile = "ArrayEquationsTest_res.mat",
//     simulationOptions = "startTime = 0.0, stopTime = 1.0, numberOfIntervals = 500, tolerance = 1e-6, method = 'dassl', fileNamePrefix = 'ArrayEquationsTest', options = '', outputFormat = 'mat', variableFilter = '.*', cflags = '', simflags = ''",
//     messages = ""
// end SimulationResult;
// ""
// 560.0
// 570.0
// 580.0
// 372.5
// endResult