Xpetra Version of the Day
Loading...
Searching...
No Matches
Xpetra_MatrixFactory.cpp
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
47#ifndef XPETRA_MATRIXFACTORY_CPP
48#define XPETRA_MATRIXFACTORY_CPP
49
51#if 0
52namespace Xpetra {
53 template<>
54 Teuchos::RCP<Matrix > MatrixFactory2<double,int,int,Xpetra::Matrix<double, int, int>::node_type>::BuildCopy(const Teuchos::RCP<const Matrix > A) {
55 RCP<const CrsMatrixWrap> oldOp = Teuchos::rcp_dynamic_cast<const CrsMatrixWrap>(A);
56 if (oldOp == Teuchos::null)
57 throw Exceptions::BadCast("Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
58
59 RCP<const CrsMatrix> oldCrsOp = oldOp->getCrsMatrix();
60
61#ifdef HAVE_XPETRA_EPETRA
62#ifndef XPETRA_EPETRA_NO_32BIT_GLOBAL_INDICES
63 RCP<const EpetraCrsMatrixT<GlobalOrdinal,Node> > oldECrsOp = Teuchos::rcp_dynamic_cast<const EpetraCrsMatrixT<GlobalOrdinal,Node> >(oldCrsOp);
64 if (oldECrsOp != Teuchos::null) {
65 // Underlying matrix is Epetra
66 RCP<CrsMatrix> newECrsOp(new EpetraCrsMatrixT<GlobalOrdinal,Node>(*oldECrsOp));
67 RCP<CrsMatrixWrap> newOp (new CrsMatrixWrap (newECrsOp));
68
69 return newOp;
70 }
71#endif
72#endif
73
74#ifdef HAVE_XPETRA_TPETRA
75 // Underlying matrix is Tpetra
76 RCP<const TpetraCrsMatrix> oldTCrsOp = Teuchos::rcp_dynamic_cast<const TpetraCrsMatrix>(oldCrsOp);
77 if (oldTCrsOp != Teuchos::null) {
78 RCP<CrsMatrix> newTCrsOp(new TpetraCrsMatrix(*oldTCrsOp));
79 RCP<CrsMatrixWrap> newOp (new CrsMatrixWrap(newTCrsOp));
80
81 return newOp;
82 }
83#else
84 throw Exceptions::BadCast("Cast from Xpetra::Matrix to Xpetra::EpetraCrsMatrix or Xpetra::TpetraCrsMatrix failed");
85#endif
86
87 return Teuchos::null; // make compiler happy
88 }
89
90#ifdef HAVE_XPETRA_INT_LONG_LONG
91 template<>
92 Teuchos::RCP<Xpetra::Matrix<double,int,long long,typename Xpetra::Matrix<double, int, long long>::node_type> > MatrixFactory2<double,int,long long,typename Xpetra::Matrix<double, int, long long>::node_type>::BuildCopy(const Teuchos::RCP<const Xpetra::Matrix<double,int,long long,typename Matrix<double, int, long long>::node_type> > A) {
93 RCP<const CrsMatrixWrap> oldOp = Teuchos::rcp_dynamic_cast<const CrsMatrixWrap>(A);
94 if (oldOp == Teuchos::null)
95 throw Exceptions::BadCast("Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
96
97 RCP<const CrsMatrix> oldCrsOp = oldOp->getCrsMatrix();
98
99#ifdef HAVE_XPETRA_EPETRA
100#ifndef XPETRA_EPETRA_NO_64BIT_GLOBAL_INDICES
101 RCP<const EpetraCrsMatrixT<GlobalOrdinal,Node> > oldECrsOp = Teuchos::rcp_dynamic_cast<const EpetraCrsMatrixT<GlobalOrdinal,Node> >(oldCrsOp);
102 if (oldECrsOp != Teuchos::null) {
103 // Underlying matrix is Epetra
104 RCP<CrsMatrix> newECrsOp(new EpetraCrsMatrixT<GlobalOrdinal,Node>(*oldECrsOp));
105 RCP<CrsMatrixWrap> newOp (new CrsMatrixWrap (newECrsOp));
106
107 return newOp;
108 }
109#endif
110#endif
111
112#ifdef HAVE_XPETRA_TPETRA
113 // Underlying matrix is Tpetra
114 RCP<const TpetraCrsMatrix> oldTCrsOp = Teuchos::rcp_dynamic_cast<const TpetraCrsMatrix>(oldCrsOp);
115 if (oldTCrsOp != Teuchos::null) {
116 RCP<CrsMatrix> newTCrsOp(new TpetraCrsMatrix(*oldTCrsOp));
117 RCP<CrsMatrixWrap> newOp (new CrsMatrixWrap(newTCrsOp));
118
119 return newOp;
120 }
121#else
122 throw Exceptions::BadCast("Cast from Xpetra::Matrix to Xpetra::EpetraCrsMatrix or Xpetra::TpetraCrsMatrix failed");
123#endif
124
125 return Teuchos::null; // make compiler happy
126 }
127
128#endif // HAVE_XPETRA_INT_LONG_LONG
129
130} // namespace Xpetra
131#endif // if 0
132#endif
Xpetra-specific matrix class.
Xpetra namespace