OK, I will try to discuss this with Onno on Thursday once he has finished being Examiner In Chief.

cheers
--cjc

On 10 October 2016 at 18:29, Mitchell, Lawrence <lawrence.mitchell@imperial.ac.uk> wrote:

> On 10 Oct 2016, at 18:19, Colin Cotter <colin.cotter@imperial.ac.uk> wrote:
>
> OK, this is the problem then - you need to find an appropriate preconditioner.
>
> You need to start by finding an appropriate preconditioner for the system without the R space part.
>
> Lawrence: once that is done, is there anything special to do to incorporate it into the system with the R space part [asking because you seem to have thought about this previously], or do we just recommend block preconditioning in this case?

I think you could use block first.  It's easier.  Otherwise I think the firedrake webpages sketch how to do sherman-morrison.

Lawrence
_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake



--
http://www.imperial.ac.uk/people/colin.cotter

www.cambridge.org/9781107663916