Ifpack2 Templated Preconditioning Package  Version 1.0
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
Ifpack2_Details_OverlappingRowGraph_def.hpp
1 /*@HEADER
2 // ***********************************************************************
3 //
4 // Ifpack2: Templated Object-Oriented Algebraic Preconditioner Package
5 // Copyright (2009) Sandia Corporation
6 //
7 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
8 // license for use of this work by or on behalf of the U.S. Government.
9 //
10 // Redistribution and use in source and binary forms, with or without
11 // modification, are permitted provided that the following conditions are
12 // met:
13 //
14 // 1. Redistributions of source code must retain the above copyright
15 // notice, this list of conditions and the following disclaimer.
16 //
17 // 2. Redistributions in binary form must reproduce the above copyright
18 // notice, this list of conditions and the following disclaimer in the
19 // documentation and/or other materials provided with the distribution.
20 //
21 // 3. Neither the name of the Corporation nor the names of the
22 // contributors may be used to endorse or promote products derived from
23 // this software without specific prior written permission.
24 //
25 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36 //
37 // Questions? Contact Michael A. Heroux (maherou@sandia.gov)
38 //
39 // ***********************************************************************
40 //@HEADER
41 */
42 
43 #ifndef IFPACK2_OVERLAPPINGROWGRAPH_DEF_HPP
44 #define IFPACK2_OVERLAPPINGROWGRAPH_DEF_HPP
45 
46 #include <Ifpack2_Details_OverlappingRowGraph_decl.hpp>
47 #include <Tpetra_Import.hpp>
48 #include <Tpetra_Export.hpp>
49 
50 
51 namespace Ifpack2 {
52 namespace Details {
53 
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),
68  rowMap_ (rowMap),
69  colMap_ (colMap),
70  numGlobalRows_ (numGlobalRows),
71  numGlobalCols_ (numGlobalCols),
72  numGlobalNonzeros_ (numGlobalNonzeros),
73  maxNumEntries_ (maxNumEntries),
74  nonoverlappingImporter_ (nonoverlappingImporter),
75  overlappingImporter_ (overlappingImporter)
76 {}
77 
78 
79 template<class GraphType>
81 
82 
83 template<class GraphType>
86 {
87  return nonoverlappingGraph_->getComm ();
88 }
89 
90 
91 #ifdef TPETRA_ENABLE_DEPRECATED_CODE
92 template<class GraphType>
93 TPETRA_DEPRECATED
96 {
97  return Teuchos::null;
98 }
99 #endif // TPETRA_ENABLE_DEPRECATED_CODE
100 
101 
102 template<class GraphType>
105 {
106  return rowMap_;
107 }
108 
109 
110 template<class GraphType>
113 {
114  return colMap_;
115 }
116 
117 
118 template<class GraphType>
121 {
122  return nonoverlappingGraph_->getDomainMap ();
123 }
124 
125 
126 template<class GraphType>
129 {
130  return nonoverlappingGraph_->getRangeMap ();
131 }
132 
133 
134 template<class GraphType>
137 {
138  TEUCHOS_TEST_FOR_EXCEPTION(true, std::logic_error, "Not implemented");
139 }
140 
141 
142 template<class GraphType>
145 {
146  TEUCHOS_TEST_FOR_EXCEPTION(true, std::logic_error, "Not implemented");
147 }
148 
149 
150 template<class GraphType>
152 {
153  return numGlobalRows_;
154 }
155 
156 
157 template<class GraphType>
159 {
160  return numGlobalCols_;
161 }
162 
163 
164 template<class GraphType>
166 {
167  return nonoverlappingGraph_->getNodeNumRows () +
168  overlappingGraph_->getNodeNumRows ();
169 }
170 
171 
172 template<class GraphType>
174 {
175  return this->getNodeNumRows ();
176 }
177 
178 
179 template<class GraphType>
180 typename GraphType::global_ordinal_type
182 {
183  return nonoverlappingGraph_->getIndexBase ();
184 }
185 
186 
187 template<class GraphType>
189 {
190  return numGlobalNonzeros_;
191 }
192 
193 
194 template<class GraphType>
196 {
197  return nonoverlappingGraph_->getNodeNumEntries () +
198  overlappingGraph_->getNodeNumEntries ();
199 }
200 
201 
202 template<class GraphType>
203 size_t
205 getNumEntriesInGlobalRow (global_ordinal_type globalRow) const
206 {
207  const local_ordinal_type localRow = rowMap_->getLocalElement (globalRow);
210  } else {
211  return getNumEntriesInLocalRow (localRow);
212  }
213 }
214 
215 
216 template<class GraphType>
217 size_t
219 getNumEntriesInLocalRow (local_ordinal_type localRow) const
220 {
221  using Teuchos::as;
222  const size_t numMyRowsA = nonoverlappingGraph_->getNodeNumRows ();
223  if (as<size_t> (localRow) < numMyRowsA) {
224  return nonoverlappingGraph_->getNumEntriesInLocalRow (localRow);
225  } else {
226  return overlappingGraph_->getNumEntriesInLocalRow (as<local_ordinal_type> (localRow - numMyRowsA));
227  }
228 }
229 
230 
231 template<class GraphType>
233 {
234  throw std::runtime_error("Ifpack2::OverlappingRowGraph::getGlobalMaxNumRowEntries() not supported.");
235 }
236 
237 
238 template<class GraphType>
240 {
241  return maxNumEntries_;
242 }
243 
244 
245 template<class GraphType>
247 {
248  return true;
249 }
250 
251 
252 template<class GraphType>
254 {
255  return true;
256 }
257 
258 
259 template<class GraphType>
261 {
262  return false;
263 }
264 
265 
266 template<class GraphType>
268 {
269  return true;
270 }
271 
272 
273 template<class GraphType>
274 void
276 getGlobalRowCopy (global_ordinal_type globalRow,
278  size_t& numIndices) const
279 {
280  const local_ordinal_type localRow = rowMap_->getLocalElement (globalRow);
283  } else {
284  if (Teuchos::as<size_t> (localRow) < nonoverlappingGraph_->getNodeNumRows ()) {
285  nonoverlappingGraph_->getGlobalRowCopy (globalRow, indices, numIndices);
286  } else {
287  overlappingGraph_->getGlobalRowCopy (globalRow, indices, numIndices);
288  }
289  }
290 }
291 
292 
293 template<class GraphType>
294 void
296 getLocalRowCopy (local_ordinal_type localRow,
298  size_t& numIndices) const
299 {
300  using Teuchos::as;
301  const size_t numMyRowsA = nonoverlappingGraph_->getNodeNumRows ();
302  if (as<size_t> (localRow) < numMyRowsA) {
303  nonoverlappingGraph_->getLocalRowCopy (localRow, indices, numIndices);
304  } else {
305  const local_ordinal_type localRowOffset =
306  localRow - as<local_ordinal_type> (numMyRowsA);
307  overlappingGraph_->getLocalRowCopy (localRowOffset, indices, numIndices);
308  }
309 }
310 
311 } // namespace Details
312 } // namespace Ifpack2
313 
314 #define IFPACK2_DETAILS_OVERLAPPINGROWGRAPH_INSTANT(LO,GO,N) \
315  template class Ifpack2::Details::OverlappingRowGraph<Tpetra::CrsGraph< LO, GO, N > >; \
316  template class Ifpack2::Details::OverlappingRowGraph<Tpetra::RowGraph< LO, GO, N > >;
317 
318 #endif // IFPACK2_OVERLAPPINGROWGRAPH_DEF_HPP
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
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
virtual size_t getGlobalMaxNumRowEntries() const
The maximum number of entries in any row on any process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:232
virtual Teuchos::RCP< const map_type > getRowMap() const
The Map that describes the distribution of rows over processes.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:104
virtual size_t getNodeMaxNumRowEntries() const
The maximum number of entries in any row on the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:239
virtual bool isGloballyIndexed() const
Whether this graph is globally indexed.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:260
virtual bool hasColMap() const
Whether this graph has a column Map.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:246
virtual Teuchos::RCP< const export_type > getExporter() const
Export object (from row Map to range Map).
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:144
virtual size_t getNodeNumRows() const
The number of rows owned by the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:165
virtual global_ordinal_type getIndexBase() const
The index base for global indices for this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:181
virtual size_t getNodeNumEntries() const
The number of entries in this graph owned by the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:195
virtual Teuchos::RCP< const map_type > getDomainMap() const
The Map that describes the domain of this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:120
Sparse graph (Tpetra::RowGraph subclass) with ghost rows.
Definition: Ifpack2_Details_OverlappingRowGraph_decl.hpp:65
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 ~OverlappingRowGraph()
Destructor.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:80
virtual bool isLocallyIndexed() const
Whether this graph is locally indexed.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:253
virtual void getGlobalRowCopy(global_ordinal_type globalRow, const Teuchos::ArrayView< global_ordinal_type > &indices, 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:276
virtual bool isFillComplete() const
true if fillComplete() has been called, else false.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:267
virtual global_size_t getGlobalNumCols() const
The global number of columns in this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:158
virtual size_t getNodeNumCols() const
The number of columns owned by the calling process.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:173
virtual void getLocalRowCopy(local_ordinal_type localRow, const Teuchos::ArrayView< local_ordinal_type > &indices, 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:296
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:219
TypeTo as(const TypeFrom &t)
virtual Teuchos::RCP< const import_type > getImporter() const
Import object (from domain Map to column Map).
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:136
virtual Teuchos::RCP< const map_type > getColMap() const
The Map that describes the distribution of columns over processes.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:112
virtual global_size_t getGlobalNumRows() const
The global number of rows in this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:151
virtual Teuchos::RCP< const map_type > getRangeMap() const
The Map that describes the range of this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:128
virtual global_size_t getGlobalNumEntries() const
The global number of entries in this graph.
Definition: Ifpack2_Details_OverlappingRowGraph_def.hpp:188
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:205