chiAB_tau.hpp

// Copyright (c) 2018--present, The Simons Foundation
// This file is part of TRIQS/ctint and is licensed under the terms of GPLv3 or later.
// SPDX-License-Identifier: GPL-3.0-or-later
// See LICENSE in the root of this distribution for details.
 
#pragma once
#include "../qmc_config.hpp"
#include "../nfft_buf.hpp"
#include "../container_set.hpp"
 
namespace triqs_ctint::measures {

  struct chiAB_tau {
 
    chiAB_tau(params_t const &params_, qmc_config_t &qmc_config_, container_set *results);

    void accumulate(mc_weight_t sign);

    void collect_results(mpi::communicator const &comm);
 
    private:
    // Capture the parameters FIXME We cannot choose const, as we call try_insert
    params_t const &params;
 
    // The Monte-Carlo configuration
    qmc_config_t &qmc_config;
 
    // Container for the accumulation
    gf_view<imtime> chiAB_tau_;
 
    // The bosonic operator vectors
    using op_term_t = std::tuple<dcomplex, std::pair<int, int>, std::pair<int, int>>;
    std::vector<std::vector<op_term_t>> A_vec;
    std::vector<std::vector<op_term_t>> B_vec;
 
    // The average sign
    mc_weight_t Z = 0.0;
 
    // The tau-mesh
    mesh::imtime tau_mesh;
  };
 
} // namespace triqs_ctint::measures