re-initialize force to zero at the begin of apply --> fix most of the parallelism problems, but there are still some discrepencies in the solution depending on the number of proc used. These discrepencies are due to the last term of the surface integrals, involving the deformation tensor. Moreover, the last contribution of these surface integrals has been corrected (some terms where missing)