matthiasm@0: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
matthiasm@0: /*
matthiasm@0:     This file is Copyright (c) 2012 Chris Cannam
matthiasm@0:   
matthiasm@0:     Permission is hereby granted, free of charge, to any person
matthiasm@0:     obtaining a copy of this software and associated documentation
matthiasm@0:     files (the "Software"), to deal in the Software without
matthiasm@0:     restriction, including without limitation the rights to use, copy,
matthiasm@0:     modify, merge, publish, distribute, sublicense, and/or sell copies
matthiasm@0:     of the Software, and to permit persons to whom the Software is
matthiasm@0:     furnished to do so, subject to the following conditions:
matthiasm@0: 
matthiasm@0:     The above copyright notice and this permission notice shall be
matthiasm@0:     included in all copies or substantial portions of the Software.
matthiasm@0: 
matthiasm@0:     THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
matthiasm@0:     EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
matthiasm@0:     MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
matthiasm@0:     NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS BE LIABLE FOR
matthiasm@0:     ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF
matthiasm@0:     CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
matthiasm@0:     WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
matthiasm@0: */
matthiasm@0: 
matthiasm@0: #ifndef _YINUTIL_H_
matthiasm@0: #define _YINUTIL_H_
matthiasm@0: 
matthiasm@0: #include "vamp-sdk/FFT.h"
matthiasm@0: #include "MeanFilter.h"
matthiasm@0: 
matthiasm@0: #include <cmath>
matthiasm@0: 
matthiasm@0: #include <iostream>
matthiasm@0: #include <vector>
matthiasm@0: #include <exception>
matthiasm@0: 
matthiasm@0: using std::vector;
matthiasm@0: 
matthiasm@0: class YinUtil
matthiasm@0: {
matthiasm@0: public:
matthiasm@0:     static double sumSquare(const double *in, const size_t startInd, const size_t endInd);
matthiasm@0:     static void difference(const double *in, double *yinBuffer, const size_t yinBufferSize);
matthiasm@0:     static void fastDifference(const double *in, double *yinBuffer, const size_t yinBufferSize);
matthiasm@0:     static void cumulativeDifference(double *yinBuffer, const size_t yinBufferSize);
matthiasm@0:     static int absoluteThreshold(const double *yinBuffer, const size_t yinBufferSize, const double thresh);
matthiasm@0:     static vector<double> yinProb(const double *yinBuffer, const size_t prior, const size_t yinBufferSize);
matthiasm@0:     static double parabolicInterpolation(const double *yinBuffer, const size_t tau,
matthiasm@0:                                          const size_t yinBufferSize);
matthiasm@0: };
matthiasm@0: 
matthiasm@0: #endif