GEANT321/matx55/xmm55.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 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.50  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE XMM55(A,B,C)
  13 C *** ---------------------------------------------                     05/07/81
  14 C *** MULTIPLICATION OF 5X5 MATRICES: A*B -> C.
  15 C *** FOR:I,J=1,5   C(I,J)=SUM:K=1,5(A(I,K)B(K,J))
  16 C *** NB: B AND C MAY WELL BE THE SAME MATRIX
  17 C
  18 C        Author: A. Haas (Freiburg University)
  19 C
  20 C *** ---------------------------------------------
  21 #if !defined(CERNLIB_SINGLE)
  22       DOUBLE PRECISION A,B,C,B1J,B2J,B3J,B4J,B5J
  23 #endif
  24       DIMENSION A(25),B(25),C(25)
  25       DO 10 J5=5,25,5
  26          B5J=B(J5)
  27          J4=J5-1
  28          B4J=B(J4)
  29          J3=J4-1
  30          B3J=B(J3)
  31          J2=J3-1
  32          B2J=B(J2)
  33          J1=J2-1
  34          B1J=B(J1)
  35          C(J1)=A( 1)*B1J+A( 6)*B2J+A(11)*B3J+A(16)*B4J+A(21)*B5J
  36          C(J2)=A( 2)*B1J+A( 7)*B2J+A(12)*B3J+A(17)*B4J+A(22)*B5J
  37          C(J3)=A( 3)*B1J+A( 8)*B2J+A(13)*B3J+A(18)*B4J+A(23)*B5J
  38          C(J4)=A( 4)*B1J+A( 9)*B2J+A(14)*B3J+A(19)*B4J+A(24)*B5J
  39          C(J5)=A( 5)*B1J+A(10)*B2J+A(15)*B3J+A(20)*B4J+A(25)*B5J
  40    10 CONTINUE
  41 *
  42       END