Ticket #38571: arpack.patch
| File arpack.patch, 2.7 KB (added by wimmer@…, 11 years ago) |
|---|
-
scipy/sparse/linalg/eigen/arpack/ARPACK/SRC/dseupd.f
subroutine dseupd (rvec , howmny, select, d , 760 760 & ldq , workl(iw+ncv), workl(ihb), 761 761 & ncv , temp , ierr) 762 762 c 763 c %-----------------------------------------------------% 764 c | Make a copy of the last row into | 765 c | workl(iw+ncv:iw+2*ncv), as it is needed again in | 766 c | the Ritz vector purification step below | 767 c %-----------------------------------------------------% 768 c 769 do 67 j = 1, nconv 770 workl(iw+ncv+j-1) = workl(ihb+j-1) 771 67 continue 772 763 773 else if (rvec .and. howmny .eq. 'S') then 764 774 c 765 775 c Not yet implemented. See remark 2 above. … … subroutine dseupd (rvec , howmny, select, d , 830 840 if (rvec .and. (type .eq. 'SHIFTI' .or. type .eq. 'CAYLEY')) then 831 841 c 832 842 do 110 k=0, nconv-1 833 workl(iw+k) = workl(i q+k*ldq+ncv-1)843 workl(iw+k) = workl(iw+ncv+k) 834 844 & / workl(iw+k) 835 845 110 continue 836 846 c 837 847 else if (rvec .and. type .eq. 'BUCKLE') then 838 848 c 839 849 do 120 k=0, nconv-1 840 workl(iw+k) = workl(i q+k*ldq+ncv-1)850 workl(iw+k) = workl(iw+ncv+k) 841 851 & / (workl(iw+k)-one) 842 852 120 continue 843 853 c -
scipy/sparse/linalg/eigen/arpack/ARPACK/SRC/sseupd.f
subroutine sseupd(rvec , howmny, select, d , 760 760 & ldq , workl(iw+ncv), workl(ihb), 761 761 & ncv , temp , ierr) 762 762 c 763 c %-----------------------------------------------------% 764 c | Make a copy of the last row into | 765 c | workl(iw+ncv:iw+2*ncv), as it is needed again in | 766 c | the Ritz vector purification step below | 767 c %-----------------------------------------------------% 768 c 769 do 67 j = 1, nconv 770 workl(iw+ncv+j-1) = workl(ihb+j-1) 771 67 continue 772 763 773 else if (rvec .and. howmny .eq. 'S') then 764 774 c 765 775 c Not yet implemented. See remark 2 above. … … subroutine sseupd(rvec , howmny, select, d , 830 840 if (rvec .and. (type .eq. 'SHIFTI' .or. type .eq. 'CAYLEY')) then 831 841 c 832 842 do 110 k=0, nconv-1 833 workl(iw+k) = workl(i q+k*ldq+ncv-1)843 workl(iw+k) = workl(iw+ncv+k) 834 844 & / workl(iw+k) 835 845 110 continue 836 846 c 837 847 else if (rvec .and. type .eq. 'BUCKLE') then 838 848 c 839 849 do 120 k=0, nconv-1 840 workl(iw+k) = workl(i q+k*ldq+ncv-1)850 workl(iw+k) = workl(iw+ncv+k) 841 851 & / (workl(iw+k)-one) 842 852 120 continue 843 853 c