I think there is an error in this line:
|
copy!(bp,y_p); mul!(bp,C,du,1.0,-1.0) # bp = C*(A^-1 y_u) - y_p |
To my understanding it should be something like:
copy!(bp,y_p); mul!(bp,C,x_u,-1.0,1.0) # bp = y_p - C*(A^-1 y_u)
I was having wrong results in my application and this modification gives the good solution.
I think there is an error in this line:
GridapSolvers.jl/src/LinearSolvers/SchurComplementSolvers.jl
Line 65 in 0e62602
To my understanding it should be something like:
I was having wrong results in my application and this modification gives the good solution.