GEANT321/matx55/dmm55.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 DMM55(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       DOUBLE PRECISION A,B,C,B1J,B2J,B3J,B4J,B5J
  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    10 C(J5)=A( 5)*B1J+A(10)*B2J+A(15)*B3J+A(20)*B4J+A(25)*B5J
  38 *
  39       END