Rewrite subroutine to use MKL - matrix multiplication

p { margin-bottom: 0.1in; line-height: 120%; }

Hi,

My first post on the forum!

I'm working with an older code written i Fortran at uni. The program consists of several subroutines and I would like to rewrite parts (subroutines) of the code so it runs faster (parallel). My knowledge of Fortran is limited/"ok". The first routine I like to look into is MT:

description:

            “Matrix multiplication

             B = AA*B*A (AA is the transpose of A)

             B = Full, symmetric matrix

             A = Unit matrix with an n by n submatrix (argument A) on the diagonal”

For code, please see attachment.

The subroutine is often called like this:

   do J = 1,10,3

      CALL MT(UM,FS,WE,12,3,J,PU,ER)

   continue

Then UM(3,3), FS(12,12) and WE(3)

Is there any way I can use MKL to rewrite the code? Using dgemm? Or maybe use openmp?

Any clues?

Thanks,

Lars