Logo Search packages:      
Sourcecode: zaptel version File versions

mkfilter.h

/* mkfilter -- given n, compute recurrence relation
   to implement Butterworth, Bessel or Chebyshev filter of order n
   A.J. Fisher, University of York   <fisher@minster.york.ac.uk>
   September 1992 */

#include <string.h>

/* Header file */

#define global
#define unless(x)   if(!(x))
#define until(x)    while(!(x))

#define VERSION       "4.6"
#undef      PI
#define PI      3.14159265358979323846  /* Microsoft C++ does not define M_PI ! */
#define TWOPI         (2.0 * PI)
#define EPS     1e-10
#define MAXORDER    10
#define MAXPZ         512         /* .ge. 2*MAXORDER, to allow for doubling of poles in BP filter;
                         high values needed for FIR filters */
#define MAXSTRING   256

typedef void (*proc)();
typedef unsigned int uint;

extern "C"
  { double atof(const char*);
    int atoi(char*);
    void exit(int);
  };

extern const char *progname;
extern void readdata(char*, double&, int&, double*, int&, double*);

inline double sqr(double x)       { return x*x;                    }
inline bool seq(char *s1, char *s2) { return strcmp(s1,s2) == 0;         }
inline bool onebit(uint m)        { return (m != 0) && ((m & m-1) == 0);     }

inline double asinh(double x)
  { /* Microsoft C++ does not define */
    return log(x + sqrt(1.0 + sqr(x)));
  }

inline double fix(double x)
  { /* nearest integer */
    return (x >= 0.0) ? floor(0.5+x) : -floor(0.5-x);
  }


Generated by  Doxygen 1.6.0   Back to index