var/external/wavelet/PrintFilter.c

   1 // Author: Aime' Fournier
   2 // E-mail: fournier@ucar.edu
   3
   4 #include "qf.h"                 // Data structs and function prototypes.
   5 #include "stdio.h"
   6
   7 void PrintFilter(int nam,int ran,int kin,pqf *f) {
   8   int k;
   9   real fs=0.,fse;
  10
  11   if( kin==0 ) fse=sqrt(2.);
  12   else if( kin==1 ) fse=0.;
  13   else abort();
  14   printf("Filter %s%02d[%d] has alpha=%d, omega=%2d, f={",nam,ran,kin,f->alpha,f->omega);
  15   for( k=f->alpha; k<=f->omega; k++ ) {
  16      printf(" %10.6f",f->f[k]);
  17      fs+=f->f[k];               // f sum
  18   }
  19   printf("} and sum-error %9.2e.\n",fs-fse);
  20 }