Add Monte-Carlo integration method and tests
Added matrix implementation, more tests for integrator, and a factory for generating current density functions for magnetic field calculations. Currently only supports arbitrary torus.
Added tests for current density function
Added tests for DoubleVector and DoubleMatrix. Also added more functionality to DoubleMatrix.