Fix cubeoid tetgen mesh input

It was generating extra nodes for the first boundary for some reason!

Needed small tweak to fpdiff tol in hlib/matrix_generation
Add a global nodal quadrature which should use J=1

Not sure if we need some other stuff to do this though, might not work
yet.