Lines Matching refs:iKs
131 iKs = 0
201 iKs = iK
203 uBrg(:,:,iKs+1) = uBrg(:,:,iKs)
211 call Au1GMR ( EGmass, uBrg(:,:,iKs+1), ilwork, iBC,iper)
215 call bc3per (iBC, uBrg(:,:,iKs+1), iper, ilwork, nflow)
224 do jK = 1, iKs+1
228 temp = uBrg(:,:,iKs+1) * uBrg(:,:,1) ! {u_{i+1}*u_1} vector
235 uBrg(:,:,iKs+1) = uBrg(:,:,iKs+1) - beta * uBrg(:,:,jK-1)
237 temp = uBrg(:,:,iKs+1) * uBrg(:,:,jK) !{u_{i+1}*u_j} vector
242 HBrg(jK,iKs) = beta ! put this in the Hessenberg Matrix
252 HBrg(iKs+1,iKs) = unorm ! this fills the 1 sub diagonal band
256 uBrg(:,:,iKs+1) = uBrg(:,:,iKs+1) / unorm ! normalize the next Krylov
270 do jK = 1, iKs-1
271 tmp = Rcos(jK) * HBrg(jK, iKs) +
272 & Rsin(jK) * HBrg(jK+1,iKs)
273 HBrg(jK+1,iKs) = -Rsin(jK) * HBrg(jK, iKs) +
274 & Rcos(jK) * HBrg(jK+1,iKs)
275 HBrg(jK, iKs) = tmp
278 tmp = sqrt(HBrg(iKs,iKs)**2 + HBrg(iKs+1,iKs)**2)
279 Rcos(iKs) = HBrg(iKs, iKs) / tmp
280 Rsin(iKs) = HBrg(iKs+1,iKs) / tmp
281 HBrg(iKs, iKs) = tmp
282 HBrg(iKs+1,iKs) = zero
286 tmp = Rcos(iKs) * eBrg(iKs) + Rsin(iKs) * eBrg(iKs+1)
287 eBrg(iKs+1) =-Rsin(iKs) * eBrg(iKs) + Rcos(iKs) * eBrg(iKs+1)
288 eBrg(iKs) = tmp
293 echeck=abs(eBrg(iKs+1))
306 do jK = iKs, 1, -1
315 do jK = 1, iKs
326 echeck=abs(eBrg(iKs+1))
510 iKs = 0
578 iKs = iK
580 uBrg(:,:,iKs+1) = uBrg(:,:,iKs)
586 & uBrg(:,:,iKs+1) )
592 call bc3per (iBC, uBrg(:,:,iKs+1), iper, ilwork, nflow)
598 do jK = 1, iKs+1
602 temp = uBrg(:,:,iKs+1) * uBrg(:,:,1) ! {u_{i+1}*u_1} vector
609 uBrg(:,:,iKs+1)=uBrg(:,:,iKs+1)-beta * uBrg(:,:,jK-1)
611 temp = uBrg(:,:,iKs+1) * uBrg(:,:,jK) !{u_{i+1}*u_j} vector
616 HBrg(jK,iKs) = beta ! put this in the Hessenberg Matrix
627 HBrg(iKs+1,iKs) = unorm ! this fills the 1 sub diagonal band
631 uBrg(:,:,iKs+1) = uBrg(:,:,iKs+1) / unorm ! normalize the next Krylov
645 do jK = 1, iKs-1
646 tmp = Rcos(jK) * HBrg(jK, iKs) +
647 & Rsin(jK) * HBrg(jK+1,iKs)
648 HBrg(jK+1,iKs) = -Rsin(jK) * HBrg(jK, iKs) +
649 & Rcos(jK) * HBrg(jK+1,iKs)
650 HBrg(jK, iKs) = tmp
653 tmp = sqrt(HBrg(iKs,iKs)**2 + HBrg(iKs+1,iKs)**2)
654 Rcos(iKs) = HBrg(iKs, iKs) / tmp
655 Rsin(iKs) = HBrg(iKs+1,iKs) / tmp
656 HBrg(iKs, iKs)= tmp
657 HBrg(iKs+1,iKs)= zero
661 tmp = Rcos(iKs) * eBrg(iKs) + Rsin(iKs) * eBrg(iKs+1)
662 eBrg(iKs+1)=-Rsin(iKs) * eBrg(iKs) + Rcos(iKs) * eBrg(iKs+1)
663 eBrg(iKs) = tmp
668 echeck=abs(eBrg(iKs+1))
669 if (echeck .le. epsnrm.and. iKs .ge. minIters) exit
681 do jK = iKs, 1, -1
690 do jK = 1, iKs
700 echeck=abs(eBrg(iKs+1))