--- trunk/finley/test/python/slip_stress.py	2006/11/08 03:26:21	891
+++ trunk/finley/test/python/slip_stress.py	2006/11/08 06:14:29	892
@@ -73,7 +73,8 @@
          return  self.__pde_u.getSolution()
 
       def solve_g(self,u,tol=1.e-8):
-         dp= (self.slip-jump(u))*lam_lmbd/FunctionOnContactZero(self.domain).getSize()**2
+         dp=-(self.slip-jump(u))*lam_lmbd/FunctionOnContactZero(self.domain).getSize()
+         print dp
          return  dp
 
 
@@ -97,6 +98,6 @@
 u0=Vector(0.,Solution(dom))
 p0=Vector(1.,FunctionOnContactZero(dom))
 prop.initialize(fixed_u_mask=mask,slip=Data(s,FunctionOnContactZero(dom)), density=rho,lmbd=lam_lmbd, mu=lam_mu)
-u,p=prop.solve(u0,p0,iter_max=50,tolerance=0.1,accepted_reduction=10.)
+u,p=prop.solve(u0,p0,iter_max=50,tolerance=0.1,accepted_reduction=0.99)
 saveVTK("dis.xml",u=u)
 saveVTK("fault.xml",sigma=p,s=jump(u))