OK, then we don't need to discuss this at the moment but keep this in our mind for future when needed. I'll remove this part from the pull request. (I think, the linear mapping is valid only for constant gradient of shape functions, e.g. linear triangles and tetrahedrals. For other element types, we need NR approach or something else. )