Changeset 2f686d3 in OpenModelica

Timestamp:

2013-12-12T10:37:34+01:00 (10 years ago)

Author:

Volker Waurich <volker3189@…>

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:

278ccf63

Parents:

b875dff

Message:

• use analyticJacobian to solve reduced linear equationsystem when using +d=partlintornsystem

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

File:

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

Compiler/BackEnd/HpcOmSimCode.mo

@@ -922 +922 @@
         BackendDAE.TORNSYSTEM(tearingvars = tvarIdcs, residualequations = resEqIdcs, otherEqnVarTpl = otherEqnVarTpl, linear = linear) = comp;
         true = linear;
-        print("handle linear torn systems of size: "+&intString(listLength(tvarIdcs)+listLength(otherEqnVarTpl))+&"\n");
+        Debug.fcall(Flags.HPCOM_DUMP,print,"handle linear torn systems of size: "+&intString(listLength(tvarIdcs)+listLength(otherEqnVarTpl))+&"\n");
            //print("handle tornsystem with compnumber:"+&intString(compIdx)+&"\n");
            //BackendDump.dumpEqSystem(systIn,"the original system");
@@ -1016 +1016 @@
   list<list<BackendDAE.Var>> xa_i_lst, xa_i_lst1, r_i_lst, r_i_lst1, a_i_lst, a_i_lst1;
   list<DAE.ComponentRef> tcrs,ovcrs;
+  Option<list<tuple<Integer, Integer, BackendDAE.Equation>>> jac;
 algorithm
    // handle torn systems for the linear case
@@ -1095 +1096 @@
    
    //// get the strongComponent for the residual equations and add it at the end of the new StrongComponents
-   rComp := BackendDAE.EQUATIONSYSTEM(residualEqs,tearingVars,NONE(),BackendDAE.JAC_NO_ANALYTIC());   
+   
+   jac := buildLinearJacobian(a_i_lst);  // TODO: check this flatten
+   //print("Jac:\n" +& BackendDump.dumpJacobianStr(jac) +& "\n");
+   rComp := BackendDAE.EQUATIONSYSTEM(residualEqs,tearingVars,jac,BackendDAE.JAC_TIME_VARYING());   
+   //rComp := BackendDAE.EQUATIONSYSTEM(residualEqs,tearingVars,NONE(),BackendDAE.JAC_NO_ANALYTIC());   
+   //rComp := BackendDAE.TORNSYSTEM(tearingVars,residualEqs,{},true);
+   //BackendDump.dumpEquationList(resEqsOut,"the equatinos of the system\n");
+   //BackendDump.dumpVarList(tvars, "the vars of the system\n");
    
    oComps := List.appendElt(rComp,compsNew);
@@ -1102 +1110 @@
    //printPartLinTornInfo(tcrs,reqns,otherEqnsLst,ovcrs,xa_i_lst,g_i_lst,r_i_lst,h_i_lst,a_i_lst,hs_i_lst,hs,compsNew);
 end reduceLinearTornSystem2; 
+
+
+protected function buildLinearJacobian "builds the jac structure out of the given jacobian-entries (all on the diagonal).
+author:Waurich TUD 2013-12"
+  input list<list<BackendDAE.Var>> inElements;  //outer list refers to the row, inner list to the column
+  output Option<list<tuple<Integer, Integer, BackendDAE.Equation>>> outJac;
+protected
+  list<tuple<Integer, Integer, BackendDAE.Equation>> jac;
+algorithm
+  jac := List.fold1(List.intRange(listLength(inElements)),buildLinearJacobian1,inElements,{});  
+  jac := listReverse(jac); 
+  outJac := SOME(jac); 
+end buildLinearJacobian;
+
+
+protected function buildLinearJacobian1 "builds the jac structure out of the given jacobian-entries (all on the diagonal).
+author:Waurich TUD 2013-12"
+  input Integer rowIdx;
+  input list<list<BackendDAE.Var>> inElements;
+  input list<tuple<Integer, Integer, BackendDAE.Equation>> inJac;
+  output list<tuple<Integer, Integer, BackendDAE.Equation>> outJac;
+protected
+  list<BackendDAE.Var> elements;
+  list<tuple<Integer, Integer, BackendDAE.Equation>> jac;
+algorithm
+  elements := listGet(inElements,rowIdx);
+  outJac := List.fold2(List.intRange(listLength(inElements)),buildLinearJacobian2,elements,rowIdx,inJac);  
+end buildLinearJacobian1;
+
+
+protected function buildLinearJacobian2 "builds the jac structure out of the given jacobian-entries (all on the diagonal).
+author:Waurich TUD 2013-12"
+  input Integer colIdx;
+  input list<BackendDAE.Var> inElements;
+  input Integer rowIdx;
+  input list<tuple<Integer, Integer, BackendDAE.Equation>> inJac;
+  output list<tuple<Integer, Integer, BackendDAE.Equation>> outJac;
+protected
+  DAE.ComponentRef cref;
+  DAE.Exp exp;
+  BackendDAE.Equation eq;
+  BackendDAE.Var elem;
+  tuple<Integer,Integer,BackendDAE.Equation> entry;
+  list<tuple<Integer, Integer, BackendDAE.Equation>> jac;
+algorithm
+  elem := listGet(inElements,colIdx);
+  cref := BackendVariable.varCref(elem);
+  exp := DAE.CREF(cref,DAE.T_REAL_DEFAULT);
+  exp := DAE.UNARY(DAE.UMINUS(DAE.T_REAL_DEFAULT),exp);
+  eq := BackendDAE.RESIDUAL_EQUATION(exp,DAE.emptyElementSource,false);
+  entry := (colIdx,rowIdx,eq);
+  outJac := entry::inJac;
+end buildLinearJacobian2;