double start[64], elapsed[64];
#ifdef _OPENMP
#pragma omp threadprivate (start, elapsed)
#endif