Testing fused discretisation
============================

Runs both standard Gauss and 'fused' Gauss
