Fix is_cellwise_constant(ReferenceGrad(x))

[pbrubeck]

Adding a special case to detect that the reference gradient of a P1 field is cellwise constant.

[pbrubeck]

Currently, we cannot iterate grad on a Coefficient more than p+1 times (where p is the polynomial degree of the Coefficient), as grad(piecewise_constant) gives a domain-less Zero