AnyAlgebra.h

Go to the documentation of this file.

#ifndef ANYALGEBRA_H_
#define ANYALGEBRA_H_
 
#include <iostream>
 
#include "MUQ/config.h"
 
#include "MUQ/Modeling/LinearAlgebra/ScalarAlgebra.h"
#include "MUQ/Modeling/LinearAlgebra/EigenVectorAlgebra.h"
#include "MUQ/Modeling/LinearAlgebra/EigenMatrixAlgebra.h"
#if MUQ_HAS_SUNDIALS==1
#include "MUQ/Modeling/LinearAlgebra/SundialsAlgebra.h"
#endif
 
namespace muq {
  namespace Modeling {
    class AnyAlgebra {
    public:
 
      AnyAlgebra();
 
 
      unsigned int Size(boost::any const& obj, int const dim=-1) const;
 
 
      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 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;
 
    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_info 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;
 
    };
  } // namespace Modeling
} // namespace muq
 
#endif