Remove broadcast in Krylov methods

[amontoison]

• SYMMLQ

src/symmlq.jl:354:   sprod .= sprod ./ sprod[(ix % window) + 1]

• TriLQR

src/trilqr.jl:167:   vₖ .= r₀ ./ βₖ              # v₁ = (b - Ax₀) / β₁
src/trilqr.jl:168:   uₖ .= s₀ ./ γₖ              # u₁ = (c - Aᴴy₀) / γ₁
src/trilqr.jl:333:   wₖ₋₁ .= vₖ₋₁ ./ conj(δₖ₋₁)
src/trilqr.jl:340:   wₖ₋₁ .= wₖ₋₁ ./ conj(δₖ₋₁)
src/trilqr.jl:348:   wₖ₋₁ .= wₖ₋₁ ./ conj(δₖ₋₁)
src/trilqr.jl:385:   vₖ .= q ./ βₖ₊₁  # βₖ₊₁vₖ₊₁ = q
src/trilqr.jl:388:   uₖ .= p ./ γₖ₊₁  # γₖ₊₁uₖ₊₁ = p

• BiLQ

src/bilq.jl:196:    vₖ .= r₀ ./ βₖ              # v₁ = (b - Ax₀) / β₁
src/bilq.jl:197:    uₖ .= c ./ conj(γₖ)         # u₁ = c / γ̄₁
src/bilq.jl:322:    vₖ .= q ./ βₖ₊₁        # βₖ₊₁vₖ₊₁ = q
src/bilq.jl:323:    uₖ .= p ./ conj(γₖ₊₁)  # γ̄ₖ₊₁uₖ₊₁ = p

• USYMLQ

src/usymlq.jl:173:    vₖ .= r₀ ./ βₖ              # v₁ = (b - Ax₀) / β₁
src/usymlq.jl:174:    uₖ .= c ./ γₖ               # u₁ = c / γ₁
src/usymlq.jl:289:    vₖ .= q ./ βₖ₊₁  # βₖ₊₁vₖ₊₁ = q
src/usymlq.jl:292:    uₖ .= p ./ γₖ₊₁  # γₖ₊₁uₖ₊₁ = p
src/usymqr.jl:177:    vₖ .= r₀ ./ βₖ               # v₁ = (b - Ax₀) / β₁
src/usymqr.jl:178:    uₖ .= c ./ γₖ                # u₁ = c / γ₁
src/usymqr.jl:266:    wₖ .= wₖ ./ δₖ
src/usymqr.jl:273:    wₖ .= wₖ ./ δₖ
src/usymqr.jl:281:    wₖ .= wₖ ./ δₖ
src/usymqr.jl:301:    vₖ .= q ./ βₖ₊₁  # βₖ₊₁vₖ₊₁ = q
src/usymqr.jl:304:    uₖ .= p ./ γₖ₊₁  # γₖ₊₁uₖ₊₁ = p

• DQGMRES

src/dqgmres.jl:186:    V[1] .= r₀ ./ rNorm
src/dqgmres.jl:232:    V[next_pos] .= w ./ Haux  # vₖ₊₁ = w / hₖ₊₁.ₖ
src/dqgmres.jl:273:    P[pos] .= P[pos] ./ H[1]

• TriCG

src/tricg.jl:352:        gx₂ₖ .= -conj(δₖ) .* gx₂ₖ₋₁
src/tricg.jl:359:        gx₂ₖ₋₁ .= conj(ηₖ) .* gx₂ₖ₋₁ .+ conj(λₖ) .* gx₂ₖ
src/tricg.jl:360:        gy₂ₖ₋₁ .= conj(ηₖ) .* gy₂ₖ₋₁ .+ conj(λₖ) .* gy₂ₖ
src/tricg.jl:362:        gx₂ₖ .= vₖ .- conj(σₖ) .* gx₂ₖ
src/tricg.jl:363:        gy₂ₖ .=    .- conj(σₖ) .* gy₂ₖ
src/tricg.jl:365:        gx₂ₖ₋₁ .=    .- gx₂ₖ₋₁ .- conj(δₖ) .* gx₂ₖ
src/tricg.jl:366:        gy₂ₖ₋₁ .= uₖ .- gy₂ₖ₋₁ .- conj(δₖ) .* gy₂ₖ

• GPMR

src/gpmr.jl:243:    V[1] .= b₀ ./ β
src/gpmr.jl:248:    U[1] .= c₀ ./ γ
src/gpmr.jl:461:    V[k+1] .= q ./ Haux  # hₖ₊₁.ₖvₖ₊₁ = q
src/gpmr.jl:469:    U[k+1] .= p ./ Faux  # fₖ₊₁.ₖuₖ₊₁ = p

• QMR

src/qmr.jl:202:    vₖ .= r₀ ./ βₖ               # v₁ = (b - Ax₀) / β₁
src/qmr.jl:203:    uₖ .= c ./ conj(γₖ)          # u₁ = c / γ̄₁
src/qmr.jl:308:    wₖ .= wₖ ./ δₖ
src/qmr.jl:315:    wₖ .= wₖ ./ δₖ
src/qmr.jl:323:    wₖ .= wₖ ./ δₖ
src/qmr.jl:335:    vₖ .= q ./ βₖ₊₁        # βₖ₊₁vₖ₊₁ = q
src/qmr.jl:336:    uₖ .= p ./ conj(γₖ₊₁)  # γ̄ₖ₊₁uₖ₊₁ = p

