f_min = np.sum(f_mM, axis=0)