-
Notifications
You must be signed in to change notification settings - Fork 32
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Saddle Operator #277
Comments
Looking forward to it! Thanks! More generally, |
@dpo It was nice to meet you at Seattle. |
It was great to (finally) meet you as well @lrsantos11 ! Now that I think of it, |
@lrsantos11 The following already works m, n = 5, 3
A = rand(m, n)
M = rand(m, m); M = M + M'
N = rand(n, n); N = N + N'
opA = LinearOperator(A)
opM = LinearOperator(M)
opN = LinearOperator(N)
K = [opM opA ;
opA' opN] but it doesn’t let us specify that |
Given operators
opA
,opM
andopN
, creates saddle point linear operator.I'll contribute on this, using
BlockDiagonalOperator
as template.LinearOperators.jl/src/special-operators.jl
Lines 229 to 240 in 9025458
I already have implemented outside LO. Its a matter of just including, testing and PRing.
The text was updated successfully, but these errors were encountered: