git.gag.com Git - debian/gnuradio/blob - gnuradio-core/src/lib/filter/gr_hilbert_fc.cc

   1 /* -*- c++ -*- */
   2 /*
   3  * Copyright 2004 Free Software Foundation, Inc.
   4  *
   5  * This file is part of GNU Radio
   6  *
   7  * GNU Radio is free software; you can redistribute it and/or modify
   8  * it under the terms of the GNU General Public License as published by
   9  * the Free Software Foundation; either version 2, or (at your option)
  10  * any later version.
  11  *
  12  * GNU Radio is distributed in the hope that it will be useful,
  13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15  * GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License
  18  * along with GNU Radio; see the file COPYING.  If not, write to
  19  * the Free Software Foundation, Inc., 51 Franklin Street,
  20  * Boston, MA 02110-1301, USA.
  21  */
  22
  23 #ifdef HAVE_CONFIG_H
  24 #include "config.h"
  25 #endif
  26
  27 #include <gr_hilbert_fc.h>
  28 #include <gr_firdes.h>
  29 #include <gr_fir_fff.h>
  30 #include <gr_fir_util.h>
  31
  32 // public constructor
  33 gr_hilbert_fc_sptr
  34 gr_make_hilbert_fc (unsigned int ntaps)
  35 {
  36   return gr_hilbert_fc_sptr (new gr_hilbert_fc (ntaps));
  37 }
  38
  39 gr_hilbert_fc::gr_hilbert_fc (unsigned int ntaps)
  40   : gr_sync_block ("hilbert_fc",
  41                    gr_make_io_signature (1, 1, sizeof (float)),
  42                    gr_make_io_signature (1, 1, sizeof (gr_complex))),
  43     d_ntaps (ntaps | 0x1),      // ensure ntaps is odd
  44     d_hilb (gr_fir_util::create_gr_fir_fff (gr_firdes::hilbert (d_ntaps)))
  45 {
  46   set_history (d_ntaps);
  47 }
  48
  49 gr_hilbert_fc::~gr_hilbert_fc ()
  50 {
  51   delete d_hilb;
  52 }
  53
  54 int
  55 gr_hilbert_fc::work (int noutput_items,
  56                           gr_vector_const_void_star &input_items,
  57                           gr_vector_void_star &output_items)
  58 {
  59   float *in = (float *) input_items[0];
  60   gr_complex *out = (gr_complex *) output_items[0];
  61
  62   for (int i = 0; i < noutput_items; i++)
  63     out[i] = gr_complex (in[i + d_ntaps/2],
  64                          d_hilb->filter (&in[i]));
  65
  66   return noutput_items;
  67 }