Added flag to use only one instance of algebraic loop solver in cpp

runtime

Added flag to use only one instance of algebraich loop solver in cpp
runtime

[cppRuntime] Simplified handling algebraic loops in cpp runtime
Use only one solver for linear and one for nonlinear system. If needed a
own solver can used for a system

removed accidentally added files

added new AlglooperSolver class for linear and nonlinear systems

adapted linear and nonlinear solvers on new interfaces ILinearAlgloopSolver and INonLinearAlgloopSolver

added restart attribute to solve methode of algloopsolver, to restart without initialisation of algloop system

Adapted cpp template for new algloop solver interface

fix for intialize algloops at first solve call

fix for init algloop solvers used by analytic jacobians

fix for start algebraic loop solver without initializing

Combined initialize and evaluate function from algebraic loops

for parallel simulation use multiple instances of algloopsolver

fix typo

added flag to use only one instance of algloop solver in cpp runtime

removed initiequation call for some kind of linear equations

added missing code for jacobian  algloopsolvers

fix call for create kinsol solver in cpp runtime

Belonging to [master]:
  - #2574

Author: niklwors

.gitignore

@@ -56,6 +56,8 @@ $RECYCLE.BIN/
 .Spotlight-V100
 .Trashes
 
+
+
 # Directories potentially created on remote AFP share
 .AppleDB
 .AppleDesktop
@@ -153,5 +155,7 @@ Parser/ParModelica_Lexer_BaseModelica_Lexer.h
 
 SimulationRuntime/**/test_files/*.jar
 SimulationRuntime/cpp/build_msvc/
-
+# Visual Studio 2015/2017 cache/options directory
+SimulationRuntime/cpp/.vs/
+SimulationRuntime/cpp/CMakeSettings.json
 tags

Compiler/Template/CodegenCpp.tpl

@@ -3551,7 +3551,7 @@ package Flags
   constant ConfigFlag LABELED_REDUCTION;
   constant ConfigFlag LOAD_MSL_MODEL;
   constant ConfigFlag Load_PACKAGE_FILE;
-
+   constant ConfigFlag SINGLE_INSTANCE_AGLSOLVER;
   function set
     input DebugFlag inFlag;
     input Boolean inValue;

Compiler/Template/SimCodeTV.mo

@@ -1459,6 +1459,9 @@ constant ConfigFlag IGNORE_REPLACEABLE = CONFIG_FLAG(117, "ignoreReplaceable",
     "evalRecursionLimit", NONE(), EXTERNAL(), INT_FLAG(256), NONE(),
     Util.gettext("The recursion limit used when evaluating constant function calls."));
 
+  constant ConfigFlag SINGLE_INSTANCE_AGLSOLVER = CONFIG_FLAG(127, "singleInstanceAglSolver",
+  NONE(), EXTERNAL(), BOOL_FLAG(false), NONE(),
+  Util.gettext("Sets to instantiate only  one algebraic loop solver all algebraic loops"));
 
 protected
 // This is a list of all configuration flags. A flag can not be used unless it's
@@ -1590,7 +1593,8 @@ constant list<ConfigFlag> allConfigFlags = {
   BUILDING_MODEL,
   POST_OPT_MODULES_DAE,
   EVAL_LOOP_LIMIT,
-  EVAL_RECURSION_LIMIT
+  EVAL_RECURSION_LIMIT,
+  SINGLE_INSTANCE_AGLSOLVER
 };
 
 public function new

Compiler/Util/Flags.mo

@@ -69,7 +69,7 @@ else()
   message(STATUS "Libs will be installed in ${LIBINSTALLEXT}")
 endif()
 
-PROJECT(CppSolverInterface)
+PROJECT(CppSimulationRuntime)
 SET(CMAKE_VERBOSE_MAKEFILE ON)
 LIST(APPEND CMAKE_MODULE_PATH "${CMAKE_CURRENT_LIST_DIR}/CMake")
 SET(CMAKE_POSITION_INDEPENDENT_CODE ON)
@@ -93,11 +93,14 @@ OPTION(USE_CPP_03 "USE_CPP_03" OFF)
 #Set Variables
 IF(NOT FMU_TARGET)
  IF(MSVC)
-  SET(LIBINSTALLEXT "lib/omc/cpp/msvc" CACHE STRING "library directory")
+  SET(LIBINSTALLEXT "lib/omc/cpp/msvc" CACHE STRINGz "library directory")
  ELSE(MSVC)
   SET(LIBINSTALLEXT "lib/omc/cpp" CACHE STRING "library directory")
 ENDIF()
 ENDIF(NOT FMU_TARGET)
+
+MESSAGE(STATUS "Using library folder extension" ${LIBINSTALLEXT})
+
 SET(MODELICA_MODEL "ModelicaSystem")
 SET(LIBPREFIX "OMCpp")
 IF(BUILD_SHARED_LIBS)
@@ -350,7 +353,7 @@ IF(WIN32)
   SET(MICO_LIB_HOME $ENV{OMDEV}/lib/mico-msys-mingw/)
   SET(MICO_INCLUDE_HOME  $ENV{OMDEV}/include/mico-msys-mingw/)
 
-  SET(INSTALL_OMDEV_LIBS ON)
+  SET(INSTALL_OMDEV_LIBS OFF)
 ENDIF(WIN32)
 
 # Find OpenMP

SimulationRuntime/cpp/CMakeLists.txt

@@ -268,7 +268,7 @@ install (FILES
   				DESTINATION include/omc/cpp)
 ENDIF("${CMAKE_CXX_COMPILER_ID}" STREQUAL "GNU" OR "${CMAKE_CXX_COMPILER_ID}" STREQUAL "Intel" OR "${CMAKE_CXX_COMPILER_ID}" MATCHES "Clang")
 
-install(FILES $<TARGET_PDB_FILE:${ModelicaName}> DESTINATION ${LIBINSTALLEXT} OPTIONAL)
+
 
 install(TARGETS ${ModelicaName} DESTINATION ${LIBINSTALLEXT})
 install (FILES  ${CMAKE_SOURCE_DIR}/Include/Core/Modelica.h ${CMAKE_SOURCE_DIR}/Include/Core/ModelicaDefine.h DESTINATION include/omc/cpp/Core)

SimulationRuntime/cpp/Core/Modelica/CMakeLists.txt

@@ -53,6 +53,8 @@ void Initialization::initializeSystem()
   if( _solver->stateSelection())
   {
     _system->initEquations();
+	continous_system->stepCompleted(0.0);
+
 
     /* report a warning about strange start values */
     if(_solver->stateSelection())

