mc_weight_type Accept() { 
 Config.DT.confirm_removeTwoOperators(); 
 det->complete_operation(); 
#ifdef DEBUG
 std::cout << "CONFIG ACCEPT: " << Config.DT << std::endl;
 for (int a = 0; a<Config.Na; ++a) print_det(Config.dets[a]);
#endif
 Config.update_Sign();
 return Config.ratioNewSign_OldSign();
}   
Пример #2
0
  mc_weight_type accept() { 
    Config.DT.confirm_insert_and_remove_One_Operator();
    det->complete_operation(); 

#ifdef DEBUG
    std::cout << "CONFIG ACCEPT: " << Config.DT << std::endl;
    std::cout << "ROLL MATRIX: " << roll_matrix << std::endl;
    for (int a = 0; a < Config.Na; ++a) print_det(Config.dets[a]);
#endif

    Config.update_Sign();
    return Config.ratioNewSign_OldSign() * det->roll_matrix(roll_matrix);
  }   
  mc_weight_type Accept() { 
   Config.DT.confirm_insertTwoOperators();
   det->complete_operation(); 
   if (Config.RecordStatisticConfigurations) {
    HISTO_Length_Kinks_Accepted << deltaTau;
   }
#ifdef DEBUG
   std::cout << "CONFIG ACCEPT: " << Config.DT << std::endl;
   for (int a = 0; a<Config.Na; ++a) print_det(Config.dets[a]);
#endif
   Config.update_Sign();
   return Config.ratioNewSign_OldSign();
  }