Vector Class Reference

#include <Vector.h>

Collaboration diagram for Vector:

Classes
struct	VectorDataTag

Public Types
typedef struct Vector::VectorDataTag	VectorData

Public Member Functions
void	prepareExtend (int new_extention)
void	alloc (int n, int ext)
	Vector (int _n=0)
	Vector (int _n, int _ext)
	Vector (int _n, double *dd, char externalData=0)
	Vector (char *filename)
	Vector (char *line, int guess_on_size)
	Vector (Vector a, Vector b)
char *	getFromLine (char *line)
void	extend ()
void	setSize (int _n)
void	exactshape ()
void	print ()
void	save (char *filename, char ascii)
void	save (FILE *f, char ascii)
void	appendToMatrixFile (char *saveFileName, char **names=NULL)
void	setExternalData (int _n, double *dd)
Vector	clone ()
void	copyFrom (Vector r, int _n=0)
	Vector (const Vector &P)
Vector &	operator= (const Vector &P)
void	destroyCurrentBuffer ()
	~Vector ()
unsigned	sz ()
int	operator== (const Vector Q)
int	equals (const Vector Q)
	operator double * () const
void	setPart (int i, Vector v, int n=0, int ii=0)
double	euclidianNorm ()
double	L1Norm ()
double	LnftyNorm ()
double	euclidianDistance (Vector v)
double	L1Distance (Vector v)
double	LnftyDistance (Vector v)
double	square ()
void	multiply (double a)
void	multiply (Vector R, double a)
void	zero (int _i=0, int _n=0)
void	set (double dd)
void	shift (int s)
double	scalarProduct (Vector v)
void	oneByOneMutiply (Vector r)
void	oneByOneInvert ()
double	mmin ()
double	mmax ()
bool	isNull ()
Vector	operator- (Vector v)
Vector	operator+ (Vector v)
Vector	operator-= (Vector v)
Vector	operator+= (Vector v)
void	addInPlace (double a, Vector v)
void	addInPlace (double a, int i, Matrix m)
void	transposeAndMultiply (Vector vR, Matrix M)
void	diagonalizeAndMultiply (Matrix M)
void	permutIn (Vector vR, VectorInt viP)
void	permutOut (Vector vR, VectorInt viP)

Public Attributes
VectorData *	d = nullptr

Static Public Attributes
static Vector	emptyVector

Detailed Description

Definition at line 37 of file Vector.h.

Member Typedef Documentation

VectorData

typedef struct Vector::VectorDataTag Vector::VectorData

Constructor & Destructor Documentation

Vector() [1/7]

Vector::Vector

(

int

_n = 0

)

Definition at line 49 of file Vector.cpp.

{
    alloc(n,n);
    zero();
}

Vector() [2/7]

Vector::Vector	(	int	_n,
		int	_ext
	)

Definition at line 55 of file Vector.cpp.

{
    alloc(n,ext); 
    zero();
}

Vector() [3/7]

Vector::Vector	(	int	_n,
		double *	dd,
		char	externalData = 0
	)

Definition at line 61 of file Vector.cpp.

{
    alloc(n,n);
    if (dd) 
    {
        if (_exte) { d->externalData=1; d->p=dd; }
        else memcpy(d->p,dd,n*sizeof(double));
    }
    else zero();
}

Vector() [4/7]

Vector::Vector

(

char *

filename

)

Definition at line 431 of file Vector.cpp.

{
    unsigned _n;
    std::ifstream ifp(filename, std::ios::in | std::ios::binary);
    ifp.read(reinterpret_cast<char*>(&_n), sizeof(int));
    alloc(_n,_n);
    ifp.read(reinterpret_cast<char*>(d->p), d->n*sizeof(double));
    ifp.close();
}

Vector() [5/7]

Vector::Vector	(	char *	line,
		int	guess_on_size
	)

Definition at line 543 of file Vector.cpp.

