add simple `solve` function
For the simplest use we want to just mumps.solve(a, b)
. Currently the alternative is
ctx = mumps.MUMPSContext()
ctx.factor(a)
ctx.solve(b)
This is a pretty trivial wrapper to write, but IMO makes the interface quite a bit more ergonomic for the basic usecase.
Edited by Joseph Weston