space_basis.hpp

 
#ifndef EXATN_NUMERICS_SPACE_BASIS_HPP_
#define EXATN_NUMERICS_SPACE_BASIS_HPP_
 
#include "tensor_basic.hpp"
#include "basis_vector.hpp"
 
#include <vector>
 
namespace exatn{
 
namespace numerics{
 
struct SymmetryRange{
 DimOffset lower;    //lower bound of the subrange
 DimOffset upper;    //upper bound of the subrange
 SymmetryId symm_id; //symmetry id of the subrange
};
 
 
class SpaceBasis{
public:
 
 SpaceBasis(DimExtent space_dim);
 SpaceBasis(DimExtent space_dim,
            const std::vector<SymmetryRange> & symmetry_subranges);
 
 SpaceBasis(const SpaceBasis & space_basis) = default;
 SpaceBasis & operator=(const SpaceBasis & space_basis) = default;
 SpaceBasis(SpaceBasis && space_basis) noexcept = default;
 SpaceBasis & operator=(SpaceBasis && space_basis) noexcept = default;
 virtual ~SpaceBasis() = default;
 
 void printIt() const;
 
 DimExtent getDimension() const;
 
 const std::vector<SymmetryRange> & getSymmetrySubranges() const;
 
 void registerSymmetrySubrange(const SymmetryRange subrange);
 
private:
 
 DimExtent basis_dim_;                        //basis dimension (number of basis vectors)
 std::vector<SymmetryRange> symmetry_ranges_; //symmetry subranges (may be none)
};
 
} //namespace numerics
 
} //namespace exatn
 
#endif //EXATN_NUMERICS_SPACES_BASIS_HPP_