{
    char *tline=line,*oldtline=NULL;

    if (gn==0) 
    {
        while ((*tline!=EOL1)&&(*tline!=EOL2))
        {
            while ((*tline==' ')|| 
                (*tline=='\t'))tline++;
            if ((*tline==EOL1)||(*tline==EOL2)||(*tline==0)) break;
            while (isANumber(*tline)) tline++;
            gn++;
            if (oldtline==tline)
            {
                alloc(gn-1,gn-1);
                return;
                //tline[10]=0;
                //printf("Error in simulation output file. The full line is:\n"
                //       "%s\n"
                //       "There is an error here:\n"
                //       "%s\n",line,tline);
            }
            oldtline=tline;
        };
    };

    if (gn==0) { alloc(0,0); return; };
    alloc(gn,gn);
    getFromLine(line);
}

Vector() [6/7]

Vector::Vector	(	Vector	a,
		Vector	b
	)

Definition at line 72 of file Vector.cpp.

{
    int n=a.sz()+b.sz();
    alloc(n,n);
    memcpy(d->p       ,a,a.sz()*sizeof(double));
    memcpy(d->p+a.sz(),b,b.sz()*sizeof(double));
}

Vector() [7/7]

Vector::Vector

(

const Vector &

P

)

Definition at line 188 of file Vector.cpp.

{
    // shallow copy
    d=A.d;
    (d->ref_count)++ ;
}

~Vector()

Vector::~Vector

(

)

Definition at line 172 of file Vector.cpp.

{
    destroyCurrentBuffer();
}

Member Function Documentation

addInPlace() [1/2]

void Vector::addInPlace	(	double	a,
		Vector	v
	)

Definition at line 416 of file Vector.cpp.

{
    int ni=sz();
    double *xp1=d->p, *xp2=v.d->p;
    if (a==1.0)  while (ni--) *(xp1++)+=     *(xp2++);
    else         while (ni--) *(xp1++)+=a * (*(xp2++));
}

addInPlace() [2/2]

void Vector::addInPlace	(	double	a,
		int	i,
		Matrix	m
	)

Definition at line 424 of file Vector.cpp.

{
    int ni=sz();
    double *xp1=d->p, *xp2=((double**)m)[i];
    while (ni--) *(xp1++)+=a * (*(xp2++));
}

alloc()

void Vector::alloc	(	int	n,
		int	ext
	)

Definition at line 36 of file Vector.cpp.

{
    d=new VectorData();
    d->n=_n;
    d->extention=_ext;
    d->ref_count=1;

    if (_ext==0) { d->p=NULL; return; };

    d->p=(double*)malloc(_ext*sizeof(double)); 
    if (d->p==NULL) { printf("memory allocation error\n"); getchar(); exit(253); }
}

appendToMatrixFile()

void Vector::appendToMatrixFile	(	char *	saveFileName,
		char **	names = NULL
	)

Definition at line 612 of file Vector.cpp.

{
    FILE *f=fopen(saveFileName,"rb+");
    if (f==NULL)
    {
        Matrix t(*this);
        if (names) t.setColNames(names);
        t.save(saveFileName,0);
        return;
    }
    int nc=sz();
    fseek(f,0,SEEK_END);
    unsigned l=ftell(f);
    if (l==0)
    {
        Matrix t(*this);
        if (names) t.setColNames(names);
        t.save(saveFileName,0);
        return;
    }
    fwrite(d->p,sizeof(double)*nc,1,f);
    unsigned nlfile;
    fseek(f,13,SEEK_SET);
    if (fread(&nlfile,sizeof(int),1,f) != sizeof(int)) {
        std::cerr << "error while appending to file " << saveFileName << std::endl; exit(255);    
    } 
    nlfile++;
    fseek(f,13,SEEK_SET);
    fwrite(&nlfile, sizeof(unsigned), 1, f);
    fclose(f);
}

clone()

Vector Vector::clone

(

)

Definition at line 207 of file Vector.cpp.

