AnyAlgebra2.h

Go to the documentation of this file.

#ifndef ANYALGEBRA2_H_
#define ANYALGEBRA2_H_
 
namespace muq {
  namespace Modeling {
 
      class AnyMat{
 
          AnyMat(std::shared_ptr<boost::any> objIn);
          
          unsigned Size(const int dim=-1);
          double Norm();
 
          AnyMat operator*(AnyMat const& otherMat) const;
 
          // Access operators
          double  operator(int i, int j=-1) const;
          double& operator(int i, int j=-1);
 
          AnyMat operator+(AnyMat const& otherMat) const;
          AnyMat operator-(AnyMat const& otherMat) const;
 
          AnyMat& operator+=(AnyMat const& otherMat);
          AnyMat& operator-=(AnyMat const& otherMat);
 
          AnyMat operator=(AnyMat const& otherMat);
 
          AnyMat Solve(AnyMat const& rhs);
 
          AnyMat Cholesky();
          
          AnyMat Zero(int rows, int cols=-1) const;
          AnyMat Ones(int rows, int cols=-1) const;
 
          bool IsZero() const;
 
          AnyMat Identity(int rows, int cols) const;
 
          double LogDeterminant() const;
 
          AnyMat Concatenate(AnyMat const& otherMat, int axis=-1) const;
          
      private:
          std::shared_ptr<boost::any> obj;
      };
 
      
    class AnyAlgebra2 {
    public:
 
      static unsigned int Size(boost::any const& obj, int const dim=-1);
      
 
 
      double Norm(boost::any const& obj) const;
 
 
      double InnerProduct(boost::any const& vec1, boost::any const& vec2) const;
 
 
      boost::any OuterProduct(boost::any const& vec1, boost::any const& vec2) const;
 
 
      boost::any AccessElement(boost::any const& obj, unsigned int const i = 0, unsigned int const j = 0) const;
 
 
      boost::any Zero(std::type_info const& type, unsigned int const rows = 0, unsigned int const cols = 0) const;
 
 
      bool IsZero(boost::any const& obj) const;
 
 
      boost::any Identity(std::type_info const& type, unsigned int const rows=0, unsigned int const cols=0) const;
 
 
      boost::any Identity(std::type_index const& type, unsigned int const rows=0, unsigned int const cols=0) const;
 
 
      boost::any Add(boost::any const& in0, boost::any const& in1) const;
 
 
      boost::any Subtract(boost::any const& in0, boost::any const& in1) const;
 
 
      boost::any Multiply(boost::any const& in0, boost::any const& in1) const;
 
 
      boost::any ApplyInverse(boost::any const& A, boost::any const& x) const;
 
 
      boost::any Apply(boost::any const& A, boost::any const& x) const;
 
 
      boost::any Inverse(boost::any const& obj) const;
 
 
      boost::any SquareRoot(boost::any const& obj) const;
 
 
      double LogDeterminate(boost::any const& obj) const;
 
 
      boost::any Concatenate(boost::any const& vec1, boost::any const& vec2) const;
 
    private:
 
 
      virtual unsigned int SizeImpl(boost::any const& obj) const;
      
 
      virtual double NormImpl(boost::any const& obj) const;
            
 
      virtual double InnerProductImpl(boost::any const& vec1, boost::any const& vec2) const;
 
 
      virtual boost::any OuterProductImpl(boost::any const& vec1, boost::any const& vec2) const;
 
 
      virtual boost::any AccessElementImpl(boost::any const& vec, unsigned int const i) const;
 
 
      virtual boost::any IdentityImpl(std::type_index const& type, unsigned int const rows, unsigned int const cols) const;
      
 
      virtual boost::any AddImpl(boost::any const& in0, boost::any const& in1) const;
      
 
      virtual boost::any SubtractImpl(boost::any const& in0, boost::any const& in1) const;
 
 
      virtual boost::any MultiplyImpl(boost::any const& in0, boost::any const& in1) const;
 
 
      virtual boost::any ApplyImpl(boost::any const& A, boost::any const& x) const;
      
 
      virtual boost::any ApplyInverseImpl(boost::any const& A, boost::any const& x) const;
 
 
      virtual boost::any ZeroImpl(std::type_info const& type, unsigned int const rows, unsigned int const cols) const;
 
 
      virtual bool IsZeroImpl(boost::any const& obj) const;
 
 
      virtual boost::any InverseImpl(boost::any const& obj) const;
 
 
      virtual boost::any SquareRootImpl(boost::any const& obj) const;
 
 
      virtual double LogDeterminateImpl(boost::any const& obj) const;
 
 
      virtual boost::any ConcatenateImpl(boost::any const& vec1, boost::any const& vec2) const;
 
    };
  } // namespace Modeling
} // namespace muq
 
#endif