//
// test floating point precision nco
//
void autotest_nco_basic() {
nco_crcf p = nco_crcf_create(LIQUID_NCO);
unsigned int i; // loop index
float s, c; // sine/cosine result
float tol=1e-4f; // error tolerance
float f=0.0f; // frequency to test
nco_crcf_set_phase( p, 0.0f);
CONTEND_DELTA( nco_crcf_cos(p), 1.0f, tol );
CONTEND_DELTA( nco_crcf_sin(p), 0.0f, tol );
nco_crcf_sincos(p, &s, &c);
CONTEND_DELTA( s, 0.0f, tol );
CONTEND_DELTA( c, 1.0f, tol );
nco_crcf_set_phase(p, M_PI/2);
CONTEND_DELTA( nco_crcf_cos(p), 0.0f, tol );
CONTEND_DELTA( nco_crcf_sin(p), 1.0f, tol );
nco_crcf_sincos(p, &s, &c);
CONTEND_DELTA( s, 1.0f, tol );
CONTEND_DELTA( c, 0.0f, tol );
// cycle through one full period in 64 steps
nco_crcf_set_phase(p, 0.0f);
f = 2.0f * M_PI / 64.0f;
nco_crcf_set_frequency(p, f);
for (i=0; i<128; i++) {
nco_crcf_sincos(p, &s, &c);
CONTEND_DELTA( s, sinf(i*f), tol );
CONTEND_DELTA( c, cosf(i*f), tol );
nco_crcf_step(p);
}
// double frequency: cycle through one full period in 32 steps
nco_crcf_set_phase(p, 0.0f);
f = 2.0f * M_PI / 32.0f;
nco_crcf_set_frequency(p, f);
for (i=0; i<128; i++) {
nco_crcf_sincos(p, &s, &c);
CONTEND_DELTA( s, sinf(i*f), tol );
CONTEND_DELTA( c, cosf(i*f), tol );
nco_crcf_step(p);
}
// destroy NCO object
nco_crcf_destroy(p);
}
// autotest helper function
// _theta : input phase
// _cos : expected output: cos(_theta)
// _sin : expected output: sin(_theta)
// _type : NCO type (e.g. LIQUID_NCO)
// _tol : error tolerance
void nco_crcf_phase_test(float _theta,
float _cos,
float _sin,
int _type,
float _tol)
{
// create object
nco_crcf nco = nco_crcf_create(_type);
// set phase
nco_crcf_set_phase(nco, _theta);
// compute cosine and sine outputs
float c = nco_crcf_cos(nco);
float s = nco_crcf_sin(nco);
// run tests
CONTEND_DELTA( c, _cos, _tol );
CONTEND_DELTA( s, _sin, _tol );
// destroy object
nco_crcf_destroy(nco);
}
void quiet_beacon_encoder_emit(quiet_beacon_encoder *e, quiet_beacon_sample_t *buf, size_t buflen) {
for (size_t i = 0; i < buflen; i++) {
buf[i] = nco_crcf_cos(e->nco) * e->gain;
nco_crcf_step(e->nco);
}
}