{
    // a deep copy
    Vector r(sz());
    r.copyFrom(*this);
    return r;
}

copyFrom()

void Vector::copyFrom	(	Vector	r,
		int	_n = 0
	)

Definition at line 215 of file Vector.cpp.

{
    if (_n==0) _n=r.sz();
    setSize(_n);
    if (_n) memcpy(d->p,r.d->p,_n*sizeof(double));
}

destroyCurrentBuffer()

void Vector::destroyCurrentBuffer

(

)

Definition at line 177 of file Vector.cpp.

{
    if (!d) return;
    (d->ref_count) --;
    if (d->ref_count==0)
    {
        if ((d->p)&&(!d->externalData)) free(d->p);
        delete d;
    };
}

diagonalizeAndMultiply()

void Vector::diagonalizeAndMultiply

(

Matrix

M

)

Definition at line 316 of file Vector.cpp.

{
    int i,j,nl=M.nLine(),nc=M.nColumn();
    if ((int)sz()!=nl)
    {
        printf("(matrix_diagonal * matrix) error");
        getchar(); exit(249);
    }
    double **p1=M,*p2=(*this);

    for (i=0; i<nl; i++) 
        for (j=0; j<nc; j++)
            p1[i][j]*=p2[i];

}

equals()

int Vector::equals

(

const Vector

Q

)

Definition at line 135 of file Vector.cpp.

{
  if (Q.d==d) return 1;
  if (Q.d==emptyVector.d)
  {
      double *cP=d->p;
      int i=sz();
      while (i--) if (*(cP++)) return 0;
      return 1;
  }

  if (sz() != Q.sz()) return 0;

  double *cP = d->p, *cQ = Q.d->p;
  int i = sz();

  while( i-- )
  {
    if (*cP!=*cQ) return 0;
    cP++; cQ++;
  }

  return 1;
}

euclidianDistance()

double Vector::euclidianDistance

(

Vector

v

)

Definition at line 244 of file Vector.cpp.

{
    Vector t=(*this)-v;
    return ::euclidianNorm(sz(), t.d->p);
/*
    double *xp1=d->p, *xp2=v.d->p, sum=0;
    int ni=sz();
    while (ni--) sum+=sqr(*(xp1++)-*(xp2++));
    return sqrt(sum);
*/
}

euclidianNorm()

double Vector::euclidianNorm

(

)

Definition at line 222 of file Vector.cpp.

{
    return ::euclidianNorm(sz(), d->p);
}

exactshape()

void Vector::exactshape

(

)

Definition at line 125 of file Vector.cpp.

{
    if (d->extention!=d->n)
    {
        d->p=(double*)realloc(d->p,d->n*sizeof(double));
        if (d->p==NULL) { printf("memory allocation error\n"); getchar(); exit(253); }
        d->extention=d->n;
    };
}

extend()

void Vector::extend

(

)

Definition at line 119 of file Vector.cpp.

{
    d->n++;
    if (d->n>d->extention) prepareExtend(d->extention+100);
}

getFromLine()

char * Vector::getFromLine

(

char *

line

)

Definition at line 575 of file Vector.cpp.

{
    double *dp=d->p;
    int n=sz(),k;
    char *tline=line, *oldtline; 
    for (k=0; k<n; k++) 
    {
        while ((*tline==' ')|| 
                (*tline=='\t'))tline++;
        if ((*tline==EOL1)||(*tline==EOL2)) 
        {
            setSize(k);
            return tline;
        }
        oldtline=tline;
        while(isANumber(*tline)) tline++;
        if (!isANumber(*(tline-1)))
        {
            setSize(k);
            return tline;
            //tline[10]=0;
            //printf( "Error in simulation output file. The full line is:\n"
            //        "%s\n"
            //        "There is an error here:\n"
            //        "%s\n",line,tline);
        }
        if (oldtline==tline)
        {
            setSize(k);
            return tline;
        };
        if (*tline) { *tline='\0'; tline++; }
        dp[k]=atof(oldtline);
    }
    return tline;
}

