Move lambdify_mpf and friends into syutil.
Refactor how the nodal basis inside of a triangle is obtained.
Avoid potential SymPy symbol collisions in the triangle basis.
Factor out the triangle orthonormal basis function.
Add support for tetrahedra.

Further testing, especially at higher orders, is still required.
Merge branch 'feature/tets' into develop.