3 * Copyright 2003 Free Software Foundation, Inc.
5 * This file is part of GNU Radio
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)
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.
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.
23 #ifndef _GR_ROTATOR_H_
24 #define _GR_ROTATOR_H_
26 #include <gr_complex.h>
30 gr_complex d_phase_incr;
33 gr_rotator () : d_phase (1), d_phase_incr (0) { }
35 void set_phase (gr_complex phase) { d_phase = phase; }
36 void set_phase_incr (gr_complex incr) { d_phase_incr = incr; }
38 gr_complex rotate (gr_complex in){
39 d_phase *= d_phase_incr; // incr our phase (complex mult == add phases)
41 d_phase /= abs(d_phase); // ensure multiplication is rotation
42 // FIXME. This is expensive. Maybe workaround using
43 // double precision complex???
45 return in * d_phase; // rotate in by phase
50 #endif /* _GR_ROTATOR_H_ */