f_max = np.sum(f_mM, axis=1)