update lodf LAPACK code

NREL-Sienna · Sep 29, 2023 · 9eea136 · 9eea136
1 parent 451381c
commit 9eea136
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/src/lodf_calculations.jl b/src/lodf_calculations.jl
@@ -123,9 +123,11 @@ function _calculate_LODF_matrix_DENSE(
             push!(m_V, 1.0 - ptdf_denominator_t[iline, iline])
         end
     end
-    lodf_t = LinearAlgebra.diagm(m_V) \ ptdf_denominator_t
-    lodf_t[LinearAlgebra.diagind(lodf_t)] .= -1.0
-    return lodf_t
+    (mV, bipiv, binfo) = getrf!(Matrix(LinearAlgebra.diagm(m_V)))
+    _binfo_check(binfo)
+    getrs!('N', mV, bipiv, ptdf_denominator_t)
+    ptdf_denominator_t[LinearAlgebra.diagind(ptdf_denominator_t)] .= -1.0
+    return ptdf_denominator_t
 end
 
 function _pardiso_sequential_LODF!(