This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / f / mxdipr.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1996/04/01 15:02:32  mclareni
6 * Mathlib gen
7 *
8 *
9 #include "gen/pilot.h"
10       SUBROUTINE MXDIPR(A,B,C,IAD,JBD,IJD,IA,KA,JB,LB)
11 C
12 C     THIS ROUTINE COMPUTES THE DIRECT, OR TENSOR, OR KRONECKER, PRODUCT
13 C     OF TWO MATRICES A AND B, AND STORES THE RESULT IN C.
14 C
15       DIMENSION A(IAD,1),B(JBD,1),C(IJD,1)
16 C
17       IAJB=IA*JB
18       KALB=KA*LB
19       IF(IAJB*KALB .EQ. 0) RETURN
20       JB1=JB-1
21       LB1=LB-1
22       DO 1 I = 1,IAJB
23       IIA=(I+JB1)/JB
24       IIB=MOD(I-1,JB)+1
25       DO 1 K = 1,KALB
26       KKA=(K+LB1)/LB
27       KKB=MOD(K-1,LB)+1
28     1 C(I,K)=A(IIA,KKA)*B(IIB,KKB)
29       RETURN
30       END