FuncMatrix.hpp

/*
 * FuncMatrix.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_MATRIXFUNC_HPP__
#define __PDS_MATRIXFUNC_HPP__
 
 
#include <Pds/Matrix>
 
 
namespace Pds{
 
    
    Pds::Matrix Peaks(unsigned int N,double L=3);
 
 
    Pds::Matrix Mountain(unsigned int N,double L=2);
 
 
    bool Meshgrid(const Pds::Vector &SpaceX,const Pds::Vector &SpaceY,Pds::Matrix &X,Pds::Matrix &Y);
    
    bool Meshgrid(const Pds::Vector &SpaceR,Pds::Matrix &X,Pds::Matrix &Y);
    
    Pds::Matrix Operate(double (*func)(double x,double y),const Pds::Matrix &X,const Pds::Matrix &Y);
   
    Pds::Matrix Operate(double (*func)(double x,double y,double z),
                        const Pds::Matrix &X,
                        const Pds::Matrix &Y,
                        const Pds::Matrix &Z);
 
    Pds::Matrix MatrixId(unsigned int N);
   
    Pds::Matrix MatrixId(const Pds::Size &S);
   
    Pds::Matrix MatrixId(unsigned int Nlin,unsigned int Ncol);
    
    Pds::Matrix RandN(unsigned int N);
   
    Pds::Matrix RandN(const Pds::Size &S);
   
    Pds::Matrix RandN(unsigned int Nlin,unsigned int Ncol);
    
    Pds::Matrix RandU(unsigned int N);
    
    Pds::Matrix RandU(const Pds::Size &S);
    
    Pds::Matrix RandU(unsigned int Nlin,unsigned int Ncol);
    
    
    Pds::Matrix RandC(unsigned int N,double p=0.5);
    
    Pds::Matrix RandC(const Pds::Size &S,double p=0.5);
    
    Pds::Matrix RandC(unsigned int Nlin,unsigned int Ncol,double p);
    
    Pds::Matrix Diag(Matrix A);
 
    Pds::Matrix Eye(unsigned int N);
    
    Pds::Matrix Ones(const Pds::Size &S);
    
    Pds::Matrix Ones(unsigned int Nlin, unsigned int Ncol);
    
    Pds::Matrix Ones(unsigned int N);
    
    
    Pds::Matrix Zeros(const Pds::Size &S);
    
    Pds::Matrix Zeros(unsigned int Nlin, unsigned int Ncol);
    
    Pds::Matrix Zeros(unsigned int N);
 
 
    Pds::Matrix CountingTable(unsigned int N,unsigned int M);
    
    Pds::Matrix CountingTable(unsigned int N, int min, int max, unsigned int from_id);
 
    Pds::Matrix MultisetIndexSum(unsigned int N,unsigned int M);
    
    Pds::Matrix SetIndexSum(unsigned int N,unsigned int M);
 
 
    std::string MultisetIndexSumToString(const Pds::Matrix &ID,std::string Separator=", ",unsigned int n=1);
 
    Pds::Matrix Monomials(const Pds::Matrix &X,unsigned int M,Pds::Matrix &ID);
 
    Pds::Matrix Monomials(const Pds::Matrix &X,unsigned int M);
 
    std::vector<unsigned int>  Find(const Pds::Matrix &A);
 
 
    Pds::Matrix  MergeVer(const std::initializer_list<Pds::Matrix> list);
 
    Pds::Matrix  MergeVer(const std::list<Pds::Matrix> &list);
    
 
    Pds::Matrix  MergeVer(const std::vector<Pds::Matrix> &vec);
 
    Pds::Matrix  MergeVer(const std::vector<Pds::Vector> &vec);
 
    Pds::Matrix  MergeHor(const std::initializer_list<Pds::Matrix> list);
 
    Pds::Matrix  MergeHor(const std::list<Pds::Matrix> &list);
 
    Pds::Matrix  MergeHor(const std::vector<Pds::Matrix> &vec);
 
    Pds::Matrix  MergeHor(const std::vector<Pds::Vector> &vec);
    
    Pds::Matrix RegressorMatrix(const std::initializer_list<Pds::Matrix> list);
 
 
    Pds::Matrix RegressorMatrix(const Pds::Matrix &B);
 
    Pds::Matrix  RegressorMatrix(const Pds::Matrix &B,unsigned int M);
 
    Pds::Matrix  RegressorMatrix(const Pds::Matrix &B,unsigned int M,Pds::Matrix &ID);
    
}
#endif