/[escript]/branches/symbolic_from_3470/doc/examples/usersguide/heatedblock.py
ViewVC logotype

Diff of /branches/symbolic_from_3470/doc/examples/usersguide/heatedblock.py

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

revision 3788 by caltinay, Wed Sep 21 04:48:06 2011 UTC revision 3789 by caltinay, Tue Jan 31 04:55:05 2012 UTC
# Line 44  C=Tensor4(0.,Function(mydomain)) Line 44  C=Tensor4(0.,Function(mydomain))
44  for i in range(mydomain.getDim()):  for i in range(mydomain.getDim()):
45    for j in range(mydomain.getDim()):    for j in range(mydomain.getDim()):
46       C[i,i,j,j]+=lam       C[i,i,j,j]+=lam
47       C[j,i,j,i]+=mu       C[i,j,i,j]+=mu
48       C[j,i,i,j]+=mu       C[i,j,j,i]+=mu
49  msk=whereZero(x[0])*[1.,0.,0.] \  msk=whereZero(x[0])*[1.,0.,0.] \
50     +whereZero(x[1])*[0.,1.,0.] \     +whereZero(x[1])*[0.,1.,0.] \
51     +whereZero(x[2])*[0.,0.,1.]     +whereZero(x[2])*[0.,0.,1.]

Legend:
Removed from v.3788  
changed lines
  Added in v.3789

  ViewVC Help
Powered by ViewVC 1.1.26