made changes to leja quadrature rule to allow for more than one point to be added at a time. This works for interpolation but not for quadrature. The leja weights assume points are added one at a time. Note there also seems to be a problem with the sparse grid when num_dims L> 2. Add some unit tests for these situations