Moved model iteration out of optimal discretize

This will allow future magnetic functions to get the magnetization from
the tesseroids without having to modify the discretization code.
Also made some aesthetic changes to the Cython code.
Test if tesseroid obeys Laplace equation
Say that tesseroid vs shell test is accurate to 1%
Separate tesseroid gravity gradients into kernels