Salve a tutti mi rivolgo a voi per trovare una soluzione ad un problema che mi snerva da alcuni giorni .
Dunque io ho una serie di chiamate in fortran a librerie blas mkl dgemm , vorrei parallelilizzare le chiamate e dunque metto

!$omp section
call DGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C ,LDC)

mettendo default (SHARED)

Problema : le chiamate restituiscono valori discordi rispetto alla versione non paralellizata .
Sapete aiutarmi ?