ENH: Add volume and discrete source space creation support

The 3% mismatch between the `src['interpolator']` matrices seems to be driven by slight differences in the transform reading in C and Python. Should be good for review and merge.