function restrict(kf,kc)

fprintf('restrict residuals from level %d to level %d \n', kf,kc)
 