Generate "fused" instationary operator
The conceptual separation of temporal and spatial operator is more of curse than a blessing in code generation (esp. once you try to go HPC). We should generate a fused operator. Marian has already worked on that and can be asked for details. Another question is how the weights from the time stepping scheme can be expressed in the DSL (or if not explicitly expressed can be represented in UFL expressions!).