diff --git a/tests/Makefile.am b/tests/Makefile.am
index a51b8eb82a17313818ff956ca3f15a232df8df65..d33853a8dc4bdd42fe8d93199009e36d0f9b5af6 100644
--- a/tests/Makefile.am
+++ b/tests/Makefile.am
@@ -22,7 +22,7 @@ AM_LDFLAGS = ../src/.libs/libswiftsim.a $(HDF5_LDFLAGS) $(HDF5_LIBS) $(FFTW_LIBS
 # List of programs and scripts to run in the test suite
 TESTS = testGreetings testMaths testReading.sh testSingle testKernel testSymmetry \
         testPair.sh testPairPerturbed.sh test27cells.sh test27cellsPerturbed.sh  \
-        testParser.sh testSPHStep test125cells.sh testFFT \
+        testParser.sh testSPHStep test125cells.sh test125cellsPerturbed.sh testFFT \
         testAdiabaticIndex testRiemannExact testRiemannTRRS testRiemannHLLC \
         testMatrixInversion testThreadpool testDump testLogger \
         testVoronoi1D testVoronoi2D testVoronoi3D
@@ -92,6 +92,6 @@ testLogger_SOURCES = testLogger.c
 # Files necessary for distribution
 EXTRA_DIST = testReading.sh makeInput.py testPair.sh testPairPerturbed.sh \
 	     test27cells.sh test27cellsPerturbed.sh testParser.sh \
-	     test125cells.sh testParserInput.yaml difffloat.py \
+	     test125cells.sh test125cells.sh testParserInput.yaml difffloat.py \
 	     tolerance_125.dat tolerance_27_normal.dat tolerance_27_perturbed.dat \
 	     tolerance_pair_normal.dat tolerance_pair_perturbed.dat
diff --git a/tests/test125cellsPerturbed.sh b/tests/test125cellsPerturbed.sh
new file mode 100755
index 0000000000000000000000000000000000000000..79005fe65559c2f699ecc84c4b0da155606d77c2
--- /dev/null
+++ b/tests/test125cellsPerturbed.sh
@@ -0,0 +1,24 @@
+#!/bin/bash
+for v in {0..3}
+do
+    for p in {0..2}
+    do
+	echo ""
+
+	rm -f brute_force_125_perturbed.dat swift_dopair_125_perturbed.dat
+
+  echo "./test125cells -n 6 -r 1 -d 0.1 -v $v -p $p -f perturbed"
+
+	./test125cells -n 6 -r 1 -d 0.1 -v $v -p $p -f perturbed
+
+	if [ -e brute_force_125_perturbed.dat ]
+	then
+	    python ./difffloat.py brute_force_125_perturbed.dat swift_dopair_125_perturbed.dat ./tolerance_125_perturbed.dat 6
+	else
+	    exit 1
+        fi
+
+    done
+done
+	
+exit $?
diff --git a/tests/test125cellsPerturbed.sh.in b/tests/test125cellsPerturbed.sh.in
new file mode 100755
index 0000000000000000000000000000000000000000..a59c77b7f649711f8dfbea4a2bc81a5e7767b0e7
--- /dev/null
+++ b/tests/test125cellsPerturbed.sh.in
@@ -0,0 +1,24 @@
+#!/bin/bash
+for v in {0..3}
+do
+    for p in {0..2}
+    do
+	echo ""
+
+	rm -f brute_force_125_perturbed.dat swift_dopair_125_perturbed.dat
+
+  echo "./test125cells -n 6 -r 1 -d 0.1 -v $v -p $p -f perturbed"
+
+	./test125cells -n 6 -r 1 -d 0.1 -v $v -p $p -f perturbed
+
+	if [ -e brute_force_125_perturbed.dat ]
+	then
+	    python @srcdir@/difffloat.py brute_force_125_perturbed.dat swift_dopair_125_perturbed.dat @srcdir@/tolerance_125_perturbed.dat 6
+	else
+	    exit 1
+        fi
+
+    done
+done
+	
+exit $?