Implement and test pseudocode in report
Add testSpd methods for debugging slow convergence.

Tests seem to pass, but will leave them in there for now.
Flatten the GBC 3D P-wave image.
Make images of conical faults in cylindrical coordinates.
Experimenting