From 7ae729b1b1777d658ee354526109faa4102717f4 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= <marek.necada@aalto.fi>
Date: Tue, 26 Jul 2016 07:41:36 +0300
Subject: [PATCH] Print lattice size to stderr along with the progress info,
 doc nmax limit

Former-commit-id: 646f95f462565702ff15ca66df0294a0568a60f8
---
 qpms/qpms_p.py | 11 +++++++++--
 1 file changed, 9 insertions(+), 2 deletions(-)

diff --git a/qpms/qpms_p.py b/qpms/qpms_p.py
index 4747be3..758accf 100644
--- a/qpms/qpms_p.py
+++ b/qpms/qpms_p.py
@@ -470,9 +470,15 @@ def Ã(m,n,μ,ν,kdlj,θlj,φlj,r_ge_d,J):
     J: 1, 2, 3 or 4
         Type of the wave in the old center.
 
-    Returns:
-    --------
+    Returns
+    -------
     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) 
                 +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):
         timec = time.time()
         print('%.4f: running scatter_plane_wave_rectarray' % timec, file = sys.stderr)
+        print('xN = %d, yN = %d' % (xN, yN), file = sys.stderr)
         sys.stderr.flush()
     nelem = TMatrices.shape[-1]
     if ((nelem != TMatrices.shape[-3]) or (2 != TMatrices.shape[-2]) or (2 != TMatrices.shape[-4])):