double_counting.hpp

Go to the documentation of this file.

// Copyright (c) 2025--present, The Simons Foundation
// This file is part of TRIQS/modest 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 <nda/nda.hpp>
#include "utils/defs.hpp"
#include "triqs/gfs.hpp"
 
using namespace triqs::gfs;
 
namespace triqs::modest {
 
  std::pair<double, double> dc_formulas(std::string const method, double const N_tot, double const N_sigma, long const n_orb, double const U,
                                        double const J);
  std::pair<nda::array<nda::matrix<double>, 2>, nda::matrix<double>> double_counting(nda::array<nda::matrix<dcomplex>, 2> const &density_matrix,
                                                                                     double U_int, double J_hund, std::string const method);
 
  class dc_solver {
 
    private:
    long n_sigma;
    std::string method;
    double U_int;
    double J_hund;
 
    nda::array<nda::matrix<double>, 2> get_density_matrix_from_gf(block_gf<imfreq, matrix_valued> const &gimp);
 
    public:
    dc_solver(long n_sigma, std::string method, double U_int, double J_hund);
 
    std::vector<nda::matrix<dcomplex>> dc_self_energy(block_gf<imfreq, matrix_valued> const &gimp);
 
    nda::matrix<double> dc_energy(block_gf<imfreq, matrix_valued> const &gimp);
  };
 
} // namespace triqs::modest