• BiLQR

src/bilqr.jl:180:   vₖ .= r₀ ./ βₖ              # v₁ = (b - Ax₀) / β₁
src/bilqr.jl:181:   uₖ .= s₀ ./ conj(γₖ)        # u₁ = (c - Aᴴy₀) / γ̄₁
src/bilqr.jl:357:   wₖ₋₁ .= uₖ₋₁ ./ conj(δₖ₋₁)
src/bilqr.jl:364:   wₖ₋₁ .= wₖ₋₁ ./ conj(δₖ₋₁)
src/bilqr.jl:372:   wₖ₋₁ .= wₖ₋₁ ./ conj(δₖ₋₁)
src/bilqr.jl:407:   vₖ .= q ./ βₖ₊₁        # βₖ₊₁vₖ₊₁ = q
src/bilqr.jl:408:   uₖ .= p ./ conj(γₖ₊₁)  # γ̄ₖ₊₁uₖ₊₁ = p

• TriMR

src/trimr.jl:429:        gx₂ₖ₋₁ .= vₖ ./ δ₂ₖ₋₁
src/trimr.jl:430:        gx₂ₖ .= -(σ₂ₖ₋₁ / δ₂ₖ) .* gx₂ₖ₋₁
src/trimr.jl:431:        gy₂ₖ .= uₖ ./ δ₂ₖ
src/trimr.jl:440:        gx₂ₖ₋₁ .= (vₖ .- η₂ₖ₋₃ .* gx₂ₖ₋₃ .- σ₂ₖ₋₂ .* gx₂ₖ₋₂                   ) ./ δ₂ₖ₋₁
src/trimr.jl:441:        gx₂ₖ   .= (   .- λ₂ₖ₋₃ .* gx₂ₖ₋₃ .- η₂ₖ₋₂ .* gx₂ₖ₋₂ .- σ₂ₖ₋₁ .* gx₂ₖ₋₁) ./ δ₂ₖ
src/trimr.jl:442:        gy₂ₖ₋₁ .= (   .- η₂ₖ₋₃ .* gy₂ₖ₋₃ .- σ₂ₖ₋₂ .* gy₂ₖ₋₂                   ) ./ δ₂ₖ₋₁
src/trimr.jl:443:        gy₂ₖ   .= (uₖ .- λ₂ₖ₋₃ .* gy₂ₖ₋₃ .- η₂ₖ₋₂ .* gy₂ₖ₋₂ .- σ₂ₖ₋₁ .* gy₂ₖ₋₁) ./ δ₂ₖ
src/trimr.jl:448:        g₂ₖ₋₁ .= (vₖ .- μ₂ₖ₋₅ .* g₂ₖ₋₅ .- λ₂ₖ₋₄ .* g₂ₖ₋₄ .- η₂ₖ₋₃ .* g₂ₖ₋₃ .- σ₂ₖ₋₂ .* g₂ₖ₋₂                  ) ./ δ₂ₖ₋₁
src/trimr.jl:449:        g₂ₖ   .= (                     .- μ₂ₖ₋₄ .* g₂ₖ₋₄ .- λ₂ₖ₋₃ .* g₂ₖ₋₃ .- η₂ₖ₋₂ .* g₂ₖ₋₂ .- σ₂ₖ₋₁ .* g₂ₖ₋₁) ./ δ₂ₖ
src/trimr.jl:455:        g₂ₖ₋₁ .= (     .- μ₂ₖ₋₅ .* g₂ₖ₋₅ .- λ₂ₖ₋₄ .* g₂ₖ₋₄ .- η₂ₖ₋₃ .* g₂ₖ₋₃ .- σ₂ₖ₋₂ .* g₂ₖ₋₂                  ) ./ δ₂ₖ₋₁
src/trimr.jl:456:        g₂ₖ   .= (uₖ                     .- μ₂ₖ₋₄ .* g₂ₖ₋₄ .- λ₂ₖ₋₃ .* g₂ₖ₋₃ .- η₂ₖ₋₂ .* g₂ₖ₋₂ .- σ₂ₖ₋₁ .* g₂ₖ₋₁) ./ δ₂ₖ

• GMRES

src/gmres.jl:217:      V[1] .= r₀ ./ rNorm
src/gmres.jl:307:      V[inner_iter+1] .= q ./ Hbis  # hₖ₊₁.ₖvₖ₊₁ = q

• FGMRES

src/fgmres.jl:223:      V[1] .= r₀ ./ rNorm
src/fgmres.jl:314:      V[inner_iter+1] .= q ./ Hbis  # hₖ₊₁.ₖvₖ₊₁ = q

• DIOM

src/diom.jl:186:    V[1] .= r₀ ./ rNorm
src/diom.jl:230:    V[next_pos] .= w ./ Haux  # vₖ₊₁ = w / hₖ₊₁.ₖ
src/diom.jl:274:    P[ppos] .= P[ppos] ./ H[1]

• FOM

src/fom.jl:215:      V[1] .= r₀ ./ rNorm
src/fom.jl:296:      V[inner_iter+1] .= q ./ Hbis  # hₖ₊₁.ₖvₖ₊₁ = q