Gz newtonraphson

Adds volume elements to sph-exa.

Marking as draft for the following reasons:

• CUDA code for volume elements is not implemented
• It would make sense to address issue #36 first
• Will break the Evrard test