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

Diff of /trunk/doc/examples/usersguide/heatedblock.py

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

revision 3346 by caltinay, Fri Nov 12 01:19:02 2010 UTC revision 3667 by gross, Tue Nov 15 23:29:23 2011 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.3346  
changed lines
  Added in v.3667
  ViewVC Help
Powered by ViewVC 1.1.26