Xpetra Version of the Day
Xpetra_CrsMatrixFactory.hpp
Go to the documentation of this file.
1// @HEADER
2//
3// ***********************************************************************
4//
5// Xpetra: A linear algebra interface package
6// Copyright 2012 Sandia Corporation
7//
8// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
9// the U.S. Government retains certain rights in this software.
10//
11// Redistribution and use in source and binary forms, with or without
12// modification, are permitted provided that the following conditions are
13// met:
14//
15// 1. Redistributions of source code must retain the above copyright
16// notice, this list of conditions and the following disclaimer.
17//
18// 2. Redistributions in binary form must reproduce the above copyright
19// notice, this list of conditions and the following disclaimer in the
20// documentation and/or other materials provided with the distribution.
21//
22// 3. Neither the name of the Corporation nor the names of the
23// contributors may be used to endorse or promote products derived from
24// this software without specific prior written permission.
25//
26// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37//
38// Questions? Contact
39// Jonathan Hu (jhu@sandia.gov)
40// Andrey Prokopenko (aprokop@sandia.gov)
41// Ray Tuminaro (rstumin@sandia.gov)
42//
43// ***********************************************************************
44//
45// @HEADER
46#ifndef XPETRA_CRSMATRIXFACTORY_HPP
47#define XPETRA_CRSMATRIXFACTORY_HPP
48
49#include "Xpetra_ConfigDefs.hpp"
50
51#include "Xpetra_CrsMatrix.hpp"
52
53#ifdef HAVE_XPETRA_TPETRA
54#include "Xpetra_TpetraCrsMatrix.hpp"
55#endif
56
57#ifdef HAVE_XPETRA_EPETRA
58#include "Xpetra_EpetraCrsMatrix.hpp"
59#endif
60
61#include "Xpetra_Exceptions.hpp"
62
63namespace Xpetra {
64
65 template <class Scalar,
66 class LocalOrdinal,
67 class GlobalOrdinal,
68 class Node = KokkosClassic::DefaultNode::DefaultNodeType>
69 class CrsMatrixFactory {
70 private:
72 CrsMatrixFactory() {}
73
74 public:
76 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
77 Build (const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap)
78 {
79 TEUCHOS_TEST_FOR_EXCEPTION(rowMap->lib() == UseEpetra, std::logic_error,
80 "Can't create Xpetra::EpetraCrsMatrix with these scalar/LO/GO types");
81#ifdef HAVE_XPETRA_TPETRA
82 if (rowMap->lib() == UseTpetra)
83 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, 0) );
84#endif
85
86 XPETRA_FACTORY_END;
87 }
88
90 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
91 Build (const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap,
92 size_t maxNumEntriesPerRow,
93 const Teuchos::RCP<Teuchos::ParameterList>& plist = Teuchos::null)
94 {
95 XPETRA_MONITOR("CrsMatrixFactory::Build");
96
97#ifdef HAVE_XPETRA_TPETRA
98 if (rowMap->lib() == UseTpetra)
99 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, maxNumEntriesPerRow, plist) );
100#endif
101
102 XPETRA_FACTORY_ERROR_IF_EPETRA(rowMap->lib());
103 XPETRA_FACTORY_END;
104 }
105
107 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
108 Build (const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > >& rowMap,
109 const ArrayRCP<const size_t>& NumEntriesPerRowToAlloc,
110 const Teuchos::RCP<Teuchos::ParameterList>& plist = Teuchos::null)
111 {
112#ifdef HAVE_XPETRA_TPETRA
113 if (rowMap->lib() == UseTpetra)
114 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, NumEntriesPerRowToAlloc, plist) );
115#endif
116
117 XPETRA_FACTORY_ERROR_IF_EPETRA(rowMap->lib());
118 XPETRA_FACTORY_END;
119 }
120
122 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
123 Build (const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& rowMap,
124 const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >& colMap,
125 size_t maxNumEntriesPerRow,
126 const Teuchos::RCP<Teuchos::ParameterList>& plist = Teuchos::null)
127 {
128 XPETRA_MONITOR("CrsMatrixFactory::Build");
129
130#ifdef HAVE_XPETRA_TPETRA
131 if (rowMap->lib() == UseTpetra)
132 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
133#endif
134
135 XPETRA_FACTORY_ERROR_IF_EPETRA(rowMap->lib());
136 XPETRA_FACTORY_END;
137 }
138
140 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
141 XPETRA_MONITOR("CrsMatrixFactory::Build");
142
143#ifdef HAVE_XPETRA_TPETRA
144 if (rowMap->lib() == UseTpetra)
145 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
146#endif
147
148 XPETRA_FACTORY_ERROR_IF_EPETRA(rowMap->lib());
149 XPETRA_FACTORY_END;
150 }
151
153 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
154 XPETRA_MONITOR("CrsMatrixFactory::Build");
155
156#ifdef HAVE_XPETRA_TPETRA
157 if (graph->getRowMap()->lib() == UseTpetra)
158 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(graph, plist) );
159#endif
160
161 XPETRA_FACTORY_ERROR_IF_EPETRA(graph->getRowMap()->lib());
162 XPETRA_FACTORY_END;
163 }
164
165
166 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
167 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
168 const Import<LocalOrdinal,GlobalOrdinal,Node> &importer,
169 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap = Teuchos::null,
170 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap = Teuchos::null,
171 const Teuchos::RCP<Teuchos::ParameterList>& params = Teuchos::null) {
172 XPETRA_MONITOR("CrsMatrixFactory::Build");
173
174#ifdef HAVE_XPETRA_TPETRA
175 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
176 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,importer,domainMap,rangeMap,params));
177#endif
178
179 XPETRA_FACTORY_ERROR_IF_EPETRA(sourceMatrix->getRowMap()->lib());
180 XPETRA_FACTORY_END;
181 }
182
183 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
184 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
185 const Export<LocalOrdinal,GlobalOrdinal,Node> &exporter,
186 const RCP<Map<LocalOrdinal,GlobalOrdinal,Scalar> > & domainMap = Teuchos::null,
187 const RCP<Map<LocalOrdinal,GlobalOrdinal,Scalar> > & rangeMap = Teuchos::null,
188 const Teuchos::RCP<Teuchos::ParameterList>& params = Teuchos::null) {
189 XPETRA_MONITOR("CrsMatrixFactory::Build");
190
191#ifdef HAVE_XPETRA_TPETRA
192 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
193 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,exporter,domainMap,rangeMap,params));
194#endif
195
196 XPETRA_FACTORY_ERROR_IF_EPETRA(sourceMatrix->getRowMap()->lib());
197 XPETRA_FACTORY_END;
198 }
199
200 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
201 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
202 const Import<LocalOrdinal,GlobalOrdinal,Node> &RowImporter,
203 const RCP<const Import<LocalOrdinal,GlobalOrdinal,Node> > DomainImporter,
204 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap,
205 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap,
206 const Teuchos::RCP<Teuchos::ParameterList>& params) {
207 XPETRA_MONITOR("CrsMatrixFactory::Build");
208
209#ifdef HAVE_XPETRA_TPETRA
210 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
211 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,RowImporter,DomainImporter,domainMap,rangeMap,params));
212#endif
213
214 XPETRA_FACTORY_ERROR_IF_EPETRA(sourceMatrix->getRowMap()->lib());
215 XPETRA_FACTORY_END;
216 }
217
218 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
219 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
220 const Export<LocalOrdinal,GlobalOrdinal,Node> &RowExporter,
221 const RCP<const Export<LocalOrdinal,GlobalOrdinal,Node> > DomainExporter,
222 const RCP<Map<LocalOrdinal,GlobalOrdinal,Scalar> > & domainMap,
223 const RCP<Map<LocalOrdinal,GlobalOrdinal,Scalar> > & rangeMap,
224 const Teuchos::RCP<Teuchos::ParameterList>& params) {
225 XPETRA_MONITOR("CrsMatrixFactory::Build");
226
227#ifdef HAVE_XPETRA_TPETRA
228 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
229 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,RowExporter,DomainExporter,domainMap,rangeMap,params));
230#endif
231
232 XPETRA_FACTORY_ERROR_IF_EPETRA(sourceMatrix->getRowMap()->lib());
233 XPETRA_FACTORY_END;
234 }
235
236#ifdef HAVE_XPETRA_KOKKOS_REFACTOR
237 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build (
238 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rowMap,
239 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& colMap,
240 const typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type& lclMatrix,
241 const Teuchos::RCP<Teuchos::ParameterList>& params = null) {
242 XPETRA_MONITOR("CrsMatrixFactory::Build");
243
244#ifdef HAVE_XPETRA_TPETRA
245 if (rowMap->lib() == UseTpetra)
246 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, lclMatrix, params));
247#endif
248
249 XPETRA_FACTORY_ERROR_IF_EPETRA(rowMap->lib());
250 XPETRA_FACTORY_END;
251 }
252
253 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build (
254 const typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type& lclMatrix,
255 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rowMap,
256 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& colMap,
257 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& domainMap = Teuchos::null,
258 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rangeMap = Teuchos::null,
259 const Teuchos::RCP<Teuchos::ParameterList>& params = null) {
260 XPETRA_MONITOR("CrsMatrixFactory::Build");
261
262#ifdef HAVE_XPETRA_TPETRA
263 if (rowMap->lib() == UseTpetra)
264 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(lclMatrix, rowMap, colMap, domainMap, rangeMap, params));
265#endif
266
267 XPETRA_FACTORY_ERROR_IF_EPETRA(rowMap->lib());
268 XPETRA_FACTORY_END;
269 }
270
271 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build (
272 const typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type& lclMatrix,
273 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rowMap,
274 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& colMap,
275 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& domainMap,
276 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rangeMap,
277 const Teuchos::RCP<const Import<LocalOrdinal,GlobalOrdinal,Node>>& importer,
278 const Teuchos::RCP<const Export<LocalOrdinal,GlobalOrdinal,Node>>& exporter,
279 const Teuchos::RCP<Teuchos::ParameterList>& params = null) {
280 XPETRA_MONITOR("CrsMatrixFactory::Build");
281
282#ifdef HAVE_XPETRA_TPETRA
283 if (rowMap->lib() == UseTpetra)
284 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(lclMatrix, rowMap, colMap, domainMap, rangeMap, importer, exporter, params));
285#endif
286
287 TEUCHOS_TEST_FOR_EXCEPTION(rowMap->lib() == UseEpetra, std::logic_error, "Epetra doesn't support this matrix constructor");
288
289 XPETRA_FACTORY_END;
290 }
291#endif
292
293 };
294
295// we need the Epetra specialization only if Epetra is enabled
296#if (defined(HAVE_XPETRA_EPETRA) && !defined(XPETRA_EPETRA_NO_32BIT_GLOBAL_INDICES))
297
298 // Specializtion for SC=double, LO=int, GO=int and Node=EpetraNode
299 // Used both for Epetra and Tpetra
300 template <>
301 class CrsMatrixFactory<double, int, int, EpetraNode> {
302 typedef double Scalar;
303 typedef int LocalOrdinal;
304 typedef int GlobalOrdinal;
305 typedef EpetraNode Node;
306
307 private:
309 CrsMatrixFactory() {}
310
311 public:
313 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
314 Build (const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap)
315 {
316 XPETRA_MONITOR("CrsMatrixFactory::Build");
317#ifdef HAVE_XPETRA_TPETRA
318 if (rowMap->lib() == UseTpetra)
319 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, 0) );
320#endif
321#ifdef HAVE_XPETRA_EPETRA
322 if(rowMap->lib() == UseEpetra)
323 return rcp( new EpetraCrsMatrixT<int,Node>(rowMap));
324#endif
325 XPETRA_FACTORY_END;
326 }
327
328 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
329 XPETRA_MONITOR("CrsMatrixFactory::Build");
330
331#ifdef HAVE_XPETRA_TPETRA
332 if (rowMap->lib() == UseTpetra)
333 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, maxNumEntriesPerRow, plist) );
334#endif
335
336 if (rowMap->lib() == UseEpetra)
337 return rcp( new EpetraCrsMatrixT<int,Node>(rowMap, maxNumEntriesPerRow, plist) );
338
339 XPETRA_FACTORY_END;
340 }
341
342 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
343 XPETRA_MONITOR("CrsMatrixFactory::Build");
344
345#ifdef HAVE_XPETRA_TPETRA
346 if (rowMap->lib() == UseTpetra)
347 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, NumEntriesPerRowToAlloc, plist) );
348#endif
349
350 if (rowMap->lib() == UseEpetra)
351 return rcp( new EpetraCrsMatrixT<int,Node>(rowMap, NumEntriesPerRowToAlloc, plist) );
352
353 XPETRA_FACTORY_END;
354 }
355
357 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
358 XPETRA_MONITOR("CrsMatrixFactory::Build");
359
360#ifdef HAVE_XPETRA_TPETRA
361 if (rowMap->lib() == UseTpetra)
362 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
363#endif
364
365 if (rowMap->lib() == UseEpetra)
366 return rcp( new EpetraCrsMatrixT<int,Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
367
368 XPETRA_FACTORY_END;
369 }
370
372 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
373 XPETRA_MONITOR("CrsMatrixFactory::Build");
374
375#ifdef HAVE_XPETRA_TPETRA
376 if (rowMap->lib() == UseTpetra)
377 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
378#endif
379
380 if (rowMap->lib() == UseEpetra)
381 return rcp( new EpetraCrsMatrixT<int,Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
382
383 XPETRA_FACTORY_END;
384 }
385
387 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
388 XPETRA_MONITOR("CrsMatrixFactory::Build");
389
390#ifdef HAVE_XPETRA_TPETRA
391 if (graph->getRowMap()->lib() == UseTpetra)
392 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(graph, plist) );
393#endif
394
395 if (graph->getRowMap()->lib() == UseEpetra)
396 return rcp( new EpetraCrsMatrixT<int,Node>(graph, plist) );
397
398 XPETRA_FACTORY_END;
399 }
400
401
403 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
404 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
405 const Import<LocalOrdinal,GlobalOrdinal,Node> &importer,
406 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap = Teuchos::null,
407 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap = Teuchos::null,
408 const Teuchos::RCP<Teuchos::ParameterList>& params = Teuchos::null) {
409 XPETRA_MONITOR("CrsMatrixFactory::Build");
410
411#ifdef HAVE_XPETRA_TPETRA
412 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
413 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,importer,domainMap,rangeMap,params) );
414#endif
415
416 if (sourceMatrix->getRowMap()->lib() == UseEpetra)
417 return rcp( new EpetraCrsMatrixT<int,Node>(sourceMatrix,importer,domainMap,rangeMap,params) );
418
419 XPETRA_FACTORY_END;
420 }
421
423 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
424 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
425 const Export<LocalOrdinal,GlobalOrdinal,Node> &exporter,
426 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap = Teuchos::null,
427 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap = Teuchos::null,
428 const Teuchos::RCP<Teuchos::ParameterList>& params = Teuchos::null) {
429 XPETRA_MONITOR("CrsMatrixFactory::Build");
430
431#ifdef HAVE_XPETRA_TPETRA
432 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
433 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,exporter,domainMap,rangeMap,params) );
434#endif
435
436 if (sourceMatrix->getRowMap()->lib() == UseEpetra)
437 return rcp( new EpetraCrsMatrixT<int,Node>(sourceMatrix,exporter,domainMap,rangeMap,params) );
438
439 XPETRA_FACTORY_END;
440 }
441
443 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
444 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
445 const Import<LocalOrdinal,GlobalOrdinal,Node> & RowImporter,
446 const RCP<const Import<LocalOrdinal,GlobalOrdinal,Node> > DomainImporter,
447 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap,
448 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap,
449 const Teuchos::RCP<Teuchos::ParameterList>& params) {
450 XPETRA_MONITOR("CrsMatrixFactory::Build");
451
452#ifdef HAVE_XPETRA_TPETRA
453 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
454 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,RowImporter,DomainImporter,domainMap,rangeMap,params) );
455#endif
456
457 if (sourceMatrix->getRowMap()->lib() == UseEpetra)
458 return rcp( new EpetraCrsMatrixT<int,Node>(sourceMatrix,RowImporter,DomainImporter,domainMap,rangeMap,params) );
459
460 XPETRA_FACTORY_END;
461 }
462
464 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
465 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
466 const Export<LocalOrdinal,GlobalOrdinal,Node> &RowExporter,
467 const RCP<const Export<LocalOrdinal,GlobalOrdinal,Node> > DomainExporter,
468 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap,
469 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap,
470 const Teuchos::RCP<Teuchos::ParameterList>& params) {
471 XPETRA_MONITOR("CrsMatrixFactory::Build");
472
473#ifdef HAVE_XPETRA_TPETRA
474 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
475 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,RowExporter,DomainExporter,domainMap,rangeMap,params) );
476#endif
477
478 if (sourceMatrix->getRowMap()->lib() == UseEpetra)
479 return rcp( new EpetraCrsMatrixT<int,Node>(sourceMatrix,RowExporter,DomainExporter,domainMap,rangeMap,params) );
480
481 XPETRA_FACTORY_END;
482 }
483
484
485#ifdef HAVE_XPETRA_KOKKOS_REFACTOR
486 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build (
487 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rowMap,
488 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& colMap,
489 const typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type& lclMatrix,
490 const Teuchos::RCP<Teuchos::ParameterList>& params = null) {
491 XPETRA_MONITOR("CrsMatrixFactory::Build");
492
493#ifdef HAVE_XPETRA_TPETRA
494 if (rowMap->lib() == UseTpetra)
495 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, lclMatrix, params));
496#endif
497
498 if (rowMap->lib() == UseEpetra)
499 return rcp( new EpetraCrsMatrixT<int,Node>(rowMap, colMap, lclMatrix, params) );
500
501 XPETRA_FACTORY_END;
502 }
503
504 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build (
505 const typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type& lclMatrix,
506 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rowMap,
507 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& colMap,
508 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& domainMap = Teuchos::null,
509 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rangeMap = Teuchos::null,
510 const Teuchos::RCP<Teuchos::ParameterList>& params = null) {
511 XPETRA_MONITOR("CrsMatrixFactory::Build");
512
513#ifdef HAVE_XPETRA_TPETRA
514 if (rowMap->lib() == UseTpetra)
515 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(lclMatrix, rowMap, colMap, domainMap, rangeMap, params));
516#endif
517
518 if (rowMap->lib() == UseEpetra)
519 return rcp( new EpetraCrsMatrixT<int,Node>(lclMatrix, rowMap, colMap, domainMap, rangeMap, params) );
520
521 XPETRA_FACTORY_END;
522 }
523
524 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build (
525 const typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type& lclMatrix,
526 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rowMap,
527 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& colMap,
528 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& domainMap,
529 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rangeMap,
530 const Teuchos::RCP<const Import<LocalOrdinal,GlobalOrdinal,Node>>& importer,
531 const Teuchos::RCP<const Export<LocalOrdinal,GlobalOrdinal,Node>>& exporter,
532 const Teuchos::RCP<Teuchos::ParameterList>& params = null) {
533 XPETRA_MONITOR("CrsMatrixFactory::Build");
534
535#ifdef HAVE_XPETRA_TPETRA
536 if (rowMap->lib() == UseTpetra)
537 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(lclMatrix, rowMap, colMap, domainMap, rangeMap, importer, exporter, params));
538#endif
539
540 TEUCHOS_TEST_FOR_EXCEPTION(rowMap->lib() == UseEpetra, std::logic_error, "Epetra doesn't support this matrix constructor");
541
542 XPETRA_FACTORY_END;
543 }
544#endif
545
546 };
547#endif
548
549// we need the Epetra specialization only if Epetra is enabled
550#if (defined(HAVE_XPETRA_EPETRA) && !defined(XPETRA_EPETRA_NO_64BIT_GLOBAL_INDICES))
551
552 template <>
553 class CrsMatrixFactory<double, int, long long, EpetraNode> {
554 typedef double Scalar;
555 typedef int LocalOrdinal;
556 typedef long long GlobalOrdinal;
557 typedef EpetraNode Node;
558
559 private:
561 CrsMatrixFactory() {}
562
563 public:
565 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
566 Build (const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap)
567 {
568 XPETRA_MONITOR("CrsMatrixFactory::Build");
569#ifdef HAVE_XPETRA_TPETRA
570 if (rowMap->lib() == UseTpetra)
571 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, 0) );
572#endif
573#ifdef HAVE_XPETRA_EPETRA
574 if(rowMap->lib() == UseEpetra)
575 return rcp( new EpetraCrsMatrixT<long long,Node>(rowMap,0));
576#endif
577 XPETRA_FACTORY_END;
578 }
579
580 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
581 XPETRA_MONITOR("CrsMatrixFactory::Build");
582
583#ifdef HAVE_XPETRA_TPETRA
584 if (rowMap->lib() == UseTpetra)
585 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, maxNumEntriesPerRow, plist) );
586#endif
587
588 if (rowMap->lib() == UseEpetra)
589 return rcp( new EpetraCrsMatrixT<long long, Node>(rowMap, maxNumEntriesPerRow, plist) );
590
591 XPETRA_FACTORY_END;
592 }
593
594 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
595 XPETRA_MONITOR("CrsMatrixFactory::Build");
596
597#ifdef HAVE_XPETRA_TPETRA
598 if (rowMap->lib() == UseTpetra)
599 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, NumEntriesPerRowToAlloc, plist) );
600#endif
601
602 if (rowMap->lib() == UseEpetra)
603 return rcp( new EpetraCrsMatrixT<long long, Node>(rowMap, NumEntriesPerRowToAlloc, plist) );
604
605 XPETRA_FACTORY_END;
606 }
607
609 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
610 XPETRA_MONITOR("CrsMatrixFactory::Build");
611
612#ifdef HAVE_XPETRA_TPETRA
613 if (rowMap->lib() == UseTpetra)
614 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
615#endif
616
617 if (rowMap->lib() == UseEpetra)
618 return rcp( new EpetraCrsMatrixT<long long, Node>(rowMap, colMap, maxNumEntriesPerRow, plist) );
619
620 XPETRA_FACTORY_END;
621 }
622
624 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
625 XPETRA_MONITOR("CrsMatrixFactory::Build");
626
627#ifdef HAVE_XPETRA_TPETRA
628 if (rowMap->lib() == UseTpetra)
629 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
630#endif
631
632 if (rowMap->lib() == UseEpetra)
633 return rcp( new EpetraCrsMatrixT<long long, Node>(rowMap, colMap, NumEntriesPerRowToAlloc, plist) );
634
635 XPETRA_FACTORY_END;
636 }
637
639 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null) {
640 XPETRA_MONITOR("CrsMatrixFactory::Build");
641
642#ifdef HAVE_XPETRA_TPETRA
643 if (graph->getRowMap()->lib() == UseTpetra)
644 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(graph, plist) );
645#endif
646
647 if (graph->getRowMap()->lib() == UseEpetra)
648 return rcp( new EpetraCrsMatrixT<long long, Node>(graph, plist) );
649
650 XPETRA_FACTORY_END;
651 }
652
653
655 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
656 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
657 const Import<LocalOrdinal,GlobalOrdinal,Node> &importer,
658 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap = Teuchos::null,
659 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap = Teuchos::null,
660 const Teuchos::RCP<Teuchos::ParameterList>& params = Teuchos::null) {
661 XPETRA_MONITOR("CrsMatrixFactory::Build");
662
663#ifdef HAVE_XPETRA_TPETRA
664 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
665 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,importer,domainMap,rangeMap,params) );
666#endif
667
668 if (sourceMatrix->getRowMap()->lib() == UseEpetra)
669 return rcp( new EpetraCrsMatrixT<long long, Node>(sourceMatrix,importer,domainMap,rangeMap,params) );
670
671 XPETRA_FACTORY_END;
672 }
673
675 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
676 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
677 const Export<LocalOrdinal,GlobalOrdinal,Node> &exporter,
678 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap = Teuchos::null,
679 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap = Teuchos::null,
680 const Teuchos::RCP<Teuchos::ParameterList>& params = Teuchos::null) {
681 XPETRA_MONITOR("CrsMatrixFactory::Build");
682
683#ifdef HAVE_XPETRA_TPETRA
684 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
685 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,exporter,domainMap,rangeMap,params) );
686#endif
687
688 if (sourceMatrix->getRowMap()->lib() == UseEpetra)
689 return rcp( new EpetraCrsMatrixT<long long, Node>(sourceMatrix,exporter,domainMap,rangeMap,params) );
690
691 XPETRA_FACTORY_END;
692 }
693
695 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
696 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
697 const Import<LocalOrdinal,GlobalOrdinal,Node> & RowImporter,
698 const RCP<const Import<LocalOrdinal,GlobalOrdinal,Node> > DomainImporter,
699 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap,
700 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap,
701 const Teuchos::RCP<Teuchos::ParameterList>& params) {
702 XPETRA_MONITOR("CrsMatrixFactory::Build");
703
704#ifdef HAVE_XPETRA_TPETRA
705 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
706 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,RowImporter,DomainImporter,domainMap,rangeMap,params) );
707#endif
708
709 if (sourceMatrix->getRowMap()->lib() == UseEpetra)
710 return rcp( new EpetraCrsMatrixT<long long,Node>(sourceMatrix,RowImporter,DomainImporter,domainMap,rangeMap,params) );
711
712 XPETRA_FACTORY_END;
713 }
714
716 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build(
717 const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix,
718 const Export<LocalOrdinal,GlobalOrdinal,Node> &RowExporter,
719 const RCP<const Export<LocalOrdinal,GlobalOrdinal,Node> > DomainExporter,
720 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & domainMap,
721 const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > & rangeMap,
722 const Teuchos::RCP<Teuchos::ParameterList>& params) {
723 XPETRA_MONITOR("CrsMatrixFactory::Build");
724
725#ifdef HAVE_XPETRA_TPETRA
726 if (sourceMatrix->getRowMap()->lib() == UseTpetra)
727 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(sourceMatrix,RowExporter,DomainExporter,domainMap,rangeMap,params) );
728#endif
729
730 if (sourceMatrix->getRowMap()->lib() == UseEpetra)
731 return rcp( new EpetraCrsMatrixT<long long,Node>(sourceMatrix,RowExporter,DomainExporter,domainMap,rangeMap,params) );
732
733 XPETRA_FACTORY_END;
734 }
735
736#ifdef HAVE_XPETRA_KOKKOS_REFACTOR
737 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build (
738 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rowMap,
739 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& colMap,
740 const typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type& lclMatrix,
741 const Teuchos::RCP<Teuchos::ParameterList>& params = null) {
742 XPETRA_MONITOR("CrsMatrixFactory::Build");
743
744#ifdef HAVE_XPETRA_TPETRA
745 if (rowMap->lib() == UseTpetra)
746 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(rowMap, colMap, lclMatrix, params));
747#endif
748
749 if (rowMap->lib() == UseEpetra)
750 return rcp( new EpetraCrsMatrixT<long long, Node>(rowMap, colMap, lclMatrix, params) );
751
752 XPETRA_FACTORY_END;
753 }
754 static RCP<CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > Build (
755 const typename Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>::local_matrix_type& lclMatrix,
756 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rowMap,
757 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& colMap,
758 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& domainMap = Teuchos::null,
759 const Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& rangeMap = Teuchos::null,
760 const Teuchos::RCP<Teuchos::ParameterList>& params = null) {
761 XPETRA_MONITOR("CrsMatrixFactory::Build");
762
763#ifdef HAVE_XPETRA_TPETRA
764 if (rowMap->lib() == UseTpetra)
765 return rcp( new TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(lclMatrix, rowMap, colMap, domainMap, rangeMap, params));
766#endif
767
768 if (rowMap->lib() == UseEpetra)
769 return rcp( new EpetraCrsMatrixT<long long, Node>(lclMatrix, rowMap, colMap, domainMap, rangeMap, params) );
770
771 XPETRA_FACTORY_END;
772 }
773#endif
774
775 };
776#endif
777
778}
779
780#define XPETRA_CRSMATRIXFACTORY_SHORT
781#endif
Xpetra_ConfigDefs.hpp
XPETRA_MONITOR
#define XPETRA_MONITOR(funcName)
Definition: Xpetra_ConfigDefs.hpp:132
Xpetra_CrsMatrix.hpp
Xpetra_EpetraCrsMatrix.hpp
Xpetra_Exceptions.hpp
XPETRA_FACTORY_ERROR_IF_EPETRA
#define XPETRA_FACTORY_ERROR_IF_EPETRA(lib)
Definition: Xpetra_Exceptions.hpp:70
XPETRA_FACTORY_END
#define XPETRA_FACTORY_END
Definition: Xpetra_Exceptions.hpp:77
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap=Teuchos::null, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor using FusedImport.
Definition: Xpetra_CrsMatrixFactory.hpp:403
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap)
Constructor for empty matrix (intended use is an import/export target - can't insert entries directly...
Definition: Xpetra_CrsMatrixFactory.hpp:314
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying column Map and fixed number of entries for each row.
Definition: Xpetra_CrsMatrixFactory.hpp:357
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Definition: Xpetra_CrsMatrixFactory.hpp:342
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::CrsMatrixFactory
CrsMatrixFactory()
Private constructor. This is a static class.
Definition: Xpetra_CrsMatrixFactory.hpp:309
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &RowImporter, const RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > DomainImporter, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Constructor using FusedImport.
Definition: Xpetra_CrsMatrixFactory.hpp:443
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying a previously constructed graph.
Definition: Xpetra_CrsMatrixFactory.hpp:387
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &RowExporter, const RCP< const Export< LocalOrdinal, GlobalOrdinal, Node > > DomainExporter, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Constructor using FusedExport.
Definition: Xpetra_CrsMatrixFactory.hpp:464
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying column Map and number of entries in each row.
Definition: Xpetra_CrsMatrixFactory.hpp:372
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap=Teuchos::null, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor using FusedExport.
Definition: Xpetra_CrsMatrixFactory.hpp:423
Xpetra::CrsMatrixFactory< double, int, int, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Definition: Xpetra_CrsMatrixFactory.hpp:328
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying column Map and fixed number of entries for each row.
Definition: Xpetra_CrsMatrixFactory.hpp:609
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap)
Constructor for empty matrix (intended use is an import/export target - can't insert entries directly...
Definition: Xpetra_CrsMatrixFactory.hpp:566
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Definition: Xpetra_CrsMatrixFactory.hpp:580
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &RowExporter, const RCP< const Export< LocalOrdinal, GlobalOrdinal, Node > > DomainExporter, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Constructor using FusedExport.
Definition: Xpetra_CrsMatrixFactory.hpp:716
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap=Teuchos::null, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor using FusedExport.
Definition: Xpetra_CrsMatrixFactory.hpp:675
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying a previously constructed graph.
Definition: Xpetra_CrsMatrixFactory.hpp:639
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::CrsMatrixFactory
CrsMatrixFactory()
Private constructor. This is a static class.
Definition: Xpetra_CrsMatrixFactory.hpp:561
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying column Map and number of entries in each row.
Definition: Xpetra_CrsMatrixFactory.hpp:624
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &RowImporter, const RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > DomainImporter, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Constructor using FusedImport.
Definition: Xpetra_CrsMatrixFactory.hpp:695
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Definition: Xpetra_CrsMatrixFactory.hpp:594
Xpetra::CrsMatrixFactory< double, int, long long, EpetraNode >::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap=Teuchos::null, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Constructor using FusedImport.
Definition: Xpetra_CrsMatrixFactory.hpp:655
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap=Teuchos::null, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Definition: Xpetra_CrsMatrixFactory.hpp:166
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap)
Constructor for empty matrix (intended use is an import/export target - can't insert entries directly...
Definition: Xpetra_CrsMatrixFactory.hpp:77
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying (possibly different) number of entries in each row.
Definition: Xpetra_CrsMatrixFactory.hpp:108
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &RowExporter, const RCP< const Export< LocalOrdinal, GlobalOrdinal, Node > > DomainExporter, const RCP< Map< LocalOrdinal, GlobalOrdinal, Scalar > > &domainMap, const RCP< Map< LocalOrdinal, GlobalOrdinal, Scalar > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Definition: Xpetra_CrsMatrixFactory.hpp:218
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, const RCP< Map< LocalOrdinal, GlobalOrdinal, Scalar > > &domainMap=Teuchos::null, const RCP< Map< LocalOrdinal, GlobalOrdinal, Scalar > > &rangeMap=Teuchos::null, const Teuchos::RCP< Teuchos::ParameterList > &params=Teuchos::null)
Definition: Xpetra_CrsMatrixFactory.hpp:183
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &sourceMatrix, const Import< LocalOrdinal, GlobalOrdinal, Node > &RowImporter, const RCP< const Import< LocalOrdinal, GlobalOrdinal, Node > > DomainImporter, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &domainMap, const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rangeMap, const Teuchos::RCP< Teuchos::ParameterList > &params)
Definition: Xpetra_CrsMatrixFactory.hpp:200
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying column Map and fixed number of entries for each row.
Definition: Xpetra_CrsMatrixFactory.hpp:123
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, size_t maxNumEntriesPerRow, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying fixed number of entries for each row.
Definition: Xpetra_CrsMatrixFactory.hpp:91
Xpetra::CrsMatrixFactory::CrsMatrixFactory
CrsMatrixFactory()
Private constructor. This is a static class.
Definition: Xpetra_CrsMatrixFactory.hpp:72
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &rowMap, const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &colMap, const ArrayRCP< const size_t > &NumEntriesPerRowToAlloc, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying column Map and number of entries in each row.
Definition: Xpetra_CrsMatrixFactory.hpp:140
Xpetra::CrsMatrixFactory::Build
static RCP< CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph, const Teuchos::RCP< Teuchos::ParameterList > &plist=Teuchos::null)
Constructor specifying a previously constructed graph.
Definition: Xpetra_CrsMatrixFactory.hpp:153
TEUCHOS_TEST_FOR_EXCEPTION
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
rcp
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Xpetra
Xpetra namespace
Definition: Xpetra_BlockedCrsMatrix.hpp:89