From e8ab2af46e1a82d05bace4b4d1a34b1734c33bdc Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= <marek@necada.org>
Date: Tue, 13 Mar 2018 16:47:03 +0000
Subject: [PATCH] =?UTF-8?q?(nefunk=C4=8Dn=C3=AD)=20test=20besselov=C3=BDch?=
 =?UTF-8?q?=20fc=C3=AD?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Former-commit-id: 3d5ce98a7e4ff78f35395a27099ea41e77fd1584
---
 .../besselDJcases.REMOVED.git-id              |  1 +
 .../besselDYcases.REMOVED.git-id              |  1 +
 .../besselJcases.REMOVED.git-id               |  1 +
 .../besselYcases.REMOVED.git-id               |  1 +
 tests/bessel_precision/compileall.sh          |  8 ++++
 tests/bessel_precision/sagebesselgen.py       | 38 +++++++++++++++++++
 6 files changed, 50 insertions(+)
 create mode 100644 tests/bessel_precision/besselDJcases.REMOVED.git-id
 create mode 100644 tests/bessel_precision/besselDYcases.REMOVED.git-id
 create mode 100644 tests/bessel_precision/besselJcases.REMOVED.git-id
 create mode 100644 tests/bessel_precision/besselYcases.REMOVED.git-id
 create mode 100755 tests/bessel_precision/compileall.sh
 create mode 100644 tests/bessel_precision/sagebesselgen.py

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)
+
+