b = np.dot(np.dot(tempp.T, f_mM[i:ntc - 1, n - j - 2::-1]),