Commit 51682803 authored by gbalduzz's avatar gbalduzz
Browse files

cleaned double insertion proposal

parent 2e6a3463

include/dca/phys/dca_step/cluster_solver/ctint/structs/solver_configuration.hpp

+11 −13
Original line number Diff line number Diff line
@@ -151,16 +151,14 @@ void SolverConfiguration::insertRandom(Rng& rng) { 
  Vertex v(aux_spin, ushort(indices.first), current_tag_++, tau);

  push_back(v);



  if (double_insertion_prob_) {

    // TODO: generalize to multiband n_bands > 2

    if (indices.second != -1 && double_insertion_prob_) {

  if (indices.second != -1) {

    const double tau2 = rng() * max_tau_;

    const bool aux_spin2 = rng() > 0.5;

    Vertex v2(aux_spin2, ushort(indices.second), current_tag_++, tau2);

    push_back(v2);

    last_insertion_size_ = 2;

  }

    else

  else {

    last_insertion_size_ = 1;

  }

  assert(2 * size() == getSector(0).size() + getSector(1).size());