* * $Id$ * * $Log$ * Revision 1.1.1.1 1996/03/06 15:37:36 mclareni * Add geane321 source directories * * * MATRIX MULTIPLICATION 5 BY 5. #include "geant321/pilot.h" *CMZ : 3.21/02 29/03/94 15.41.50 by S.Giani *-- Author : SUBROUTINE RMM55(A,B,C) C *** --------------------------------------------- 05/07/81 C *** MULTIPLICATION OF 5X5 MATRICES: A*B -> C. C *** FOR:I,J=1,5 C(I,J)=SUM:K=1,5(A(I,K)B(K,J)) C *** NB: B AND C MAY WELL BE THE SAME MATRIX C C Author: A. Haas (Freiburg University) C C *** --------------------------------------------- DIMENSION A(25),B(25),C(25) DO 10 J5=5,25,5 B5J=B(J5) J4=J5-1 B4J=B(J4) J3=J4-1 B3J=B(J3) J2=J3-1 B2J=B(J2) J1=J2-1 B1J=B(J1) C(J1)=A( 1)*B1J+A( 6)*B2J+A(11)*B3J+A(16)*B4J+A(21)*B5J C(J2)=A( 2)*B1J+A( 7)*B2J+A(12)*B3J+A(17)*B4J+A(22)*B5J C(J3)=A( 3)*B1J+A( 8)*B2J+A(13)*B3J+A(18)*B4J+A(23)*B5J C(J4)=A( 4)*B1J+A( 9)*B2J+A(14)*B3J+A(19)*B4J+A(24)*B5J C(J5)=A( 5)*B1J+A(10)*B2J+A(15)*B3J+A(20)*B4J+A(25)*B5J 10 CONTINUE * END