diff --git a/tests/bessel_precision/besselDJcases.REMOVED.git-id b/tests/bessel_precision/besselDJcases.REMOVED.git-id
new file mode 100644
index 0000000..05d7222
--- /dev/null
+++ b/tests/bessel_precision/besselDJcases.REMOVED.git-id
@@ -0,0 +1 @@
+12584d0a62798e0d69d43fc030ba6fad056477e4
\ No newline at end of file
diff --git a/tests/bessel_precision/besselDYcases.REMOVED.git-id b/tests/bessel_precision/besselDYcases.REMOVED.git-id
new file mode 100644
index 0000000..05d7222
--- /dev/null
+++ b/tests/bessel_precision/besselDYcases.REMOVED.git-id
@@ -0,0 +1 @@
+12584d0a62798e0d69d43fc030ba6fad056477e4
\ No newline at end of file
diff --git a/tests/bessel_precision/besselJcases.REMOVED.git-id b/tests/bessel_precision/besselJcases.REMOVED.git-id
new file mode 100644
index 0000000..9b87a7e
--- /dev/null
+++ b/tests/bessel_precision/besselJcases.REMOVED.git-id
@@ -0,0 +1 @@
+a8fbb5a38d9893cec329d03ffe9d8426203b9043
\ No newline at end of file
diff --git a/tests/bessel_precision/besselYcases.REMOVED.git-id b/tests/bessel_precision/besselYcases.REMOVED.git-id
new file mode 100644
index 0000000..03b21e5
--- /dev/null
+++ b/tests/bessel_precision/besselYcases.REMOVED.git-id
@@ -0,0 +1 @@
+5d2526ef87d9a1b5f5660fea4e9c5ccf1ce9394d
\ No newline at end of file
diff --git a/tests/bessel_precision/compileall.sh b/tests/bessel_precision/compileall.sh
new file mode 100755
index 0000000..50bf761
--- /dev/null
+++ b/tests/bessel_precision/compileall.sh
@@ -0,0 +1,8 @@
+#!/bin/bash
+
+c99 -ggdb -o jtest -DJTEST besseltest.c
+c99 -ggdb -o ytest -DYTEST besseltest.c
+c99 -ggdb -o djtest -DDJTEST besseltest.c
+c99 -ggdb -o dytest -DDYTEST besseltest.c
+c99 -ggdb -o djtest_steed -DJTEST_STEED besseltest.c
+
diff --git a/tests/bessel_precision/sagebesselgen.py b/tests/bessel_precision/sagebesselgen.py
new file mode 100644
index 0000000..6e90aaa
--- /dev/null
+++ b/tests/bessel_precision/sagebesselgen.py
@@ -0,0 +1,38 @@
+from __future__ import print_function
+
+def printhankelrow(lMax, x, file=sys.stdout):
+    print(lMax, N(x), end=' ', file=file);
+    for l in range(lMax+1):
+        print(N(spherical_hankel1(l,x)), end = ' ', file=file)
+    print('', file=file)
+
+
+def printbesselJrow(lMax, x, file=sys.stdout):
+    print(lMax, N(x), end=' ', file=file);
+    for l in range(lMax+1):
+        print(N(spherical_bessel_J(l,x)), end = ' ', file=file)
+    print('', file=file)
+
+
+def printbesselYrow(lMax, x, file=sys.stdout):
+    print(lMax, N(x), end=' ', file=file);
+    for l in range(lMax+1):
+        print(N(spherical_bessel_Y(l,x)), end = ' ', file=file)
+    print('', file=file)
+
+def printbesselDJrow(lMax, x, file=sys.stdout):
+    print(lMax, N(x), end=' ', file=file);
+    for l in range(lMax+1):
+        print(N(-spherical_bessel_J(l+1,x)
+            + l/x*spherical_bessel_J(l,x)), end = ' ', file=file)
+    print('', file=file)
+
+
+def printbesselDYrow(lMax, x, file=sys.stdout):
+    print(lMax, N(x), end=' ', file=file);
+    for l in range(lMax+1):
+        print(N(-spherical_bessel_J(l+1,x)
+            + l/x*spherical_bessel_J(l,x)), end = ' ', file=file)
+    print('', file=file)
+
+