Commit e4e15c01 authored by gbalduzz's avatar gbalduzz
Browse files

Use correct domain for inversion of points in k-space.

parent 4ffcef97

include/dca/phys/dca_step/cluster_solver/shared_tools/accumulation/tp/tp_accumulator.hpp

+4 −4
Original line number Diff line number Diff line
@@ -312,8 +312,8 @@ std::complex<typename TpAccumulator<Parameters, linalg::CPU>::Real> TpAccumulato 
  auto minus_w2 = [=](const int w) { return 2 * n_pos_frqs_ - 1 - w; };

  auto plus_w1 = [=](const int w) { return w - n_pos_frqs_; };

  auto minus_k = [=](const int k) {

    const static int k0 = RDmn::parameter_type::origin_index();

    return RDmn::parameter_type::subtract(k, k0);

    const static int k0 = KDmn::parameter_type::origin_index();

    return KDmn::parameter_type::subtract(k, k0);



  };
@@ -333,8 +333,8 @@ void TpAccumulator<Parameters, linalg::CPU>::getGMultiband(int s, int k1, int k2 
  auto plus_w1 = [=](const int w) { return w - n_pos_frqs_; };



  auto minus_k = [=](const int k) {

    const static int k0 = RDmn::parameter_type::origin_index();

    return RDmn::parameter_type::subtract(k, k0);

    const static int k0 = KDmn::parameter_type::origin_index();

    return KDmn::parameter_type::subtract(k, k0);



  };