From 82259b77f5e046021f06cd6f2be02ba678ae2f95 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= Date: Thu, 15 Dec 2016 08:02:47 +0200 Subject: [PATCH] Fix constant TMatrix bug in constmultipole rectarray Former-commit-id: b767b7178b2c0c7ed8c438f17470d86fda48bea4 --- qpms/qpms_p.py | 1 + 1 file changed, 1 insertion(+) diff --git a/qpms/qpms_p.py b/qpms/qpms_p.py index 4f4ab2a..511ef4a 100644 --- a/qpms/qpms_p.py +++ b/qpms/qpms_p.py @@ -1582,6 +1582,7 @@ def scatter_constmultipole_rectarray(omega, epsilon_b, xN, yN, xd, yd, TMatrices pq_0 = np.broadcast_to(pq_0, (N, 2, nelem)) MP_0 = np.empty((N,2,nelem),dtype=np.complex_) for j in range(N): # I wonder how this can be done without this loop... + xij, yij = xyind[j] MP_0[j] = np.tensordot(TMatrices[xij, yij],pq_0[j],axes=([-2,-1],[-2,-1])) MP_0.shape = (N*2*nelem,)