Last change
on this file since 5847 was
5847,
checked in by steve, 15 years ago

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

File size:
640 bytes

Line  

1  #!/usr/bin/env python 

2  

3  """Sequential program computing the Mandelbrot set. 

4  

5  Ole Nielsen, SUT 2003 

6  """ 

7  

8  from mandelbrot import calculate_region 

9  from mandelplot import plot 

10  import time 

11  

12  # User definable parameters 

13  kmax = 2**15 # Maximal number of iterations (=number of colors) 

14  M = N = 200 # width = height = N (200, 400, 600, 700 are good) 

15  

16  # Region in complex plane 

17  real_min = 2.0 

18  real_max = 1.0 

19  imag_min = 1.5 

20  imag_max = 1.5 

21  

22  # Compute Mandelbrot set 

23  t0 = time.time() 

24  A = calculate_region(real_min, real_max, imag_min, imag_max, kmax, M, N) 

25  print 'Computed region in %.2f seconds' %(time.time()t0) 

26  

27  # Plot result 

28  plot(A, kmax) 

29  

30  

31  

32  

Note: See
TracBrowser
for help on using the repository browser.