Update tutorials and improve tutorial testing framework

doc/tutorials/05-raspberry_electrophoresis/05-raspberry_electrophoresis.ipynb

@@ -500,7 +500,8 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "Efield = np.array((0.1, 0, 0)) # an electric field of 0.1 is the upper limit of the linear response regime for this model\n",
+    "E = 0.1 # an electric field of 0.1 is the upper limit of the linear response regime for this model\n",
+    "Efield = np.array([E, 0, 0])\n",
     "for p in system.part:\n",
     "    p.ext_force = p.q * Efield"
    ]

testsuite/scripts/tutorials/test_05-raspberry_electrophoresis.py

@@ -17,15 +17,26 @@
 
 import unittest as ut
 import importlib_wrapper
+import numpy as np
 
 tutorial, skipIfMissingFeatures = importlib_wrapper.configure_and_import(
     "@TUTORIALS_DIR@/05-raspberry_electrophoresis/05-raspberry_electrophoresis.py",
-    gpu=True, num_iterations=10, num_steps_per_iteration=10)
+    gpu=True, box_l=20., E=0.2, num_iterations=150, num_steps_per_iteration=400)
 
 
 @skipIfMissingFeatures
 class Tutorial(ut.TestCase):
     system = tutorial.system
+    def test_trajectory(self):
+        trajectory = np.loadtxt('posVsTime.dat')[:,1:4]
+        # the raspberry should have traveled mostly on the x-axis
+        dist = np.abs(trajectory[-1, :] - trajectory[0, :])
+        self.assertGreater(dist[0], 2 * dist[1])
+        self.assertGreater(dist[0], 2 * dist[2])
+        # the velocity should be highest on the x-axis
+        vel = np.abs(np.mean(trajectory[1:, :] - trajectory[:-1, :], axis=0))
+        self.assertGreater(vel[0], 2 * vel[1])
+        self.assertGreater(vel[0], 2 * vel[2])
 
 
 if __name__ == "__main__":