probabilistic-numerics / probnum

Probabilistic Numerics in Python.
http://probnum.org
MIT License
435 stars 56 forks source link

Encoding prior knowledge for the solution of linear systems #526

Open JonathanWenger opened 3 years ago

JonathanWenger commented 3 years ago

Is your feature request related to a problem? Please describe. Currently it is not always clear how a user can supply prior knowledge to the probabilistic linear solver via a prior belief over the inverse / solution and matrix.

Give an example use case. Say the user has knowledge about the spectrum, smallest eigenvalue of A or a preconditioner. Currently it is not clear how this should be encoded into a prior belief.

# Known spectrum of A

# Prior belief
belief = LinearSystemBelief.from_spectrum(spectrum)
belief = LinearSystemBelief.from_preconditioner(spectrum)
belief = LinearSystemBelief.from_min_eigenvalue(spectrum)

Describe the solution you'd like. Provide specific prior belief classes with constructors from_* which cover the most common use cases of encoding prior knowledge.

mmahsereci commented 3 years ago

Thanks, Jonathan, I think this is an important issue. I have two questions: