Added RL test files and simulation scripts.

Author: joaquinffernandez

test/causalize/RLTest/TestRL1.mo

@@ -0,0 +1,41 @@
+/* Model for causalization testing =============================================
+ * Description: model with N algebraic loops.
+   Ua[i] and Uc[i] must be simultaneously solved.
+==============================================================================*/
+
+model TestRL1
+  constant Integer N = 100;
+  parameter Real U0 = 1.0;
+  Real iL[N], Ua[N], Uc[N];
+  parameter Real Ra = 1, Rb = 1, Rc = 1, L = 1;
+equation
+  L*der(iL[1]) = U0 - Ua[1];
+  for i in 2:N loop
+    L*der(iL[i]) = Uc[i-1] - Ua[i];
+  end for;
+  for i in 1:N-1 loop
+    Ua[i] = Rb*iL[i+1] + Uc[i]*Rb/Rc;
+  end for;
+  for i in 1:N-1 loop
+    iL[i] = Ua[i]/Ra + Uc[i]/Rc + iL[i+1];
+  end for;
+  iL[N] = Ua[N]/Ra + Uc[N]/Rc;
+  Ua[N] = Uc[N]*Rb/Rc;
+end TestRL1;
+
+//Algebraic Loops:
+//der(iL[1])
+//der(iL[1]) = U0-Ua[1]
+//
+//der(iL[i])
+//L*der(iL[i]) = Uc[i-1]-Ua[i]
+//
+//Ua[i]
+//Uc[i]
+//Ua[i] = Rb*iL[i+1]+Uc[i]*Rb/Rc
+//iL[i] = Ua[i]/Ra+Uc[i]/Rc+iL[i+1]
+//
+//Uc[N]
+//Ua[N]
+//iL[N] = Ua[N]/Ra+Uc[N]/Rc
+//Ua[N] = Uc[N]*Rb/Rc

test/causalize/RLTest/TestRL2.mo

@@ -0,0 +1,35 @@
+/* Model for causalization testing =============================================
+ * Description: model without algebraic loops.
+==============================================================================*/
+
+model TestRL2
+  constant Integer N=1000000;
+  Real iL[N],iR[N],uL[N];
+  parameter Real L=1,R=1,L1=1,I=1,R0=1;
+equation
+  for i in 1:N loop
+    L*der(iL[i])=uL[i];
+  end for;  
+  for i in 1:N-1 loop
+    iR[i]-iR[i+1]-iL[i]=0;
+    uL[i]-uL[i+1]-R*iR[i+1]=0;
+  end for;  
+  iR[N]-iL[N]+I=0;
+  uL[1]+(R+R0)*iR[1]=0;
+end TestRL2;
+
+//Algebraic Loops:
+//der(iL[i])
+//L*der(iL[i]) = uL[i]
+//
+//iR[i]
+//iR[i]-iR[i+1]-iL[i] = 0
+//
+//uL[i+1]
+//uL[i]-uL[i+1]-R*iR[i+1] = 0
+//
+//iR[N]
+//iR[N]-iL[N]+I = 0
+//
+//uL[1]
+//uL[1]+(R+R0)*iR[1] = 0

test/causalize/RLTest/TestRL3.mo

@@ -0,0 +1,34 @@
+/* Model for causalization testing =============================================
+ * Description: model with a large algebraic loop of size k*N.
+   iR[i] and iR[i+1] can be computed out of each other.
+   uL[i] and uL[i+1] can be computed out of each other.
+==============================================================================*/
+
+model TestRL3
+  constant Integer N=100;
+  Real iL[N],iR[N],uL[N];
+  parameter Real L=1,R=1,L1=1,U=1,R0=1;
+equation
+  for i in 1:N loop
+    L*der(iL[i])=uL[i];
+  end for;  
+  for i in 1:N-1 loop
+    iR[i]-iR[i+1]-iL[i]=0;
+    uL[i]-uL[i+1]-R*iR[i+1]=0;
+  end for;
+  U-uL[1]-R*iR[1]=0;
+  uL[N]-(iR[100]-iL[N])*R0=0;
+end TestRL3;
+
+//Algebraic Loops:
+//der(iL[i])
+//L*der(iL[i]) = uL[i]
+//
+//iR[i]
+//uL[i+1]
+//uL[1]
+//iR[N]
+//iR[i]-iR[i+1]-iL[i] = 0
+//uL[i]-uL[i+1]-R*iR[i+1] = 0
+//U-uL[1]-R*iR[1] = 0
+//uL[N]-(iR[N]-iL[N])*R0 = 0