Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
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
5 changes: 4 additions & 1 deletion CHANGES
Original file line number Diff line number Diff line change
@@ -1,5 +1,8 @@
arpack-ng - 3.9.0


[ John Doe ]
* Avoid calling [c|z]dotc for better portability on macOS

[ Dima Pasechnik ]
* [BUG FIX] autotools: replace obsolete AC_TRY_COMPILE macros.
* Support for NAG's nagfor Fortran compiler
Expand Down
8 changes: 4 additions & 4 deletions PARPACK/SRC/BLACS/pcgetv0.f
Original file line number Diff line number Diff line change
Expand Up @@ -197,8 +197,8 @@ subroutine pcgetv0
Real
& pscnorm2, slapy2
Complex
& cdotc
external cdotc, pscnorm2, slapy2
& ccdotc
external ccdotc, pscnorm2, slapy2
c
c %-----------------%
c | Data Statements |
Expand Down Expand Up @@ -335,7 +335,7 @@ subroutine pcgetv0
c
first = .FALSE.
if (bmat .eq. 'G') then
cnorm = cdotc (n, resid, 1, workd, 1)
cnorm = ccdotc (n, resid, 1, workd, 1)
call cgsum2d( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
rnorm0 = sqrt(slapy2(real (cnorm),aimag(cnorm)))
else if (bmat .eq. 'I') then
Expand Down Expand Up @@ -394,7 +394,7 @@ subroutine pcgetv0
end if
c
if (bmat .eq. 'G') then
cnorm = cdotc (n, resid, 1, workd, 1)
cnorm = ccdotc (n, resid, 1, workd, 1)
call cgsum2d( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
rnorm = sqrt(slapy2(real (cnorm),aimag(cnorm)))
else if (bmat .eq. 'I') then
Expand Down
10 changes: 5 additions & 5 deletions PARPACK/SRC/BLACS/pcnaitr.f
Original file line number Diff line number Diff line change
Expand Up @@ -303,10 +303,10 @@ subroutine pcnaitr
c %--------------------%
c
Complex
& cdotc
& ccdotc
Real
& pslamch, pscnorm2, clanhs, slapy2
external cdotc, pscnorm2, clanhs, pslamch, slapy2
external ccdotc, pscnorm2, clanhs, pslamch, slapy2
c
c %---------------------%
c | Intrinsic Functions |
Expand Down Expand Up @@ -573,7 +573,7 @@ subroutine pcnaitr
c %-------------------------------------%
c
if (bmat .eq. 'G') then
cnorm = cdotc (n, resid, 1, workd(ipj), 1)
cnorm = ccdotc (n, resid, 1, workd(ipj), 1)
call cgsum2d( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
wnorm = sqrt( slapy2(real(cnorm),aimag(cnorm)) )
else if (bmat .eq. 'I') then
Expand Down Expand Up @@ -647,7 +647,7 @@ subroutine pcnaitr
c %------------------------------%
c
if (bmat .eq. 'G') then
cnorm = cdotc (n, resid, 1, workd(ipj), 1)
cnorm = ccdotc (n, resid, 1, workd(ipj), 1)
call cgsum2d( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
rnorm = sqrt( slapy2(real(cnorm),aimag(cnorm)) )
else if (bmat .eq. 'I') then
Expand Down Expand Up @@ -749,7 +749,7 @@ subroutine pcnaitr
c %-----------------------------------------------------%
c
if (bmat .eq. 'G') then
cnorm = cdotc (n, resid, 1, workd(ipj), 1)
cnorm = ccdotc (n, resid, 1, workd(ipj), 1)
call cgsum2d( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
rnorm1 = sqrt( slapy2(real(cnorm),aimag(cnorm)) )
else if (bmat .eq. 'I') then
Expand Down
6 changes: 3 additions & 3 deletions PARPACK/SRC/BLACS/pcnaup2.f
Original file line number Diff line number Diff line change
Expand Up @@ -254,10 +254,10 @@ subroutine pcnaup2
c %--------------------%
c
Complex
& cdotc
& ccdotc
Real
& pscnorm2, pslamch, slapy2
external cdotc, pscnorm2, pslamch, slapy2
external ccdotc, pscnorm2, pslamch, slapy2
c
c %---------------------%
c | Intrinsic Functions |
Expand Down Expand Up @@ -767,7 +767,7 @@ subroutine pcnaup2
end if
c
if (bmat .eq. 'G') then
cmpnorm = cdotc (n, resid, 1, workd, 1)
cmpnorm = ccdotc (n, resid, 1, workd, 1)
call cgsum2d( comm, 'All', ' ', 1, 1, cmpnorm, 1, -1, -1 )
rnorm = sqrt(slapy2(real(cmpnorm),aimag(cmpnorm)))
else if (bmat .eq. 'I') then
Expand Down
6 changes: 3 additions & 3 deletions PARPACK/SRC/BLACS/pcneupd.f
Original file line number Diff line number Diff line change
Expand Up @@ -340,8 +340,8 @@ subroutine pcneupd
external scnrm2,pslamch,slapy2
c
Complex
& cdotc
external cdotc
& ccdotc
external ccdotc
c
c %---------------------%
c | Intrinsic Functions |
Expand Down Expand Up @@ -743,7 +743,7 @@ subroutine pcneupd
c | inner product can be set to j. |
c %------------------------------------------%
c
workev(j) = cdotc(j, workl(ihbds), 1,
workev(j) = ccdotc(j, workl(ihbds), 1,
& workl(invsub+(j-1)*ldq), 1)
40 continue
c
Expand Down
8 changes: 4 additions & 4 deletions PARPACK/SRC/BLACS/pzgetv0.f
Original file line number Diff line number Diff line change
Expand Up @@ -197,8 +197,8 @@ subroutine pzgetv0
Double precision
& pdznorm2 , dlapy2
Complex*16
& zdotc
external zdotc , pdznorm2 , dlapy2
& zzdotc
external zzdotc , pdznorm2 , dlapy2
c
c %-----------------%
c | Data Statements |
Expand Down Expand Up @@ -335,7 +335,7 @@ subroutine pzgetv0
c
first = .FALSE.
if (bmat .eq. 'G') then
cnorm = zdotc (n, resid, 1, workd, 1)
cnorm = zzdotc (n, resid, 1, workd, 1)
call zgsum2d ( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
rnorm0 = sqrt(dlapy2 (dble (cnorm),dimag (cnorm)))
else if (bmat .eq. 'I') then
Expand Down Expand Up @@ -394,7 +394,7 @@ subroutine pzgetv0
end if
c
if (bmat .eq. 'G') then
cnorm = zdotc (n, resid, 1, workd, 1)
cnorm = zzdotc (n, resid, 1, workd, 1)
call zgsum2d ( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
rnorm = sqrt(dlapy2 (dble (cnorm),dimag (cnorm)))
else if (bmat .eq. 'I') then
Expand Down
10 changes: 5 additions & 5 deletions PARPACK/SRC/BLACS/pznaitr.f
Original file line number Diff line number Diff line change
Expand Up @@ -303,10 +303,10 @@ subroutine pznaitr
c %--------------------%
c
Complex*16
& zdotc
& zzdotc
Double precision
& pdlamch, pdznorm2, zlanhs, dlapy2
external zdotc, pdznorm2, zlanhs, pdlamch, dlapy2
external zzdotc, pdznorm2, zlanhs, pdlamch, dlapy2
c
c %---------------------%
c | Intrinsic Functions |
Expand Down Expand Up @@ -573,7 +573,7 @@ subroutine pznaitr
c %-------------------------------------%
c
if (bmat .eq. 'G') then
cnorm = zdotc (n, resid, 1, workd(ipj), 1)
cnorm = zzdotc (n, resid, 1, workd(ipj), 1)
call zgsum2d( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
wnorm = sqrt( dlapy2(dble(cnorm),dimag(cnorm)) )
else if (bmat .eq. 'I') then
Expand Down Expand Up @@ -647,7 +647,7 @@ subroutine pznaitr
c %------------------------------%
c
if (bmat .eq. 'G') then
cnorm = zdotc (n, resid, 1, workd(ipj), 1)
cnorm = zzdotc (n, resid, 1, workd(ipj), 1)
call zgsum2d( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
rnorm = sqrt( dlapy2(dble(cnorm),dimag(cnorm)) )
else if (bmat .eq. 'I') then
Expand Down Expand Up @@ -749,7 +749,7 @@ subroutine pznaitr
c %-----------------------------------------------------%
c
if (bmat .eq. 'G') then
cnorm = zdotc (n, resid, 1, workd(ipj), 1)
cnorm = zzdotc (n, resid, 1, workd(ipj), 1)
call zgsum2d( comm, 'All', ' ', 1, 1, cnorm, 1, -1, -1 )
rnorm1 = sqrt( dlapy2(dble(cnorm),dimag(cnorm)) )
else if (bmat .eq. 'I') then
Expand Down
6 changes: 3 additions & 3 deletions PARPACK/SRC/BLACS/pznaup2.f
Original file line number Diff line number Diff line change
Expand Up @@ -254,10 +254,10 @@ subroutine pznaup2
c %--------------------%
c
Complex*16
& zdotc
& zzdotc
Double precision
& pdznorm2, pdlamch, dlapy2
external zdotc, pdznorm2, pdlamch, dlapy2
external zzdotc, pdznorm2, pdlamch, dlapy2
c
c %---------------------%
c | Intrinsic Functions |
Expand Down Expand Up @@ -767,7 +767,7 @@ subroutine pznaup2
end if
c
if (bmat .eq. 'G') then
cmpnorm = zdotc (n, resid, 1, workd, 1)
cmpnorm = zzdotc (n, resid, 1, workd, 1)
call zgsum2d( comm, 'All', ' ', 1, 1, cmpnorm, 1, -1, -1 )
rnorm = sqrt(dlapy2(dble(cmpnorm),dimag(cmpnorm)))
else if (bmat .eq. 'I') then
Expand Down
6 changes: 3 additions & 3 deletions PARPACK/SRC/BLACS/pzneupd.f
Original file line number Diff line number Diff line change
Expand Up @@ -340,8 +340,8 @@ subroutine pzneupd
external dznrm2,pdlamch,dlapy2
c
Complex*16
& zdotc
external zdotc
& zzdotc
external zzdotc
c
c %---------------------%
c | Intrinsic Functions |
Expand Down Expand Up @@ -743,7 +743,7 @@ subroutine pzneupd
c | inner product can be set to j. |
c %------------------------------------------%
c
workev(j) = zdotc(j, workl(ihbds), 1,
workev(j) = zzdotc(j, workl(ihbds), 1,
& workl(invsub+(j-1)*ldq), 1)
40 continue
c
Expand Down
8 changes: 4 additions & 4 deletions PARPACK/SRC/MPI/pcgetv0.f
Original file line number Diff line number Diff line change
Expand Up @@ -200,8 +200,8 @@ subroutine pcgetv0
Real
& pscnorm2, slapy2
Complex
& cdotc
external cdotc, pscnorm2, slapy2
& ccdotc
external ccdotc, pscnorm2, slapy2
c
c %-----------------%
c | Data Statements |
Expand Down Expand Up @@ -331,7 +331,7 @@ subroutine pcgetv0
c
first = .FALSE.
if (bmat .eq. 'G') then
cnorm_buf = cdotc (n, resid, 1, workd, 1)
cnorm_buf = ccdotc (n, resid, 1, workd, 1)
call MPI_ALLREDUCE( [cnorm_buf], buf2, 1,
& MPI_COMPLEX, MPI_SUM, comm, ierr )
cnorm = buf2(1)
Expand Down Expand Up @@ -393,7 +393,7 @@ subroutine pcgetv0
end if
c
if (bmat .eq. 'G') then
cnorm_buf = cdotc (n, resid, 1, workd, 1)
cnorm_buf = ccdotc (n, resid, 1, workd, 1)
call MPI_ALLREDUCE( [cnorm_buf], buf2, 1,
& MPI_COMPLEX, MPI_SUM, comm, ierr )
cnorm = buf2(1)
Expand Down
10 changes: 5 additions & 5 deletions PARPACK/SRC/MPI/pcnaitr.f
Original file line number Diff line number Diff line change
Expand Up @@ -307,10 +307,10 @@ subroutine pcnaitr
c %--------------------%
c
Complex
& cdotc
& ccdotc
Real
& pslamch10, pscnorm2, clanhs, slapy2
external cdotc, pscnorm2, clanhs, pslamch10, slapy2
external ccdotc, pscnorm2, clanhs, pslamch10, slapy2
c
c %---------------------%
c | Intrinsic Functions |
Expand Down Expand Up @@ -576,7 +576,7 @@ subroutine pcnaitr
c %-------------------------------------%
c
if (bmat .eq. 'G') then
cnorm_buf = cdotc (n, resid, 1, workd(ipj), 1)
cnorm_buf = ccdotc (n, resid, 1, workd(ipj), 1)
call MPI_ALLREDUCE( [cnorm_buf], buf2, 1,
& MPI_COMPLEX, MPI_SUM, comm, ierr )
cnorm = buf2(1)
Expand Down Expand Up @@ -653,7 +653,7 @@ subroutine pcnaitr
c %------------------------------%
c
if (bmat .eq. 'G') then
cnorm_buf = cdotc (n, resid, 1, workd(ipj), 1)
cnorm_buf = ccdotc (n, resid, 1, workd(ipj), 1)
call MPI_ALLREDUCE( [cnorm_buf], buf2, 1,
& MPI_COMPLEX, MPI_SUM, comm, ierr )
cnorm = buf2(1)
Expand Down Expand Up @@ -758,7 +758,7 @@ subroutine pcnaitr
c %-----------------------------------------------------%
c
if (bmat .eq. 'G') then
cnorm_buf = cdotc (n, resid, 1, workd(ipj), 1)
cnorm_buf = ccdotc (n, resid, 1, workd(ipj), 1)
call MPI_ALLREDUCE( [cnorm_buf], buf2, 1,
& MPI_COMPLEX, MPI_SUM, comm, ierr )
cnorm = buf2(1)
Expand Down
6 changes: 3 additions & 3 deletions PARPACK/SRC/MPI/pcnaup2.f
Original file line number Diff line number Diff line change
Expand Up @@ -257,10 +257,10 @@ subroutine pcnaup2
c %--------------------%
c
Complex
& cdotc
& ccdotc
Real
& pscnorm2, pslamch10, slapy2
external cdotc, pscnorm2, pslamch10, slapy2
external ccdotc, pscnorm2, pslamch10, slapy2
c
c %---------------------%
c | Intrinsic Functions |
Expand Down Expand Up @@ -770,7 +770,7 @@ subroutine pcnaup2
end if
c
if (bmat .eq. 'G') then
cmpnorm_buf = cdotc (n, resid, 1, workd, 1)
cmpnorm_buf = ccdotc (n, resid, 1, workd, 1)
call MPI_ALLREDUCE( [cmpnorm_buf], buf2, 1,
& MPI_COMPLEX, MPI_SUM, comm, ierr )
cmpnorm = buf2(1)
Expand Down
6 changes: 3 additions & 3 deletions PARPACK/SRC/MPI/pcneupd.f
Original file line number Diff line number Diff line change
Expand Up @@ -340,8 +340,8 @@ subroutine pcneupd
external scnrm2,pslamch10,slapy2
c
Complex
& cdotc
external cdotc
& ccdotc
external ccdotc
c
c %---------------------%
c | Intrinsic Functions |
Expand Down Expand Up @@ -743,7 +743,7 @@ subroutine pcneupd
c | inner product can be set to j. |
c %------------------------------------------%
c
workev(j) = cdotc(j, workl(ihbds), 1,
workev(j) = ccdotc(j, workl(ihbds), 1,
& workl(invsub+(j-1)*ldq), 1)
40 continue
c
Expand Down
8 changes: 4 additions & 4 deletions PARPACK/SRC/MPI/pzgetv0.f
Original file line number Diff line number Diff line change
Expand Up @@ -200,8 +200,8 @@ subroutine pzgetv0
Double precision
& pdznorm2 , dlapy2
Complex*16
& zdotc
external zdotc , pdznorm2 , dlapy2
& zzdotc
external zzdotc , pdznorm2 , dlapy2
c
c %-----------------%
c | Data Statements |
Expand Down Expand Up @@ -331,7 +331,7 @@ subroutine pzgetv0
c
first = .FALSE.
if (bmat .eq. 'G') then
cnorm_buf = zdotc (n, resid, 1, workd, 1)
cnorm_buf = zzdotc (n, resid, 1, workd, 1)
call MPI_ALLREDUCE( [cnorm_buf], buf2, 1,
& MPI_DOUBLE_COMPLEX , MPI_SUM, comm, ierr )
cnorm = buf2(1)
Expand Down Expand Up @@ -393,7 +393,7 @@ subroutine pzgetv0
end if
c
if (bmat .eq. 'G') then
cnorm_buf = zdotc (n, resid, 1, workd, 1)
cnorm_buf = zzdotc (n, resid, 1, workd, 1)
call MPI_ALLREDUCE( [cnorm_buf], buf2, 1,
& MPI_DOUBLE_COMPLEX , MPI_SUM, comm, ierr )
cnorm = buf2(1)
Expand Down
Loading