omit analyzePattern for refactorization

Signed-off-by: Jan Dinkelbach <jdinkelbach@eonerc.rwth-aachen.de>

Include/dpsim/Simulation.h

@@ -56,6 +56,13 @@ namespace DPsim {
 		/// System list
 		CPS::SystemTopology mSystem;
 
+		/// Start time point to measure calculation time
+		std::chrono::time_point<std::chrono::steady_clock> mSimulationStartTimePoint;
+		/// End time point to measure calculation time
+		std::chrono::time_point<std::chrono::steady_clock> mSimulationEndTimePoint;
+		/// Measured calculation time for simulation
+		std::chrono::duration<double> mSimulationCalculationTime;
+
 		// #### Logging ####
 		/// Simulation log level
 		CPS::Logger::Level mLogLevel;

Source/MNASolverEigenSparse.cpp

@@ -119,8 +119,8 @@ void MnaSolverEigenSparse<VarType>::recomputeSystemMatrix(Real time) {
 	for (auto comp : mMNAIntfVariableComps)
 		comp->mnaApplySystemMatrixStamp(mVariableSystemMatrix);
 
-	// Calculate factorization of current matrix
-	mLuFactorizationVariableSystemMatrix.analyzePattern(mVariableSystemMatrix);
+	// Refactorization of matrix assuming that structure remained
+	// constant by omitting analyzePattern
 	mLuFactorizationVariableSystemMatrix.factorize(mVariableSystemMatrix);
 }
 

Source/Simulation.cpp

@@ -357,9 +357,16 @@ void Simulation::start() {
 	mLog->info("Start simulation: {}", mName);
 	mLog->info("Time step: {:e}", mTimeStep);
 	mLog->info("Final time: {:e}", mFinalTime);
+
+	mSimulationStartTimePoint = std::chrono::steady_clock::now();
 }
 
 void Simulation::stop() {
+
+	mSimulationEndTimePoint = std::chrono::steady_clock::now();
+	mSimulationCalculationTime = mSimulationEndTimePoint-mSimulationStartTimePoint;
+	mLog->info("Simulation calculation time: {:.6f}", mSimulationCalculationTime.count());
+
 	mScheduler->stop();
 
 	for (auto ifm : mInterfaces)