Indeed that would be possible, but in connection with OpenMP you would have to decide whether to put the OpenMP parallelization inside, or not. At the moment the library does not do any OpenMP at all so that the user can decide how they want to parallelize the operations.
1 Like
Another small question: Are you sure the line above is correct? I think it should be
tid = omp_get_thread_num()
Thanks!
It is, that’s an openmp “Sentinel”, part of the openmp standard for Fortran. Meaning that the line will take effect only if openmp is activated.
2 Likes