annotate maths/CosineDistance.cpp @ 74:769da847732b

* A few fixes prompted by vamp-plugin-tester
author cannam
date Mon, 08 Jun 2009 12:30:20 +0000
parents dfe38135e4c7
children e5907ae6de17
rev   line source
cannam@31 1 /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*- vi:set ts=8 sts=4 sw=4: */
cannam@31 2
cannam@31 3 /*
cannam@31 4 QM DSP Library
cannam@31 5
cannam@31 6 Centre for Digital Music, Queen Mary, University of London.
cannam@31 7 This file copyright 2008 Kurt Jacobson.
cannam@31 8 All rights reserved.
cannam@31 9 */
cannam@31 10
cannam@31 11 #include "CosineDistance.h"
cannam@31 12
cannam@31 13 #include <iostream>
cannam@31 14 #include <limits>
cannam@31 15
cannam@31 16 using std::cerr;
cannam@31 17
cannam@31 18 double CosineDistance::distance(const vector<double> &v1,
cannam@31 19 const vector<double> &v2)
cannam@31 20 {
cannam@31 21 dist = 1.0; dDenTot = 0; dDen1 = 0; dDen2 = 0; dSum1 =0;
cannam@74 22 double small = 1e-20;
cannam@31 23
cannam@31 24 //check if v1, v2 same size
cannam@31 25 if (v1.size() != v2.size())
cannam@31 26 {
cannam@31 27 cerr << "CosineDistance::distance: ERROR: vectors not the same size\n";
cannam@31 28 return 1.0;
cannam@31 29 }
cannam@31 30 else
cannam@31 31 {
cannam@31 32 for(int i=0; i<v1.size(); i++)
cannam@31 33 {
cannam@31 34 dSum1 += v1[i]*v2[i];
cannam@31 35 dDen1 += v1[i]*v1[i];
cannam@31 36 dDen2 += v2[i]*v2[i];
cannam@31 37 }
cannam@74 38 dDenTot = sqrt(fabs(dDen1*dDen2)) + small;
cannam@31 39 dist = 1-((dSum1)/dDenTot);
cannam@31 40 return dist;
cannam@31 41 }
cannam@31 42 }