void evaluate(
const ConstBasisDataSlice<ValueType>& basisData,
const ConstGeometricalDataSlice<CoordinateType>& geomData,
_1dSliceOf3dArray<ValueType>& result) const {
assert(basisData.componentCount() == 1);
assert(geomData.dimWorld() == 3);
// result := gradient of the shape function extended outside
// the surface so that its normal derivative on the surf. is zero
result(0) = basisData.derivatives(0, 0) * geomData.jacobianInverseTransposed(0, 0) +
basisData.derivatives(0, 1) * geomData.jacobianInverseTransposed(0, 1);
result(1) = basisData.derivatives(0, 0) * geomData.jacobianInverseTransposed(1, 0) +
basisData.derivatives(0, 1) * geomData.jacobianInverseTransposed(1, 1);
result(2) = basisData.derivatives(0, 0) * geomData.jacobianInverseTransposed(2, 0) +
basisData.derivatives(0, 1) * geomData.jacobianInverseTransposed(2, 1);
}
