Vector.hpp

/*
 * Vector.hpp
 * 
 * Copyright 2018 Fernando Pujaico Rivera <fernando.pujaico.rivera@gmail.com>
 * 
 * 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., 51 Franklin Street, Fifth Floor, Boston,
 * MA 02110-1301, USA.
 * 
 */
 
#ifndef __PDS_VECTOR_HPP__
#define __PDS_VECTOR_HPP__
 
 
#include <Pds/Matrix>
#include <Pds/Size>
 
#include <sstream>      // std::ostringstream
 
namespace Pds{
 
class Vector: public Matrix
{
 
public:
 
    Vector(void);
    
    Vector(unsigned int N);
 
 
    Vector(const Pds::Size &S);
    
    Vector(const Pds::Matrix &B);
 
 
    template <class Datum>
    Vector(const Pds::Array<Datum> &B);
    
    Vector(const char *str);
 
 
    Vector(const std::initializer_list<double> list);
    
    Vector(unsigned int N,double val);
 
    Vector(const Pds::Matrix &B,unsigned int col);
    
    Vector(double (*func)(double),const Pds::Vector &B );
 
    Vector(double (*func)(double),const Pds::Matrix &B );
    
 
    Vector(Pds::Ra::FormatType Type,std::string filepath);
    
    ~Vector();
 
    Pds::Vector Conv(const Pds::Vector &B, bool Same=false) const;
 
    Pds::Vector XCorr(const Pds::Vector &B, bool Same=false) const;
    
    Pds::Matrix MulTComp(double b,const Pds::Vector &B) const;
    
    Pds::Vector Unit(void) const;
    
    
    bool Reshape(unsigned int Nlin,unsigned int Ncol);
 
    bool Reshape(unsigned int Nel);
 
 
    bool FusionVer(std::list<Pds::Matrix> &list);
 
    
private:
    void SortWith_merge(Pds::Matrix &A, 
                        unsigned int init, 
                        unsigned int mid, 
                        unsigned int end, 
                        double **tmp1, 
                        double **tmp2);
                        
    void SortWith_mergesort(Pds::Matrix &A, 
                            unsigned int init,
                            unsigned int end,
                            double **tmp1,
                            double **tmp2);
public:
    bool SortWith(Pds::Matrix &A);
 
    const double &operator[](const unsigned int &id) const
    {
        return this->array[id][0];
    }
 
 
    const double &GetRaw(unsigned int lin) const
    {
        return this->array[lin][0];
    }
 
 
    void SetRaw(unsigned int lin,const double &val)
    {
        this->array[lin][0]=val;
    }
    bool Copy(const Pds::Matrix &B);
 
    bool Vectorize(const Pds::Matrix &B)
    {   
        return this->Copy(B);
    }
    
    
    bool CopyRow(unsigned int lin, const Pds::Matrix &B);
 
    void ElementAddAssigRaw(unsigned int lin,const double &val)
    {
        this->array[lin][0]+=val;
    }
}; // Class Vector
 
} // namespace Pds
 
 
#endif