## File: r3mmt.f90

package info (click to toggle)
elkcode 5.4.24-2
 `12345678910111213141516171819202122232425262728293031323334353637` `````` ! Copyright (C) 2003-2004 J. K. Dewhurst, S. Sharma and C. Ambrosch-Draxl. ! This file is distributed under the terms of the GNU Lesser General Public ! License. See the file COPYING for license details. !BOP ! !ROUTINE: r3mmt ! !INTERFACE: subroutine r3mmt(a,b,c) ! !INPUT/OUTPUT PARAMETERS: ! a : input matrix 1 (in,real(3,3)) ! b : input matrix 2 (in,real(3,3)) ! c : output matrix (out,real(3,3)) ! !DESCRIPTION: ! Multiplies a real matrix with the transpose of another. ! ! !REVISION HISTORY: ! Created January 2003 (JKD) !EOP !BOC implicit none ! arguments real(8), intent(in) :: a(3,3),b(3,3) real(8), intent(out) :: c(3,3) c(1,1)=a(1,1)*b(1,1)+a(1,2)*b(1,2)+a(1,3)*b(1,3) c(2,1)=a(2,1)*b(1,1)+a(2,2)*b(1,2)+a(2,3)*b(1,3) c(3,1)=a(3,1)*b(1,1)+a(3,2)*b(1,2)+a(3,3)*b(1,3) c(1,2)=a(1,1)*b(2,1)+a(1,2)*b(2,2)+a(1,3)*b(2,3) c(2,2)=a(2,1)*b(2,1)+a(2,2)*b(2,2)+a(2,3)*b(2,3) c(3,2)=a(3,1)*b(2,1)+a(3,2)*b(2,2)+a(3,3)*b(2,3) c(1,3)=a(1,1)*b(3,1)+a(1,2)*b(3,2)+a(1,3)*b(3,3) c(2,3)=a(2,1)*b(3,1)+a(2,2)*b(3,2)+a(2,3)*b(3,3) c(3,3)=a(3,1)*b(3,1)+a(3,2)*b(3,2)+a(3,3)*b(3,3) return end subroutine !EOC ``````