results.hpp

// Copyright (c) 2022--present, The Simons Foundation
// Copyright (c) 2022--present, Max Planck Institute for Polymer Research, Mainz, Germany
// This file is part of TRIQS/ctseg 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 <optional>
#include <triqs/stat/histograms.hpp>
#include <triqs/gfs.hpp>
 
using namespace triqs::gfs;
 
namespace triqs_ctseg {

  struct results_t {

    block_gf<imtime> G_tau;

    std::optional<block_gf<imtime>> F_tau;

    std::optional<block2_gf<imtime>> nn_tau;

    std::optional<block2_gf<dlr_imfreq>> nn_nu_dlr;

    std::optional<gf<imtime>> Sperp_tau;

    std::optional<std::map<std::pair<std::string, std::string>, nda::matrix<double>>> nn_static;

    std::optional<std::map<std::string, nda::array<double, 1>>> densities;

    std::optional<std::vector<double>> pert_order_Delta;

    std::optional<double> average_order_Delta;

    std::optional<std::vector<double>> pert_order_Jperp;

    std::optional<double> average_order_Jperp;

    std::optional<nda::vector<double>> state_hist;

    std::optional<block2_gf<prod<imfreq, imfreq>, tensor_valued<4>>> g2w;

    std::optional<block2_gf<prod<imfreq, imfreq, imfreq>, tensor_valued<4>>> g3w;

    double average_sign;
  };

  void h5_write(h5::group h5group, std::string subgroup_name, results_t const &c);

  void h5_read(h5::group h5group, std::string subgroup_name, results_t &c);
 
} // namespace triqs_ctseg