retime.hpp

Go to the documentation of this file.

// Copyright (c) 2014-2017 Commissariat à l'énergie atomique et aux énergies alternatives (CEA)
// Copyright (c) 2014-2017 Centre national de la recherche scientifique (CNRS)
// Copyright (c) 2019-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 "../utility/macros.hpp"
 
#include <fmt/format.h>
#include <h5/h5.hpp>
 
#include <iostream>
#include <string>
#include <utility>
 
namespace triqs::mesh {


  class C2PY_RENAME(MeshReTime) retime : public detail::linear<retime, double> {
    public:
    using mesh_point_t = detail::linear<retime, double>::mesh_point_t;

    C2PY_DEPRECATED_PARAMETER_NAME(n_max : n_t)
    retime(double t_min = 0.0, double t_max = 0.0, long n_t = 0) : linear(t_min, t_max, n_t) {}

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

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

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

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

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

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

    friend std::ostream &operator<<(std::ostream &sout, retime const &m) {
      return sout << fmt::format("Real time mesh with t_min = {}, t_max = {}, N = {}", m.a_, m.b_, m.N_);
    }
  };

  auto evaluate(retime const &m, auto const &f, double t) { return m.evaluate(f, t); }

 
} // namespace triqs::mesh