def test_run_short(Nmax=NMAX, Max_iter=1000, tol=1):
    Phi = electrostatics.calculate_potential(
        Nmax=Nmax, Max_iter=Max_iter, tol=tol)
    assert_almost_equal(Phi.min(), -100)
    assert_almost_equal(Phi.max(), 100)
    assert_almost_equal(np.trace(Phi), 539.21017867137766,
                        decimal=0)
def test_run_long(Nmax=NMAX, Max_iter=30000, tol=TOL):
    Phi = electrostatics.calculate_potential(
        Nmax=Nmax, Max_iter=Max_iter, tol=tol)
    assert_almost_equal(Phi.min(), -100)
    assert_almost_equal(Phi.max(), 100)
    assert_almost_equal(np.trace(Phi), 535.78614641692775,
                        decimal=5)
 def setUp(self, Nmax=NMAX, Max_iter=30000, tol=TOL):
     self.Phi = electrostatics.calculate_potential(
         Nmax=Nmax, Max_iter=Max_iter, tol=tol)
     self.Delta = 1.0