SimulationRuntime/cpp/Core/SimController/Initialization.cpp

@@ -15,7 +15,8 @@ install(FILES $<TARGET_PDB_FILE:${SolverName}> DESTINATION ${LIBINSTALLEXT} OPTI
 
 install(TARGETS ${SolverName} DESTINATION ${LIBINSTALLEXT})
 install(FILES
-  ${CMAKE_SOURCE_DIR}/Include/Core/Solver/IAlgLoopSolver.h
+  ${CMAKE_SOURCE_DIR}/Include/Core/Solver/ILinearAlgLoopSolver.h
+  ${CMAKE_SOURCE_DIR}/Include/Core/Solver/INonLinearAlgLoopSolver.h
   ${CMAKE_SOURCE_DIR}/Include/Core/Solver/ILinSolverSettings.h
   ${CMAKE_SOURCE_DIR}/Include/Core/Solver/INonLinSolverSettings.h
   ${CMAKE_SOURCE_DIR}/Include/Core/Solver/ISolver.h

SimulationRuntime/cpp/Core/Solver/CMakeLists.txt

@@ -20,15 +20,15 @@ AlgLoopSolverFactory::~AlgLoopSolverFactory()
 {
 }
 
-shared_ptr<IAlgLoopSolver> AlgLoopSolverFactory::createLinearAlgLoopSolver(ILinearAlgLoop* algLoop)
+shared_ptr<ILinearAlgLoopSolver> AlgLoopSolverFactory::createLinearAlgLoopSolver(shared_ptr<ILinearAlgLoop> algLoop)
 {
       try
       {
         string linsolver_name = _global_settings->getSelectedLinSolver();
 		shared_ptr<ILinSolverSettings> algsolversetting= createLinSolverSettings(linsolver_name);
 		_linalgsolversettings.push_back(algsolversetting);
-        shared_ptr<IAlgLoopSolver> algsolver= createLinSolver(algLoop,linsolver_name,algsolversetting);
-        _algsolvers.push_back(algsolver);
+        shared_ptr<ILinearAlgLoopSolver> algsolver= createLinSolver(linsolver_name,algsolversetting,algLoop);
+        _linear_algsolvers.push_back(algsolver);
         return algsolver;
       }
       catch(std::exception &arg)
@@ -39,24 +39,23 @@ shared_ptr<IAlgLoopSolver> AlgLoopSolverFactory::createLinearAlgLoopSolver(ILine
 }
 
 /// Creates a nonlinear solver according to given system of equations of type algebraic loop
-shared_ptr<IAlgLoopSolver> AlgLoopSolverFactory::createNonLinearAlgLoopSolver(INonLinearAlgLoop* algLoop)
+shared_ptr<INonLinearAlgLoopSolver> AlgLoopSolverFactory::createNonLinearAlgLoopSolver(shared_ptr<INonLinearAlgLoop> algLoop)
 {
-  if(algLoop->getDimReal() > 0)
-  {
-
-    string nonlinsolver_name = _global_settings->getSelectedNonLinSolver();
-    shared_ptr<INonLinSolverSettings> algsolversetting= createNonLinSolverSettings(nonlinsolver_name);
-    algsolversetting->setContinueOnError(_global_settings->getNonLinearSolverContinueOnError());
-    _algsolversettings.push_back(algsolversetting);
-
-    shared_ptr<IAlgLoopSolver> algsolver= createNonLinSolver(algLoop,nonlinsolver_name,algsolversetting);
-    _algsolvers.push_back(algsolver);
-    return algsolver;
-  }
-  else
-  {
-    // TODO: Throw an error message here.
-    throw ModelicaSimulationError(MODEL_FACTORY,"AlgLoop solver is not available");
-  }
+	try
+    {
+		string nonlinsolver_name = _global_settings->getSelectedNonLinSolver();
+		shared_ptr<INonLinSolverSettings> algsolversetting= createNonLinSolverSettings(nonlinsolver_name);
+		algsolversetting->setContinueOnError(_global_settings->getNonLinearSolverContinueOnError());
+		_algsolversettings.push_back(algsolversetting);
+
+		shared_ptr<INonLinearAlgLoopSolver> algsolver= createNonLinSolver(nonlinsolver_name,algsolversetting,algLoop);
+		_non_linear_algsolvers.push_back(algsolver);
+		return algsolver;
+	}
+    catch(std::exception &arg)
+    {
+      throw ModelicaSimulationError(MODEL_FACTORY,"Linear AlgLoop solver is not available");
+    }
+
 }
 /** @} */ // end of coreSystem

SimulationRuntime/cpp/Core/System/AlgLoopSolverFactory.cpp

@@ -12,13 +12,18 @@ LinearAlgLoopDefaultImplementation::LinearAlgLoopDefaultImplementation()
   ,_b(NULL)
   ,_AData(NULL)
   ,_Ax(NULL)
+  ,_x0(NULL)
+ , _firstcall(true)
 {
+
 }
 
 LinearAlgLoopDefaultImplementation::~LinearAlgLoopDefaultImplementation()
 {
   if(_b)
     delete [] _b;
+ if (_x0)
+	 delete [] _x0;
 }
 
 /// Provide number (dimension) of variables according to data type
@@ -37,6 +42,9 @@ void LinearAlgLoopDefaultImplementation::initialize()
     delete [] _b;
   _b     = new double[_dimAEq];
   memset(_b,0,_dimAEq*sizeof(double));
+  if(_x0)
+	  delete [] _x0;
+  _x0 = new double[_dimAEq];
 };
 
 void LinearAlgLoopDefaultImplementation::getb(double* res) const
@@ -52,7 +60,10 @@ bool LinearAlgLoopDefaultImplementation::getUseSparseFormat(){
 void LinearAlgLoopDefaultImplementation::setUseSparseFormat(bool value){
   _useSparseFormat = value;
 }
-
+void LinearAlgLoopDefaultImplementation::getRealStartValues(double* vars) const
+{
+    memcpy(vars, _x0, sizeof(double) * _dimAEq);
+}
 
 
 //void LinearAlgLoopDefaultImplementation::getSparseAdata(double* data, int nonzeros)