Local solver kernels about firedrake HOT 9 CLOSED

I also have some advection schemes where it would be great if we could assemble matrices from bilinear forms that contain operators equivalent to local solver kernels.

from firedrake.

This would also be a way to do hybridised systems.

from firedrake.

It's ok, I've figured out that we just need to assemble a block inverse to get a hybridised equation.

from firedrake.

@colinjcotter: What is the state of this? Is it implemented already?

from firedrake.

It is not, we are discussing to implement something like this to use
hybridisation. The main thing is that it needs to be able to operate on
mixed function spaces.

from firedrake.

My comment about on 23rd Dec 2014 only applies to the Helmholtz case, as it turns out, where pressure can be eliminated analytically. For Poisson, and more complicated things, it is necessary to assemble a local mixed function space block and invert it. We'll need to sketch out a plan for that.

from firedrake.

Coming as part of @thomasgibson's work.

from firedrake.

Related to #975

from firedrake.

Basic functionality is there, please raise specific slate issues.

from firedrake.