Added Counter for Pivot Faults in KLUAdapter

Signed-off-by: Philipp Fensch <[email protected]>
m-mirz · Sep 26, 2023 · 1b6d19f · 1b6d19f
1 parent 279f1b5
commit 1b6d19f
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 1 deletion.
diff --git a/dpsim/include/dpsim/KLUAdapter.h b/dpsim/include/dpsim/KLUAdapter.h
@@ -38,7 +38,7 @@ namespace DPsim
 		std::vector<Int> mVaryingRows;
 
 		/// KLU-specific structs
-        klu_common mCommon;
+		klu_common mCommon;
 		klu_numeric* mNumeric = nullptr;
 		klu_symbolic* mSymbolic = nullptr;
 
@@ -47,6 +47,9 @@ namespace DPsim
 		/// AMD_ORDERING is defined in SuiteSparse/AMD
 		int mPreordering = AMD_ORDERING;
 
+		/// Count Pivot faults
+		int mPivotFaults = 0;
+
 		PARTIAL_REFACTORIZATION_METHOD mPartialRefactorizationMethod = PARTIAL_REFACTORIZATION_METHOD::FACTORIZATION_PATH;
 
 		/// Temporary value to store the number of nonzeros

diff --git a/dpsim/src/KLUAdapter.cpp b/dpsim/src/KLUAdapter.cpp
@@ -18,6 +18,7 @@ KLUAdapter::~KLUAdapter()
         klu_free_symbolic(&mSymbolic, &mCommon);
     if (mNumeric)
         klu_free_numeric(&mNumeric, &mCommon);
+    mSLog->info("Number of Pivot Faults: {}", mPivotFaults);
 }
 
 KLUAdapter::KLUAdapter()
@@ -150,6 +151,7 @@ void KLUAdapter::partialRefactorize(SparseMatrix &systemMatrix,
         if (mCommon.status == KLU_PIVOT_FAULT)
         {
             /* pivot became too small => fully factorize again */
+            mPivotFaults++;
             factorize(systemMatrix);
         }
     }