Multiplies two complex rectangular matrices, AB.

Required Arguments

A — Complex NRA by NCA rectangular matrix.   (Input)

B — Complex NRB by NCB rectangular matrix.   (Input)

C — Complex NRC by NCC rectangular matrix containing the product A * B.   (Output)

Optional Arguments

NRA — Number of rows of A.   (Input)
Default: NRA = size (A,1).

NCA — Number of columns of A.   (Input)
Default: NCA = size (A,2).

LDA — Leading dimension of A exactly as specified in the dimension statement of the calling program.   (Input)
Default: LDA = size (A,1).

NRB — Number of rows of B.   (Input)
NRB must be equal to NCA.
Default: NRB = size (B,1).

NCB — Number of columns of B.   (Input)
Default: NCB = size (B,2).

LDB — Leading dimension of B exactly as specified in the dimension statement of the calling program.   (Input)
Default: LDB = size (B,1).

NRC — Number of rows of C.   (Input)
NRC must be equal to NRA.
Default: NRC = size (C,1).

NCC — Number of columns of C.   (Input)
NCC must be equal to NCB.
Default: NCC = size (C,2).

LDC — Leading dimension of C exactly as specified in the dimension statement of the calling program.   (Input)
Default: LDC = size (C,1).

FORTRAN 90 Interface

Generic:          CALL MCRCR (A, B, C [,…])

Specific:          The specific interface names are S_MCRCR and D_MCRCR.

FORTRAN 77 Interface

Single:            CALL MCRCR (NRA, NCA, A, LDA, NRB, NCB, B, LDB, NRC, NCC, C, LDC)

Double:          The double precision name is DMCRCR.

Description

Given the complex rectangular matrices A and B, MCRCR computes the complex rectangular matrix
C = AB.

Example

Multiply a 3 ´ 4 complex matrix by a 4 ´ 3 complex matrix. The output matrix will be a 3 ´ 3 complex matrix.

 

      USE MCRCR_INT

      USE WRCRN_INT

 

      IMPLICIT   NONE

!                                 Declare variables

      INTEGER    NCA, NCB, NCC, NRA, NRB, NRC

      PARAMETER  (NCA=4, NCB=3, NCC=3, NRA=3, NRB=4, NRC=3)

!

      COMPLEX    A(NRA,NCA), B(NRB,NCB), C(NRC,NCC)

!                                 Set values for A

!            A = ( 1.0 + 1.0i  -1.0+ 2.0i  0.0 + 1.0i  0.0 - 2.0i )

!                ( 3.0 + 7.0i  6.0 - 4.0i  2.0 - 1.0i  0.0 + 1.0i )

!                ( 1.0 + 0.0i  1.0 - 2.0i  -2.0+ 0.0i  0.0 + 0.0i )

!

!                                 Set values for B

!            B = ( 2.0 + 1.0i  3.0 + 2.0i  3.0 + 1.0i )

!                ( 2.0 - 1.0i  4.0 - 2.0i  5.0 - 3.0i )

!                ( 1.0 + 0.0i  0.0 - 1.0i  0.0 + 1.0i )

!                ( 2.0 + 1.0i  1.0 + 2.0i  0.0 - 1.0i )

!

      DATA A/(1.0,1.0), (3.0,7.0), (1.0,0.0), (-1.0,2.0), (6.0,-4.0), &

          (1.0,-2.0), (0.0,1.0), (2.0,-1.0), (-2.0,0.0), (0.0,-2.0), &

          (0.0,1.0), (0.0,0.0)/

      DATA B/(2.0,1.0), (2.0,-1.0), (1.0,0.0), (2.0,1.0), (3.0,2.0), &

          (4.0,-2.0), (0.0,-1.0), (1.0,2.0), (3.0,1.0), (5.0,-3.0), &

          (0.0,1.0), (0.0,-1.0)/

!                                 Compute C = A*B

      CALL MCRCR (A, B, C)

!                                 Print results

      CALL WRCRN ('C = A*B', C)

      END

Output

 

                     C = A*B
                 1                2                3
1  (  3.00,  5.00)  (  6.00, 13.00)  (  0.00, 17.00)
2  (  8.00,  4.00)  (  8.00, -2.00)  ( 22.00,-12.00)
3  (  0.00, -4.00)  (  3.00, -6.00)  (  2.00,-14.00)

Visual Numerics, Inc. Visual Numerics - Developers of IMSL and PV-WAVE http://www.vni.com/ PHONE: 713.784.3131 FAX:713.781.9260