Mercurial > octave-nkf

comparison libcruft/arpack/src/zneupd.f @ 12277:3d38b4916cbf release-3-4-x

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
avoid memory overrun in ARPACK
author David Bateman <dbateman@free.fr>
date Fri, 28 Jan 2011 14:05:34 -0500
parents 9f5d2ef078e8
children
comparison
equal deleted inserted replaced
12276:e4ed10d59f18 12277:3d38b4916cbf
299 character type*6 299 character type*6
300 integer bounds, ierr , ih , ihbds, iheig , nconv , 300 integer bounds, ierr , ih , ihbds, iheig , nconv ,
301 & invsub, iuptri, iwev , j , ldh , ldq , 301 & invsub, iuptri, iwev , j , ldh , ldq ,
302 & mode , msglvl, ritz , wr , k , irz , 302 & mode , msglvl, ritz , wr , k , irz ,
303 & ibd , outncv, iq , np , numcnv, jj , 303 & ibd , outncv, iq , np , numcnv, jj ,
304 & ishift 304 & ishift, nconv2
305 Complex*16 305 Complex*16
306 & rnorm, temp, vl(1) 306 & rnorm, temp, vl(1)
307 Double precision 307 Double precision
308 & conds, sep, rtemp, eps23 308 & conds, sep, rtemp, eps23
309 logical reord 309 logical reord
590 c %-----------------------------------------------% 590 c %-----------------------------------------------%
591 c 591 c
592 call ztrsen ('None' , 'V' , select , 592 call ztrsen ('None' , 'V' , select ,
593 & ncv , workl(iuptri), ldh , 593 & ncv , workl(iuptri), ldh ,
594 & workl(invsub), ldq , workl(iheig), 594 & workl(invsub), ldq , workl(iheig),
595 & nconv , conds , sep , 595 & nconv2 , conds , sep ,
596 & workev , ncv , ierr) 596 & workev , ncv , ierr)
597 c 597 c
598 if (nconv2 .lt. nconv) then
599 nconv = nconv2
600 end if
601
598 if (ierr .eq. 1) then 602 if (ierr .eq. 1) then
599 info = 1 603 info = 1
600 go to 9000 604 go to 9000
601 end if 605 end if
602 c 606 c