diff --git a/src/lodf_calculations.jl b/src/lodf_calculations.jl
index cc44d64a..575d69ed 100644
--- a/src/lodf_calculations.jl
+++ b/src/lodf_calculations.jl
@@ -139,6 +139,19 @@ function _pardiso_sequential_LODF!(
     @assert LinearAlgebra.ishermitian(A)
     ps = Pardiso.MKLPardisoSolver()
     Pardiso.set_matrixtype!(ps, Pardiso.REAL_SYM)
+    Pardiso.pardisoinit(ps)
+    # Pardiso.set_msglvl!(ps, Pardiso.MESSAGE_LEVEL_ON)
+    defaults = Pardiso.get_iparms(ps)
+    Pardiso.set_iparm!(ps, 1, 1)
+    for (ix, v) in enumerate(defaults[2:end])
+        Pardiso.set_iparm!(ps, ix + 1, v)
+    end
+    Pardiso.set_iparm!(ps, 2, 2)
+    Pardiso.set_iparm!(ps, 59, 2)
+    Pardiso.set_iparm!(ps, 12, 1)
+    Pardiso.set_iparm!(ps, 11, 0)
+    Pardiso.set_iparm!(ps, 13, 0)
+    Pardiso.set_iparm!(ps, 32, 1)
     #Pardiso.set_msglvl!(ps, Pardiso.MESSAGE_LEVEL_ON)
     Pardiso.set_phase!(ps, Pardiso.ANALYSIS)
     Pardiso.pardiso(