[furystep]

Hi all,

I have been using cfMesh for a while and I noticed a weird problem about the mesh results.

I run cfMesh with cartesianMesh on two identical testcases (same input on same machine) but I found that the points' coordinates of the results appear to be different (last 2/3 decimals can be different). I know that cfMesh sets dynamically the number of threads to use.

I am not an expert in parallelization and omp library but if i force it to use only one thread the results are identical, while with more than one thread (set from bash env variable or using omp.h) the are not.

I tried setting omp_set_dynamic equal to 0 at beginning of the code (as suggested in this https://stackoverflow.com/questions/...is-not-working thread) but it didn't work.

Could someone tell me if there is any way to make the results to be perfectly identical even with more than one thread (by using idk some barriers or commands)?