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

Print lattice size to stderr along with the progress info, doc nmax 

limit


Former-commit-id: 646f95f462565702ff15ca66df0294a0568a60f8
This commit is contained in:
Marek Nečada 2016-07-26 07:41:36 +03:00
parent 3171c07c3e
commit 7ae729b1b1
1 changed files with 9 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
qpms/qpms_p.py
View File

@ -470,9 +470,15 @@ def Ã(m,n,μ,ν,kdlj,θlj,φlj,r_ge_d,J):
J: 1, 2, 3 or 4 J: 1, 2, 3 or 4
Type of the wave in the old center. Type of the wave in the old center.
Returns: Returns
-------- -------
TODO TODO
Bugs
----
gevero's gaunt coefficient implementation fails for large m, n (the unsafe territory
is somewhere around -72, 80)
""" """
exponent=(math.lgamma(2*n+1)-math.lgamma(n+2)+math.lgamma(2*ν+3)-math.lgamma(ν+2) exponent=(math.lgamma(2*n+1)-math.lgamma(n+2)+math.lgamma(2*ν+3)-math.lgamma(ν+2)
+math.lgamma(n+ν+m-μ+1)-math.lgamma(n-m+1)-math.lgamma(ν+μ+1) +math.lgamma(n+ν+m-μ+1)-math.lgamma(n-m+1)-math.lgamma(ν+μ+1)
@ -1121,6 +1127,7 @@ def scatter_plane_wave_rectarray(omega, epsilon_b, xN, yN, xd, yd, TMatrices, k_
if (watch_time): if (watch_time):
timec = time.time() timec = time.time()
print('%.4f: running scatter_plane_wave_rectarray' % timec, file = sys.stderr) print('%.4f: running scatter_plane_wave_rectarray' % timec, file = sys.stderr)
print('xN = %d, yN = %d' % (xN, yN), file = sys.stderr)
sys.stderr.flush() sys.stderr.flush()
nelem = TMatrices.shape[-1] nelem = TMatrices.shape[-1]
if ((nelem != TMatrices.shape[-3]) or (2 != TMatrices.shape[-2]) or (2 != TMatrices.shape[-4])): if ((nelem != TMatrices.shape[-3]) or (2 != TMatrices.shape[-2]) or (2 != TMatrices.shape[-4])):