Skip to content

Commit

Permalink
ga_diag_seq(): make all integer arguments 32-bit in the call to dsygv
Browse files Browse the repository at this point in the history
  • Loading branch information
edoapra committed Oct 30, 2023
1 parent adccda8 commit c833499
Showing 1 changed file with 6 additions and 3 deletions.
9 changes: 6 additions & 3 deletions global/src/ga_diag_seq.F
Original file line number Diff line number Diff line change
Expand Up @@ -28,11 +28,12 @@ subroutine gai_diag_seq(g_a, g_s, g_v, evals)
integer l_wrk
MA_ACCESS_INDEX_TYPE k_wrk
integer n2
INTGR4 n2_4
#endif
integer l_a, l_s
MA_ACCESS_INDEX_TYPE k_a, k_s
integer dim1, dim2, type, me
INTGR4 n4,ierr
INTGR4 n4,ierr,one4
logical status
c
c
Expand Down Expand Up @@ -94,8 +95,10 @@ subroutine gai_diag_seq(g_a, g_s, g_v, evals)
call rsg(n, n, dbl_mb(k_a), dbl_mb(k_s), evals, 1,
$ dbl_mb(k_v), dbl_mb(k_fv1), dbl_mb(k_fv2), ierr)
#else
call dsygv(1,'V','U',n4,dbl_mb(k_a),n,dbl_mb(k_s),n4,
$ evals,dbl_mb(k_wrk),n2, ierr)
one4=1
n2_4=n2
call dsygv(one4,'V','U',n4,dbl_mb(k_a),n4,dbl_mb(k_s),n4,
$ evals,dbl_mb(k_wrk),n2_4, ierr)
if (ierr.ne.0)
$ call ga_error('ga_diag_seq: dsygv failed',ierr)
c We used to copy to preserve code symmetry with EISPACK
Expand Down

0 comments on commit c833499

Please sign in to comment.