Changeset 92c84b1 in OpenModelica

Timestamp:

2012-11-15T22:46:18+01:00 (11 years ago)

Author:

Jens Frenkel <jens.frenkel@…>

Branches:

Added-citation-metadata, maintenance/v1.14, maintenance/v1.15, maintenance/v1.16, maintenance/v1.17, maintenance/v1.18, maintenance/v1.19, maintenance/v1.20, maintenance/v1.21, maintenance/v1.22, maintenance/v1.23, master, omlib-staging

Children:

c572316

Parents:

9339461

Message:

• bugfix ExpressionSolve

git-svn-id: ​https://openmodelica.org/svn/OpenModelica/trunk@13910 f25d12d1-65f4-0310-ae8a-bbce733d8d8e

File:

• Compiler/BackEnd/ExpressionSolve.mo (modified) (4 diffs)

Compiler/BackEnd/ExpressionSolve.mo

@@ -308 +308 @@
   (outExp,outAsserts) := matchcontinue (inExp1,inExp2,inExp3,linearExps)
     local
-      DAE.Exp lhs,dere,zeroe,rhs,e,a,z;
+      DAE.Exp lhs,dere,zeroe,rhs,e,z,a;
       DAE.ComponentRef cr,cr1;
       DAE.Exp invCr;
+      DAE.Type tp;
       list<DAE.Exp> factors;
-      DAE.Type tp;
       list<DAE.Statement> asserts;
-      String estr,se1,se2,sa;
     
      // cr = (e1(0)-e2(0))/(der(e1-e2,cr)) 
@@ -356 +355 @@
         true = Expression.isZero(inExp1);
         (e,a) = solve3(inExp2,inExp3);
-        (rhs,asserts) = solve(e,a,inExp3);
-        tp = Expression.typeof(a);
+        tp = Expression.typeof(e);
         (z,_) = Expression.makeZeroExpression(Expression.arrayDimension(tp));
-        se1 = ExpressionDump.printExpStr(inExp1);
-        se2 = ExpressionDump.printExpStr(inExp2);
-        sa = ExpressionDump.printExpStr(a);
-        estr = stringAppendList({"Singular expression ",se1," = ",se2," because ",sa," is Zero!"});
-      then
-        (rhs,DAE.STMT_ASSERT(DAE.RELATION(a,DAE.NEQUAL(tp),z,-1,NONE()),DAE.SCONST(estr),DAE.ASSERTIONLEVEL_ERROR,DAE.emptyElementSource)::asserts);
+        (rhs,asserts) = solve(e,z,inExp3);
+        asserts = generateAssert(inExp1,inExp2,inExp3,a,asserts);
+      then
+        (rhs,asserts);
        
     // swapped args: a*(b-c) = 0  solve for b     
@@ -371 +367 @@
         true = Expression.isZero(inExp2);
         (e,a) = solve3(inExp1,inExp3);
-        (rhs,asserts) = solve(e,a,inExp3);
-        tp = Expression.typeof(a);
+        tp = Expression.typeof(e);
         (z,_) = Expression.makeZeroExpression(Expression.arrayDimension(tp));
-        se1 = ExpressionDump.printExpStr(inExp2);
-        se2 = ExpressionDump.printExpStr(inExp1);
-        sa = ExpressionDump.printExpStr(a);
-        estr = stringAppendList({"Singular expression ",se1," = ",se2," because ",sa," is Zero!"});
-      then
-        (rhs,DAE.STMT_ASSERT(DAE.RELATION(a,DAE.NEQUAL(tp),z,-1,NONE()),DAE.SCONST(estr),DAE.ASSERTIONLEVEL_ERROR,DAE.emptyElementSource)::asserts);
+        (rhs,asserts) = solve(e,z,inExp3);
+        asserts = generateAssert(inExp1,inExp2,inExp3,a,asserts);
+      then
+        (rhs,asserts);
 /*
     case (_,_,DAE.CREF(componentRef = cr),_)
@@ -422 +415 @@
 end solve2;
 
+protected function generateAssert
+  input DAE.Exp inExp1;
+  input DAE.Exp inExp2;
+  input DAE.Exp inExp3;
+  input DAE.Exp a;
+  input list<DAE.Statement> inAsserts;
+  output list<DAE.Statement> outAsserts;
+algorithm
+  (outExp,outAsserts) := matchcontinue (inExp1,inExp2,inExp3,a,inAsserts)
+    local
+      DAE.Exp a,z;
+      DAE.Type tp;
+      list<DAE.Statement> asserts;
+      String estr,se1,se2,sa;
+    case (_,_,_,_,_)
+      equation
+        // zero check already done
+        true = Expression.isConst(a);
+      then
+        inAsserts;
+    else
+      equation    
+        tp = Expression.typeof(a);
+        (z,_) = Expression.makeZeroExpression(Expression.arrayDimension(tp));
+        se1 = ExpressionDump.printExpStr(inExp2);
+        se2 = ExpressionDump.printExpStr(inExp1);
+        sa = ExpressionDump.printExpStr(a);
+        estr = stringAppendList({"Singular expression ",se1," = ",se2," because ",sa," is Zero!"});
+      then
+        DAE.STMT_ASSERT(DAE.RELATION(a,DAE.NEQUAL(tp),z,-1,NONE()),DAE.SCONST(estr),DAE.ASSERTIONLEVEL_ERROR,DAE.emptyElementSource)::inAsserts;
+  end matchcontinue;
+end generateAssert;
+
 protected function solve3
 "function: solve3