Commit f1e22413 authored by Doak, Peter W.'s avatar Doak, Peter W.
Browse files

possible sigma update bugfix

parent 0938a196

include/dca/phys/dca_data/dca_data.hpp

+3 −4
Original line number Diff line number Diff line
@@ -587,14 +587,13 @@ void DcaData<Parameters, DT>::initialize_G0() { 


    util::Timer("G_0 initialization", concurrency_.id() == concurrency_.first());



    Real chemical_potential = parameters_.get_chemical_potential();

    Real beta = parameters_.get_beta();

    // Compute G0_k_w.

    compute_G0_k_w(H_DCA, chemical_potential, parameters_.get_coarsegraining_threads(), G0_k_w);

    compute_G0_k_w(H_DCA, parameters_.get_chemical_potential(),

                   parameters_.get_coarsegraining_threads(), G0_k_w);

    Symmetrize<Parameters>::execute(G0_k_w, H_symmetry, true);



    // Compute G0_k_t.

    compute_G0_k_t(H_DCA, chemical_potential, beta, G0_k_t);

    compute_G0_k_t(H_DCA, parameters_.get_chemical_potential(), parameters_.get_beta(), G0_k_t);

    Symmetrize<Parameters>::execute(G0_k_t, H_symmetry, true);



    // Compute G0_r_w.

include/dca/phys/dca_step/cluster_solver/ctaux/ctaux_cluster_solver.hpp

+1 −1
Original line number Diff line number Diff line
@@ -286,7 +286,7 @@ template <dca::linalg::DeviceType device_t, class Parameters, class Data, DistTy 
template <typename dca_info_struct_t>

double CtauxClusterSolver<device_t, Parameters, Data, DIST>::finalize(

    dca_info_struct_t& dca_info_struct) {

  //collect_measurements();

  collect_measurements();

  symmetrize_measurements();



  // Compute new Sigma.

test/unit/math/function_transform/function_transform_test.cpp

+4 −2
Original line number Diff line number Diff line
@@ -31,8 +31,9 @@ using McOptions = MockMcOptions<double>; 
#include "dca/parallel/no_concurrency/no_concurrency.hpp"

#include "dca/parallel/no_threading/no_threading.hpp"

#include "dca/profiling/null_profiler.hpp"



#ifdef DCA_HAVE_ADIOS2

adios2::ADIOS* adios_ptr;

#endif

#ifdef DCA_HAVE_MPI

#include "dca/parallel/mpi_concurrency/mpi_concurrency.hpp"

dca::parallel::MPIConcurrency* concurrency_ptr;
@@ -167,10 +168,11 @@ int main(int argc, char** argv) { 
  concurrency_ptr = &concurrency;

#endif



  #ifdef DCA_HAVE_ADIOS2

  //ADIOS expects MPI_COMM pointer or nullptr

  adios2::ADIOS adios("", concurrency_ptr->get(), "C++");

  adios_ptr = &adios;



#endif

  ::testing::InitGoogleTest(&argc, argv);



  // ::testing::TestEventListeners& listeners = ::testing::UnitTest::GetInstance()->listeners();

test/unit/phys/dca_step/cluster_solver/shared_tools/accumulation/tp/CMakeLists.txt

+3 −2
Original line number Diff line number Diff line
@@ -13,14 +13,15 @@ if (DCA_HAVE_ADIOS2) 
    INCLUDE_DIRS ${DCA_INCLUDE_DIRS};${PROJECT_SOURCE_DIR}

    LIBS     ${DCA_LIBS}

    )

endif()



dca_add_gtest(tp_accumulator_test

  FAST

  INCLUDE_DIRS ${DCA_INCLUDE_DIRS};${PROJECT_SOURCE_DIR}

  LIBS     ${DCA_LIBS} ${DCA_GPU_LIBS}

  )



  endif()





dca_add_gtest(tp_accumulator_test_complex

  FAST

  GTEST_MAIN