The wrong code is this line m_vCntGenes[m_ninner] = (X[m_vIDX_cl.nonzero(), m_vIDX_ge.nonzero() > 0).sum() A comprehensive bracket is missing here, and there is no echo between the right bracket and the left bracket (X[m_vIDX_cl.nonzero()