Reduce_timestep() [5/7]

template<typename T >

void cntr::Reduce_timestep	(	int	tstp,
		int	root,
		herm_matrix_timestep< T > &	Gred,
		herm_matrix< T > &	G
	)

MPI reduce for the herm_matrix to a herm_matrix_timestep

Purpose

MPI reduce for the herm_matrix to a herm_matrix_timestep on rank root. Works for scalar or square-matrix contour objects.

Parameters

tstp	time step
root	Index of root
Gred	The reduced herm_matrix_timestep on rank root.
G	The herm_matrix on the individual ranks.

Definition at line 230 of file cntr_mpitools_impl.hpp.

References cntr::herm_matrix< T >::get_timestep(), cntr::herm_matrix< T >::nt(), cntr::herm_matrix_timestep< T >::ntau(), cntr::herm_matrix< T >::ntau(), Reduce_timestep(), cntr::herm_matrix_timestep< T >::resize(), cntr::herm_matrix_timestep< T >::size1(), cntr::herm_matrix< T >::size1(), cntr::herm_matrix_timestep< T >::size2(), cntr::herm_matrix< T >::size2(), and cntr::herm_matrix_timestep< T >::tstp().

                      {
    assert(tstp <= G.nt());
    int taskid;
    MPI_Comm_rank(MPI_COMM_WORLD, &taskid);
    if (taskid == root) {
        assert(tstp == Gred.tstp());
        assert(G.ntau() == Gred.ntau());
        assert(G.size1() == Gred.size1());
        assert(G.size2() == Gred.size2());
    }

    herm_matrix_timestep<T> Gtemp;
    Gtemp.resize(tstp, G.ntau(), G.size1());
    G.get_timestep(tstp, Gtemp);

    Reduce_timestep(tstp, root, Gred, Gtemp);

}

Here is the call graph for this function: