/[escript]/trunk/finley/test/python/axisymm-splitB.py
ViewVC logotype

Diff of /trunk/finley/test/python/axisymm-splitB.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2468 by gross, Thu Jun 11 02:55:51 2009 UTC revision 2534 by caltinay, Thu Jul 16 06:49:19 2009 UTC
# Line 102  while istep < nstep: Line 102  while istep < nstep:
102                             X=-dt*r*(dev_stress-teta3*p*kronecker(dom)),                             X=-dt*r*(dev_stress-teta3*p*kronecker(dom)),
103                             y=sigma_d*face_mask*r_b)                             y=sigma_d*face_mask*r_b)
104      U_star=momentumStep1.getSolution()      U_star=momentumStep1.getSolution()
105      saveVTK("u.xml",u=U_star,u0=U)      saveVTK("u.vtu",u=U_star,u0=U)
106      #      #
107      #  step 2:      #  step 2:
108      #      #
# Line 156  while istep < nstep: Line 156  while istep < nstep:
156      #      #
157      dom.setX(dom.getX()+U*dt)      dom.setX(dom.getX()+U*dt)
158      t=t+dt      t=t+dt
159      if (istep-1)%w_step==0:saveVTK("u.%d.xml"%((istep-1)/w_step),p=p,eta=eta_d,U=U_star,U_star=U_star,gamma=gamma)      if (istep-1)%w_step==0:saveVTK("u.%d.vtu"%((istep-1)/w_step),p=p,eta=eta_d,U=U_star,U_star=U_star,gamma=gamma)
160      if istep == 3: 1/0      if istep == 3: 1/0

Legend:
Removed from v.2468  
changed lines
  Added in v.2534

  ViewVC Help
Powered by ViewVC 1.1.26