spaces.hpp

 
#ifndef EXATN_NUMERICS_SPACES_HPP_
#define EXATN_NUMERICS_SPACES_HPP_
 
#include "tensor_basic.hpp"
#include "space_basis.hpp"
 
#include <utility>
#include <string>
#include <vector>
 
namespace exatn{
 
namespace numerics{
 
class VectorSpace{
public:
 
 VectorSpace(DimExtent space_dim);
 VectorSpace(DimExtent space_dim,
             const std::string & space_name);
 VectorSpace(DimExtent space_dim,
             const std::string & space_name,
             const std::vector<SymmetryRange> & symmetry_subranges);
 
 VectorSpace(const VectorSpace & vector_space) = default;
 VectorSpace & operator=(const VectorSpace & vector_space) = default;
 VectorSpace(VectorSpace && vector_space) noexcept = default;
 VectorSpace & operator=(VectorSpace && vector_space) noexcept = default;
 virtual ~VectorSpace() = default;
 
 void printIt() const;
 
 DimExtent getDimension() const;
 
 const std::string & getName() const;
 
 const std::vector<SymmetryRange> & getSymmetrySubranges() const;
 
 void registerSymmetrySubrange(const SymmetryRange subrange);
 
 SpaceId getRegisteredId() const;
 
 friend class SpaceRegister;
 
private:
 
 void resetRegisteredId(SpaceId id);
 
 SpaceBasis basis_;       //basis defining the vector space
 std::string space_name_; //optional space name
 SpaceId id_;             //registered space id (defaults to SOME_SPACE)
 
};
 
 
class Subspace{
public:
 
 Subspace(const VectorSpace * vector_space,
          DimOffset lower_bound,
          DimOffset upper_bound);
 Subspace(const VectorSpace * vector_space,
          std::pair<DimOffset,DimOffset> bounds);
 Subspace(const VectorSpace * vector_space,
          DimOffset lower_bound,
          DimOffset upper_bound,
          const std::string & subspace_name);
 Subspace(const VectorSpace * vector_space,
          std::pair<DimOffset,DimOffset> bounds,
          const std::string & subspace_name);
 
 Subspace(const Subspace & subspace) = default;
 Subspace & operator=(const Subspace & subspace) = default;
 Subspace(Subspace && subspace) noexcept = default;
 Subspace & operator=(Subspace && subspace) noexcept = default;
 virtual ~Subspace() = default;
 
 void printIt() const;
 
 DimExtent getDimension() const;
 
 DimOffset getLowerBound() const;
 
 DimOffset getUpperBound() const;
 
 std::pair<DimOffset,DimOffset> getBounds() const;
 
 const VectorSpace * getVectorSpace() const;
 
 const std::string & getName() const;
 
 SubspaceId getRegisteredId() const;
 
 friend class SubspaceRegister;
 
private:
 
 void resetRegisteredId(SubspaceId id);
 
 const VectorSpace * vector_space_; //non-owning pointer to the vector space
 DimOffset lower_bound_;            //lower bound defining the subspace of the vector space
 DimOffset upper_bound_;            //upper bound defining the subspace of the vector space
 std::string subspace_name_;        //optional subspace name
 SubspaceId id_;                    //registered subspace id (defaults to UNREG_SUBSPACE)
 
};
 
} //namespace numerics
 
} //namespace exatn
 
#endif //EXATN_NUMERICS_SPACES_HPP_