QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

(nefunkční) test besselových fcí 

Former-commit-id: 3d5ce98a7e4ff78f35395a27099ea41e77fd1584
This commit is contained in:
Marek Nečada 2018-03-13 16:47:03 +00:00
parent 5b5e77f985
commit e8ab2af46e
6 changed files with 50 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
tests/bessel_precision/besselDJcases.REMOVED.git-id Normal file
View File

@ -0,0 +1 @@
12584d0a62798e0d69d43fc030ba6fad056477e4

1
tests/bessel_precision/besselDYcases.REMOVED.git-id Normal file
View File

@ -0,0 +1 @@
12584d0a62798e0d69d43fc030ba6fad056477e4

1
tests/bessel_precision/besselJcases.REMOVED.git-id Normal file
View File

@ -0,0 +1 @@
a8fbb5a38d9893cec329d03ffe9d8426203b9043

1
tests/bessel_precision/besselYcases.REMOVED.git-id Normal file
View File

@ -0,0 +1 @@
5d2526ef87d9a1b5f5660fea4e9c5ccf1ce9394d

8
tests/bessel_precision/compileall.sh Executable file
View File

@ -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

38
tests/bessel_precision/sagebesselgen.py Normal file
View File

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