solver_core.hpp

#pragma once
#include <memory>
#include <triqs/gfs.hpp>
#include <triqs/mesh.hpp>
#include <triqs/utility/macros.hpp>
#include <nda/nda.hpp>
#include <h5/h5.hpp>
#include "Nevanlinna_parameters_t.hpp"
#include "kernel.hpp"
 
namespace triqs_Nevanlinna {

  class solver_core {
 
    public:
    ~solver_core() = default;

    solver_core(Nevanlinna_parameters_t const &p);
 
    // Copy/Move construction
    solver_core(solver_core const &) = delete;
    solver_core(solver_core &&)      = default;

    solver_core &operator=(solver_core const &) = delete;
    solver_core &operator=(solver_core &&)      = default;

    void solve(triqs::gfs::gf_const_view<triqs::mesh::imfreq> g_iw);

    [[nodiscard]] triqs::gfs::gf<triqs::mesh::refreq> evaluate(const triqs::mesh::refreq &grid, double eta);

    [[nodiscard]] triqs::gfs::gf<triqs::mesh::refreq> evaluate(const triqs::mesh::refreq &grid, double eta,
                                                               nda::array_const_view<std::complex<double>, 3> theta);

    [[nodiscard]] C2PY_PROPERTY_GET(Pick_eigenvalues) nda::vector<double> get_Pick_eigenvalues() const { return _kernel->get_Pick_eigenvalues(); };

    [[nodiscard]] C2PY_PROPERTY_GET(size) size_t size() const { return _kernel->size(); };
 
    private:
    // vector of Nevanlinna factorization kernels for multi-orbital factorization
    std::unique_ptr<kernel> _kernel{};
  };
} // namespace triqs_Nevanlinna