1 May
2019
1 May
'19
11:43 a.m.
On Mon, Apr 29, 2019 at 4:30 PM Cotter, Colin J <colin.cotter@imperial.ac.uk> wrote: Dear Firedrakers, Greetings from Montreal. I'm thinking about a double splitting method to solve the system
(M + D + L)u = f,
which as a preconditioner works out as
(M + D)^{-1}M(M + L)^{-1} r,
where r is a residual.
In the absence of a better idea, I'll implement this through a Python preconditioner. But is there a smarter way to do this?
I think writing a Python PC is the obvious thing to do, given that you don't have a natural algebraic splitting of M, D, and L. Lawrence