From 7cfe888ec50ed608585f5fb2e16b8cc26d0ee62f Mon Sep 17 00:00:00 2001 From: marcosvanella Date: Wed, 27 Nov 2024 11:52:41 -0500 Subject: [PATCH] FDS Source : Make explicitly IPARM(8)=2 in MKL ULMAT/UGLMAT. --- Source/pres.f90 | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) diff --git a/Source/pres.f90 b/Source/pres.f90 index e753f95952..b5b123bab6 100644 --- a/Source/pres.f90 +++ b/Source/pres.f90 @@ -1690,7 +1690,6 @@ SUBROUTINE ULMAT_SOLVE_ZONE(NM,IPZ) CASE(MKL_PARDISO_FLAG) LIBRARY_SELECT #ifdef WITH_MKL !.. Back substitution and iterative refinement - IPARM(8) = 0 ! max numbers of iterative refinement steps PHASE = 33 ! only solving CALL PARDISO(ZM%PT_H, MAXFCT, MNUM, ZM%MTYPE, PHASE, ZM%NUNKH, & ZM%A_H, ZM%IA_H, ZM%JA_H, PERM, NRHS, IPARM, MSGLVL, ZM%F_H, ZM%X_H, ERROR) @@ -2685,7 +2684,7 @@ SUBROUTINE ULMAT_DEFINE_IPARM IPARM(4) = 0 ! no iterative-direct algorithm IPARM(5) = 0 ! no user fill-in reducing permutation IPARM(6) = 0 ! =0 solution on the first n components of x -IPARM(8) = 0 ! numbers of iterative refinement steps +IPARM(8) = 2 ! numbers of iterative refinement steps IPARM(10) =13 ! perturb the pivot elements with 1E-13 IPARM(11) = 1 ! use nonsymmetric permutation and scaling MPS IPARM(13) = 1 ! maximum weighted matching algorithm is switched-off (default for symmetric). @@ -3156,7 +3155,6 @@ SUBROUTINE GLMAT_SOLVER(T,DT) IPARM(41) = ZSL%LOWER_ROW IPARM(42) = ZSL%UPPER_ROW !.. Back substitution and iterative refinement - IPARM(8) = 0 ! max numbers of iterative refinement steps PHASE = 33 ! only solving #ifdef SINGLE_PRECISION_PSN_SOLVE ZSL%F_H_FB(1:ZSL%NUNKH_LOCAL) = REAL(ZSL%F_H(1:ZSL%NUNKH_LOCAL),FB)