CTITLESAFXCC5 -- CORRELATIONS FOR SPATIAL PREDICTION
C***********************************************************************
C                 COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991            *
C***********************************************************************
CA    AUTHOR            D CORRIGAN
CA    DESIGNER          D CORRIGAN
CA    LANGUAGE          FORTRAN 77
CA    SYSTEM            IBM (SEE CRAY)
CA    WRITTEN           AUGUST 1988
C     REVISED           11-14-88  REPLACED CALL FROM SNDOT TO RNDOTP
C                                 OF MATH ADVANTAGE LIBRARY.
C     REVISED           01-15-89  JJC FOR SPARC PRODUCTION.
C     REVISED           02-05-89  CHANGED BACK TO SNDOT DUE TO THE VERY
C                                 SLOW PERFORMANCE OF MATHADV ROUTINES.
CA
CA
CA   THIS ROUTINE PERFORMS CORRELATIONS FOR SPATIAL PREDICTION
CA
CA
C    SUBROUTINES CALLED :
C                             ARMVE
C                             ARSET
C                             RNDOTP(MATH ADVANTAGE    )
C
C
C     CALL SAFXCC5 (D,NX,LX,R,C,A,B)
C
C              CORRELATIONS FOR SPATIAL PREDICTION
C
C     ARGUMENTS :
C
C               D   INPUT         NUMBER OF INPUT TRACES
C               NX  INPUT         NUMBER OF TRACES IN DESIGN GATE
C               LX  INPUT         NUMBER OF TRACES IN SPATIAL FILTER
C               R   OUTPUT        AUTOCORRELATION RESULTS
C               C   OUTPUT        CROSS CORRELATION RESULTS
C               A   INPUT/OUTPUT  AUXILLARY ARRAY
C               B   INPUT/OUTPUT  AUXILLARY ARRAY
C
C  ---------------------------------------------------------------------
C
      SUBROUTINE SAFXCC5( D,NX,LX,R,C,A,B )
C
      IMPLICIT INTEGER (A-Z)
      REAL A(1),B(1)
      COMPLEX C(1),D(1),R(1)
C
C
C
C       MOVE INPUT ARRAY TO B AND ZERO PAD
C
      IB  = 2*NX + 1
      CALL ARMVE(D,B,2*NX)
      CALL ARSET( B(IB),2*LX,0. )
C
C       COMPLEX CORRELATION
C
      IR  = 1
      II  = 2
      CALL SNDOT( LX,NX,A(IR),2,1,B(IR),2,0,B(IR),2,2 )
      CALL SNDOT( LX,NX,A(IR),2,3,B(II),2,0,B(II),2,2 )
      CALL SNDOT( LX,NX,A(II),2,2,B(IR),2,0,B(II),2,2 )
      CALL SNDOT( LX,NX,A(II),2,3,B(II),2,0,B(IR),2,2 )
C
CJJ   CALL RNDOTP( B(IR),2,0,B(IR),2,2,A(IR),2,NX,LX,0 )
CJJ   CALL RNDOTP( B(II),2,0,B(II),2,2,A(IR),2,NX,LX,2 )
CJJ   CALL RNDOTP( B(IR),2,0,B(II),2,2,A(II),2,NX,LX,1 )
CJJ   CALL RNDOTP( B(II),2,0,B(IR),2,2,A(II),2,NX,LX,2 )
C       MOVE OUTPUT TO R
C
      MX    = LX - 1
      MH    = MX/2
      DO 200 IX = 1,LX
      IR    = 2*IX - 1
      R(IX) = CMPLX( A(IR),A(IR+1) )
      C(IX) = CMPLX( 0.,0. )
200   CONTINUE
C
      IX    = MH + 1
      IR    = 2*IX - 1
      C(IX) = CMPLX( A(IR),A(IR+1) )
C
C
C
      RETURN
      END