Export MKLPardisoKKTSolver, fix printing

- Add paper to citing section

Author: migarstka

docs/src/citing.md

@@ -31,6 +31,18 @@ The following paper describes the algorithm used for our clique graph based merg
   issn    = {2405-8963},
   journal = {IFAC-PapersOnLine},
 }
-
 ```
 A preprint can be downloaded [here](https://arxiv.org/abs/1911.05615).
+
+The acceleration strategy is described in this paper:
+```
+@InProceedings{Garstka_2022,
+  title={Safeguarded Anderson acceleration for parametric nonexpansive operators},
+  author={Garstka, Michael and Cannon, Mark and Goulart, Paul},
+  booktitle={2022 European Control Conference (ECC)},
+  pages={435--440},
+  year={2022},
+  organization={IEEE}
+}
+```
+A preprint can be downloaded [here](https://arxiv.org/pdf/2208.02847).

src/linear_solver/kktsolver_pardiso.jl

@@ -1,6 +1,6 @@
 using .Pardiso: Pardiso, set_solver!, pardisoinit, set_phase!, set_iparm!, set_msglvl!, get_iparm, pardiso, PardisoSolver, MKLPardisoSolver, get_nprocs, set_nprocs!, set_matrixtype!, set_perm!
 using AMD
-export PardisoDirectKKTSolver, PardisoIndirectKKTSolver, MKLPardisoSolver
+export PardisoDirectKKTSolver, PardisoIndirectKKTSolver, MKLPardisoKKTSolver
 
 abstract type AbstractPardisoKKTSolver <: AbstractKKTSolver end
 

src/printing.jl

@@ -25,7 +25,7 @@ function print_header(ws::COSMO.Workspace{T}) where {T <: AbstractFloat}
 	if ws.ci.decompose && ws.ci.num_decomposable > 0
 		println("Decomp:   Num of original PSD cones: $(ws.ci.num_psd_cones)\n" * " "^10*"Num decomposable PSD cones: $(ws.ci.num_decomposable)\n" * " "^10*"Num PSD cones after decomposition: $(ws.ci.num_decom_psd_cones)\n" * " "^10*"Merge Strategy: $(stringify(ws.settings.merge_strategy))")
 	end
-	println("Settings: ϵ_abs = $(@sprintf("%.1e", settings.eps_abs)), ϵ_rel = $(@sprintf("%.1e", settings.eps_rel)),\n" * " "^10 * "ϵ_prim_inf = $(@sprintf("%.1e", settings.eps_prim_inf)), ϵ_dual_inf = $(@sprintf("%.1e", settings.eps_dual_inf)),\n" * " "^10 * "ρ = $(@sprintf("%.3g", settings.rho)), σ = $(@sprintf("%.3g", settings.sigma)), α = $(@sprintf("%.3g", settings.alpha)),\n" * " "^10 * "max_iter = $(settings.max_iter),\n" * " "^10 * "scaling iter = $(settings.scaling) ($(scaling_status)),\n" * " "^10 * "check termination every $(settings.check_termination) iter,\n" * " "^10 * "check infeasibility every $(settings.check_infeasibility) iter,\n" *	 " "^10 * "KKT system solver: $(print_lin_sys(settings.kkt_solver.ObjectType))")
+	println("Settings: ϵ_abs = $(@sprintf("%.1e", settings.eps_abs)), ϵ_rel = $(@sprintf("%.1e", settings.eps_rel)),\n" * " "^10 * "ϵ_prim_inf = $(@sprintf("%.1e", settings.eps_prim_inf)), ϵ_dual_inf = $(@sprintf("%.1e", settings.eps_dual_inf)),\n" * " "^10 * "ρ = $(@sprintf("%.3g", settings.rho)), σ = $(@sprintf("%.3g", settings.sigma)), α = $(@sprintf("%.3g", settings.alpha)),\n" * " "^10 * "max_iter = $(settings.max_iter),\n" * " "^10 * "scaling iter = $(settings.scaling) ($(scaling_status)),\n" * " "^10 * "check termination every $(settings.check_termination) iter,\n" * " "^10 * "check infeasibility every $(settings.check_infeasibility) iter,\n" *	 " "^10 * "KKT system solver: $(print_lin_sys(settings.kkt_solver))")
 
 	print_accelerator(ws.accelerator, ws.settings.safeguard, ws.settings.safeguard_tol, tab = 10)
 
@@ -75,7 +75,7 @@ function print_result(status::Symbol, iter::Int, safeguarding_iter::Int, cost::T
 	println("Optimal objective: $(@sprintf("%.4g", cost))\nRuntime: $(round.(rt; digits = 3))s ($(round.(rt * 1000; digits = 2))ms)\n")
 	nothing
 end
-
+print_lin_sys(s::OptionsFactory) = print_lin_sys(s.ObjectType)
 function print_lin_sys(s::Type{<:AbstractKKTSolver})
 	lin_sys_dict = Dict("QdldlKKTSolver" => "QDLDL", "CholmodKKTSolver" => "CHOLMOD", "PardisoDirectKKTSolver" => "Pardiso (direct)", "PardisoIndirectKKTSolver" => "Pardiso (indirect)",  "MKLPardisoKKTSolver" => "MKL Pardiso")
 	return get(lin_sys_dict, string(s), string(s))