Lines Matching refs:iKs
152 iKs = 0
222 iKs = iK
226 temp = uBrg(:,:,iKs)
236 uBrg(:,:,iKs+1) = temp ! u_{i+1}= J u_i In Johan Thesis p 15c
247 do jK = 1, iKs+1
251 temp = uBrg(:,:,iKs+1) * uBrg(:,:,1) ! {u_{i+1}*u_1} vector
258 uBrg(:,:,iKs+1) = uBrg(:,:,iKs+1) - beta * uBrg(:,:,jK-1)
260 temp = uBrg(:,:,iKs+1) * uBrg(:,:,jK) !{u_{i+1}*u_j} vector
265 HBrg(jK,iKs) = beta ! put this in the Hessenberg Matrix
275 HBrg(iKs+1,iKs) = unorm ! this fills the 1 sub diagonal band
279 uBrg(:,:,iKs+1) = uBrg(:,:,iKs+1) / unorm ! normalize the next Krylov
293 do jK = 1, iKs-1
294 tmp = Rcos(jK) * HBrg(jK, iKs) +
295 & Rsin(jK) * HBrg(jK+1,iKs)
296 HBrg(jK+1,iKs) = -Rsin(jK) * HBrg(jK, iKs) +
297 & Rcos(jK) * HBrg(jK+1,iKs)
298 HBrg(jK, iKs) = tmp
301 tmp = sqrt(HBrg(iKs,iKs)**2 + HBrg(iKs+1,iKs)**2)
302 Rcos(iKs) = HBrg(iKs, iKs) / tmp
303 Rsin(iKs) = HBrg(iKs+1,iKs) / tmp
304 HBrg(iKs, iKs) = tmp
305 HBrg(iKs+1,iKs) = zero
309 tmp = Rcos(iKs) * eBrg(iKs) + Rsin(iKs) * eBrg(iKs+1)
310 eBrg(iKs+1) = -Rsin(iKs) * eBrg(iKs) + Rcos(iKs) * eBrg(iKs+1)
311 eBrg(iKs) = tmp
316 echeck=abs(eBrg(iKs+1))
317 if (echeck .le. epsnrm.and. iKs .ge. minIters) exit
329 do jK = iKs, 1, -1
338 do jK = 1, iKs
348 if (abs(eBrg(iKs+1)) .le. epsnrm) exit