Enable the Jacobian of the domain constraint. What a useless default,
why on earth is it turned off?? Especially when the solves necessary
for the Jacobian are done either way!