*-- Author : S.Burke / J.V. Morris SUBROUTINE FKMLT(C1,C2,C3) ********************************************************************** * * * Multiplies two 5*5 symmetric matrices * * * * Output is C3 = C1.C2 (note that C3 is NOT symmetric!!!) * * * ********************************************************************** IMPLICIT DOUBLE PRECISION (A-H,O-Z) DIMENSION C1(5,5),C2(5,5),C3(5,5) ********************************************************************** C3(1,1) = C1(1,1)*C2(1,1) + C1(2,1)*C2(2,1) + C1(3,1)*C2(3,1) & + C1(4,1)*C2(4,1) + C1(5,1)*C2(5,1) C3(2,1) = C1(2,1)*C2(1,1) + C1(2,2)*C2(2,1) + C1(3,2)*C2(3,1) & + C1(4,2)*C2(4,1) + C1(5,2)*C2(5,1) C3(3,1) = C1(3,1)*C2(1,1) + C1(3,2)*C2(2,1) + C1(3,3)*C2(3,1) & + C1(4,3)*C2(4,1) + C1(5,3)*C2(5,1) C3(4,1) = C1(4,1)*C2(1,1) + C1(4,2)*C2(2,1) + C1(4,3)*C2(3,1) & + C1(4,4)*C2(4,1) + C1(5,4)*C2(5,1) C3(5,1) = C1(5,1)*C2(1,1) + C1(5,2)*C2(2,1) + C1(5,3)*C2(3,1) & + C1(5,4)*C2(4,1) + C1(5,5)*C2(5,1) C3(1,2) = C1(1,1)*C2(2,1) + C1(2,1)*C2(2,2) + C1(3,1)*C2(3,2) & + C1(4,1)*C2(4,2) + C1(5,1)*C2(5,2) C3(2,2) = C1(2,1)*C2(2,1) + C1(2,2)*C2(2,2) + C1(3,2)*C2(3,2) & + C1(4,2)*C2(4,2) + C1(5,2)*C2(5,2) C3(3,2) = C1(3,1)*C2(2,1) + C1(3,2)*C2(2,2) + C1(3,3)*C2(3,2) & + C1(4,3)*C2(4,2) + C1(5,3)*C2(5,2) C3(4,2) = C1(4,1)*C2(2,1) + C1(4,2)*C2(2,2) + C1(4,3)*C2(3,2) & + C1(4,4)*C2(4,2) + C1(5,4)*C2(5,2) C3(5,2) = C1(5,1)*C2(2,1) + C1(5,2)*C2(2,2) + C1(5,3)*C2(3,2) & + C1(5,4)*C2(4,2) + C1(5,5)*C2(5,2) C3(1,3) = C1(1,1)*C2(3,1) + C1(2,1)*C2(3,2) + C1(3,1)*C2(3,3) & + C1(4,1)*C2(4,3) + C1(5,1)*C2(5,3) C3(2,3) = C1(2,1)*C2(3,1) + C1(2,2)*C2(3,2) + C1(3,2)*C2(3,3) & + C1(4,2)*C2(4,3) + C1(5,2)*C2(5,3) C3(3,3) = C1(3,1)*C2(3,1) + C1(3,2)*C2(3,2) + C1(3,3)*C2(3,3) & + C1(4,3)*C2(4,3) + C1(5,3)*C2(5,3) C3(4,3) = C1(4,1)*C2(3,1) + C1(4,2)*C2(3,2) + C1(4,3)*C2(3,3) & + C1(4,4)*C2(4,3) + C1(5,4)*C2(5,3) C3(5,3) = C1(5,1)*C2(3,1) + C1(5,2)*C2(3,2) + C1(5,3)*C2(3,3) & + C1(5,4)*C2(4,3) + C1(5,5)*C2(5,3) C3(1,4) = C1(1,1)*C2(4,1) + C1(2,1)*C2(4,2) + C1(3,1)*C2(4,3) & + C1(4,1)*C2(4,4) + C1(5,1)*C2(5,4) C3(2,4) = C1(2,1)*C2(4,1) + C1(2,2)*C2(4,2) + C1(3,2)*C2(4,3) & + C1(4,2)*C2(4,4) + C1(5,2)*C2(5,4) C3(3,4) = C1(3,1)*C2(4,1) + C1(3,2)*C2(4,2) + C1(3,3)*C2(4,3) & + C1(4,3)*C2(4,4) + C1(5,3)*C2(5,4) C3(4,4) = C1(4,1)*C2(4,1) + C1(4,2)*C2(4,2) + C1(4,3)*C2(4,3) & + C1(4,4)*C2(4,4) + C1(5,4)*C2(5,4) C3(5,4) = C1(5,1)*C2(4,1) + C1(5,2)*C2(4,2) + C1(5,3)*C2(4,3) & + C1(5,4)*C2(4,4) + C1(5,5)*C2(5,4) C3(1,5) = C1(1,1)*C2(5,1) + C1(2,1)*C2(5,2) + C1(3,1)*C2(5,3) & + C1(4,1)*C2(5,4) + C1(5,1)*C2(5,5) C3(2,5) = C1(2,1)*C2(5,1) + C1(2,2)*C2(5,2) + C1(3,2)*C2(5,3) & + C1(4,2)*C2(5,4) + C1(5,2)*C2(5,5) C3(3,5) = C1(3,1)*C2(5,1) + C1(3,2)*C2(5,2) + C1(3,3)*C2(5,3) & + C1(4,3)*C2(5,4) + C1(5,3)*C2(5,5) C3(4,5) = C1(4,1)*C2(5,1) + C1(4,2)*C2(5,2) + C1(4,3)*C2(5,3) & + C1(4,4)*C2(5,4) + C1(5,4)*C2(5,5) C3(5,5) = C1(5,1)*C2(5,1) + C1(5,2)*C2(5,2) + C1(5,3)*C2(5,3) & + C1(5,4)*C2(5,4) + C1(5,5)*C2(5,5) RETURN END *