search

JuliaLinearAlgebra / Preconditioners.jl

A few preconditioners for iterative solvers.
https://julialinearalgebra.github.io/Preconditioners.jl/
Other
50 stars 11 forks source link

Bug in AMGPreconditioner #20

Closed pmartorell closed 2 years ago

pmartorell commented 2 years ago

Hi,

since AlgebraicMultigrid v0.4.1 the code breaks in the following snipped:

using IterativeSolvers
using Preconditioners
using AlgebraicMultigrid

A = poisson(1000)
b = A*ones(1000)
p = AMGPreconditioner{SmoothedAggregation}(A)
cg(A,b,Pl=p)

However, the raw preconditioner can still be used from AlgebraicMultigrid package:

p = aspreconditioner(smoothed_aggregation(A))
mohamed82008 commented 2 years ago

Thanks for the report. I opened a PR to fix this.