implemented face sampling, not verifeid yet
fixed bug in compute face residual, was setting res_x instead of res_mu
fixed the way face residuals are computed