JuliaIBPM / ViscousFlow.jl

A framework for simulating viscous incompressible flows about arbitrary body shapes.
https://JuliaIBPM.github.io/ViscousFlow.jl/latest
Other
100 stars 28 forks source link

Saddle point solver mutates the rhs argument #12

Closed jdeldre closed 4 years ago

jdeldre commented 4 years ago

The saddle point solver accepts the right-hand side tuple as an input, but the ldiv! changes the constraint part of that tuple. This was not a big issue, since that rhs is not usually needed after the call to ldiv!. But it has come up in recent applications.