diff options
Diffstat (limited to 'gnuradio-runtime/lib/math/qa_fxpt_nco.cc')
| -rw-r--r-- | gnuradio-runtime/lib/math/qa_fxpt_nco.cc | 104 |
1 files changed, 54 insertions, 50 deletions
diff --git a/gnuradio-runtime/lib/math/qa_fxpt_nco.cc b/gnuradio-runtime/lib/math/qa_fxpt_nco.cc index da10f49c92..2c4a8c327a 100644 --- a/gnuradio-runtime/lib/math/qa_fxpt_nco.cc +++ b/gnuradio-runtime/lib/math/qa_fxpt_nco.cc @@ -34,76 +34,80 @@ static const float SIN_COS_TOLERANCE = 1e-5; -//static const float SIN_COS_FREQ = 5003; +// static const float SIN_COS_FREQ = 5003; static const float SIN_COS_FREQ = 4096; static const int SIN_COS_BLOCK_SIZE = 100000; -static double max_d(double a, double b) -{ - return fabs(a) > fabs(b) ? a : b; -} +static double max_d(double a, double b) { return fabs(a) > fabs(b) ? a : b; } -BOOST_AUTO_TEST_CASE(t0) { - gr::nco<float,float> ref_nco; - gr::fxpt_nco new_nco; - double max_error = 0, max_phase_error = 0; +BOOST_AUTO_TEST_CASE(t0) +{ + gr::nco<float, float> ref_nco; + gr::fxpt_nco new_nco; + double max_error = 0, max_phase_error = 0; - ref_nco.set_freq((float)(2 * GR_M_PI / SIN_COS_FREQ)); - new_nco.set_freq((float)(2 * GR_M_PI / SIN_COS_FREQ)); + ref_nco.set_freq((float)(2 * GR_M_PI / SIN_COS_FREQ)); + new_nco.set_freq((float)(2 * GR_M_PI / SIN_COS_FREQ)); - BOOST_CHECK(std::abs(ref_nco.get_freq() - new_nco.get_freq()) <= SIN_COS_TOLERANCE); + BOOST_CHECK(std::abs(ref_nco.get_freq() - new_nco.get_freq()) <= SIN_COS_TOLERANCE); - for(int i = 0; i < SIN_COS_BLOCK_SIZE; i++) { - float ref_sin = ref_nco.sin(); - float new_sin = new_nco.sin(); - //printf ("i = %6d\n", i); - BOOST_CHECK(std::abs(ref_sin - new_sin) <= SIN_COS_TOLERANCE); + for (int i = 0; i < SIN_COS_BLOCK_SIZE; i++) { + float ref_sin = ref_nco.sin(); + float new_sin = new_nco.sin(); + // printf ("i = %6d\n", i); + BOOST_CHECK(std::abs(ref_sin - new_sin) <= SIN_COS_TOLERANCE); - max_error = max_d(max_error, ref_sin-new_sin); + max_error = max_d(max_error, ref_sin - new_sin); - float ref_cos = ref_nco.cos(); - float new_cos = new_nco.cos(); - BOOST_CHECK(std::abs(ref_cos - new_cos) <= SIN_COS_TOLERANCE); + float ref_cos = ref_nco.cos(); + float new_cos = new_nco.cos(); + BOOST_CHECK(std::abs(ref_cos - new_cos) <= SIN_COS_TOLERANCE); - max_error = max_d(max_error, ref_cos-new_cos); + max_error = max_d(max_error, ref_cos - new_cos); - ref_nco.step(); - new_nco.step(); + ref_nco.step(); + new_nco.step(); - BOOST_CHECK(std::abs(ref_nco.get_phase() - new_nco.get_phase()) <= SIN_COS_TOLERANCE); + BOOST_CHECK(std::abs(ref_nco.get_phase() - new_nco.get_phase()) <= + SIN_COS_TOLERANCE); - max_phase_error = max_d(max_phase_error, ref_nco.get_phase()-new_nco.get_phase()); - } - // printf ("Fxpt max error %.9f, max phase error %.9f\n", max_error, max_phase_error); + max_phase_error = + max_d(max_phase_error, ref_nco.get_phase() - new_nco.get_phase()); + } + // printf ("Fxpt max error %.9f, max phase error %.9f\n", max_error, + // max_phase_error); } -BOOST_AUTO_TEST_CASE(t1) { - gr::nco<float,float> ref_nco; - gr::fxpt_nco new_nco; - gr_complex* ref_block = new gr_complex[SIN_COS_BLOCK_SIZE]; - gr_complex* new_block = new gr_complex[SIN_COS_BLOCK_SIZE]; - double max_error = 0; +BOOST_AUTO_TEST_CASE(t1) +{ + gr::nco<float, float> ref_nco; + gr::fxpt_nco new_nco; + gr_complex* ref_block = new gr_complex[SIN_COS_BLOCK_SIZE]; + gr_complex* new_block = new gr_complex[SIN_COS_BLOCK_SIZE]; + double max_error = 0; - ref_nco.set_freq((float)(2 * GR_M_PI / SIN_COS_FREQ)); - new_nco.set_freq((float)(2 * GR_M_PI / SIN_COS_FREQ)); + ref_nco.set_freq((float)(2 * GR_M_PI / SIN_COS_FREQ)); + new_nco.set_freq((float)(2 * GR_M_PI / SIN_COS_FREQ)); - BOOST_CHECK(std::abs(ref_nco.get_freq() - new_nco.get_freq()) <= SIN_COS_TOLERANCE); + BOOST_CHECK(std::abs(ref_nco.get_freq() - new_nco.get_freq()) <= SIN_COS_TOLERANCE); - ref_nco.sincos((gr_complex*)ref_block, SIN_COS_BLOCK_SIZE); - new_nco.sincos((gr_complex*)new_block, SIN_COS_BLOCK_SIZE); + ref_nco.sincos((gr_complex*)ref_block, SIN_COS_BLOCK_SIZE); + new_nco.sincos((gr_complex*)new_block, SIN_COS_BLOCK_SIZE); - for(int i = 0; i < SIN_COS_BLOCK_SIZE; i++) { - BOOST_CHECK(std::abs(ref_block[i].real() - new_block[i].real()) <= SIN_COS_TOLERANCE); - max_error = max_d (max_error, ref_block[i].real()-new_block[i].real()); + for (int i = 0; i < SIN_COS_BLOCK_SIZE; i++) { + BOOST_CHECK(std::abs(ref_block[i].real() - new_block[i].real()) <= + SIN_COS_TOLERANCE); + max_error = max_d(max_error, ref_block[i].real() - new_block[i].real()); - BOOST_CHECK(std::abs(ref_block[i].imag() - new_block[i].imag()) <= SIN_COS_TOLERANCE); - max_error = max_d (max_error, ref_block[i].imag()-new_block[i].imag()); - } - BOOST_CHECK(std::abs(ref_nco.get_phase() - new_nco.get_phase()) <= SIN_COS_TOLERANCE); - // printf ("Fxpt max error %.9f, max phase error %.9f\n", max_error, max_phase_error); - delete[] ref_block; - delete[] new_block; + BOOST_CHECK(std::abs(ref_block[i].imag() - new_block[i].imag()) <= + SIN_COS_TOLERANCE); + max_error = max_d(max_error, ref_block[i].imag() - new_block[i].imag()); + } + BOOST_CHECK(std::abs(ref_nco.get_phase() - new_nco.get_phase()) <= SIN_COS_TOLERANCE); + // printf ("Fxpt max error %.9f, max phase error %.9f\n", max_error, + // max_phase_error); + delete[] ref_block; + delete[] new_block; } - |