libalgebra-demo/libalgebra/lie.h

/* *************************************************************

Copyright 2010 Terry Lyons, Stephen Buckley, Djalil Chafai, 
Greg Gyurkó and Arend Janssen. 

Distributed under the terms of the GNU General Public License, 
Version 3. (See accompanying file License.txt)

************************************************************* */




//  lie.h


// Include once wrapper
#ifndef DJC_COROPA_LIBALGEBRA_LIEH_SEEN
#define DJC_COROPA_LIBALGEBRA_LIEH_SEEN


template<typename SCA, typename RAT, DEG n_letters, DEG max_degree>
class lie : public algebra<lie_basis<SCA, RAT, n_letters, max_degree> >
{
public:
        typedef lie_basis<SCA, RAT, n_letters, max_degree> BASIS;
        typedef typename BASIS::KEY KEY;
        typedef sparse_vector<BASIS> VECT;
        typedef algebra<BASIS> ALG;
        typedef typename ALG::iterator iterator;
        typedef typename ALG::const_iterator const_iterator;
public:
        lie(void) {}
        lie(const lie& l)
                : ALG(l) {}
        lie(const ALG& a)
                : ALG(a) {}
        lie(const VECT& v)
                : ALG(v) {}
        explicit lie(const KEY& k)
                : ALG(k) {}
        explicit lie(LET letter, const SCA& s)
                : ALG(VECT::basis.keyofletter(letter), s) {}
public:
        inline friend lie replace(const lie& src, const std::vector<LET>& s, const std::vector<lie*>& v)        
        {
                lie result;
                std::map<KEY, lie> table;
                const_iterator i;
                for (i = src.begin(); i != src.end(); ++i)
                        result.add_scal_prod(VECT::basis.replace(i->first, s, v, table), i->second);
                return result;
        }
};

// Include once wrapper
#endif // DJC_COROPA_LIBALGEBRA_LIEH_SEEN

//EOF.