wannier_loader.hpp

#pragma once
 
#include "./tb_hk.hpp"
 
#include <nda/nda.hpp>
 
#include <array>
#include <complex>
#include <string>
#include <tuple>
#include <utility>
#include <vector>
 
namespace triqs::experimental::lattice {
 
  using dcomplex = std::complex<double>;


  using r_vector_list_t = std::vector<std::array<long, 3>>;

  using hopping_list_t = std::vector<nda::array<dcomplex, 2>>;

  using position_op_list_t = std::vector<nda::array<dcomplex, 3>>;

  using w90_tb_data_t = std::tuple<r_vector_list_t, hopping_list_t, position_op_list_t, nda::matrix<double>>;

  using w90_hr_data_t = std::tuple<r_vector_list_t, hopping_list_t>;

  w90_tb_data_t read_wannier90_tb_data(std::string const &w90_path_and_seedname);

  inline tb_hk make_tb_hk_from_w90_tb_file(std::string const &w90_path_and_seedname) {
    auto [Rs, hoppings, _unused1, _unused2] = read_wannier90_tb_data(w90_path_and_seedname);
    return {std::move(Rs), std::move(hoppings)};
  };

  w90_hr_data_t read_wannier90_hr_data(std::string const &w90_path_and_seedname);

  inline tb_hk make_tb_hk_from_w90_hr_file(std::string const &w90_path_and_seedname) {
    auto [Rs, hoppings] = read_wannier90_hr_data(w90_path_and_seedname);
    return {std::move(Rs), std::move(hoppings)};
  };

 
} // namespace triqs::experimental::lattice