/****************************************************************File Name:   vector.CAuthor: Tian Zhang, CS Dept., Univ. of Wisconsin-Madison, 1995               Copyright(c) 1995 by Tian Zhang                   All Rights ReservedPermission to use, copy and modify this software must be grantedby the author and provided that the above copyright notice appear in all relevant copies and that both that copyright notice and this permission notice appear in all relevant supporting documentations. Comments and additions may be sent the author at zhang@cs.wisc.edu.******************************************************************/#include "global.h"#include "util.h"#include "vector.h"	Vector::Vector() {dim = 0; value=NULL;}void Vector::Init(short d) {  dim = d;  if (0 == d)    value = NULL;  else    {      value = new double[dim];      memset(value, 0, dim*sizeof(double));    }}void Vector::Reset() {	memset(value, 0, dim*sizeof(double));	}Vector::Vector(const Vector& r) {	dim = r.dim;	memcpy(value,r.value,dim*sizeof(double));	}Vector::~Vector() {if (value!=NULL) delete [] value;}double Vector::Value(short i) const {return value[i];}void Vector::operator=(const Vector& r) {	if (this != &r) {		dim = r.dim;		memcpy(value,r.value,dim*sizeof(double));		}	}	void Vector::operator=(const int val) {	for (short i=0; i<dim; i++)		value[i] = val;}void Vector::operator=(const float val) {	for (short i=0; i<dim; i++) 		value[i]=val;	}void Vector::operator=(const double val) {	for (short i=0; i<dim; i++)		value[i]=val;	}// euclidian distance : sum of squares on each dimensiondouble Vector::operator||(const Vector& v2) const {	double tmp = 0;	for (short i=0; i<dim; i++) 		tmp += (value[i]-v2.value[i])*(value[i]-v2.value[i]);	return tmp;	}	// manhatan distance : sum of absolute values on each dimensiondouble Vector::operator^(const Vector& v2) const {	double tmp = 0;	for (short i=0; i<dim; i++) 		tmp += fabs(value[i]-v2.value[i]);	return tmp;	}void Vector::operator+=(const Vector &v2) {	for (short i=0; i<dim; i++) 		value[i] += v2.value[i];	}void Vector::operator-=(const Vector &v2) {	for (short i=0; i<dim; i++) 		value[i] -= v2.value[i];	}void Vector::operator*=(double cnt) {	for (short i=0; i<dim; i++) 		value[i] *= cnt;	}void Vector::operator/=(double cnt) {	for (short i=0; i<dim; i++) 		value[i] /= cnt;	}// dot productdouble Vector::operator&&(const Vector& v2) const {	double tmp = 0;	for (short i=0; i<dim; i++) 		tmp += value[i]*v2.value[i];	return tmp;	}// time by Nvoid Vector::Mul(const Vector& v1, double N) {	for (short i=0;i<dim;i++)		value[i]=v1.value[i]*N;	}	// divide by Nvoid Vector::Div(const Vector& v1, double N) {	for (short i=0;i<dim;i++)		value[i]=v1.value[i]/N;	}// plusvoid Vector::Add(const Vector& v1, const Vector& v2) {	for (short i=0; i<dim; i++)		value[i]=v1.value[i]+v2.value[i];	}// minusvoid Vector::Sub(const Vector& v1, const Vector& v2) {	for (short i=0; i<dim; i++)		value[i]=v1.value[i]-v2.value[i];	}void Vector::AddSqr(const Vector& v) {	for (short i=0;i<dim;i++)		value[i]+=v.value[i]*v.value[i];	}// weighting and movingvoid Vector::Transform(const Vector &W,const Vector &M) {	for (short i=0; i<dim; i++)        	value[i]=(value[i]-M.value[i])*W.value[i];	}// reverse weighting and reverse movingvoid Vector::Reverse_Transform(const Vector &W, const Vector &M) {	for (short i=0; i<dim; i++)		value[i]=value[i]/W.value[i]+M.value[i];	}// input a vectoristream &operator>>(istream &fi, Vector &r) {	for (short i=0;i<r.dim;i++) 		fi >> r.value[i];	return fi;	}ifstream &operator>>(ifstream &fi, Vector &r) {	for (short i=0;i<r.dim;i++) 		fi >> r.value[i];	return fi;	}	// output a vectorostream &operator<<(ostream &fo, const Vector &r) {	for (short i=0; i<r.dim; i++)		fo << r.value[i] << ' '; 	return fo;	}ofstream &operator<<(ofstream &fo, const Vector &r) {	for (short i=0; i<r.dim; i++)		fo << r.value[i] << ' '; 	return fo;	}