3 * Copyright 2006 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.
27 #include <gr_pa_2x2_phase_combiner.h>
28 #include <gr_io_signature.h>
30 gr_pa_2x2_phase_combiner_sptr
31 gr_make_pa_2x2_phase_combiner()
33 return gr_pa_2x2_phase_combiner_sptr(new gr_pa_2x2_phase_combiner());
36 gr_pa_2x2_phase_combiner::gr_pa_2x2_phase_combiner ()
37 : gr_sync_block ("pa_2x2_phase_combiner",
38 gr_make_io_signature (1, 1, NM * sizeof (gr_complex)),
39 gr_make_io_signature (1, 1, sizeof (float)))
45 gr_pa_2x2_phase_combiner::set_theta(float theta)
48 gr_complex j = gr_complex(0,1);
49 d_phase[0] = exp(j * (float) (M_PI * (sin(theta) + cos(theta))));
50 d_phase[1] = exp(j * (float) (M_PI * cos(theta)));
51 d_phase[2] = exp(j * (float) (M_PI * sin(theta)));
52 d_phase[3] = exp(j * (float) 0.0);
56 gr_pa_2x2_phase_combiner::work (int noutput_items,
57 gr_vector_const_void_star &input_items,
58 gr_vector_void_star &output_items)
60 const gr_complex *in = (const gr_complex *) input_items[0];
61 gr_complex *out = (gr_complex *) output_items[0];
63 for (int i = 0; i < noutput_items; i++){
65 acc += in[0] * d_phase[0];
66 acc += in[1] * d_phase[1];
67 acc += in[2] * d_phase[2];
68 acc += in[3] * d_phase[3];