Col_amp.h

// -*- C++ -*-
 
/*
 * Col_amp.h
 *      Contains the declarations of the class Col_amp, related types and operators
 *  Created on: Jul 7, 2010
 *  Author: Malin Sjodahl
 */
 
#ifndef COLORFULL_Col_amp_h
#define COLORFULL_Col_amp_h
 
#include "Col_str.h"
 
namespace ColorFull {
 
typedef std::vector<Col_str> col_amp;
 
 
class Col_amp {
 
public:
 
        Col_amp() {Scalar = Scalar * 0;}
 
        Col_amp( const std::string str ){Col_amp_of_str( str );}
 
        Col_amp( Col_str Cs ) {
                Scalar = Scalar * 0;
                ca.push_back(Cs);
        }
 
        col_amp ca;
 
        Polynomial Scalar;
 
        void read_in_Col_amp( std::string filename );
 
        void write_out_Col_amp( std::string filename ) const;
 
        const Col_str & at( int i ) const{ return ca.at(i); }
 
        Col_str & at( int i ) {return ca.at(i);}
 
        uint size() const {return ca.size();}
 
        bool empty() const { return ca.empty(); }
 
        void clear() { ca.clear(); }
 
        void erase( int i );
 
        void append( Col_str Cs ) {ca.push_back( Cs );}
 
        void append( col_amp ca_in );
 
 
        // Functions for probing the Col_amp
 
        bool gluons_only() const;
 
        int n_gluon() const;
 
        int n_quark() const;
 
        int n_quark_check() const;
 
        int n_gluon_check() const;
 
        int longest_quark_line() const;
 
        // Functions for manipulating the Col_amp
 
        void remove_1_rings();
 
        void remove_0_rings();
 
        void remove_empty_Col_strs();
 
        void collect_col_strs();
 
        void normal_order_col_strs();
 
        void normal_order();
 
        void simplify();
 
        void conjugate();
 
        void contract_next_neighboring_gluons( );
 
        void contract_2_rings( );
 
        void contract_Quark_line_gluons( );
 
        void contract_a_gluon( );
 
        void contract_all_gluons( );
 
        void contract_quarks( const Col_amp & Ca1, const Col_amp & Ca2 );
 
private:
 
 
        void contract_a_gluon( Col_str & Cs );
 
        void contract_Quark_line_gluons( Quark_line & Ql );
 
        void contract_Quark_line_gluons( Col_str & Cs );
 
        void contract_all_gluons( Col_str & Cs );
 
        void Col_amp_of_str( const std::string str );
 
};
 
// Define operators involving Col_amp
 
std::ostream& operator<<( std::ostream& out, const col_amp & ca );
 
std::ostream& operator<<( std::ostream& out, const Col_amp & Ca );
 
bool operator==( const Col_amp & Ca1, const Col_amp & Ca2 );
 
bool operator!=( const Col_amp & Ca1, const Col_amp & Ca2 );
 
Col_amp operator+( const Col_amp & Ca, const Col_str & Cs );
 
Col_amp operator+( const Col_str & Cs, const Col_amp & Ca );
 
Col_amp operator+( const Col_amp & Ca1, const Col_amp & Ca2 );
 
Col_amp operator+=( Col_amp & Ca, const Col_str & Cs );
 
Col_amp operator+=( Col_amp & Ca1, const Col_amp & Ca2 );
 
Col_amp operator-( const Col_amp & Ca1, const Col_amp & Ca2 );
 
Col_amp operator-( const Col_amp & Ca, const Col_str & Cs );
 
Col_amp operator-( const Col_amp & Ca, const Col_str & Cs );
 
Col_amp operator*( const Col_amp & Ca, const int i );
 
Col_amp operator*( const int i, const Col_amp & Ca );
 
Col_amp operator*( const Col_amp & Ca, const cnum c );
 
Col_amp operator*( const cnum c, const Col_amp & Ca );
 
Col_amp operator*( const Col_amp & Ca, const double d );
 
Col_amp operator*( const double d, const Col_amp & Ca );
 
Col_amp operator*( const Col_amp & Ca, const Monomial & Mon );
 
Col_amp operator*( const Monomial & Mon, const Col_amp & Ca );
 
Col_amp operator*( const Col_amp & Ca, const Polynomial & Poly );
 
Col_amp operator*( const Polynomial & Poly, const  Col_amp & Ca );
 
Col_amp operator*( const Col_amp & Ca, const Quark_line & Ql );
 
Col_amp operator*( const Quark_line & Ql, const Col_amp & Ca );
 
Col_amp operator*( const Col_amp & Ca, const Col_str & Cs );
 
Col_amp operator*( const Col_str & Cs, const Col_amp & Ca );
 
Col_amp operator*( const Col_amp & Ca1, const Col_amp & Ca2 );
 
Col_amp operator*=( Col_amp & Ca1, const Col_amp & Ca2 );
 
}// end namespace ColorFull
 
#endif /* COLORFULL_Col_amp_h */