Basic HTML version of Foils prepared August 7 98

Foil 17 PURE for Mandelbrot Sets

From DoD HPF Training -- 3. Parallel Constructs in HPF DoD Training and Others -- 1995-98. by Chuck Koelbel -- Rice University
! The caller (Explicit interface not shown)
FORALL ( i=1:n, j=1:m )
  • k(i,j) = mandelbrot( CMPLX((i-1)*1.0/(n-1), &
    • (j-1)*1.0/(m-1)), 1000 )
END FORALL
! The callee
PURE INTEGER FUNCTION mandelbrot(x, itol)
COMPLEX, INTENT(IN) :: x
INTEGER, INTENT(IN) :: itol
COMPLEX xtmp
INTEGER k
  • k = 0
  • xtmp = -x
  • DO WHILE (ABS(xtmp)<2.0 .AND. k<itol)
    • xtmp = xtmp*xtmp - x
    • k = k + 1
  • END DO
  • mandelbrot = k
END FUNCTION mandelbrot

© Northeast Parallel Architectures Center, Syracuse University, npac@npac.syr.edu

If you have any comments about this server, send e-mail to webmaster@npac.syr.edu.

Page produced by wwwfoil on Sun Aug 9 1998