Un produit matriciel en somme que tu peux paralléliser de plusieurs façons :



!$ OMP PARALLEL DEFAULT(SHARED)
!$ OMP DO
do i=1,nar
do j=1,nbc
...