Wrong dimensions in LocalFunctionBasicGeometry
The dimensions used to evaluate the jacobian
in this line of code are wrong if mydimension != coorddimension
.
(I think) the code should be:
FieldMatrix< ctype, mydimension, coorddimension > gradPhi;
localFunction().jacobian( local, gradPhi );