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

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

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

revision 3320 by jfenwick, Thu Jan 28 02:03:15 2010 UTC revision 3321 by gross, Thu Oct 28 05:02:41 2010 UTC
# Line 49  l40=Line(p4,p0) Line 49  l40=Line(p4,p0)
49  l37=Line(p3,p7)  l37=Line(p3,p7)
50  l62=Line(p6,p2)  l62=Line(p6,p2)
51    
52  bottom=PlaneSurface(CurveLoop(l01,l13,l32,l20))  bottom=PlaneSurface(-CurveLoop(l01,l13,l32,l20))
53  top=PlaneSurface(CurveLoop(l45,l57,l76,l64))  top=PlaneSurface(CurveLoop(l45,l57,l76,l64))
54  front=PlaneSurface(CurveLoop(l01,l15,-l45,l40))  front=PlaneSurface(CurveLoop(l01,l15,-l45,l40))
55  back=PlaneSurface(CurveLoop(l32,-l62,-l76,-l37))  back=PlaneSurface(CurveLoop(l32,-l62,-l76,-l37))
56  left=PlaneSurface(CurveLoop(-l40,-l64,l62,l20))  left=PlaneSurface(CurveLoop(-l40,-l64,l62,l20))
57  right=PlaneSurface(CurveLoop(-l15,l13,l37,-l57))  right=PlaneSurface(CurveLoop(-l15,l13,l37,-l57))
58  v=Volume(SurfaceLoop(top,-bottom,front,back,left,right))  v=Volume(SurfaceLoop(top,bottom,front,back,left,right))
59    
60  des=Design(dim=3, order=2, element_size = 0.1, keep_files=True)  des=Design(dim=3, order=2, element_size = 0.1, keep_files=True)
61  des.setScriptFileName("brick.geo")  des.setScriptFileName("brick.geo")
62  des.addItems(v, top, bottom, back, front, left , right)  des.addItems(v, top, bottom,back, front, left , right)
63    
64  dom=MakeDomain(des)  dom=MakeDomain(des)
65  dom.write("brick.fly")  dom.write("brick.fly")

Legend:
Removed from v.3320  
changed lines
  Added in v.3321

  ViewVC Help
Powered by ViewVC 1.1.26