wick.h

Go to the documentation of this file.

// This file is part of MARTY.
// 
// MARTY is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
// 
// MARTY is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
// 
// You should have received a copy of the GNU General Public License
// along with MARTY. If not, see <https://www.gnu.org/licenses/>.

#ifndef WICK_H_INCLUDED
#define WICK_H_INCLUDED

#include "quantumField.h"

namespace mty {

class Wick: public csl::Operator<csl::AbstractFunc> {

    public:

    Wick();

    explicit
    Wick(const csl::Expr& expr);

    ~Wick(){};

    csl::Type getType() const override;

    void print(int mode=0,
               std::ostream& out=std::cout,
               bool lib = false) const override;

    void printCode(int mode = 0, std::ostream &out = std::cout) const override;

    std::string printLaTeX(int mode=0) const override;

    std::optional<csl::Expr> evaluate(
            csl::eval::mode user_mode = csl::eval::base
            ) const override;

    bool operatorAppliesOn(csl::Expr_info expr) const override;

    csl::Expr getOperand() const override;

    void setOperand(const csl::Expr& operand) override;

    bool operator==(csl::Expr_info other) const override;

    csl::unique_Expr copy_unique() const override;

    csl::Expr deepCopy() const override;

    csl::Expr refresh() const override;

    csl::Expr deepRefresh() const override;

    private:

    static int numberOfFields(const csl::Expr& expr);

    bool checkEvenNumberOfFields() const;

    void checkRightForm();

    void expandToRightForm();

    protected:

    mutable bool rightForm = false;
};

csl::Expr wick_s(const csl::Expr& expr);

} // End of namespace mty

#endif /* WICK_H_INCLUDED */