AnasaziOutputStreamTraits.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//                 Anasazi: Block Eigensolvers Package
//
// Copyright 2004 NTESS and the Anasazi contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef ANASAZI_OUTPUT_STREAM_TRAITS_HPP
#define ANASAZI_OUTPUT_STREAM_TRAITS_HPP
 
#include "AnasaziConfigDefs.hpp"
#include "AnasaziTypes.hpp"
 
#include "Teuchos_FancyOStream.hpp"
#include "Teuchos_RCP.hpp"
 
#ifdef HAVE_ANASAZI_MPI
#include <mpi.h>
#include "AnasaziGlobalComm.hpp"
#endif
 
namespace Anasazi {
 
template<class OperatorType>
struct OutputStreamTraits {
 
  // The input argument, op, is presumed to be a valid object that can be queried to
  // determine the correct output stream.
  static Teuchos::RCP<Teuchos::FancyOStream>
  getOutputStream (const OperatorType& /* op */, int rootRank = 0)
  {
    Teuchos::RCP<Teuchos::FancyOStream> fos = Teuchos::getFancyOStream(Teuchos::rcpFromRef(std::cout));
 
#ifdef HAVE_ANASAZI_MPI
    // The default implementation will output on processor 0, if parallel.
    int myRank = 0;
    int numProcs = 1;
    // Initialize MPI
    int mpiStarted = 0;
    MPI_Initialized(&mpiStarted);
    if (mpiStarted)
    {
      MPI_Comm_rank(get_global_comm(), &myRank);
      MPI_Comm_size(get_global_comm(), &numProcs);
    }
    fos->setProcRankAndSize(myRank, numProcs);
    fos->setOutputToRootOnly(rootRank);
#endif
 
    return fos;
  }
};
 
 
} // end Anasazi namespace
 
#endif
 
// end of file AnasaziOutputStreamTraits.hpp