行列の積だと１００％になったよ。
ただし、行列×ベクトルだとコンパイルした時点で並列化してくれない。
もっとも多い計算が行列×ベクトルなのに・・・。
たとえば、これだと１００％になった。


PROGRAM REI
IMPLICIT NONE
INTEGER ::I,J
REAL(8),DIMENSION (5000,5000) :: A,B,C
DO I=1,5000
DO J=1,5000
A(I,J) = 1.1d0
B(I,J) = 2.1d0
END DO
END DO
C = MATMUL(A,B)
STOP
END PROGRAM REI