GEANT321/matx55/rmm55.F

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