[firedrake] operations on matrices

[Colin Cotter]

On 07/01/15 09:36, Lawrence Mitchell wrote:

> So one hacky way to do this is as follows (but, it's hacky, there are
> probably better ways):
>
> 1. Generate the kernel for (1)
>
> 2. Generate the kernel that takes the output of (1), an nxm element
> tensor and reads field values to modify the element tensor.
>
> 3. Construct a "super-kernel" that takes all the argument from (1) and
> the additional field arguments from (2) which calls first 1., then 2.
> to generate the correctly modified element tensor.
>
> 4. Use the kernel from 3. as your assembly kernel in a par_loop
> (rather than the kernel from 1.)

Thanks. So, how do I do that in practice, given that I want to write the 
kernel for 2 in a c string?

cheers
--cjc