Generate delayed Smolyak cubatures using C routines modified from
smolyak.c
found in Knut Petras' SMOLPACK.
smolyak.quad(d, k)smolyak.size(d, k)
smolyak.quad
returns a list object with the following
elements.
Quadrature nodes in rows of matrix.
Quadrature weights.
smolyak.size
returns an integer.
Dimension of unit cube.
Depth of algorithm.