gsl_vector_integer.h

Go to the documentation of this file.

//

//  This vector integer class has been adapted from the gslwrap project.
//  It is a C++ wrapper for the GNU Scientific Library
//  Copyright (C) 2001 Ramin Nakisa

//  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.

#ifndef _GSL_VECTOR_INTEGER_H
#define _GSL_VECTOR_INTEGER_H

#include <iostream>
#include <assert.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_vector_int.h>
#include <gsl/gsl_blas.h>
#include "gsl_vector_double.h"

//

class vector_int
{
private:

//

gsl_vector_int *gsldata;

//
//  

void alloc(size_t n);

//

void calloc(size_t n);

public:

//
//

    vector_int();

//
//
    
    vector_int(const size_t& n,bool clear=true);

//
//

    ~vector_int();

//    
    vector_int( const vector_int &other );

// 
//

        void copy(const vector_int& other);

// 
//

        void resize(size_t n);

// 

        bool is_set() const;

// 

    size_t size() const;

//

    gsl_vector_int       *gslobj()       {assert(gsldata);return gsldata;}

//

        const gsl_vector_int *gslobj() const {assert(gsldata);return gsldata;}

//
//

    int get(size_t i) const;

//
//

    void  set(size_t i,int x);

//
//

        int       &operator[](size_t i);

//
//
        
        const int &operator[](size_t i) const;

//
//

        int &operator()(size_t i);

//
//

        const int &operator()(size_t i) const;

// 
//
        
        int fprintf (FILE * stream, const char * format) const;

//
//

    int fscanf (FILE * stream);

//
//

        void set_all(int x);

//
//
        
        void set_zero();

//
//

        int set_basis (size_t i);
    
//

    vector_int& operator=(const vector_int& other){copy(other);return (*this);}

//
//

    int swap_elements (size_t i, size_t j);

//

    int reverse ();

//
//

        int operator+=(const vector_int &other) {return gsl_vector_int_add (gsldata, other.gsldata);}

//
//

        int operator-=(const vector_int &other) {return gsl_vector_int_sub (gsldata, other.gsldata);}

//

        int operator*=(const vector_int &other) {return gsl_vector_int_mul (gsldata, other.gsldata);}

//

        int operator/=(const vector_int &other) {return gsl_vector_int_div (gsldata, other.gsldata);}

//
//

        int operator*=(int x) {return gsl_vector_int_scale (gsldata, x);}

//
//

        int operator+=(int x) {return gsl_vector_int_add_constant (gsldata,x);}

//
//

        int operator/=(int x) {return gsl_vector_int_scale (gsldata, 1/x);}

//
//      

        bool operator==(const vector_int& other) const;

//
//      
        
        bool operator!=(const vector_int& other) const { return (!((*this)==other));}

//

        int sum() const;

//

    double max() const;

//

    double min() const;

//

        size_t max_index();
        
//

        bool isnull();

//
//

        void load( const char *filename );

//
//
                
        void save( const char *filename ) const;

//
//

        vector_int glex() const; 

//
//

        vector_int advance_cartesian_product(size_t num_points) const; 
        
        };

#endif // GSL_VECTOR_INTEGER_H