/*******************************************************************************
VectorMax */
float
VectorMax(const float* src, unsigned length)
{
float res = FLT_MIN;
#ifdef __APPLE__
// Use the Accelerate framework if we have it
vDSP_maxv(src, 1, &res, length);
#else
for (unsigned i = 0; i < length; ++i)
{
if (src[i] > res)
{
res = src[i];
}
}
#endif
return res;
}
smpl_t
fvec_max (fvec_t * s)
{
#ifndef HAVE_ACCELERATE
uint_t j;
smpl_t tmp = 0.0;
for (j = 0; j < s->length; j++) {
tmp = (tmp > s->data[j]) ? tmp : s->data[j];
}
#else
smpl_t tmp = 0.;
#if !HAVE_AUBIO_DOUBLE
vDSP_maxv(s->data, 1, &tmp, s->length);
#else
vDSP_maxvD(s->data, 1, &tmp, s->length);
#endif
#endif
return tmp;
}
bool ofxAudioUnitFftNode::getAmplitude(std::vector<float> &outAmplitude)
{
getSamplesFromChannel(_sampleBuffer, 0);
// return empty if we don't have enough samples yet
if(_sampleBuffer.size() < _N) {
outAmplitude.clear();
return false;
}
// normalize input waveform
if(_outputSettings.normalizeInput) {
float timeDomainMax;
vDSP_maxv(&_sampleBuffer[0], 1, &timeDomainMax, _N);
vDSP_vsdiv(&_sampleBuffer[0], 1, &timeDomainMax, &_sampleBuffer[0], 1, _N);
}
PerformFFT(&_sampleBuffer[0], _window, _fftData, _fftSetup, _N);
// get amplitude
vDSP_zvmags(&_fftData, 1, _fftData.realp, 1, _N/2);
// normalize magnitudes
float two = 2.0;
vDSP_vsdiv(_fftData.realp, 1, &two, _fftData.realp, 1, _N/2);
// scale output according to requested settings
if(_outputSettings.scale == OFXAU_SCALE_LOG10) {
for(int i = 0; i < (_N / 2); i++) {
_fftData.realp[i] = log10f(_fftData.realp[i] + 1);
}
} else if(_outputSettings.scale == OFXAU_SCALE_DECIBEL) {
float ref = 1.0;
vDSP_vdbcon(_fftData.realp, 1, &ref, _fftData.realp, 1, _N / 2, 1);
float dbCorrectionFactor = 0;
switch (_outputSettings.window) {
case OFXAU_WINDOW_HAMMING:
dbCorrectionFactor = DB_CORRECTION_HAMMING;
break;
case OFXAU_WINDOW_HANNING:
dbCorrectionFactor = DB_CORRECTION_HAMMING;
break;
case OFXAU_WINDOW_BLACKMAN:
dbCorrectionFactor = DB_CORRECTION_HAMMING;
break;
}
vDSP_vsadd(_fftData.realp, 1, &dbCorrectionFactor, _fftData.realp, 1, _N / 2);
}
// restrict minimum to 0
if(_outputSettings.clampMinToZero) {
float min = 0.0;
float max = INFINITY;
vDSP_vclip(_fftData.realp, 1, &min, &max, _fftData.realp, 1, _N / 2);
}
// normalize output between 0 and 1
if(_outputSettings.normalizeOutput) {
float max;
vDSP_maxv(_fftData.realp, 1, &max, _N / 2);
if(max > 0) {
vDSP_vsdiv(_fftData.realp, 1, &max, _fftData.realp, 1, _N / 2);
}
}
outAmplitude.assign(_fftData.realp, _fftData.realp + _N/2);
return true;
}