isNull()

bool Vector::isNull

(

)

Definition at line 358 of file Vector.cpp.

{
    double *xp=d->p;
    int ni=sz();
    while (ni--) if (*(xp++)!=0) return false;
    return true;
}

L1Distance()

double Vector::L1Distance

(

Vector

v

)

Definition at line 272 of file Vector.cpp.

{
    if (sz()==0) return 0;
    double *xp1=d->p, *xp2=v.d->p, sum=0;
    int ni=sz();
    while (ni--) sum+=condorAbs(*(xp1++)-*(xp2++));
    return sum;
}

L1Norm()

double Vector::L1Norm

(

)

Definition at line 227 of file Vector.cpp.

{
    if (sz()==0) return 0;
    double *x=d->p, sum=0;
    int ni=sz();
    while (ni--) sum+=condorAbs(*(x++));
    return sum;
}

LnftyDistance()

double Vector::LnftyDistance

(

Vector

v

)

Definition at line 256 of file Vector.cpp.

{
    double *xp1=d->p, *xp2=v.d->p, sum=-1.0;
    int ni=sz();
    while (ni--) sum=::mmax(sum, condorAbs(*(xp1++)-*(xp2++)));
    return sum;
}

LnftyNorm()

double Vector::LnftyNorm

(

)

Definition at line 264 of file Vector.cpp.

{
    double *xp1=d->p, sum=-1.0;
    int ni=sz();
    while (ni--) sum=::mmax(sum, condorAbs(*(xp1++)));
    return sum;
}

mmax()

double Vector::mmax

(

)

Definition at line 349 of file Vector.cpp.

{
    if (sz()==0) return 0;
    double *xp=d->p, m=-INF;
    int ni=sz();
    while (ni--) m=::mmax(m,*(xp++));
    return m;
}

mmin()

double Vector::mmin

(

)

Definition at line 340 of file Vector.cpp.

{
    if (sz()==0) return 0;
    double *xp=d->p, m=INF;
    int ni=sz();
    while (ni--) m=::mmin(m,*(xp++));
    return m;
}

multiply() [1/2]

void Vector::multiply

(

double

a

)

Definition at line 281 of file Vector.cpp.

{
    double *xp=d->p;
    int ni=sz();
    while (ni--) *(xp++)*=a;
}

multiply() [2/2]

void Vector::multiply	(	Vector	R,
		double	a
	)

Definition at line 288 of file Vector.cpp.

{
    int ni=sz();
    R.setSize(ni);
    double *xs=d->p, *xd=R;
    while (ni--) *(xd++)=a * (*(xs++));
}

oneByOneInvert()

void Vector::oneByOneInvert

(

)

Definition at line 383 of file Vector.cpp.

{
    int i=sz();
    double *xb=(*this);
    while (i--) xb[i]=1/xb[i];
}

oneByOneMutiply()

void Vector::oneByOneMutiply

(

Vector

r

)

Definition at line 376 of file Vector.cpp.

{
    int i=::mmin(sz(),rescaling.sz());
    double *xb=(*this), *r=rescaling;
    while (i--) xb[i]*=r[i];
}

operator double *()

Vector::operator double * ( ) const

inline

Definition at line 83 of file Vector.h.

{ if (d) return d->p; else return NULL; };

operator+()

Vector Vector::operator+

(

Vector

v

)

Definition at line 390 of file Vector.cpp.

{
    int ni=sz();
    Vector r(sz());
    double *xp1=r.d->p, *xp2=v.d->p, *xp3=d->p;
    while (ni--) 
        *(xp1++)+=*(xp3++)+*(xp2++);
    return r;
}

operator+=()

Vector Vector::operator+=

(

Vector

v

)

Definition at line 408 of file Vector.cpp.

{
    int ni=sz();
    double *xp1=d->p, *xp2=v.d->p;
    while (ni--) *(xp1++)+=*(xp2++);
    return *this;
}

