/data/Documents/eclipse/workspace/ColorFull/Col_amp.h

Go to the documentation of this file.

// -*- 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 */