Add extractEigenvalues() function to Simulation

Signed-off-by: Georgii Tishenin <[email protected]>

dpsim/examples/cxx/Circuits/EMT_SinglePhaseRLC.cpp

@@ -48,7 +48,7 @@ int main(int argc, char* argv[]) {
 	sim.setFinalTime(finalTime);
 	sim.setDomain(Domain::EMT);
 	sim.addLogger(logger);
-	// sim.extractEigenvalues();
+	sim.extractEigenvalues();
 
 	sim.run();
 

dpsim/include/dpsim/Simulation.h

@@ -224,6 +224,10 @@ namespace DPsim {
 		/// Create the schedule for the independent tasks
 		void schedule();
 
+		// ### Eigenvalue extraction ###
+		/// Extract eigenvalues from power system conductance matrix
+		void extractEigenvalues();
+
 		/// Schedule an event in the simulation
 		void addEvent(Event::Ptr e) {
 			mEvents.addEvent(e);

dpsim/src/Simulation.cpp

@@ -403,6 +403,23 @@ Real Simulation::step() {
 	return mTime;
 }
 
+void Simulation::extractEigenvalues()
+{
+	if (!mInitialized)
+	{
+		initialize();
+	}
+	// TODO: [Georgii] throw exceptions for multiple solvers and for non-supported solver types
+	if (mSolvers.size() == 1)
+	{
+		auto mnaSolver = std::dynamic_pointer_cast<DPsim::MnaSolver<Real>>(mSolvers[0]);
+		if (mnaSolver)
+		{
+			mnaSolver->extractEigenvalues();
+		}
+	}
+}
+
 void Simulation::logStepTimes(String logName) {
 	auto stepTimeLog = Logger::get(logName, Logger::Level::info);
 	Logger::setLogPattern(stepTimeLog, "%v");