c = m_now[:, None] * c + (1.0 - m_now)[:, None] * c_pre