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

Jdu do práce 

Former-commit-id: ea635a016207ad637a0647127dea650ea650d6d2
This commit is contained in:
Marek Nečada 2018-03-14 05:46:46 +00:00
parent e8ab2af46e
commit 3b407ab40a
2 changed files with 29 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
tests/bessel_precision/compileall.sh
View File

@ -1,8 +1,8 @@
#!/bin/bash #!/bin/bash
c99 -ggdb -o jtest -DJTEST besseltest.c c99 -ggdb -o jtest -DJTEST besseltest.c -lgsl -lblas -lm
c99 -ggdb -o ytest -DYTEST besseltest.c c99 -ggdb -o ytest -DYTEST besseltest.c -lgsl -lblas -lm
c99 -ggdb -o djtest -DDJTEST besseltest.c c99 -ggdb -o djtest -DDJTEST besseltest.c -lgsl -lblas -lm
c99 -ggdb -o dytest -DDYTEST besseltest.c c99 -ggdb -o dytest -DDYTEST besseltest.c -lgsl -lblas -lm
c99 -ggdb -o djtest_steed -DJTEST_STEED besseltest.c c99 -ggdb -o djtest_steed -DJTEST_STEED besseltest.c -lgsl -lblas -lm

26
tests/bessel_precision/sagebesselgen.py
View File

@ -31,8 +31,30 @@ def printbesselDJrow(lMax, x, file=sys.stdout):
def printbesselDYrow(lMax, x, file=sys.stdout): def printbesselDYrow(lMax, x, file=sys.stdout):
print(lMax, N(x), end=' ', file=file); print(lMax, N(x), end=' ', file=file);
for l in range(lMax+1): for l in range(lMax+1):
print(N(-spherical_bessel_J(l+1,x) print(N(-spherical_bessel_Y(l+1,x)
+ l/x*spherical_bessel_J(l,x)), end = ' ', file=file) + l/x*spherical_bessel_Y(l,x)), end = ' ', file=file)
print('', file=file) print('', file=file)
def genall(lMax):
f = open('besselDJcases', 'w')
for o in IntegerRange(1,100):
printbesselDJrow(lMax, o, file=f)
printbesselDJrow(lMax, 1/o, file=f)
printbesselDJrow(lMax, o/sqrt(3), file=f)
f = open('besselDYcases', 'w')
for o in IntegerRange(1,100):
printbesselDYrow(lMax, o, file=f)
printbesselDYrow(lMax, 1/o, file=f)
printbesselDYrow(lMax, o/sqrt(3), file=f)
f = open('besselJcases', 'w')
for o in IntegerRange(1,100):
printbesselJrow(lMax, o, file=f)
printbesselJrow(lMax, 1/o, file=f)
printbesselJrow(lMax, o/sqrt(3), file=f)
f = open('besselYcases', 'w')
for o in IntegerRange(1,100):
printbesselYrow(lMax, o, file=f)
printbesselYrow(lMax, 1/o, file=f)
printbesselYrow(lMax, o/sqrt(3), file=f)