'omp_set_num_threads' and 'omp parallel for' in cgl_step_nonlin()
basic 'omp parallel for' optimisations in etd_step.c