function subtract(kf,kc)

fprintf('restrict and subtract solution on level %d from the solution on level %d \n',...
    kf, kc)
