matrixvector.h

一种聚类算法,名字是cocluster

/*
  MatrixVector.h
    Subroutines to handle matrix and vector

    Copyright (c) 2005, 2006
              by Hyuk Cho
    Copyright (c) 2003, 2004
    	      by Hyuk Cho, Yuqiang Guan, and Suvrit Sra
                {hyukcho, yguan, suvrit}@cs.utexas.edu
*/


#if !defined(_MATRIX_VECTOR_H_)
#define _MATRIX_VECTOR_H_


#include "Matrix.h"

//#include "RandomGenerator.h"
// document encoding 
#define NORM_TERM_FREQ			1 // normalized term frequency
#define NORM_TERM_FREQ_INV_DOC_FREQ	2 // normalized term frequency-inverse
					  // document frequency
#define Power_Method_Epsilon            0.0001
/*

// Encode documents using the specified encoding scheme.
// The elements of the input matrix contains the number of occurrences of a word
// in a document.
void encode_mat(SparseMatrix *mat, int scheme = NORM_TERM_FREQ);
void encode_mat(DenseMatrix *mat, int scheme = NORM_TERM_FREQ);
void dmatvec(int m, int n, double **a, double *x, double *y);
void dmatvecat(int m, int n, double **a, double *x, double *y);
void dqrbasis( int m, int n, double **a, double **q , double *work);
double dvec_l2normsq( int dim, double *v );
void dvec_l2normalize( int dim, double *v );
void dvec_scale( double alpha, int n, double *v );
*/

void average_vec(double vec[], int n, int num);
double norm_2 (double vec[], int n);
double norm_1(double *x, int n);
double KL_norm (double vec[], int n);
double euclidian_distance(double *v1, double *v2, int n);
double Kullback_leibler(double *x, double *y, int n);
double normalize_vec(double *vec, int n);
double normalize_vec_1(double *vec, int n);
double dot_mult(double *v1, double *v2, int n);
void Ax (double **A, double *x, int dim1, int dim2, double *result);
void Ax(double *vals, int *rowinds, int *colptrs, int dim1, int dim2, double *x, double *result);
double x_dot_y(double *x, double *y, int dim1);
void power_method(double **A_t_A, int dim, double * CV, double *init, double & Lamda);
void power_method(double *vals, int *rowinds, int *colptrs, int n_col, double * CV, double *init, double & Lamda);
double mutual_info(double ** matrix, int r, int c);
double mutual_info(double ** matrix, int r, int c, double*);


#endif // !defined(_MATRIX_VECTOR_H_)