Teko_BlockingEpetra.hpp

// @HEADER
// *****************************************************************************
//      Teko: A package for block and physics based preconditioning
//
// Copyright 2010 NTESS and the Teko contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef __Teko_BlockingEpetra_hpp__
#define __Teko_BlockingEpetra_hpp__
 
#include "Teuchos_RCP.hpp"
 
// include basic Epetra information
#include "Epetra_Map.h"
#include "Epetra_Comm.h"
#include "Epetra_CrsMatrix.h"
#include "Epetra_Vector.h"
#include "Epetra_MultiVector.h"
#include "Epetra_Export.h"
#include "Epetra_Import.h"
#include "Epetra_IntVector.h"
 
#include <vector>
 
namespace Teko {
namespace Epetra {
namespace Blocking {
 
typedef std::pair<Teuchos::RCP<Epetra_Map>, Teuchos::RCP<Epetra_Map> > MapPair;
typedef std::pair<Teuchos::RCP<Epetra_Import>, Teuchos::RCP<Epetra_Export> > ImExPair;

const MapPair buildSubMap(const std::vector<int> &gid, const Epetra_Comm &comm);

const ImExPair buildExportImport(const Epetra_Map &baseMap, const MapPair &maps);

void many2one(Epetra_MultiVector &one,
              const std::vector<Teuchos::RCP<const Epetra_MultiVector> > &many,
              const std::vector<Teuchos::RCP<Epetra_Export> > &subExport);

void one2many(std::vector<Teuchos::RCP<Epetra_MultiVector> > &many,
              const Epetra_MultiVector &single,
              const std::vector<Teuchos::RCP<Epetra_Import> > &subImport);

void buildSubVectors(const std::vector<MapPair> &maps,
                     std::vector<Teuchos::RCP<Epetra_MultiVector> > &vectors, int count);

Teuchos::RCP<Epetra_IntVector> getSubBlockColumnGIDs(const Epetra_CrsMatrix &A,
                                                     const MapPair &mapPair);

Teuchos::RCP<Epetra_CrsMatrix> buildSubBlock(int i, int j, const Epetra_CrsMatrix &A,
                                             const std::vector<MapPair> &subMaps);

void rebuildSubBlock(int i, int j, const Epetra_CrsMatrix &A, const std::vector<MapPair> &subMaps,
                     Epetra_CrsMatrix &mat);
 
}  // namespace Blocking
}  // namespace Epetra
}  // namespace Teko
 
#endif