CTITLEMTOEPI --SET UP FOR SOLVING COMPLEX BLOCK TOEPLITZ EQUATION SYSTEM
C***********************************************************************
C                 COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991            *
C***********************************************************************
C***********************************************************************
CABS  MTOEPI - SET UP FOR SOLVING COMPLEX BLOCK TOEPLITZ EQUATION SYSTEM
C
C     COPYRIGHT ATLANTIC RICHFIELD COMPANY 1989.
C
C     ALL RIGHTS RESERVED.  NO PART OF THIS PROGRAM MAY BE PHOTOCOPIED,
C     REPRODUCED, OR  TRANSLATED TO ANOTHER PROGRAM LANGUAGE WITHOUT THE
C     PRIOR CONSENT OF ATLANTIC RICHFIELD COMPANY.
C
C
CA    DESIGNER             D CORRIGAN
CA    AUTHOR               D CORRIGAN
CA    LANGUAGE             VS FORTRAN
CA    SYSTEM               IBM
CA    WRITTEN              12-13-88
CA
CA    REVISED      10-05-90 D.CORRIGAN
CA                               TO AVOID ASSUMPTIONS ABOUT
CA                               THE OVERALL MAGNITUDE OF THE
CA                               CROSS CORRELATIONS BY SETTING
CA                               ZTOL = ETOL * MAX(ABS(R))
CA    REVISED      10-19-90 CLJ  RENAME FROM BCTEP, ADD TO THE CRAY
CA
CA    PURPOSE OF PROGRAM:
CA
CA    SET UP FOR SOLVING
CA    A LINEAR SYSTEM WHOSE MATRIX IS
CA    COMPLEX AND BLOCK TOEPLITZ
CA
CA    CALLING PROCEDURE:
CA       SUBROUTINE MTOEPI( N,L,R,RX,G,F,W,* )
CA
CA    CALLING ARGUMENTS:
CA
CA    ARGUMENTS ( INPUT):
CA
CA    N   - NUMBER OF CHANNELS
CA    L   - LENGTH OF EACH N CHANNEL FILTER
CA    R   - L ORDER N TOEPLITZ MATRICES (N,N,L)
CA    RX  - MAGNITUDE OF THE MAXIMUM DIAGONAL ELEMENT
CA    G   - L N-COMPONENT VECTORS - THE RHS
CA
CA    ARGUMENTS (OUTPUT):
CA
CA    F   - L N-COMPONENT VECTORS - THE FILTERS
CA
CA    ARGUMENTS (WORKSPACE)
CA
CA    W - ALLOCATED TO HOLD:
CA
CA    A   - L - N X N MATRICES
CA    B   - L - N X N MATRICES
CA    AP  - L - N X N MATRICES
CA    BP  - L - N X N MATRICES
CA    VA  - N X N MATRIX
CA    VB  - N X N MATRIX
CA    DA  - N X N MATRIX
CA    DB  - N X N MATRIX
CA    CA  - N X N MATRIX
CA    CB  - N X N MATRIX
CA    CF  - N-COMPONENT VECTOR
CA    CI  - N-COMPONENT VECTOR
CA    GA  - N-COMPONENT VECTOR
CA    TM  - N X N MATRIX
CA
CA    LENGTH OF W : ( 4*L*N*N + 7*N*N + 3*N ) COMPLEX WORDS
CA
CA
C     SUBROUTINES CALLED:
C
C                              BCLEV
C
CEND
C***********************************************************************
C
      SUBROUTINE MTOEPI( N,L,R,RX,G,F,W,* )
C
      IMPLICIT   INTEGER(A-Z)
C
C  ---------------------------------------------------------------------
C
      COMPLEX R(N,N,L),G(N,L)
      COMPLEX F(N,L)
      COMPLEX W(1)
C
      REAL RX
C
C  ---------------------------------------------------------------------
C
C       ALLOCATE WORK AREAS IN W
C
      NN  = N*N
      LNN = L*NN
C
      IAA = 1
      IAP = IAA + LNN
      IBB = IAP + LNN
      IBP = IBB + LNN
      IVA = IBP + LNN
      IVB = IVA +  NN
      IDA = IVB +  NN
      IDB = IDA +  NN
      ICA = IDB +  NN
      ICB = ICA +  NN
      ICF = ICB +  NN
      ICI = ICF +   N
      IGA = ICI +   N
      ITM = IGA +   N
C
C       CALL MTOEPL
C
      CALL MTOEPL( N,L,R,RX,G,F,
     *             W(IAA),W(IAP),W(IBB),W(IBP),
     *             W(IVA),W(IVB),W(IDA),W(IDB),
     *             W(ICA),W(ICB),W(ICF),W(ICI),
     *             W(IGA),W(ITM),*200 )
C
      RETURN
C
C  ---------------------------------------------------------------------
C
C
200   RETURN 1
      END