- Temporarily use omcTearing for strong components bigger than 200 va…

…riables until problem with Cellier Tearing is fixed git-svn-id: https://openmodelica.org/svn/OpenModelica/trunk@24116 f25d12d1-65f4-0310-ae8a-bbce733d8d8e
OpenModelica · Jan 20, 2015 · 77598f1 · 77598f1
1 parent cf31ae1
commit 77598f1
Showing 1 changed file with 13 additions and 2 deletions.
diff --git a/Compiler/BackEnd/Tearing.mo b/Compiler/BackEnd/Tearing.mo
@@ -241,6 +241,7 @@ algorithm
       Option<list<tuple<Integer, Integer, BackendDAE.Equation>>> ojac;
       BackendDAE.JacobianType jacType;
       Boolean mixedSystem;
+	  TearingMethod inMethod2;
 
     case ((BackendDAE.EQUATIONSYSTEM(eqns=eindex, vars=vindx, jac=BackendDAE.FULL_JACOBIAN(ojac), jacType=jacType, mixedSystem=mixedSystem)), _, _, _) equation
       equality(jacType = BackendDAE.JAC_LINEAR());
@@ -256,7 +257,12 @@ algorithm
       if Flags.isSet(Flags.TEARING_DUMPVERBOSE) then
         print("Jacobian:\n" + BackendDump.dumpJacobianStr(ojac) + "\n\n");
       end if;
-      (comp1, true) = callTearingMethod(inMethod, isyst, ishared, eindex, vindx, ojac, jacType, mixedSystem);
+      if listLength(vindx) > 200 then
+        inMethod2 = OMC_TEARING();
+      else
+	    inMethod2 = inMethod;
+      end if;
+      (comp1, true) = callTearingMethod(inMethod2, isyst, ishared, eindex, vindx, ojac, jacType, mixedSystem);
     then (comp1, true);
 
     // tearing of non-linear systems
@@ -268,7 +274,12 @@ algorithm
       if Flags.isSet(Flags.TEARING_DUMPVERBOSE) then
         print("Jacobian:\n" + BackendDump.dumpJacobianStr(ojac) + "\n\n");
       end if;
-      (comp1, true) = callTearingMethod(inMethod, isyst, ishared, eindex, vindx, ojac, jacType, mixedSystem);
+      if listLength(vindx) > 200 then
+        inMethod2 = OMC_TEARING();
+      else
+	    inMethod2 = inMethod;
+      end if;
+      (comp1, true) = callTearingMethod(inMethod2, isyst, ishared, eindex, vindx, ojac, jacType, mixedSystem);
     then (comp1, true);
 
     // no component for tearing