Skip to content

Preconditioners do not work correctly #839

@termi-official

Description

@termi-official

Describe the example

Taking the example from the docs here and putting NaN into the preconditioner does not change the solver behavior at all.

Minimal Reproducible Example 👇

import LinearSolve as LS
import LinearAlgebra as LA

n = 25
A = rand(n, n)
b = rand(n)

precs = Returns((LS.InvPreconditioner(NaN*ones(n)), LA.Diagonal(NaN*ones(n))))

prob = LS.LinearProblem(A, b)
sol1 = LS.solve(prob, LS.KrylovJL_GMRES(verbose=1,precs));
sol2 = LS.solve(prob, LS.KrylovJL_GMRES(verbose=1));

Error & Stacktrace ⚠️

None. Which is the problem.

Not Working Environment (please complete the following information):

[7ed4a6bd] LinearSolve v3.48.0

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions