DiagonalNonConstantOperator

Hi,

the merge will add a new operator to the master. So far we could only generate diagonal operators, i.e. ones that only work with the diagonal of the associated matrix respectively the central stencil weight in the case of constant stencils. The new operator extends this to the case of non-constant stencils.

Cheers
Marcus