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