[NumLib] Add a restriction to the mapping the local dNdx to the global coordinate system

The discretised weak form of an inclined element can be computed numerically only in the local coordinate system \mathbf x defined in the element (note: it is not reference coordinate for the numerical integration). Therefore the rotation of \frac{\partial \mathbf N}{\mathbf x} to the global coordinate system \mathbf X is not needed. The only computations related to the coordinate rotation for an inclined element are 1) to project the global body force vector to the local coordinate system, and 2) to rotate the obtained the local velocity vector to the global system.

This MR. introduces a restriction in the computation of \frac{\partial \mathbf N}{\mathbf x} to avoid the rotation of \frac{\partial \mathbf N}{\mathbf x} to the global coordinate system, which is used by LIE process only. This restriction together with the rotation of \frac{\partial \mathbf N}{\mathbf x} to the global coordinate system can be removed once the velocity calculation in LIE process is changed.

This MR also adjusts the size of shape function matrix of element.

This MR is a part of MR !3632 (merged).

1. Feature description was added to the changelog
2. Tests covering your feature were added?
3. Any new feature or behavior change was documented?