Skip to content

Conversation

@timholy
Copy link
Member

@timholy timholy commented Jan 9, 2026

When building a Woodbury or SymWoodbury matrix, we precompute

(C⁻¹ + V * (A\U))⁻¹

However, if the final W is rank-deficient, this inverse does not
exist. To handle this case, we now support an optional keyword argument
use_pinv::Bool=false to use the pseudo-inverse instead.

When building a `Woodbury` or `SymWoodbury` matrix, we precompute

    (C⁻¹ + V * (A\U))⁻¹

However, if the final `W` is rank-deficient, this inverse does not
exist. To handle this case, we now support an optional keyword argument
`use_pinv::Bool=false` to use the pseudo-inverse instead.
@timholy timholy merged commit 1faa137 into master Jan 9, 2026
2 checks passed
@timholy timholy deleted the teh/pinv branch January 9, 2026 16:07
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant