Mercurial > octave-nkf

comparison libcruft/arpack/src/sneupd.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
351 & invsub, iuptri, iwev , iwork(1), 351 & invsub, iuptri, iwev , iwork(1),
352 & j , k , ldh , ldq , 352 & j , k , ldh , ldq ,
353 & mode , msglvl, outncv, ritzr , 353 & mode , msglvl, outncv, ritzr ,
354 & ritzi , wri , wrr , irr , 354 & ritzi , wri , wrr , irr ,
355 & iri , ibd , ishift, numcnv , 355 & iri , ibd , ishift, numcnv ,
356 & np , jj 356 & np , jj , nconv2
357 logical reord 357 logical reord
358 Real 358 Real
359 & conds , rnorm, sep , temp, 359 & conds , rnorm, sep , temp,
360 & vl(1,1), temp1, eps23 360 & vl(1,1), temp1, eps23
361 c 361 c
659 call strsen('None' , 'V' , 659 call strsen('None' , 'V' ,
660 & select , ncv , 660 & select , ncv ,
661 & workl(iuptri), ldh , 661 & workl(iuptri), ldh ,
662 & workl(invsub), ldq , 662 & workl(invsub), ldq ,
663 & workl(iheigr), workl(iheigi), 663 & workl(iheigr), workl(iheigi),
664 & nconv , conds , 664 & nconv2 , conds ,
665 & sep , workl(ihbds) , 665 & sep , workl(ihbds) ,
666 & ncv , iwork , 666 & ncv , iwork ,
667 & 1 , ierr) 667 & 1 , ierr)
668 c 668 c
669 if (nconv2 .lt. nconv) then
670 nconv = nconv2
671 end if
672
669 if (ierr .eq. 1) then 673 if (ierr .eq. 1) then
670 info = 1 674 info = 1
671 go to 9000 675 go to 9000
672 end if 676 end if
673 c 677 c