Lines Matching refs:iKss
861 iKss = 0
923 iKss = iK
925 uBrgt(:,iKss+1) = uBrgt(:,iKss)
929 call Au1GMRSclr ( EGmasst, uBrgt(:,iKss+1), ilwork, iper )
934 call bc3perSclr (iBC, uBrgt(:,iKss+1), iper)
938 do jK = 1, iKss+1
942 tempt = uBrgt(:,iKss+1) * uBrgt(:,1) ! {u_{i+1}*u_1} vector
949 uBrgt(:,iKss+1) = uBrgt(:,iKss+1) - beta * uBrgt(:,jK-1)
951 tempt = uBrgt(:,iKss+1) * uBrgt(:,jK) !{u_{i+1}*u_j} vector
956 HBrg(jK,iKss) = beta ! put this in the Hessenberg Matrix
966 HBrg(iKss+1,iKss) = unorm ! this fills the 1 sub diagonal band
970 uBrgt(:,iKss+1) = uBrgt(:,iKss+1) / unorm ! normalize the next Krylov
984 do jK = 1, iKss-1
985 tmp = Rcos(jK) * HBrg(jK, iKss) +
986 & Rsin(jK) * HBrg(jK+1,iKss)
987 HBrg(jK+1,iKss) = -Rsin(jK) * HBrg(jK, iKss) +
988 & Rcos(jK) * HBrg(jK+1,iKss)
989 HBrg(jK, iKss) = tmp
992 tmp = sqrt(HBrg(iKss,iKss)**2 + HBrg(iKss+1,iKss)**2)
993 Rcos(iKss) = HBrg(iKss, iKss) / tmp
994 Rsin(iKss) = HBrg(iKss+1,iKss) / tmp
995 HBrg(iKss, iKss) = tmp
996 HBrg(iKss+1,iKss) = zero
1000 tmp = Rcos(iKss)*eBrg(iKss) + Rsin(iKss)*eBrg(iKss+1)
1001 eBrg(iKss+1)=-Rsin(iKss)*eBrg(iKss) + Rcos(iKss)*eBrg(iKss+1)
1002 eBrg(iKss) = tmp
1006 ercheck=eBrg(iKss+1)
1008 if (abs(eBrg(iKss+1)) .le. epsnrm) exit
1020 do jK = iKss, 1, -1
1029 do jK = 1, iKss
1039 if (abs(eBrg(iKss+1)) .le. epsnrm) exit