Reduce_timestep()

template<typename T >

void cntr::herm_matrix_timestep< T >::Reduce_timestep	(	int	tstp,
		int	root
	)

MPI reduce for the herm_matrix_timestep

Purpose

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

Parameters

tstp	time step
root	Index of root

Definition at line 2476 of file cntr_herm_matrix_timestep_impl.hpp.

                                                                {
   assert(tstp == tstp_);
   int taskid;
   int len = 2 * (2 * (tstp_ + 1) + ntau_ + 1) * element_size_;
   MPI_Comm_rank(MPI_COMM_WORLD, &taskid);
   if (sizeof(T) == sizeof(double)) {
      if (taskid == root) {
         MPI_Reduce(MPI_IN_PLACE, (double *)this->data_, len,
            MPI_DOUBLE_PRECISION, MPI_SUM, root, MPI_COMM_WORLD);
      } else {
         MPI_Reduce((double *)this->data_,
            (double *)this->data_, len, MPI_DOUBLE_PRECISION,
            MPI_SUM, root, MPI_COMM_WORLD);
      }
   } else {
      std::cerr << "herm_matrix_timestep<T>::MPI_Reduce only for double "
      << std::endl;
      exit(0);
   }
}