Error in SchurComplementSolvers.jl

[oriolcg]

I think there is an error in this line:

GridapSolvers.jl/src/LinearSolvers/SchurComplementSolvers.jl

Line 65 in 0e62602

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.

[JordiManyer]

Yes, that looks like a bug. I'll have a look. In any case that solver is due a rather big revision soon.