remove probing

src/ptdf_calculations.jl

@@ -315,9 +315,8 @@ function _calculate_PTDF_matrix_MKLPardiso(
             "Distibuted slack is not supported for systems with multiple reference buses.",
         )
     elseif isempty(dist_slack) && length(ref_bus_positions) != buscount
-        v1 = @elapsed Pardiso.pardiso(ps, PTDFm_t[valid_ix, :], ABA, full_BA)
-        v2 = @elapsed PTDFm_t[valid_ix, :] .= full_BA
-        @error v1 v2
+        Pardiso.pardiso(ps, PTDFm_t[valid_ix, :], ABA, full_BA)
+        PTDFm_t[valid_ix, :] .= full_BA
         return PTDFm_t
     elseif length(dist_slack) == buscount
         @info "Distributed bus"
@@ -329,7 +328,7 @@ function _calculate_PTDF_matrix_MKLPardiso(
     else
         error("Distributed bus specification doesn't match the number of buses.")
     end
-    Pardiso.set_phase!(ps, Pardiso.RELEASE_ALL)
+    # Pardiso.set_phase!(ps, Pardiso.RELEASE_ALL)
     return
 end
 
@@ -352,10 +351,8 @@ function calculate_PTDF_matrix_MKLPardiso(
     bus_lookup::Dict{Int, Int},
     dist_slack::Vector{Float64})
     A, ref_bus_positions = calculate_A_matrix(branches, buses)
-    BA, vals1 = @timed calculate_BA_matrix(branches, bus_lookup)
-    PTDFm, vals2 = @timed _calculate_PTDF_matrix_MKLPardiso(A, BA, ref_bus_positions, dist_slack)
-    @error "BA time $vals1"
-    @error "MKL time $vals2"
+    BA = calculate_BA_matrix(branches, bus_lookup)
+    PTDFm = _calculate_PTDF_matrix_MKLPardiso(A, BA, ref_bus_positions, dist_slack)
     return PTDFm, A
 end