Unverified Commit 23e03873 authored by gbalduzz's avatar gbalduzz Committed by GitHub
Browse files

Merge branch 'master' into optional_cuda_aware_mpi

parents 9cb55a7b 329e2f2e

include/dca/phys/dca_step/symmetrization/symmetrize_single_particle_function.hpp

+1 −1
Original line number Diff line number Diff line
@@ -678,7 +678,7 @@ void symmetrize_single_particle_function::executeCluster( 
        for (int s_ind = 0; s_ind < sym_super_cell_dmn_t::dmn_size(); ++s_ind) {

          int k_new = k_symmetry_matrix(k_ind, b0, s_ind).first;  // FIXME: b0 -> b1



          int b0_new = k_symmetry_matrix(k_new, b0, s_ind).second;

          int b0_new = k_symmetry_matrix(k_ind, b0, s_ind).second;

          int b1_new = k_symmetry_matrix(k_ind, b1, s_ind).second;



          const double sign = Lattice::transformationSignOfK(b0, b1, s_ind);

include/dca/phys/models/analytic_hamiltonians/twoband_Cu.hpp

+4 −7
Original line number Diff line number Diff line
@@ -97,18 +97,15 @@ int TwoBandCu<PointGroup>::transformationSign(int b1, int b2, int s) { 
  if (b1 == b2)

    return 1;



  //  using List = typename PointGroup::PointGroup_list;

  using dca::util::IndexOf;



  constexpr bool symmetrize_off_diagonal = true;



  if (symmetrize_off_diagonal) {

    //    const bool is_odd_rotation = IndexOf<domains::Cn_2D<1, 4>, List>::value == s ||

    //                                 IndexOf<domains::Cn_2D<3, 4>, List>::value == s;

    // TODO: generalize.

    const bool is_odd_rotation = (s % 2) == 1;



    return is_odd_rotation ? -1 : 1;

    const bool is_odd_rotation = s == 1 || s == 3;

    const bool is_odd_reflection = s == 4 || s == 6;



    return is_odd_reflection || is_odd_rotation ? -1 : 1;

  }



  else {