Solve for inputs (or a combination of inputs and outputs) of an algor…

…ithm section if necessary

Therefore, the algorithm is treated as a nonlinear loop. This should be improved for discrete or hybrid systems (mixed continuous and discrete systems).
See #3081
TODO: Handle discrete algorithm outputs properly.

Compiler/BackEnd/HpcOmScheduler.mo

@@ -3549,22 +3549,22 @@ algorithm
     case(SimCode.SES_NONLINEAR(SimCode.NONLINEARSYSTEM(eqs=eqs,crefs=crefs,indexNonLinearSystem=indexNonLinearSystem,jacobianMatrix=jacobianMatrix,linearTearing=linearTearing,homotopySupport=homotopySupport,mixedSystem=mixedSystem)),_)
       equation
         eqs = List.map1(eqs,TDS_replaceSimEqSysIndex,assIn);
-        oldIdx = SimCodeUtil.eqIndex(simEqIn);
+        oldIdx = SimCodeUtil.simEqSystemIndex(simEqIn);
         newIdx = arrayGet(assIn,oldIdx);
         jacobianMatrix = TDS_replaceSimEqSysIdxInJacobianMatrix(jacobianMatrix,assIn);
         simEqSys = SimCode.SES_NONLINEAR(SimCode.NONLINEARSYSTEM(newIdx,eqs,crefs,indexNonLinearSystem,jacobianMatrix,linearTearing,homotopySupport,mixedSystem), NONE());
    then simEqSys;
     case(SimCode.SES_LINEAR(SimCode.LINEARSYSTEM(partOfMixed=partOfMixed,vars=vars,beqs=beqs,simJac=simJac,residual=eqs,jacobianMatrix=jacobianMatrix,sources=sources,indexLinearSystem=indexLinearSystem)),ass)
       equation
         eqs = List.map1(eqs,TDS_replaceSimEqSysIndex,ass);
-        oldIdx = SimCodeUtil.eqIndex(simEqIn);
+        oldIdx = SimCodeUtil.simEqSystemIndex(simEqIn);
         newIdx = arrayGet(ass,oldIdx);
         jacobianMatrix = TDS_replaceSimEqSysIdxInJacobianMatrix(jacobianMatrix,ass);
         simEqSys = SimCode.SES_LINEAR(SimCode.LINEARSYSTEM(newIdx,partOfMixed,vars,beqs,simJac,eqs,jacobianMatrix,sources,indexLinearSystem), NONE());
    then simEqSys;
     case(_,ass)
       equation
-        oldIdx = SimCodeUtil.eqIndex(simEqIn);
+        oldIdx = SimCodeUtil.simEqSystemIndex(simEqIn);
         newIdx = arrayGet(ass,oldIdx);
         simEqSys = SimCodeUtil.replaceSimEqSysIndex(simEqIn,newIdx);
    then simEqSys;
@@ -3616,7 +3616,7 @@ algorithm
    then (simEqSys,(newIdx+1,ass));
     case(_,(newIdx,ass))
       equation
-        oldIdx = SimCodeUtil.eqIndex(simEqIn);
+        oldIdx = SimCodeUtil.simEqSystemIndex(simEqIn);
         ass = arrayUpdate(ass,oldIdx,newIdx);
         simEqSys = SimCodeUtil.replaceSimEqSysIndex(simEqIn,newIdx);
    then (simEqSys,(newIdx+1,ass));
@@ -4154,7 +4154,7 @@ algorithm
       list<SimCode.SimEqSystem> eqSysLst;
     case(_,_)
       equation
-      _ = SimCodeUtil.eqIndex(eqSysIn);
+      _ = SimCodeUtil.simEqSystemIndex(eqSysIn);
       pos = List.position1OnTrue(eqSysLstIn,SimCodeUtil.equationIndexEqual,eqSysIn);
       eqSysLst = List.replaceAt(eqSysIn,pos,eqSysLstIn);
     then eqSysLst;

Compiler/SimCode/SimCode.mo

@@ -390,6 +390,14 @@ uniontype SimEqSystem
     list<DAE.Statement> statements;
   end SES_ALGORITHM;
 
+  record SES_INVERSE_ALGORITHM
+    "this should only occure inside SES_NONLINEAR"
+    Integer index;
+    list<DAE.Statement> statements;
+    list<DAE.ComponentRef> outputCrefs "this are the output crefs of the original algorithm - not of the inverted one";
+    list<DAE.ComponentRef> knownOutputCrefs "this is a subset of outputCrefs that contain only those that are already known";
+  end SES_INVERSE_ALGORITHM;
+
   record SES_LINEAR
     LinearSystem lSystem;
     Option<LinearSystem> alternativeTearing;