43 #ifndef IFPACK2_OVERLAPPINGROWGRAPH_DEF_HPP 44 #define IFPACK2_OVERLAPPINGROWGRAPH_DEF_HPP 46 #include <Ifpack2_Details_OverlappingRowGraph_decl.hpp> 47 #include <Tpetra_Import.hpp> 48 #include <Tpetra_Export.hpp> 54 template<
class GraphType>
57 const Teuchos::RCP<const row_graph_type>& overlappingGraph,
58 const Teuchos::RCP<const map_type>& rowMap,
59 const Teuchos::RCP<const map_type>& colMap,
60 const Tpetra::global_size_t numGlobalRows,
61 const Tpetra::global_size_t numGlobalCols,
62 const Tpetra::global_size_t numGlobalNonzeros,
63 const size_t maxNumEntries,
64 const Teuchos::RCP<const import_type>& nonoverlappingImporter,
65 const Teuchos::RCP<const import_type>& overlappingImporter) :
66 nonoverlappingGraph_ (nonoverlappingGraph),
67 overlappingGraph_ (overlappingGraph),
70 numGlobalRows_ (numGlobalRows),
71 numGlobalCols_ (numGlobalCols),
72 numGlobalNonzeros_ (numGlobalNonzeros),
73 maxNumEntries_ (maxNumEntries),
74 nonoverlappingImporter_ (nonoverlappingImporter),
75 overlappingImporter_ (overlappingImporter)
79 template<
class GraphType>
83 template<
class GraphType>
84 Teuchos::RCP<const Teuchos::Comm<int> >
87 return nonoverlappingGraph_->getComm ();
93 template<
class GraphType>
94 Teuchos::RCP<const Tpetra::Map<typename GraphType::local_ordinal_type, typename GraphType::global_ordinal_type, typename GraphType::node_type> >
101 template<
class GraphType>
102 Teuchos::RCP<const Tpetra::Map<typename GraphType::local_ordinal_type, typename GraphType::global_ordinal_type, typename GraphType::node_type> >
109 template<
class GraphType>
110 Teuchos::RCP<const Tpetra::Map<typename GraphType::local_ordinal_type, typename GraphType::global_ordinal_type, typename GraphType::node_type> >
113 return nonoverlappingGraph_->getDomainMap ();
117 template<
class GraphType>
118 Teuchos::RCP<const Tpetra::Map<typename GraphType::local_ordinal_type, typename GraphType::global_ordinal_type, typename GraphType::node_type> >
121 return nonoverlappingGraph_->getRangeMap ();
125 template<
class GraphType>
126 Teuchos::RCP<const Tpetra::Import<typename GraphType::local_ordinal_type, typename GraphType::global_ordinal_type, typename GraphType::node_type> >
129 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
"Not implemented");
133 template<
class GraphType>
134 Teuchos::RCP<const Tpetra::Export<typename GraphType::local_ordinal_type, typename GraphType::global_ordinal_type, typename GraphType::node_type> >
137 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
"Not implemented");
141 template<
class GraphType>
144 return numGlobalRows_;
148 template<
class GraphType>
151 return numGlobalCols_;
155 template<
class GraphType>
158 return nonoverlappingGraph_->getNodeNumRows () +
159 overlappingGraph_->getNodeNumRows ();
163 template<
class GraphType>
166 return this->getNodeNumRows ();
170 template<
class GraphType>
171 typename GraphType::global_ordinal_type
174 return nonoverlappingGraph_->getIndexBase ();
178 template<
class GraphType>
181 return numGlobalNonzeros_;
185 template<
class GraphType>
188 return nonoverlappingGraph_->getNodeNumEntries () +
189 overlappingGraph_->getNodeNumEntries ();
193 template<
class GraphType>
198 const local_ordinal_type localRow = rowMap_->getLocalElement (globalRow);
199 if (localRow == Teuchos::OrdinalTraits<local_ordinal_type>::invalid ()) {
200 return Teuchos::OrdinalTraits<size_t>::invalid();
202 return getNumEntriesInLocalRow (localRow);
207 template<
class GraphType>
213 const size_t numMyRowsA = nonoverlappingGraph_->getNodeNumRows ();
214 if (as<size_t> (localRow) < numMyRowsA) {
215 return nonoverlappingGraph_->getNumEntriesInLocalRow (localRow);
217 return overlappingGraph_->getNumEntriesInLocalRow (as<local_ordinal_type> (localRow - numMyRowsA));
222 template<
class GraphType>
225 throw std::runtime_error(
"Ifpack2::OverlappingRowGraph::getGlobalMaxNumRowEntries() not supported.");
229 template<
class GraphType>
232 return maxNumEntries_;
236 template<
class GraphType>
243 template<
class GraphType>
250 template<
class GraphType>
257 template<
class GraphType>
263 template<
class GraphType>
267 nonconst_global_inds_host_view_type& indices,
268 size_t& numIndices)
const 270 const local_ordinal_type localRow = rowMap_->getLocalElement (globalRow);
271 if (localRow == Teuchos::OrdinalTraits<local_ordinal_type>::invalid ()) {
272 numIndices = Teuchos::OrdinalTraits<size_t>::invalid ();
274 if (Teuchos::as<size_t> (localRow) < nonoverlappingGraph_->getNodeNumRows ()) {
275 nonoverlappingGraph_->getGlobalRowCopy (globalRow, indices, numIndices);
277 overlappingGraph_->getGlobalRowCopy (globalRow, indices, numIndices);
282 #ifdef TPETRA_ENABLE_DEPRECATED_CODE 283 template<
class GraphType>
287 const Teuchos::ArrayView<global_ordinal_type>& indices,
288 size_t& numIndices)
const 290 const local_ordinal_type localRow = rowMap_->getLocalElement (globalRow);
291 if (localRow == Teuchos::OrdinalTraits<local_ordinal_type>::invalid ()) {
292 numIndices = Teuchos::OrdinalTraits<size_t>::invalid ();
294 if (Teuchos::as<size_t> (localRow) < nonoverlappingGraph_->getNodeNumRows ()) {
295 nonoverlappingGraph_->getGlobalRowCopy (globalRow, indices, numIndices);
297 overlappingGraph_->getGlobalRowCopy (globalRow, indices, numIndices);
303 template<
class GraphType>
307 nonconst_local_inds_host_view_type& indices,
308 size_t& numIndices)
const 311 const size_t numMyRowsA = nonoverlappingGraph_->getNodeNumRows ();
312 if (as<size_t> (localRow) < numMyRowsA) {
313 nonoverlappingGraph_->getLocalRowCopy (localRow, indices, numIndices);
315 const local_ordinal_type localRowOffset =
316 localRow - as<local_ordinal_type> (numMyRowsA);
317 overlappingGraph_->getLocalRowCopy (localRowOffset, indices, numIndices);
322 #ifdef TPETRA_ENABLE_DEPRECATED_CODE 323 template<
class GraphType>
327 const Teuchos::ArrayView<local_ordinal_type>& indices,
328 size_t& numIndices)
const 331 const size_t numMyRowsA = nonoverlappingGraph_->getNodeNumRows ();
332 if (as<size_t> (localRow) < numMyRowsA) {
333 nonoverlappingGraph_->getLocalRowCopy (localRow, indices, numIndices);
335 const local_ordinal_type localRowOffset =
336 localRow - as<local_ordinal_type> (numMyRowsA);
337 overlappingGraph_->getLocalRowCopy (localRowOffset, indices, numIndices);
342 template<
class GraphType>
346 global_inds_host_view_type &indices)
const {
347 const local_ordinal_type LocalRow = rowMap_->getLocalElement (GlobalRow);
348 if (LocalRow == Teuchos::OrdinalTraits<local_ordinal_type>::invalid()) {
349 indices = global_inds_host_view_type();
351 if (Teuchos::as<size_t> (LocalRow) < nonoverlappingGraph_->getNodeNumRows ()) {
352 nonoverlappingGraph_->getGlobalRowView (GlobalRow, indices);
354 overlappingGraph_->getGlobalRowView (GlobalRow, indices);
359 #ifdef TPETRA_ENABLE_DEPRECATED_CODE 360 template<
class GraphType>
364 Teuchos::ArrayView<const global_ordinal_type>& indices)
const 366 const local_ordinal_type LocalRow = rowMap_->getLocalElement (GlobalRow);
367 if (LocalRow == Teuchos::OrdinalTraits<local_ordinal_type>::invalid()) {
368 indices = Teuchos::null;
370 if (Teuchos::as<size_t> (LocalRow) < nonoverlappingGraph_->getNodeNumRows ()) {
371 nonoverlappingGraph_->getGlobalRowView (GlobalRow, indices);
373 overlappingGraph_->getGlobalRowView (GlobalRow, indices);
380 template<
class GraphType>
384 local_inds_host_view_type & indices)
const {
386 const size_t numMyRowsA = nonoverlappingGraph_->getNodeNumRows ();
387 if (as<size_t> (LocalRow) < numMyRowsA) {
388 nonoverlappingGraph_->getLocalRowView (LocalRow, indices);
390 overlappingGraph_->getLocalRowView (LocalRow - as<local_ordinal_type> (numMyRowsA),
396 #ifdef TPETRA_ENABLE_DEPRECATED_CODE 397 template<
class GraphType>
401 Teuchos::ArrayView<const local_ordinal_type>& indices)
const 404 const size_t numMyRowsA = nonoverlappingGraph_->getNodeNumRows ();
405 if (as<size_t> (LocalRow) < numMyRowsA) {
406 nonoverlappingGraph_->getLocalRowView (LocalRow, indices);
408 overlappingGraph_->getLocalRowView (LocalRow - as<local_ordinal_type> (numMyRowsA),
418 #define IFPACK2_DETAILS_OVERLAPPINGROWGRAPH_INSTANT(LO,GO,N) \ 419 template class Ifpack2::Details::OverlappingRowGraph<Tpetra::CrsGraph< LO, GO, N > >; \ 420 template class Ifpack2::Details::OverlappingRowGraph<Tpetra::RowGraph< LO, GO, N > >; 422 #endif // IFPACK2_OVERLAPPINGROWGRAPH_DEF_HPP virtual size_t getNodeNumCols() const
The number of columns owned by the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:164
virtual bool isGloballyIndexed() const
Whether this graph is globally indexed.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:251
virtual size_t getNodeNumEntries() const
The number of entries in this graph owned by the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:186
virtual void getLocalRowCopy(local_ordinal_type localRow, nonconst_local_inds_host_view_type &gblColInds, size_t &numIndices) const
Copy out a list of local column indices in the given local row that are owned by the calling process...
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:306
OverlappingRowGraph(const Teuchos::RCP< const row_graph_type > &nonoverlappingGraph, const Teuchos::RCP< const row_graph_type > &overlappingGraph, const Teuchos::RCP< const map_type > &rowMap, const Teuchos::RCP< const map_type > &colMap, const Tpetra::global_size_t numGlobalRows, const Tpetra::global_size_t numGlobalCols, const Tpetra::global_size_t numGlobalNonzeros, const size_t maxNumEntries, const Teuchos::RCP< const import_type > &nonoverlappingImporter, const Teuchos::RCP< const import_type > &overlappingImporter)
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:56
virtual bool isLocallyIndexed() const
Whether this graph is locally indexed.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:244
virtual size_t getNodeNumRows() const
The number of rows owned by the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:156
virtual global_ordinal_type getIndexBase() const
The index base for global indices for this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:172
virtual Teuchos::RCP< const map_type > getRangeMap() const
The Map that describes the range of this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:119
virtual Teuchos::RCP< const import_type > getImporter() const
Import object (from domain Map to column Map).
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:127
virtual bool hasColMap() const
Whether this graph has a column Map.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:237
Ifpack2 implementation details.
virtual bool isFillComplete() const
true if fillComplete() has been called, else false.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:258
virtual Teuchos::RCP< const map_type > getRowMap() const
The Map that describes the distribution of rows over processes.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:95
virtual size_t getNumEntriesInLocalRow(local_ordinal_type localRow) const
The number of entries in the given local row that are owned by the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:210
virtual Teuchos::RCP< const map_type > getColMap() const
The Map that describes the distribution of columns over processes.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:103
virtual void getGlobalRowCopy(global_ordinal_type globalRow, nonconst_global_inds_host_view_type &gblColInds, size_t &numIndices) const
Copy out a list of column indices in the given global row that are owned by the calling process...
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:266
virtual Teuchos::RCP< const Teuchos::Comm< int > > getComm() const
The communicator over which the graph is distributed.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:85
virtual global_size_t getGlobalNumEntries() const
The global number of entries in this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:179
virtual ~OverlappingRowGraph()
Destructor.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:80
virtual Teuchos::RCP< const export_type > getExporter() const
Export object (from row Map to range Map).
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:135
virtual void getGlobalRowView(const global_ordinal_type gblRow, global_inds_host_view_type &gblColInds) const
Get a const, non-persisting view of the given global row's global column indices, as a Teuchos::Array...
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:345
virtual global_size_t getGlobalNumCols() const
The global number of columns in this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:149
virtual Teuchos::RCP< const map_type > getDomainMap() const
The Map that describes the domain of this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:111
virtual global_size_t getGlobalNumRows() const
The global number of rows in this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:142
virtual size_t getNumEntriesInGlobalRow(global_ordinal_type globalRow) const
The number of entries in the given global row that are owned by the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:196
Preconditioners and smoothers for Tpetra sparse matrices.
Definition: Ifpack2_AdditiveSchwarz_decl.hpp:73
virtual size_t getGlobalMaxNumRowEntries() const
The maximum number of entries in any row on any process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:223
virtual void getLocalRowView(const local_ordinal_type lclRow, local_inds_host_view_type &lclColInds) const
Get a constant, nonpersisting, locally indexed view of the given row of the graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:383
virtual size_t getNodeMaxNumRowEntries() const
The maximum number of entries in any row on the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:230