Non-allocating version of `calcnodes`

Currently `SymCubatures.calcnodes` allocates its own output array.  When `PumiInterface` gets the coordinates of all the elements, it would be faster to have a method that allows the user to supply the output array.
