Xpetra_ThyraUtils.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 // Tobias Wiesner (tawiesn@sandia.gov)
43 //
44 // ***********************************************************************
45 //
46 // @HEADER
47 
48 #include "Xpetra_ConfigDefs.hpp"
49 #ifdef HAVE_XPETRA_THYRA
50 
52 
53 #include "Xpetra_ThyraUtils.hpp"
54 
55 
56 namespace Xpetra {
57 
58 #ifdef HAVE_XPETRA_EPETRA
59  // implementation of "toThyra" for full specialization on SC=double, LO=GO=int and NO=EpetraNode
60  // We need the specialization in the cpp file due to a circle dependency in the .hpp files for BlockedCrsMatrix
62  ThyraUtils<double, int, int, EpetraNode>::toThyra(const Teuchos::RCP<Xpetra::BlockedCrsMatrix<double, int, int, EpetraNode> >& mat) {
63 
64  int nRows = mat->Rows();
65  int nCols = mat->Cols();
66 
68 
69  bool bTpetra = false;
70  bool bEpetra = false;
71 #ifdef HAVE_XPETRA_TPETRA
72  // Note: Epetra is enabled
73 #if ((defined(EPETRA_HAVE_OMP) && defined(HAVE_TPETRA_INST_OPENMP) && defined(HAVE_TPETRA_INST_INT_INT) && defined(HAVE_TPETRA_INST_DOUBLE)) || \
74  (!defined(EPETRA_HAVE_OMP) && defined(HAVE_TPETRA_INST_SERIAL) && defined(HAVE_TPETRA_INST_INT_INT) && defined(HAVE_TPETRA_INST_DOUBLE)))
76  if(tpetraMat!=Teuchos::null) bTpetra = true;
77 #else
78  bTpetra = false;
79 #endif
80 #endif
81 
82 #ifdef HAVE_XPETRA_EPETRA
84  if(epetraMat!=Teuchos::null) bEpetra = true;
85 #endif
86 
87  TEUCHOS_TEST_FOR_EXCEPT(bTpetra == bEpetra); // we only allow Epetra OR Tpetra
88 
89  // create new Thyra blocked operator
91  Thyra::defaultBlockedLinearOp<Scalar>();
92 
93  blockMat->beginBlockFill(nRows,nCols);
94 
95  for (int r=0; r<nRows; ++r) {
96  for (int c=0; c<nCols; ++c) {
98  Xpetra::ThyraUtils<Scalar,LocalOrdinal,GlobalOrdinal,Node>::toThyra(mat->getMatrix(r,c));
99  std::stringstream label; label << "A" << r << c;
100  blockMat->setBlock(r,c,thBlock);
101  }
102  }
103 
104  blockMat->endBlockFill();
105 
106  return blockMat;
107  }
108 #endif
109 
110 } // namespace Xpetra
111 
112 #endif // HAVE_XPETRA_THYRA
Xpetra namespace
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)