Commit 4ea33960 authored by gbalduzz's avatar gbalduzz
Browse files

more robust detection of ndd interaction.

parent 75ab4784

src/phys/dca_step/cluster_solver/ctint/structs/ctint_matrix_configuration.cpp

+7 −5
Original line number Diff line number Diff line
@@ -39,10 +39,16 @@ void MatrixConfiguration::addVertex(Vertex& v, unsigned config_id, 
  auto spin = [=](const int nu) { return nu >= n_bands_; };

  auto band = [=](const int nu) -> ushort { return nu - n_bands_ * spin(nu); };



  std::array<unsigned, 2> indices;

  const auto& nu = (*H_int_)[v.interaction_id].nu;

  const auto& r = (*H_int_)[v.interaction_id].r;



  const bool is_ndd = band(nu[0]) != band(nu[1]) || band(nu[2]) != band(nu[3]);



  auto field_type = [&](const Vertex& v, const int leg) -> short {

    const short sign = v.aux_spin ? 1 : -1;

    const InteractionElement& elem = (*H_int_)[v.interaction_id];

    if (elem.partners_id.size())

    if (is_ndd)

      return leg == 1 ? -3 * sign : 3 * sign;  // non density-density.

    else if (elem.w > 0)

      return leg == 1 ? -1 * sign : 1 * sign;  // positive dd interaction.
@@ -50,10 +56,6 @@ void MatrixConfiguration::addVertex(Vertex& v, unsigned config_id, 
      return 2 * sign;  // negative dd interaction.

  };



  std::array<unsigned, 2> indices;

  const auto& nu = (*H_int_)[v.interaction_id].nu;

  const auto& r = (*H_int_)[v.interaction_id].r;



  for (ushort leg = 0; leg < 2; ++leg) {

    assert(spin(nu[0 + 2 * leg]) == spin(nu[1 + 2 * leg]));

    const short s = spin(nu[0 + 2 * leg]);