We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 5d5a9f2 commit de4a408Copy full SHA for de4a408
src/Poisson.jl
@@ -129,7 +129,7 @@ function pcg!(p::Poisson{T};it=6) where T
129
for i in 1:it
130
perBC!(ϵ,p.perdir)
131
@inside z[I] = mult(I,p.L,p.D,ϵ)
132
- alpha = rho/perdot(z,ϵ,perdir)
+ alpha = rho/perdot(z,ϵ,p.perdir)
133
(abs(alpha)<1e-2 || abs(alpha)>1e2) && return # alpha should be O(1)
134
@loop (x[I] += alpha*ϵ[I];
135
r[I] -= alpha*z[I]) over I ∈ inside(x)
0 commit comments