MathMatrix.hpp

/*
 * MathMatrix.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_MATRIXMATH_HPP__
#define __PDS_MATRIXMATH_HPP__
 
 
#include <Pds/Matrix>
 
namespace Pds{
    
 
    Pds::Matrix Sin(const Pds::Matrix &A);
    
    Pds::Matrix Cos(const Pds::Matrix &A);
    
    Pds::Matrix Tan(const Pds::Matrix &A);
    Pds::Matrix Exp(const Pds::Matrix &A);
    
    Pds::Matrix TwoExp(const Pds::Matrix &A);
    
    Pds::Matrix Log(const Pds::Matrix &A);
 
    Pds::Matrix Ln(const Pds::Matrix &A);
    
    Pds::Matrix Log2(const Pds::Matrix &A);
 
    Pds::Matrix Log10(const Pds::Matrix &A);
 
    Pds::Matrix Sqrt(const Pds::Matrix &A);
 
    Pds::Matrix Pow(const Pds::Matrix &A,double var);
 
    Pds::Matrix Pow(const Pds::Matrix &A,const Pds::Matrix &B);
 
    Pds::Matrix Atan2(const Pds::Matrix &A,const Pds::Matrix &B);
    double Var(const Pds::Matrix &A, double *mean);
 
    double Var(const Pds::Matrix &A);
 
    double Mean(const Pds::Matrix &A);
 
    double Sum(const Pds::Matrix &A);
    
 
    double Accuracy(const Pds::Matrix &A,const Pds::Matrix &B,double Umbral=0.5);
 
 
    
    double MeanAbsoluteError(const Pds::Matrix &A,const Pds::Matrix &B);
    double Det(const Pds::Matrix &A);
 
    double Rms(const Pds::Matrix &A);
 
    double MeanAbsolute(const Pds::Matrix &A);
 
    double MeanSquare(const Pds::Matrix &A);
 
 
    double SumSquare(const Pds::Matrix &A);
 
    double Norm(const Pds::Matrix &A);
    
    double SQNL(double x);
 
    Pds::Matrix SQNL(const Pds::Matrix &A);
    
 
    double Gaussian(double x);
 
    Pds::Matrix Gaussian(const Pds::Matrix &A);
 
    Pds::Matrix Tanh(const Pds::Matrix &A);
 
    double Tanh(double x);
    
    double DTanh(double x);
 
    double Sigmoid(double x);
    
    Pds::Matrix Sigmoid(const Pds::Matrix &A);
 
    double DSigmoid(double x);
    
    unsigned int Factorial(unsigned int n);
 
 
    unsigned int NchooseK(unsigned int n,unsigned int k);
 
    unsigned int NmultichooseK(unsigned int n,unsigned int k);
 
    double UnitStep(double x);
 
    Pds::Matrix UnitStep(const Pds::Matrix &A);
 
    double UnitRamp(double x);
 
    Pds::Matrix UnitRamp(const Pds::Matrix &A);
 
    
    double Sign(double x);
 
 
 
    Pds::Matrix Sign(const Pds::Matrix &A);
    
    Pds::Matrix Round(const Pds::Matrix &A);
 
 
    double InformationGain(const Pds::Matrix &B, const std::vector<Pds::Matrix> &A,double Umbral=0.5);
    
    double Hb(double x);
 
    Pds::Matrix Hb(const Pds::Matrix &A);
 
    double HbInv(double h);
 
    Pds::Matrix HbInv(const Pds::Matrix &A);
 
    Pds::Matrix Abs(const Pds::Matrix &A);
    
    double Logit(double x);
 
    Pds::Matrix Logit(const Pds::Matrix &A);
 
    double Sinc(double x);
 
    Pds::Matrix Sinc(const Pds::Matrix &A);
    double SimpsonIntegration(double (*f)(double), double a,double b,unsigned int n);
    
    double SimpsonIntegration(double (*f)(double,double),double r, double a,double b,unsigned int n);
    
    double ImproperIntegration(double (*f)(double), double a,unsigned int n);
 
 
    double ImproperIntegration(double (*f)(double,double),double r, double a,unsigned int n);
} //namespace Pds
#endif