      subroutine su3mult_mh_mh_m(am1,am2,rm) ! Bernd Berg Mmar 02 2008.
c Alexei Bazavov Oct 20 2005. Multiply SU(3) conjugate matrix (mh) by 
c                  conjugate matrix (mh) and put result in matrix (m).
c am1 -- matrix 1
c am2 -- matrix 2
c rm -- resulting matrix
      include '../../Libs/Fortran/implicit.sta'
      include '../../Libs/Fortran/constants.par'
      dimension am1(1:18),am2(1:18),rm(1:18),am3(1:18),am4(1:18)

c TEMPORARY SOLUTION

c make hermitian conjugates
      am3(1)=am1(1)
      am3(2)=am1(4)
      am3(3)=am1(7)
      am3(4)=am1(2)
      am3(5)=am1(5)
      am3(6)=am1(8)
      am3(7)=am1(3)
      am3(8)=am1(6)
      am3(9)=am1(9)
      am3(10)=-am1(10)
      am3(11)=-am1(13)
      am3(12)=-am1(16)
      am3(13)=-am1(11)
      am3(14)=-am1(14)
      am3(15)=-am1(17)
      am3(16)=-am1(12)
      am3(17)=-am1(15)
      am3(18)=-am1(18)

      am4(1)=am2(1)
      am4(2)=am2(4)
      am4(3)=am2(7)
      am4(4)=am2(2)
      am4(5)=am2(5)
      am4(6)=am2(8)
      am4(7)=am2(3)
      am4(8)=am2(6)
      am4(9)=am2(9)
      am4(10)=-am2(10)
      am4(11)=-am2(13)
      am4(12)=-am2(16)
      am4(13)=-am2(11)
      am4(14)=-am2(14)
      am4(15)=-am2(17)
      am4(16)=-am2(12)
      am4(17)=-am2(15)
      am4(18)=-am2(18)

c multiply
      call su3mult_m_m_m(am3,am4,rm)

      return
      end
