Testing fused discretisation ============================ Runs both standard Gauss and 'fused' Gauss