refreq.hpp

Go to the documentation of this file.

// Copyright (c) 2014-2018 Commissariat à l'énergie atomique et aux énergies alternatives (CEA)
// Copyright (c) 2014-2018 Centre national de la recherche scientifique (CNRS)
// Copyright (c) 2018-2023 Simons Foundation
//
// This program 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.
//
// This program 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 may obtain a copy of the License at
//     https://www.gnu.org/licenses/gpl-3.0.txt
//
// Authors: Philipp Dumitrescu, Olivier Parcollet, Nils Wentzell

 
#pragma once
 
#include "./bases/linear.hpp"
#include "./tail_fitter.hpp"
#include "../utility/macros.hpp"
 
#include <fmt/format.h>
#include <h5/h5.hpp>
 
#include <iostream>
#include <string>
#include <utility>
 
namespace triqs::mesh {


  class C2PY_RENAME(MeshReFreq) refreq : public detail::linear<refreq, double>, public tail_fitter_handle {
    public:
    using mesh_point_t = detail::linear<refreq, double>::mesh_point_t;

    C2PY_DEPRECATED_PARAMETER_NAME(n_max : n_w)
    refreq(double w_min = 0.0, double w_max = 0.0, long n_w = 0) : linear(w_min, w_max, n_w) {}

    refreq(std::pair<double, double> window, int n_w) : refreq(window.first, window.second, n_w) {}

    static constexpr bool positive_only() { return false; }

    [[nodiscard]] C2PY_PROPERTY_GET(w_min) double w_min() const { return a_; }

    [[nodiscard]] C2PY_PROPERTY_GET(w_max) double w_max() const { return b_; }

    [[nodiscard]] static std::string hdf5_format() { return "MeshReFreq"; }

    friend void h5_write(h5::group g, std::string const &name, refreq const &m) { m.h5_write_impl(g, name, "MeshReFreq"); }

    friend void h5_read(h5::group g, std::string const &name, refreq &m) { m.h5_read_impl(g, name, "MeshReFreq"); }

    friend std::ostream &operator<<(std::ostream &sout, refreq const &m) {
      return sout << fmt::format("Real frequency mesh with w_min = {}, w_max = {}, N = {}", m.a_, m.b_, m.N_);
    }
  };

  auto evaluate(refreq const &m, auto const &f, double w) { return m.evaluate(f, w); }

 
} // namespace triqs::mesh