modifyDistanceMatrix.cc File Reference

#include <iterator>
#include <fstream>
#include <iostream>
#include <cmath>
#include <cassert>
#include <string>
#include <map>
#include <functional>

Go to the source code of this file.

Typedefs
typedef map< double, int >	CMap
Functions
void	calculateSums (float inFloat)
void	vectorToMap (float *inBegin, float *inEnd, CMap &outMap)
double	rankifyMatrix (float *inMatrix, int inSize)
double	rescale (float inFloat)
void	changeSign (float *inMatrix, int inSize)
void	symmetrify (float *inMatrix, int inSize)
	main (int argc, char **argv)
Variables
int	gCount = 0
double	gSum = 0
double	gSumOfSquares = 0

---

Typedef Documentation

typedef map<double,int> CMap

Definition at line 31 of file modifyDistanceMatrix.cc.

---

Function Documentation

void calculateSums

(

float

inFloat

)

Definition at line 25 of file modifyDistanceMatrix.cc.

References gCount, gSum, and gSumOfSquares.

Referenced by main().

                                 {
 gCount       ++;
 gSum         += inFloat;
 gSumOfSquares+= inFloat * inFloat;
}

void changeSign	(	float *	inMatrix,
		int	inSize
	)

Definition at line 100 of file modifyDistanceMatrix.cc.

Referenced by main().

               {
  inSize=inSize*inSize;

  for(float* p=inMatrix;p!=inMatrix+inSize;p++){
    *p = 1 - *p;
  }
}

main	(	int	argc,
		char **	argv
	)

Definition at line 128 of file modifyDistanceMatrix.cc.

References calculateSums(), changeSign(), rankifyMatrix(), and symmetrify().

                          {

  if(argc!=4){
    cout  << "usage: "
    << argv[0]
    << " <INFILE> <OUTFILE> <SIZE>"
    << endl;
    exit(1);
  }
  
  string lInName=argv[1];
  string lOutName=argv[2];
  int lSize=atoi(argv[3]);

  float* lMatrix=new float[lSize*lSize];
  assert(lMatrix);
  {
    ifstream lIn(lInName.c_str());
    
    while(lIn){
      float lBuffer;
      lIn.read((char*)&lBuffer,sizeof(lBuffer));
      calculateSums(lBuffer);
    }
  }
  cout << "until here "
       << flush 
       << endl;
  {
    ifstream lIn(lInName.c_str());
    ofstream lOut(lOutName.c_str());
    
    assert(lOut);

    int lCount=0;
    while(lIn){
      float lBuffer;
      lIn.read((char*)&lBuffer,
         sizeof(lBuffer));
      lMatrix[lCount++]=lBuffer;
    }

#ifdef VISTEX
     rankifyMatrix(lMatrix,
       lSize);
#endif

     //symmetrify the matrix
     symmetrify(lMatrix,
    lSize);
     
     //get from score to distance
     changeSign(lMatrix,
    lSize);


    lOut.write((char*)lMatrix,
         sizeof(float)
         *(lSize*lSize));
      
  }
}

double rankifyMatrix	(	float *	inMatrix,
		int	inSize
	)

Take the matrix make pairs of score and ID

Definition at line 45 of file modifyDistanceMatrix.cc.

References vectorToMap().

Referenced by main().

                    {
  for(int lLine=0;
      lLine<inSize;
      lLine++){
    //this map is sorted by rank inascending order
    CMap lMap;
    vectorToMap(inMatrix + lLine*inSize,
    inMatrix + (1+lLine)*inSize,
    lMap);
    //cout << endl;

    int lCount=inSize;
    for(CMap::const_iterator i=lMap.begin();
  i!=lMap.end();
  i++){
      //in the matrix now are the query ranks
      float l= float(--lCount)/inSize;
      if(l>1)
  l=1;
      if(l<0)
  l=0;
      inMatrix[(lLine*inSize)+i->second]=1-l;
    }
  }
}

double rescale

(

float

inFloat

)

Definition at line 75 of file modifyDistanceMatrix.cc.

References gCount, gSum, and gSumOfSquares.

                             {
  double lAverage=(gSum/gCount);
  double lAverageSquare=
    gSumOfSquares/gCount;
  double lSDev=
    sqrt(lAverageSquare
   -
   lAverage*lAverage);
  //...so from score to pseudo distance
  double lReturnValue=0.5-((inFloat-lAverage)/lSDev/3);

  if(lReturnValue>.9999)
    lReturnValue=1;
  if(lReturnValue<0.0001)
    lReturnValue=0;

  //  cout << "[" << lReturnValue << "]" <<flush;
  assert((lReturnValue<=1) && (lReturnValue>=0));

  return 
    lReturnValue;
}

void symmetrify	(	float *	inMatrix,
		int	inSize
	)

Definition at line 110 of file modifyDistanceMatrix.cc.

Referenced by main().

               {
  for(int i=0;
      i<inSize;
      i++){
    for(int j=0;
  j<inSize;
  j++){
      float l=min(inMatrix[i*inSize+j],
      inMatrix[j*inSize+i]);
      inMatrix[i*inSize+j]=l;
      inMatrix[j*inSize+i]=l;
    }

  }
};

void vectorToMap	(	float *	inBegin,
		float *	inEnd,
		CMap &	outMap
	)

Definition at line 33 of file modifyDistanceMatrix.cc.

Referenced by rankifyMatrix().

                  {
  int lCount=0;
  for(float* i=inBegin;
      i!=inEnd;
      i++){
    outMap.insert(make_pair(*i,lCount++));
    //    cout << "-" << flush;
  }
}

---

Variable Documentation

int gCount = 0

Definition at line 21 of file modifyDistanceMatrix.cc.

Referenced by calculateSums(), and rescale().

double gSum = 0

Definition at line 22 of file modifyDistanceMatrix.cc.

Referenced by calculateSums(), and rescale().

double gSumOfSquares = 0

Definition at line 23 of file modifyDistanceMatrix.cc.

Referenced by calculateSums(), and rescale().