## File: r3mm.f90

package info (click to toggle)
elkcode 5.4.24-2
• area: main
• in suites: buster
• size: 12,840 kB
• sloc: f90: 48,415; fortran: 22,457; perl: 965; makefile: 384; sh: 369; python: 105; ansic: 67
 file content (37 lines) | stat: -rw-r--r-- 1,094 bytes parent folder | download | duplicates (2)
 12345678910111213141516171819202122232425262728293031323334353637 ! Copyright (C) 2002-2005 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: r3mm ! !INTERFACE: subroutine r3mm(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 two real $3\times 3$ matrices. ! ! !REVISION HISTORY: ! Created April 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(2,1)+a(1,3)*b(3,1) c(2,1)=a(2,1)*b(1,1)+a(2,2)*b(2,1)+a(2,3)*b(3,1) c(3,1)=a(3,1)*b(1,1)+a(3,2)*b(2,1)+a(3,3)*b(3,1) c(1,2)=a(1,1)*b(1,2)+a(1,2)*b(2,2)+a(1,3)*b(3,2) c(2,2)=a(2,1)*b(1,2)+a(2,2)*b(2,2)+a(2,3)*b(3,2) c(3,2)=a(3,1)*b(1,2)+a(3,2)*b(2,2)+a(3,3)*b(3,2) c(1,3)=a(1,1)*b(1,3)+a(1,2)*b(2,3)+a(1,3)*b(3,3) c(2,3)=a(2,1)*b(1,3)+a(2,2)*b(2,3)+a(2,3)*b(3,3) c(3,3)=a(3,1)*b(1,3)+a(3,2)*b(2,3)+a(3,3)*b(3,3) return end subroutine !EOC