unpropagated change

NonLinearSolver/contact/contactTerms.cpp

@@ -84,7 +84,7 @@ template<> void defoDefoContactBilinearTermByPerturbation<double>::get(contactEl
     disp(j)+=_pert;
     for(int i=0; i<nbdofs; i++)
     {
-      m(i,j)=(Fplus(i)-Fminus(i))/(-2.*_pert); //minus because stiffness related toof internal force
+      m(i,j)=(Fplus(i)-Fminus(i))/(2.*_pert); //
     }
   }
   defoDefoTerm->get(cE,npts,GP,Fminus,disp);// put back the right _nodesInContact