Uniformisation of the weights of the Gauss Quadrature points (sum = volume of the reference element)
Adding Gauss Quadrature rule for a line
Modifications by Gaetan Bricteux

- fix bb computation to avoid model jumping around when adding points in z=0 plane
- pass interpol matrices as ref

move integration routines to Numeric/