/[escript]/trunk/doc/examples/pyvisi/diffusion_with_pyvisi.py
ViewVC logotype

Diff of /trunk/doc/examples/pyvisi/diffusion_with_pyvisi.py

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

revision 1147 by ksteube, Wed May 16 06:39:11 2007 UTC revision 1148 by jongui, Wed May 16 22:45:33 2007 UTC
# Line 4  from esys.escript.linearPDEs import Line Line 4  from esys.escript.linearPDEs import Line
4  from esys.finley import Rectangle  from esys.finley import Rectangle
5  from esys.pyvisi import Scene, DataCollector, Map, Camera  from esys.pyvisi import Scene, DataCollector, Map, Camera
6  from esys.pyvisi.constant import *  from esys.pyvisi.constant import *
7    import os
8    
9  PYVISI_EXAMPLE_IMAGES_PATH = "data_sample_images/"  PYVISI_EXAMPLE_IMAGES_PATH = "data_sample_images"
10  X_SIZE = 400  X_SIZE = 400
11  Y_SIZE = 400  Y_SIZE = 400
12  JPG_RENDERER = Renderer.DISPLAY  JPG_RENDERER = Renderer.ONLINE_JPG
13    
14  #... set some parameters ...  #... set some parameters ...
15  xc=[0.02,0.002]  xc=[0.02,0.002]
# Line 37  T=Tref Line 38  T=Tref
38    
39  # Create a Scene.  # Create a Scene.
40  s = Scene(renderer = JPG_RENDERER, x_size = X_SIZE, y_size = Y_SIZE)  s = Scene(renderer = JPG_RENDERER, x_size = X_SIZE, y_size = Y_SIZE)
41    
42  # Create a DataCollector reading directly from escript objects.  # Create a DataCollector reading directly from escript objects.
43  dc = DataCollector(source = Source.ESCRIPT)  dc = DataCollector(source = Source.ESCRIPT)
44    
45  Map(scene = s, data_collector = dc,  # Create a Map.
46    m = Map(scene = s, data_collector = dc,
47          viewport = Viewport.SOUTH_WEST, lut = Lut.COLOR,          viewport = Viewport.SOUTH_WEST, lut = Lut.COLOR,
48          cell_to_point = False, outline = True)          cell_to_point = False, outline = True)
49    
50    # Create a Camera.
51    c = Camera(scene = s, viewport = Viewport.SOUTH_WEST)
52    
53  # ... start iteration:  # ... start iteration:
54  while t<0.4*100:  while t<0.4:
55        i+=1        i+=1
56        t+=h        t+=h
       print "time step :",t  
57        mypde.setValue(Y=qH+rhocp/h*T)        mypde.setValue(Y=qH+rhocp/h*T)
58        T=mypde.getSolution()        T=mypde.getSolution()
59    
60        dc.setData(temp = T)        dc.setData(temp = T)
       # Create a Map.  
       #Map(scene = s, data_collector = dc,  
               #viewport = Viewport.SOUTH_WEST, lut = Lut.COLOR,  
               #cell_to_point = False, outline = True)  
   
       # Create a Camera.  
       #c= Camera(scene = s, data_collector = dc,  
       #        viewport = Viewport.SOUTH_WEST)  
61                
62        # Render the object.        # Render the object.
63        s.render(image_name = PYVISI_EXAMPLE_IMAGES_PATH +          s.render(image_name = os.path.join(PYVISI_EXAMPLE_IMAGES_PATH,  
64                "diffusion%02d.jpg" % i)                "diffusion%02d.jpg") % i)
65    
Legend:
Removed from v.1147  
changed lines
  Added in v.1148
  ViewVC Help
Powered by ViewVC 1.1.26