Frank Barknecht wrote: >> ( I replaced 2*M_PI in the code with TWO_PI to avoid a multiply >> (compiler optimisation maybe does that already)). > > Did you define that yourself? My math.h doesn't have TWO_PI ... depending on compiler options, my math.h doesn't have M_PI ... http://c-faq.com/fp/mpi.html Claude