cannam@31: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
cannam@31: 
cannam@31: /*
cannam@31:     QM DSP Library
cannam@31: 
cannam@31:     Centre for Digital Music, Queen Mary, University of London.
cannam@31:     This file copyright 2008 Kurt Jacobson.
cannam@31:     All rights reserved.
cannam@31: */
cannam@31: 
cannam@31: #include "CosineDistance.h"
cannam@31: 
cannam@31: #include <iostream>
cannam@31: #include <limits>
cannam@31: 
cannam@31: using std::cerr;
cannam@31: 
cannam@31: double CosineDistance::distance(const vector<double> &v1,
cannam@31:                                 const vector<double> &v2)
cannam@31: {
cannam@31:     dist = 1.0; dDenTot = 0; dDen1 = 0; dDen2 = 0; dSum1 =0;
cannam@31: 
cannam@31:     //check if v1, v2 same size
cannam@31:     if (v1.size() != v2.size())
cannam@31:     {
cannam@31:         cerr << "CosineDistance::distance: ERROR: vectors not the same size\n";
cannam@31:         return 1.0;
cannam@31:     }
cannam@31:     else
cannam@31:     {
cannam@31:         for(int i=0; i<v1.size(); i++)
cannam@31:         {
cannam@31:             dSum1 += v1[i]*v2[i];
cannam@31:             dDen1 += v1[i]*v1[i];
cannam@31:             dDen2 += v2[i]*v2[i];
cannam@31:         }
cannam@31:         dDenTot = sqrt(fabs(dDen1*dDen2));
cannam@31:         if(dDenTot == 0)
cannam@31:         {
cannam@31:             cerr << "CosineDistance::distance: WARNING: dividing by zero in cosine dist\n";
cannam@31:             return 1.0;
cannam@31:         }
cannam@31: 	
cannam@31:         dist = 1-((dSum1)/dDenTot);
cannam@31:         return dist;
cannam@31:     }
cannam@31: }