fix segfault with rectangular matrixes

ntessore · Oct 28, 2021 · c3c3106 · c3c3106
1 parent 3151b2c
commit c3c3106
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/convolvecl.pyx b/convolvecl.pyx
@@ -3,7 +3,7 @@
 # author: Nicolas Tessore <[email protected]>
 # license: MIT
 #
-# cython: language_level=3, boundscheck=False
+# cython: language_level=3, boundscheck=False, embedsignature=True
 # distutils: extra_compile_args=['-fopenmp', '-Ofast']
 # distutils: extra_link_args=['-fopenmp']
 '''
@@ -185,7 +185,7 @@ def mixmat(cl2, lmax=None, l1max=None, l2max=None, spins1=(0, 0),
                 for i in range(n):
                     mll1 = mll1 + twol2p1cl2_[l2min_int+i]*thr[i]*thr_[i]
                 m_[l, l1] = (2*l1 + 1)/FOUR_PI*mll1
-                if symmetric:
+                if symmetric and l1 <= lmax_ and l <= l1max_:
                     m_[l1, l] = (2*l + 1)/FOUR_PI*mll1
 
         # free local buffers

diff --git a/setup.cfg b/setup.cfg
@@ -1,6 +1,6 @@
 [metadata]
 name = convolvecl
-version = 2021.10.21
+version = 2021.10.28
 maintainer = Nicolas Tessore
 maintainer_email = [email protected]
 description = convolution of angular power spectra