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(const std::vector<double> &vec);
    
    Vector(const std::vector<int> &vec);
    
    Vector(const std::vector<bool> &vec);
    
    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();
 
public:
 
    Pds::Vector Conv(const Pds::Vector &B, Pds::Ra::Padding  Padding=Pds::Ra::Full) const;
    
    
    static Pds::Size Conv1DOutputSize(  unsigned int Anlin,
                                        unsigned int Bnlin, 
                                        unsigned int Stride,
                                        Pds::Ra::Padding  Padding);
    
    Pds::Vector Conv1D( const Pds::Vector &B, 
                        unsigned int Stride=1,
                        Pds::Ra::Padding  Padding=Pds::Ra::Full) const;
 
    Pds::Vector XCorr(const Pds::Vector &B, bool Same=false) const;
    
    Pds::Matrix MulTComp(double b,const Pds::Vector &B) const;
 
public:
 
    static bool IsVectorBlock(const std::vector<Pds::Vector> &Block);
    
    
    static Pds::Vector FlattenVectorBlock(const std::vector<Pds::Vector> &Block);
    
    
    static std::vector<Pds::Vector> RandNVectorBlock(unsigned int N,unsigned int Sz);
    
    
    static std::vector<Pds::Vector> MaxPoolingVectorBlock(  const std::vector<Pds::Vector> &A,
                                                            unsigned int LinPool);
                                                            
    
 
public:
    
    static bool ExportVectorBlockXmlToStream(std::ofstream &myfile,const std::vector<Pds::Vector> &Block);
    
    static std::string ExportVectorBlockXmlToString(const std::vector<Pds::Vector> &Block);
    
    
    static std::vector<Pds::Vector> ImportVectorBlockXmlFromString(const std::string &str);
    
public:
 
    static bool IsTensorVectorBlock(const std::vector<std::vector<Pds::Vector>> &Tensor);
    
    static std::vector<std::vector<Pds::Vector>> RandNTensorVectorBlock(unsigned int L, unsigned int N, unsigned int Sz);
    
    static std::string ExportTensorVectorBlockXmlToString(const std::vector<std::vector<Pds::Vector>> &Block);
    
    static std::vector<std::vector<Pds::Vector>> ImportTensorVectorBlockXmlFromString(const std::string &str);
    
    
    static bool ExportTensorVectorBlockXmlToStream( std::ofstream &myfile,
                                                    const std::vector<std::vector<Pds::Vector>> &Block);
    
public:
 
    Pds::Vector Unit(void) const;
    
    bool Normalize(void);
    
    
    Pds::Vector Reverse(void);
    
    Pds::Vector Cumulative(void);
    
    Pds::Vector CumulativeNorm(void);
    
    
public:
    bool Reshape(unsigned int Nlin,unsigned int Ncol);
 
    bool Reshape(unsigned int Nel);
 
 
    bool FusionVer(std::list<Pds::Matrix> &list);
 
public:
    std::vector<unsigned int> 
    NearestK(   unsigned int K,
                const std::vector<Pds::Vector> &C,
                std::vector<double> &D2) const;
    
    
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);
                               
    void SortWith_merge(unsigned int init, 
                        unsigned int mid, 
                        unsigned int end, 
                        double **tmp1);
                            
    void SortWith_mergesort(   unsigned int init,
                               unsigned int end, 
                               double **tmp1);
 
    void SortWith_mergeid(  unsigned int init, 
                            unsigned int mid, 
                            unsigned int end, 
                            double **tmp1,
                            std::vector<unsigned int> &id,
                            std::vector<unsigned int> &idtmp);
                                    
    void SortWith_mergesortid(  unsigned int init,
                                unsigned int end, 
                                double **tmp1,
                                std::vector<unsigned int> &id,
                                std::vector<unsigned int> &idtmp);
         
public:
    bool SortWith(Pds::Matrix &A);
 
    bool SortMe(void);
    
    bool SortMeReverse(void);
    
 
    std::vector<unsigned int> SortMeAndGetIds( void);
    
 
    bool SortFollowing(const std::vector<unsigned int> &id);
public:
    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;
    }
    
public:
 
    Pds::Json ToJson(void) const;
    
    
    bool FromJson(const Pds::Json &J) ;
    
public:
    
    bool Copy(const Pds::Matrix &B);
 
    bool Vectorize(const Pds::Matrix &B)
    {   
        return this->Copy(B);
    }
    
    
    bool CopyFromRow(unsigned int lin, const Pds::Matrix &B);
 
    void AddRawAssigAt(unsigned int lin,const double &val)
    {
        this->array[lin][0]+=val;
    }
}; // Class Vector
 
} // namespace Pds
 
 
#endif