Contents
     dwiener - perform Wiener deconvolution of two signals
     SUBROUTINE DWIENER(N_POINTS, ACOR, XCOR, FLTR, EROP, ISW, IERR)
     INTEGER N_POINTS, ISW, IERR
     DOUBLE PRECISION ACOR(*), XCOR(*), FLTR(*), EROP(*)
     SUBROUTINE DWIENER_64(N_POINTS, ACOR, XCOR, FLTR, EROP, ISW, IERR)
     INTEGER*8 N_POINTS, ISW, IERR
     DOUBLE PRECISION ACOR(*), XCOR(*), FLTR(*), EROP(*)
  F95 INTERFACE
     SUBROUTINE WIENER(N_POINTS, ACOR, XCOR, FLTR, EROP, ISW, IERR)
     INTEGER :: N_POINTS, ISW, IERR
     REAL(8), DIMENSION(:) :: ACOR, XCOR, FLTR, EROP
     SUBROUTINE WIENER_64(N_POINTS, ACOR, XCOR, FLTR, EROP, ISW, IERR)
     INTEGER(8) :: N_POINTS, ISW, IERR
     REAL(8), DIMENSION(:) :: ACOR, XCOR, FLTR, EROP
  C INTERFACE
     #include <sunperf.h>
     void dwiener(int n_points, double *acor, double *xcor,  dou-
               ble *fltr, double *erop, int *isw, int *ierr);
     void dwiener_64(long n_points, double *acor,  double  *xcor,
               double   *fltr,  double  *erop,  long  *isw,  long
               *ierr);
     dwiener performs Wiener deconvolution of two signals.
     N_POINTS (input)
               On entry, the number of points in the input corre-
               lations.  Unchanged on exit.
     ACOR (input)
               On entry, autocorrelation coefficients.  Unchanged
               on exit.
     XCOR (input)
               On    entry,    cross-correlation    coefficients.
               Unchanged on exit.
     FLTR (output)
               On exit, filter coefficients.
     EROP (output)
               On exit, the prediction error.
     ISW (input)
               On entry, if  ISW  .EQ.  0  then  perform  spiking
               deconvolution,  otherwise perform general deconvo-
               lution.  Unchanged on exit.
     IERR (output)
               On exit, the deconvolution was successful iff IERR
               .EQ. 0, otherwise there was an error.