Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

FDS Source : Make explicit IPARM(8)=2 in MKL ULMAT/UGLMAT. #13810

Merged
merged 1 commit into from
Nov 27, 2024
Merged

FDS Source : Make explicit IPARM(8)=2 in MKL ULMAT/UGLMAT. #13810

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 1 addition & 3 deletions Source/pres.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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)
Expand Down Expand Up @@ -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).
Expand Down Expand Up @@ -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)
Expand Down
Loading