c@170: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
c@170: /*
c@170:     Constant-Q library
c@170:     Copyright (c) 2013-2014 Queen Mary, University of London
c@170: 
c@170:     Permission is hereby granted, free of charge, to any person
c@170:     obtaining a copy of this software and associated documentation
c@170:     files (the "Software"), to deal in the Software without
c@170:     restriction, including without limitation the rights to use, copy,
c@170:     modify, merge, publish, distribute, sublicense, and/or sell copies
c@170:     of the Software, and to permit persons to whom the Software is
c@170:     furnished to do so, subject to the following conditions:
c@170: 
c@170:     The above copyright notice and this permission notice shall be
c@170:     included in all copies or substantial portions of the Software.
c@170: 
c@170:     THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
c@170:     EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
c@170:     MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
c@170:     NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY
c@170:     CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF
c@170:     CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
c@170:     WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
c@170: 
c@170:     Except as contained in this notice, the names of the Centre for
c@170:     Digital Music; Queen Mary, University of London; and Chris Cannam
c@170:     shall not be used in advertising or otherwise to promote the sale,
c@170:     use or other dealings in this Software without prior written
c@170:     authorization.
c@170: */
c@170: 
c@170: #include "Pitch.h"
c@170: 
c@170: #include <math.h>
c@170: 
c@170: float
c@170: Pitch::getFrequencyForPitch(int midiPitch,
c@170: 			    float centsOffset,
c@170: 			    float concertA)
c@170: {
c@170:     float p = float(midiPitch) + (centsOffset / 100);
c@170:     return concertA * powf(2.0, (p - 69.0) / 12.0);
c@170: }
c@170: 
c@170: int
c@170: Pitch::getPitchForFrequency(float frequency,
c@170: 			    float *centsOffsetReturn,
c@170: 			    float concertA)
c@170: {
c@170:     float p = 12.0 * (log(frequency / (concertA / 2.0)) / log(2.0)) + 57.0;
c@170: 
c@170:     int midiPitch = int(p + 0.00001);
c@170:     float centsOffset = (p - midiPitch) * 100.0;
c@170: 
c@170:     if (centsOffset >= 50.0) {
c@170: 	midiPitch = midiPitch + 1;
c@170: 	centsOffset = -(100.0 - centsOffset);
c@170:     }
c@170:     
c@170:     if (centsOffsetReturn) *centsOffsetReturn = centsOffset;
c@170:     return midiPitch;
c@170: }
c@170: