ParityOperator.hpp

Go to the documentation of this file.

#ifndef OBSERVABLE_PARITYOPERATOR_HPP
#define OBSERVABLE_PARITYOPERATOR_HPP

#include "Observable.hpp"

namespace solid
{

template <template <typename> class T1, typename T2>
class ParityOperator : public Observable<T1, T2>, public IOperatorParameters<T1, T2>
{
    using Observable<T1, T2>::_operator;
    using Observable<T1, T2>::parameters;

public:
    ParityOperator<T1, T2>(int _L)
    {
        L = _L;
        SelectTerms();
        FillParameters();
    }

private:
    int L;

    void SelectTerms()
    {
        _operator.termsEnabled.terms = {TermsEnum::FermionOneBodyInteraction};
    }

    void FillParameters() override
    {
        arma::SpMat<T2> mu = arma::eye<arma::SpMat<T2>>(L, L);
        parameters["mu"] = mu;
    }

public:
    void PreprocessingVirtual(T1<T2> &matrixElements) override { Preprocessing(matrixElements); }
    static void Preprocessing(T1<T2> &matrixElements)
    {
        std::cout << "parity pre-process..." << std::endl;
        matrixElements.diag() = arma::Col<T2>(matrixElements.diag()).transform([](T2 val) { return std::pow(-1, val); });
    }
};

} // namespace solid

#endif