operator-()

Vector Vector::operator-

(

Vector

v

)

Definition at line 366 of file Vector.cpp.

{
    int ni=sz();
    Vector r(sz());
    double *xp1=r.d->p, *xp2=v.d->p, *xp3=d->p;
    while (ni--) 
        *(xp1++)+=*(xp3++)-*(xp2++);
    return r;
}

operator-=()

Vector Vector::operator-=

(

Vector

v

)

Definition at line 400 of file Vector.cpp.

{
    int ni=sz();
    double *xp1=d->p, *xp2=v.d->p;
    while (ni--) *(xp1++)-=*(xp2++);
    return *this;
}

operator=()

Vector & Vector::operator=

(

const Vector &

P

)

Definition at line 195 of file Vector.cpp.

{
    // shallow copy
    if (this != &A)
    {
        destroyCurrentBuffer();
        d=A.d;
        (d->ref_count) ++ ;
    }
    return *this;
}

operator==()

int Vector::operator== ( const Vector  Q )

inline

Definition at line 81 of file Vector.h.

{ return d==Q.d; };

permutIn()

void Vector::permutIn	(	Vector	vR,
		VectorInt	viP
	)

Definition at line 501 of file Vector.cpp.

{
    int i,n=sz(), *ii=viP;
    if (!n) return;
    if (n!=viP.sz())
    {
        printf("error in permutation IN: sizes don't agree.\n"); getchar(); exit(255);
    }
    vR.setSize(n);
    double *ps=(*this), *pd=vR; // pointer source / destination
    for (i=0; i<n; i++)
//        *(pd++)=ps[ii[i]];
        pd[ii[i]]=*(ps++);
}

permutOut()

void Vector::permutOut	(	Vector	vR,
		VectorInt	viP
	)

Definition at line 516 of file Vector.cpp.

{
    int i,n=sz(), *ii=viP;
    if (!n) return;
    if (n!=viP.sz())
    {
        printf("error in permutation IN: sizes don't agree.\n"); getchar(); exit(255);
    }
    vR.setSize(n);
    double *ps=(*this), *pd=vR; // pointer source / destination
    for (i=0; i<n; i++)
//        pd[ii[i]]=*(ps++);
        *(pd++)=ps[ii[i]];
}

prepareExtend()

void Vector::prepareExtend

(

int

new_extention

)

Definition at line 99 of file Vector.cpp.

{
    if (d->extention<new_extention)
    {
        d->p=(double*)realloc(d->p,new_extention*sizeof(double));
        if (d->p==NULL) { printf("memory allocation error\n"); getchar(); exit(253); }

        // not really necessary (fill with zero's):
        memset(d->p+d->extention,0,(new_extention-d->extention)*sizeof(double));
        d->extention=new_extention;
    };
}

print()

void Vector::print

(

)

Definition at line 161 of file Vector.cpp.

{
    int N=sz();
    printf("[");
    if (!N || !d->p) { printf("]\n"); return; }

    double *up=d->p;
    while (--N) printf("%f,",*(up++));
    printf("%f]\n",*up);
}

save() [1/2]

void Vector::save	(	char *	filename,
		char	ascii
	)

Definition at line 441 of file Vector.cpp.

{
    FILE *f;
    if (ascii) f=fopen(filename,"w"); else f=fopen(filename,"wb");
    if (f==NULL)
    {
        printf("Cannot write to '%s'\n",filename);
        exit(255);
    }
    save(f,ascii);
    fclose(f);
}

save() [2/2]

void Vector::save	(	FILE *	f,
		char	ascii
	)

Definition at line 454 of file Vector.cpp.

