Basis.hpp

Go to the documentation of this file.

#ifndef BASIS_H_
#define BASIS_H_
 
#include "yas/serialize.hpp"
#include "yas/std_types.hpp"
 
namespace Xped {
 
// forward declaration
template <typename Symmetry, std::size_t depth, typename AllocationPolicy>
class Qbasis;
 
class Basis
{
    template <typename Symmetry_, std::size_t depth, typename AllocationPolicy>
    friend class Qbasis;
 
public:
    Basis(){};
 
    Basis(std::size_t dim_in)
        : dim_(dim_in){};
 
    inline const std::size_t dim() const { return dim_; }
 
    Basis add(const Basis& other) const;
 
    Basis combine(const Basis& other) const;
 
    auto print() const;
 
    auto printHistory() const;
 
    inline bool operator==(const Basis& other) const { return (this->dim() == other.dim()); }
 
    friend std::ostream& operator<<(std::ostream& os, const Basis& basis);
 
    template <typename Ar>
    void serialize(Ar& ar)
    {
        ar& YAS_OBJECT_NVP("Basis", ("dim", dim_), ("history", history));
    }
 
private:
    std::size_t dim_;
 
    struct fuseData
    {
        std::size_t dim1 = 0, dim2 = 0;
        std::array<std::size_t, 2> source(std::size_t combined_num) const
        {
            std::size_t tmp = combined_num / dim1;
            return {{combined_num - dim1 * tmp, tmp}};
        };
        template <typename Ar>
        void serialize(Ar& ar)
        {
            ar& YAS_OBJECT_NVP("fuseData", ("dim1", dim1), ("dim2", dim2));
        }
    };
 
    Basis::fuseData history;
};
 
} // namespace Xped
 
#ifndef XPED_COMPILED_LIB
#    include "Core/Basis.cpp"
#endif
 
#endif