Allow control of number of threads with OMP_NUM_THREADS