decompose_product.hh

Go to the documentation of this file.

 
#pragma once
 
#include "algorithms/young_project.hh"
#include "properties/TableauBase.hh"
#include "properties/Indices.hh"
 
namespace cadabra {
 
 
    class decompose_product : public Algorithm {
        public:
            decompose_product(const Kernel&, Ex& tr);
 
            virtual bool     can_apply(iterator) override;
            virtual result_t apply(iterator&) override;
 
        private:
            typedef young_project::name_tab_t  sibtab_t;
            typedef yngtab::tableaux<sibtab_t> sibtabs_t;
            typedef young_project::pos_tab_t   numtab_t;
            typedef yngtab::tableaux<numtab_t> numtabs_t;
 
            const Indices *indices_equivalent(iterator it) const;
 
            void fill_asym_ranges(TableauBase::tab_t& tab, int offset, combin::range_vector_t&);
            void project_onto_initial_symmetries(Ex& rep, iterator rr, young_project& yp,
                                                 const TableauBase *tt, iterator ff,
                                                 int offset, const TableauBase::tab_t& thetab,
                                                 bool remove_traces);
 
            iterator               f1, f2;
            const TableauBase     *t1, *t2;
            TableauBase::tab_t     t1tab, t2tab;
            const Indices         *ind1, *ind2;
            unsigned int dim;
            yngtab::filled_tableau<iterator> nt1, nt2;
 
            combin::range_vector_t asym_ranges;
        };
 
    }