---

vector_double.cc

Go to the documentation of this file.

//  This matrix class is a C++ wrapper for the GNU Scientific Library

//  This program is 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.

//  This program 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 this program; if not, write to the Free Software
//  Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.

#include<gslwrap/vector_double.h>
#include<gslwrap/vector_double.h>

#define type_is
#ifdef type_is
#define type_is_double
#endif

namespace gsl
{

//vector::create_vector_view( const gsl_vector_view &other )
vector_view 
vector::create_vector_view( const gsl_vector_view &other )
{
        vector view;
        view.gsldata = (gsl_vector*)malloc(sizeof(gsl_vector));
        *(view.gslobj()) = other.vector;
        view.gslobj()->owner = 0;
        return view;
}

void
vector::resize(size_t n)
{
        if (gsldata)
        {
                if (n==size())
                        return;
                if (!gsldata->owner)
                {
                        cout << "vector::resize ERROR can't resize a vector view" << endl;
                        exit(-1);
//                      GSL_ERROR("You can't resize a vector view", GSL_EINVAL);
                }
                free();
        }
        alloc(n);
}

void 
vector::copy(const vector& other)
{
                resize(other.size());
                gsl_vector_memcpy (gsldata,other.gsldata);
}

bool 
vector::operator==(const vector& other) const 
{
        if (size() != other.size())
                return false;
        for (int i=0;i<size(); i++)
        {
                if (this->operator[](i) != other[i])
                        return false;
        }
        return true;
}

vector_view 
vector::subvector (size_t offset, size_t n)
{
        gsl_vector_view view = gsl_vector_subvector (gsldata, offset, n);
        return vector_view::create_vector_view(view);
}

//vector_const_view 
const vector_view 
vector::subvector (size_t offset, size_t n) const 
{
        gsl_vector_view view = gsl_vector_subvector (gsldata, offset, n);
        return vector_view::create_vector_view(view);
}

// returns sum of all the elements.
double vector::sum() const 
{
        int i;
        double sum = 0;

        for ( i = 0; i < size(); i++ ) 
        {
                sum += gsl_vector_get(gsldata, i);
        }

        return( sum );
}

double 
vector::norm2() const
{
        vector t=*this;
        return gsl_blas_dnrm2(t.gslobj());
}

ostream& 
operator<< ( ostream& os, const vector & vect )
{
        os.setf( ios::fixed);
        for (int i=0;i<vect.size();i++)
        {
                os << vect[i] << endl;
        }
        return os;
}

//**************************************************************************
// Implementation of the vector_view class :
//**************************************************************************


void 
vector_view::init(const vector& other)
{
        free();
        gsldata = (gsl_vector*)malloc(sizeof(gsl_vector));
        *gsldata = *(other.gslobj());
        gsldata->owner = 0;
}

}

---