Index: CMakeLists.txt =================================================================== --- CMakeLists.txt (revision 11543) +++ CMakeLists.txt (working copy) @@ -33,6 +33,10 @@ message(STATUS "Platform: Linux-PGI") set(CMAKE_Fortran_FLAGS_DEBUG "-fPIC -g -traceback -Mbackslash") set(CMAKE_Fortran_FLAGS_RELEASE "-fPIC -O3 -traceback -Mbackslash") + elseif(CMAKE_Fortran_COMPILER_ID STREQUAL "Flang") + message(STATUS "Platform: Linux-Flang") + set(CMAKE_Fortran_FLAGS_DEBUG "-fPIC -g") + set(CMAKE_Fortran_FLAGS_RELEASE "-fPIC -O3") else(CMAKE_Fortran_COMPILER_ID STREQUAL "GNU") message(STATUS "Platform: Linux-Intel") set(CMAKE_Fortran_FLAGS_DEBUG "-fPIC -WB -traceback -g") Index: SRC/cnaitr.f =================================================================== --- SRC/cnaitr.f (revision 11543) +++ SRC/cnaitr.f (working copy) @@ -380,7 +380,7 @@ if (msglvl .gt. 1) then call ivout (logfil, 1, j, ndigit, & '_naitr: generating Arnoldi vector number') - call svout (logfil, 1, rnorm, ndigit, + call svout (logfil, 1, (/ rnorm /), ndigit, & '_naitr: B-norm of the current residual is') end if c Index: SRC/cnaup2.f =================================================================== --- SRC/cnaup2.f (revision 11543) +++ SRC/cnaup2.f (working copy) @@ -389,7 +389,7 @@ iter = iter + 1 c if (msglvl .gt. 0) then - call ivout (logfil, 1, iter, ndigit, + call ivout (logfil, 1, (/ iter /), ndigit, & '_naup2: **** Start of major iteration number ****') end if c @@ -402,9 +402,9 @@ np = kplusp - nev c if (msglvl .gt. 1) then - call ivout (logfil, 1, nev, ndigit, + call ivout (logfil, 1, (/ nev /), ndigit, & '_naup2: The length of the current Arnoldi factorization') - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_naup2: Extend the Arnoldi factorization by') end if c @@ -658,7 +658,7 @@ end if c if (msglvl .gt. 0) then - call ivout (logfil, 1, nconv, ndigit, + call ivout (logfil, 1, (/ nconv /), ndigit, & '_naup2: no. of "converged" Ritz values at this iter.') if (msglvl .gt. 1) then kp(1) = nev @@ -698,7 +698,7 @@ end if c if (msglvl .gt. 2) then - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_naup2: The number of shifts to apply ') call cvout (logfil, np, ritz, ndigit, & '_naup2: values of the shifts') Index: SRC/dgetv0.f =================================================================== --- SRC/dgetv0.f (revision 11543) +++ SRC/dgetv0.f (working copy) @@ -364,9 +364,9 @@ c %--------------------------------------% c if (msglvl .gt. 2) then - call dvout (logfil, 1, rnorm0, ndigit, + call dvout (logfil, 1, (/ rnorm0 /), ndigit, & '_getv0: re-orthonalization ; rnorm0 is') - call dvout (logfil, 1, rnorm, ndigit, + call dvout (logfil, 1, (/ rnorm /), ndigit, & '_getv0: re-orthonalization ; rnorm is') end if c @@ -397,7 +397,7 @@ 50 continue c if (msglvl .gt. 0) then - call dvout (logfil, 1, rnorm, ndigit, + call dvout (logfil, 1, (/ rnorm /), ndigit, & '_getv0: B-norm of initial / restarted starting vector') end if if (msglvl .gt. 3) then Index: SRC/dnaitr.f =================================================================== --- SRC/dnaitr.f (revision 11543) +++ SRC/dnaitr.f (working copy) @@ -373,7 +373,7 @@ if (msglvl .gt. 1) then call ivout (logfil, 1, j, ndigit, & '_naitr: generating Arnoldi vector number') - call dvout (logfil, 1, rnorm, ndigit, + call dvout (logfil, 1, (/ rnorm /), ndigit, & '_naitr: B-norm of the current residual is') end if c Index: SRC/dnaup2.f =================================================================== --- SRC/dnaup2.f (revision 11543) +++ SRC/dnaup2.f (working copy) @@ -388,7 +388,7 @@ iter = iter + 1 c if (msglvl .gt. 0) then - call ivout (logfil, 1, iter, ndigit, + call ivout (logfil, 1, (/ iter /), ndigit, & '_naup2: **** Start of major iteration number ****') end if c @@ -401,9 +401,9 @@ np = kplusp - nev c if (msglvl .gt. 1) then - call ivout (logfil, 1, nev, ndigit, + call ivout (logfil, 1, (/ nev /), ndigit, & '_naup2: The length of the current Arnoldi factorization') - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_naup2: Extend the Arnoldi factorization by') end if c @@ -435,7 +435,7 @@ update = .false. c if (msglvl .gt. 1) then - call dvout (logfil, 1, rnorm, ndigit, + call dvout (logfil, 1, (/ rnorm /), ndigit, & '_naup2: Corresponding B-norm of the residual') end if c @@ -678,7 +678,7 @@ end if c if (msglvl .gt. 0) then - call ivout (logfil, 1, nconv, ndigit, + call ivout (logfil, 1, (/ nconv /), ndigit, & '_naup2: no. of "converged" Ritz values at this iter.') if (msglvl .gt. 1) then kp(1) = nev @@ -730,7 +730,7 @@ end if c if (msglvl .gt. 2) then - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_naup2: The number of shifts to apply ') call dvout (logfil, np, ritzr, ndigit, & '_naup2: Real part of the shifts') @@ -796,7 +796,7 @@ cnorm = .false. c if (msglvl .gt. 2) then - call dvout (logfil, 1, rnorm, ndigit, + call dvout (logfil, 1, (/ rnorm /), ndigit, & '_naup2: B-norm of residual for compressed factorization') call dmout (logfil, nev, nev, h, ldh, ndigit, & '_naup2: Compressed upper Hessenberg matrix H') Index: SRC/dsaitr.f =================================================================== --- SRC/dsaitr.f (revision 11543) +++ SRC/dsaitr.f (working copy) @@ -366,7 +366,7 @@ if (msglvl .gt. 2) then call ivout (logfil, 1, j, ndigit, & '_saitr: generating Arnoldi vector no.') - call dvout (logfil, 1, rnorm, ndigit, + call dvout (logfil, 1, (/ rnorm /), ndigit, & '_saitr: B-norm of the current residual =') end if c Index: SRC/dsaup2.f =================================================================== --- SRC/dsaup2.f (revision 11543) +++ SRC/dsaup2.f (working copy) @@ -402,13 +402,13 @@ iter = iter + 1 c if (msglvl .gt. 0) then - call ivout (logfil, 1, iter, ndigit, + call ivout (logfil, 1, (/ iter /), ndigit, & '_saup2: **** Start of major iteration number ****') end if if (msglvl .gt. 1) then - call ivout (logfil, 1, nev, ndigit, + call ivout (logfil, 1, (/ nev /), ndigit, & '_saup2: The length of the current Lanczos factorization') - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_saup2: Extend the Lanczos factorization by') end if c @@ -446,7 +446,7 @@ update = .false. c if (msglvl .gt. 1) then - call dvout (logfil, 1, rnorm, ndigit, + call dvout (logfil, 1, (/ rnorm /), ndigit, & '_saup2: Current B-norm of residual for factorization') end if c @@ -694,7 +694,7 @@ end if c if (msglvl .gt. 0) then - call ivout (logfil, 1, nconv, ndigit, + call ivout (logfil, 1, (/ nconv /), ndigit, & '_saup2: no. of "converged" Ritz values at this iter.') if (msglvl .gt. 1) then kp(1) = nev @@ -742,7 +742,7 @@ if (ishift .eq. 0) call dcopy (np, workl, 1, ritz, 1) c if (msglvl .gt. 2) then - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_saup2: The number of shifts to apply ') call dvout (logfil, np, workl, ndigit, & '_saup2: shifts selected') @@ -809,7 +809,7 @@ 130 continue c if (msglvl .gt. 2) then - call dvout (logfil, 1, rnorm, ndigit, + call dvout (logfil, 1, (/ rnorm /), ndigit, & '_saup2: B-norm of residual for NEV factorization') call dvout (logfil, nev, h(1,2), ndigit, & '_saup2: main diagonal of compressed H matrix') Index: SRC/sgetv0.f =================================================================== --- SRC/sgetv0.f (revision 11543) +++ SRC/sgetv0.f (working copy) @@ -364,9 +364,9 @@ c %--------------------------------------% c if (msglvl .gt. 2) then - call svout (logfil, 1, rnorm0, ndigit, + call svout (logfil, 1, (/ rnorm0 /), ndigit, & '_getv0: re-orthonalization ; rnorm0 is') - call svout (logfil, 1, rnorm, ndigit, + call svout (logfil, 1, (/ rnorm /), ndigit, & '_getv0: re-orthonalization ; rnorm is') end if c @@ -397,7 +397,7 @@ 50 continue c if (msglvl .gt. 0) then - call svout (logfil, 1, rnorm, ndigit, + call svout (logfil, 1, (/ rnorm /), ndigit, & '_getv0: B-norm of initial / restarted starting vector') end if if (msglvl .gt. 3) then Index: SRC/snaitr.f =================================================================== --- SRC/snaitr.f (revision 11543) +++ SRC/snaitr.f (working copy) @@ -373,7 +373,7 @@ if (msglvl .gt. 1) then call ivout (logfil, 1, j, ndigit, & '_naitr: generating Arnoldi vector number') - call svout (logfil, 1, rnorm, ndigit, + call svout (logfil, 1, (/ rnorm /), ndigit, & '_naitr: B-norm of the current residual is') end if c Index: SRC/snaup2.f =================================================================== --- SRC/snaup2.f (revision 11543) +++ SRC/snaup2.f (working copy) @@ -388,7 +388,7 @@ iter = iter + 1 c if (msglvl .gt. 0) then - call ivout (logfil, 1, iter, ndigit, + call ivout (logfil, 1, (/ iter /), ndigit, & '_naup2: **** Start of major iteration number ****') end if c @@ -401,9 +401,9 @@ np = kplusp - nev c if (msglvl .gt. 1) then - call ivout (logfil, 1, nev, ndigit, + call ivout (logfil, 1, (/ nev /), ndigit, & '_naup2: The length of the current Arnoldi factorization') - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_naup2: Extend the Arnoldi factorization by') end if c @@ -435,7 +435,7 @@ update = .false. c if (msglvl .gt. 1) then - call svout (logfil, 1, rnorm, ndigit, + call svout (logfil, 1, (/ rnorm /), ndigit, & '_naup2: Corresponding B-norm of the residual') end if c @@ -678,7 +678,7 @@ end if c if (msglvl .gt. 0) then - call ivout (logfil, 1, nconv, ndigit, + call ivout (logfil, 1, (/ nconv /), ndigit, & '_naup2: no. of "converged" Ritz values at this iter.') if (msglvl .gt. 1) then kp(1) = nev @@ -730,7 +730,7 @@ end if c if (msglvl .gt. 2) then - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_naup2: The number of shifts to apply ') call svout (logfil, np, ritzr, ndigit, & '_naup2: Real part of the shifts') @@ -796,7 +796,7 @@ cnorm = .false. c if (msglvl .gt. 2) then - call svout (logfil, 1, rnorm, ndigit, + call svout (logfil, 1, (/ rnorm /), ndigit, & '_naup2: B-norm of residual for compressed factorization') call smout (logfil, nev, nev, h, ldh, ndigit, & '_naup2: Compressed upper Hessenberg matrix H') Index: SRC/ssaitr.f =================================================================== --- SRC/ssaitr.f (revision 11543) +++ SRC/ssaitr.f (working copy) @@ -366,7 +366,7 @@ if (msglvl .gt. 2) then call ivout (logfil, 1, j, ndigit, & '_saitr: generating Arnoldi vector no.') - call svout (logfil, 1, rnorm, ndigit, + call svout (logfil, 1, (/ rnorm /), ndigit, & '_saitr: B-norm of the current residual =') end if c Index: SRC/ssaup2.f =================================================================== --- SRC/ssaup2.f (revision 11543) +++ SRC/ssaup2.f (working copy) @@ -402,13 +402,13 @@ iter = iter + 1 c if (msglvl .gt. 0) then - call ivout (logfil, 1, iter, ndigit, + call ivout (logfil, 1, (/ iter /), ndigit, & '_saup2: **** Start of major iteration number ****') end if if (msglvl .gt. 1) then - call ivout (logfil, 1, nev, ndigit, + call ivout (logfil, 1, (/ nev /), ndigit, & '_saup2: The length of the current Lanczos factorization') - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_saup2: Extend the Lanczos factorization by') end if c @@ -446,7 +446,7 @@ update = .false. c if (msglvl .gt. 1) then - call svout (logfil, 1, rnorm, ndigit, + call svout (logfil, 1, (/ rnorm /), ndigit, & '_saup2: Current B-norm of residual for factorization') end if c @@ -694,7 +694,7 @@ end if c if (msglvl .gt. 0) then - call ivout (logfil, 1, nconv, ndigit, + call ivout (logfil, 1, (/ nconv /), ndigit, & '_saup2: no. of "converged" Ritz values at this iter.') if (msglvl .gt. 1) then kp(1) = nev @@ -742,7 +742,7 @@ if (ishift .eq. 0) call scopy (np, workl, 1, ritz, 1) c if (msglvl .gt. 2) then - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_saup2: The number of shifts to apply ') call svout (logfil, np, workl, ndigit, & '_saup2: shifts selected') @@ -809,7 +809,7 @@ 130 continue c if (msglvl .gt. 2) then - call svout (logfil, 1, rnorm, ndigit, + call svout (logfil, 1, (/ rnorm /), ndigit, & '_saup2: B-norm of residual for NEV factorization') call svout (logfil, nev, h(1,2), ndigit, & '_saup2: main diagonal of compressed H matrix') Index: SRC/znaitr.f =================================================================== --- SRC/znaitr.f (revision 11543) +++ SRC/znaitr.f (working copy) @@ -380,7 +380,7 @@ if (msglvl .gt. 1) then call ivout (logfil, 1, j, ndigit, & '_naitr: generating Arnoldi vector number') - call dvout (logfil, 1, rnorm, ndigit, + call dvout (logfil, 1, (/ rnorm /), ndigit, & '_naitr: B-norm of the current residual is') end if c Index: SRC/znaup2.f =================================================================== --- SRC/znaup2.f (revision 11543) +++ SRC/znaup2.f (working copy) @@ -389,7 +389,7 @@ iter = iter + 1 c if (msglvl .gt. 0) then - call ivout (logfil, 1, iter, ndigit, + call ivout (logfil, 1, (/ iter /), ndigit, & '_naup2: **** Start of major iteration number ****') end if c @@ -402,9 +402,9 @@ np = kplusp - nev c if (msglvl .gt. 1) then - call ivout (logfil, 1, nev, ndigit, + call ivout (logfil, 1, (/ nev /), ndigit, & '_naup2: The length of the current Arnoldi factorization') - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_naup2: Extend the Arnoldi factorization by') end if c @@ -658,7 +658,7 @@ end if c if (msglvl .gt. 0) then - call ivout (logfil, 1, nconv, ndigit, + call ivout (logfil, 1, (/ nconv /), ndigit, & '_naup2: no. of "converged" Ritz values at this iter.') if (msglvl .gt. 1) then kp(1) = nev @@ -698,7 +698,7 @@ end if c if (msglvl .gt. 2) then - call ivout (logfil, 1, np, ndigit, + call ivout (logfil, 1, (/ np /), ndigit, & '_naup2: The number of shifts to apply ') call zvout (logfil, np, ritz, ndigit, & '_naup2: values of the shifts')