doublevector.cpp

/*
 * Copyright (C) 2004 - 2005 by
 *     Hieu Xuan Phan & Minh Le Nguyen {hieuxuan, nguyenml}@jaist.ac.jp
 *     Graduate School of Information Science,
 *     Japan Advanced Institute of Science and Technology (JAIST)
 *
 * doublevector.cpp - this file is part of FlexCRFs.
 *
 * Begin:	Dec. 15, 2004
 * Last change:	Sep. 09, 2005
 *
 * FlexCRFs is a free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published
 * by the Free Software Foundation; either version 2 of the License,
 * or (at your option) any later version.
 *
 * FlexCRFs is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with FlexCRFs; if not, write to the Free Software Foundation,
 * Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
 */

#include "doublevector.h"
#include <cstring>


using namespace std;

// constructor with length
doublevector::doublevector(int len) {
    this->len = len;
    vect = new double[len];
}
    
// constructor with length and content
doublevector::doublevector(int len, double * vect) {
    this->len = len;
    vect = new double[len];
    for (int i = 0; i < len; i++) {
        this->vect[i] = vect[i];
    }	    
}
    
// copy constructor
doublevector::doublevector(doublevector & dv) {
    len = dv.len;
    vect = new double[len];
    for (int i = 0; i < len; i++) {
        vect[i] = dv.vect[i];
    }
}
    
// the size of vector
int doublevector::size() {
    return len;
}
    
// overloading assignment operator
void doublevector::operator=(double val) {
	//memset(vect,0,sizeof(double)*len);
    for (int i = 0; i < len; i++) {
        vect[i] = val;
    }
}
    
// overloading assignment operator
void doublevector::operator=(doublevector & dv) {
    if (len != dv.len) {
        if (vect) {
	    delete vect;
	}
	len = dv.len;
	vect = new double[len];	    
    }

	//
    	memcpy((void*)vect,(void*)dv.vect,sizeof(double)*len);	
    /*for (int i = 0; i < len; i++) {
        vect[i] = dv.vect[i];
    }*/
}    
    
// assign the same value to all vector elements
void doublevector::assign(double val) {
    *this = val;
}
    
// assign values for all elements from another doublevector
void doublevector::assign(doublevector & dv) {
    *this = dv;
}
    
// reference to an element of index "idx"
double & doublevector::operator[](int idx) {
    return vect[idx];
}
    
// sum of all vector elements
double doublevector::sum() {
    double res = 0.0;
    for (int i = 0; i < len; i++) {
        res += vect[i];
    }
    return res;
}

// component multiplication
void doublevector::comp_mult(double val) {
    for (int i = 0; i < len; i++) {
	vect[i] *= val;
    }
}

// component multiplication
void doublevector::comp_mult(doublevector * db) {
	//double sum = 0.0;
    for (int i = 0; i < len; i++) {
	vect[i] *= db->vect[i];
	//sum += vect[i];
    }
	//return sum;
}