I suspect that calling the density through LocalDensity::operator()
doesn't work yet, because the GFE function derivatives are not
transformed correctly from the reference to the grid element.
Therefore that new method still aborts with a deliberate error
message, until I have the time to check.