specify modes for computation

[jwjeffr]

The fully vectorized calculation is $O(n^2)$ in memory. Implement some slower methods that don't take as much memory, and let the user specify which method to choose.

[jwjeffr]

serial method added in 19d57f8