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

include complex gflop cost of build_g_matrix

parent a4f54704

include/dca/phys/dca_step/cluster_solver/ctaux/walker/tools/g_tools.hpp

+9 −3
Original line number Diff line number Diff line
@@ -173,9 +173,14 @@ void G_TOOLS<device_t, Parameters>::build_G_matrix(configuration_type& full_conf 
                                                 G0.ptr(0, vertex_index), LD_G0, Scalar(0.),

                                                 G.ptr(0, 0), LD_G, thread_id, stream_id);



    if constexpr (dca::util::IsComplex_t<Scalar>::value) {

      GFLOP += (8.0 * G.nrCols() * G.nrRows() * N.nrCols() + 12.0 * (G.nrCols() * G.nrRows())) * 1.0e-9;

    }

    else {

      GFLOP += 2. * G.nrCols() * G.nrRows() * N.nrCols() * 1.e-9;

    }

  }

}



/*!

 *  \f{eqnarray}{
@@ -259,7 +264,8 @@ inline void G_TOOLS<device_t, Parameters>::compute_col_on_Gamma_matrix( 
template <dca::linalg::DeviceType device_t, class Parameters>

inline double G_TOOLS<device_t, Parameters>::compute_G_vertex_to_old_vertex(

    int configuration_e_spin_index_i, int configuration_e_spin_index_j,

    dca::linalg::Matrix<Scalar, device_t>& N, std::vector<vertex_singleton_type>& configuration_e_spin) {

    dca::linalg::Matrix<Scalar, device_t>& N,

    std::vector<vertex_singleton_type>& configuration_e_spin) {

  Scalar delta = (configuration_e_spin_index_i == configuration_e_spin_index_j) ? 1. : 0.;



  vertex_singleton_type& v_j = configuration_e_spin[configuration_e_spin_index_j];