search

JuliaLinearAlgebra / AlgebraicMultigrid.jl

Algebraic Multigrid in Julia
Other
116 stars 23 forks source link

Workaround for error in corner case #76

Closed omalled closed 2 years ago

omalled commented 3 years ago

I ran into a case where Pj as in

Px, Pj, Pp = rs_direct_interpolation_pass2(At, T, splitting, Pp)

was empty. This causes an error when maximum(Pj) is called here. I don't know enough about algebraic multigrid methods to know if this is the right way to fix this issue. Still, this change seems to work in my case and all the tests pass.

ViralBShah commented 2 years ago

Can you rebase to master?

ranjanan commented 2 years ago

I've created #83 which supersedes this, where I've rebased with master and added a test. Closing this in favour of that. Thanks for the PR @omalled!