var/external/wavelet/qf_w.c

   1 // Author: Aime' Fournier
   2 // E-mail: fournier@ucar.edu
   3
   4 #include <stdio.h>
   5 #include "qf.h"         // for pqf ...
   6 #include "real.h"
   7
   8 // See 2009/8/4 e-mail from john@michalakes.us
   9 #ifndef CRAY
  10 #   ifdef NOUNDERSCORE
  11 #      define QF_W qf_w
  12 #   else
  13 #      ifdef F2CSTYLE
  14 #         define QF_W qf_w__
  15 #      else
  16 #         define QF_W qf_w_
  17 #      endif
  18 #   endif
  19 #endif
  20
  21 void QF_W (char *nam,   // in: filter name
  22            int *ran,    // in: filter length
  23            int *kin,    // in: 0 (lpf) or 1 (hpf)
  24            real *f) {   // out: filter
  25   int n;
  26   pqf *p;               // struct from qf.h
  27
  28   p=qf(nam,*ran,*kin);  // assign p
  29 //printf("qf_w::qf_w(%s,%d,%d,{",nam,*ran,*kin);
  30   for( n=0; n<*ran; n++ ) {
  31      f[n]=p->f[n];
  32 //   printf("%9.6g,",f[n]);
  33   }
  34 //printf("\b})\n");
  35
  36   free(p->fp);free(p);
  37 }