/[escript]/trunk/modellib/py_src/mechanics.py
ViewVC logotype

Diff of /trunk/modellib/py_src/mechanics.py

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

revision 825 by gross, Tue Aug 29 06:59:27 2006 UTC revision 826 by gross, Tue Aug 29 09:12:33 2006 UTC
# Line 248  class DruckerPrager(Mechanics): Line 248  class DruckerPrager(Mechanics):
248                 + (K-2./3*G)*k3Xk3 \                 + (K-2./3*G)*k3Xk3 \
249                 + sXk3-swap_axes(sXk3,1,3) \                 + sXk3-swap_axes(sXk3,1,3) \
250                 + 1./2*(swap_axes(sXk3,0,3)-swap_axes(sXk3,1,2) \                 + 1./2*(swap_axes(sXk3,0,3)-swap_axes(sXk3,1,2) \
251                        -swap_axes(sXk3,1,3)+swap_axes(sXk3,0,2))                        -swap_axes(sXk3,1,3)+swap_axes(sXk3,0,2)) \
252                 - chi/(h+G+alpha*beta*K)*outer(tmp+beta*K*k3,tmp+alpha*K*k3)\                 - chi/(h+G+alpha*beta*K)*outer(tmp+beta*K*k3,tmp+alpha*K*k3)\
253             # print S             # print S
254             return S             return S

Legend:
Removed from v.825  
changed lines
  Added in v.826

  ViewVC Help
Powered by ViewVC 1.1.26