{
    char header[32]="CONDORVBv1.0";
    if (ascii)
    {
        unsigned i;
        double *pp=d->p;
        if (ascii!=2) fprintf(f,"CONDORVAv1.0\n");
        if (sz())
        {
            for (i=0; i<sz()-1; i++) fprintf(f,"%.16e\t",pp[i]);
            fprintf(f,"%.16e\n",pp[i]);
        }
        return;
    }
    fwrite(header,sizeof(char),13,f);
    fwrite(&d->n, sizeof(int),1, f);
    fwrite(d->p, d->n*sizeof(double),1, f);
}

scalarProduct()

double Vector::scalarProduct

(

Vector

v

)

Definition at line 332 of file Vector.cpp.

{
    double *xp1=d->p, *xp2=v.d->p, sum=0;
    int ni=sz();
    while (ni--) { sum+=*(xp1++) * *(xp2++); };
    return sum;
}

set()

void Vector::set

(

double

dd

)

Definition at line 482 of file Vector.cpp.

{
    double *p=(*this);
    if (!p) return;
    int n=sz();
    while (n--) *(p++)=dd;
}

setExternalData()

void Vector::setExternalData	(	int	_n,
		double *	dd
	)

Definition at line 80 of file Vector.cpp.

{
    if ((d->extention==_n)||(!d->extention))
    {
        d->n=_n; d->extention=_n; d->externalData=1; d->p=dd; 
    } else
    {
        printf("do not use this function ('setExternalData'): it's too dangerous.\n");
        getchar(); exit(255);
    }
}

setPart()

void Vector::setPart	(	int	i,
		Vector	v,
		int	n = 0,
		int	ii = 0
	)

Definition at line 475 of file Vector.cpp.

{
    if (n==0) n=v.sz()-ii;
    n=::mmin((int)n,(int)sz()-i);
    memcpy(d->p+i, ((double*)v)+ii, n*sizeof(double));
}

setSize()

void Vector::setSize

(

int

_n

)

Definition at line 112 of file Vector.cpp.

{
    d->n=_n;
    if (_n==0) { if (d->p) free(d->p); d->p=NULL; d->extention=0; return; }
    prepareExtend(_n);
}

shift()

void Vector::shift

(

int

s

)

Definition at line 490 of file Vector.cpp.

{
    int n=sz();
    if (!n) return;
    double *ps=(*this), *pd=ps; // pointer source / destination
    if (s==0) return;
    if (s>0) { n-=s; pd+=s; } 
    else { n+=s; ps+=s; }
    memmove(pd,ps,n*sizeof(double));
}

square()

double Vector::square

(

)

Definition at line 236 of file Vector.cpp.

{
    double *xp=d->p, sum=0;
    int ni=sz();
    while (ni--) sum+=sqr(*(xp++));
    return sum;
}

sz()

unsigned Vector::sz ( )

inline

Definition at line 79 of file Vector.h.

{return d->n;};

transposeAndMultiply()

void Vector::transposeAndMultiply	(	Vector	vR,
		Matrix	M
	)

Definition at line 297 of file Vector.cpp.

{
    if ((int)sz()!=M.nLine())
    {
        printf("error in V^t * M.\n"); getchar(); exit(254);
    }
    int n=sz(), szr=M.nColumn(), i;
    vR.setSize(szr);
    double sum, *dv=(*this), **dm=M, *dd=vR;

    while (szr--)
    {
        sum=0.0;
        i=n;
        while (i--) sum+=dv[i]*dm[i][szr];
        dd[szr]=sum;
    }
}

zero()

void Vector::zero	(	int	_i = 0,
		int	_n = 0
	)

Definition at line 93 of file Vector.cpp.

{
    if (_n==0) _n=d->n-i;
    if (d->p) memset(d->p+i,0,_n*sizeof(double));
}

Member Data Documentation

d

VectorData* Vector::d = nullptr

Definition at line 50 of file Vector.h.

emptyVector

Vector Vector::emptyVector

static

Definition at line 119 of file Vector.h.

---

The documentation for this class was generated from the following files:

• xmipp/external/condor/Vector.h
• xmipp/external/condor/Vector.cpp