Changeset 5847 for anuga_core/source/pypar-numeric/mandelbrot_example

Timestamp:

Oct 18, 2008, 9:13:18 PM (16 years ago)

Author:

steve

Message:

Changed parallel_api so that global mesh only needs to
be constructed on processor 0

Location:

anuga_core/source/pypar-numeric/mandelbrot_example

Files:

• compile.py (modified) (1 diff)
• mandel_sequential.py (modified) (1 diff)
• mandelplot.py (modified) (5 diffs)

anuga_core/source/pypar-numeric/mandelbrot_example/compile.py

@@ -149 +149 @@
     if os.name == 'posix' and os.uname()[4] == 'x86_64':
       #Extra flags for 64 bit architectures
-      #s += ' -fPIC -m64' #gcc
-      s += ' -fPIC -tp amd64' #pgcc
+      s += ' -fPIC -m64' #gcc
+      #s += ' -fPIC -tp amd64' #pgcc
       
     

anuga_core/source/pypar-numeric/mandelbrot_example/mandel_sequential.py

@@ -12 +12 @@
 # User definable parameters
 kmax = 2**15  # Maximal number of iterations (=number of colors)
-M = N = 700   # width = height = N (200, 400, 600, 700 are good)
+M = N = 200   # width = height = N (200, 400, 600, 700 are good)
 
 # Region in complex plane

anuga_core/source/pypar-numeric/mandelbrot_example/mandelplot.py

@@ -18 +18 @@
     """
     
-    from Tkinter import Frame, Canvas, TOP, NW, PhotoImage
+    #from Tkinter import Frame, Canvas, TOP, NW, PhotoImage
+    from pylab import imshow, show
     from Numeric import transpose, Float
-    from Image import new             # PIL    
+    #from Image import new             # PIL    
     import time
 
@@ -27 +28 @@
     # User definable parameters
 
-    imtype = 'ppm'      # Image format (e.g 'ppm', 'bmp', 'tiff')
-    filename ='mandel'  # Base filename
+    #imtype = 'ppm'      # Image format (e.g 'ppm', 'bmp', 'tiff')
+    #filename ='mandel'  # Base filename
     
     exponent = 0.998    # Exponent for morphing colors, good with kmax = 2**15
@@ -47 +48 @@
     A = transpose(A).astype(Float)          
 
-    im = new("RGB", A.shape)
+    #im = new("RGB", A.shape)
 
     
@@ -78 +79 @@
     
     # Save image to file        
-    im.putdata(L)
-    im.save(filename + '.' + imtype, imtype)
+    imshow(A)
+    show()
+    
+    #im.save(filename + '.' + imtype, imtype)
     print 'Computed plot in %.2f seconds: ' %(time.time()-t0)
 
@@ -88 +91 @@
     #   sys.exit()
         
-    import os
-    os.system('xv %s' %(filename + '.' + imtype))
+    #import os
+    #os.system('xv %s' %(filename + '.' + imtype))