!---------------------------------------------------------------------------
SUBROUTINE trho_apply(nchir,ndim_tc,phi3,y3,ifirst)
USE Numeric_Kinds_Module
USE CommonInfo_Module
USE APH_Module
USE Aph_Ham_Module
IMPLICIT NONE
!==============================================================================
! Apply T_rho using DGEMM
! DGEMM(transa,transb,l,n,m,alpha,a,lda,b,ldb,beta,c,ldc)
!==============================================================================
! I N P U T / O U T P U T
INTEGER,INTENT(IN) :: nchir, ndim_tc,ifirst

REAL(dp),INTENT(IN) :: phi3(nchir,nrho,ntheta)
REAL(dp),INTENT(INOUT) :: y3(nchir,nrho,ntheta)

!==============================================================================
! I N T E R N A L S
INTEGER :: itheta

!==============================================================================
 DO itheta=1,ntheta
   Call DGEMM( 'n', 't', nchir, nrho, nrho, 1.d0, phi3(:,:,itheta), nchir,  &
               t_rho, nrho, 0.d0,y3(:,:,itheta),nchir)
 ENDDO

RETURN
END SUBROUTINE trho_apply