join_gamma.hh

Go to the documentation of this file.

#pragma once
 
#include "Algorithm.hh"
#include "properties/GammaMatrix.hh"
 
namespace cadabra {
 
    class join_gamma : public Algorithm {
        public:
            join_gamma(const Kernel&, Ex&, bool expand, bool use_gendelta);
 
            virtual bool     can_apply(iterator);
            virtual result_t apply(iterator&);
 
            bool expand;
            std::vector<int>     only_expand;
            const GammaMatrix   *gm1, *gm2;
        private:
            void regroup_indices_(sibling_iterator, sibling_iterator, unsigned int,
                                  std::vector<Ex>&, std::vector<Ex>& );
            void append_prod_(const std::vector<Ex>& r1, const std::vector<Ex>& r2,
                              unsigned int num1, unsigned int num2, unsigned int i, multiplier_t mult,
                              Ex& rep, iterator loc);
 
            bool                use_generalised_delta_;
            Ex::iterator   gamma_name_;
            str_node::bracket_t gamma_bracket_;
        };
 
    }