Belos Package Browser (Single Doxygen Collection)  Development
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
test_gcrodr_complex.cpp
Go to the documentation of this file.
1 //@HEADER
2 // ************************************************************************
3 //
4 // Belos: Block Linear Solvers Package
5 // Copyright 2004 Sandia Corporation
6 //
7 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
8 // the U.S. Government retains certain rights in this software.
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 #include <BelosConfigDefs.hpp>
45 #include <BelosLinearProblem.hpp>
46 #include <BelosGCRODRSolMgr.hpp>
47 #include <complex>
48 
49 #include "MyMultiVec.hpp"
50 #include "MyBetterOperator.hpp"
51 #include "MyOperator.hpp"
52 
53 //
54 // mfh 20 Jan 2014: This test ensures the following:
55 //
56 // 1. Belos::GCRODRSolMgr can compile whether its ScalarType (first)
57 // template parameter is real or complex.
58 // 2. Belos::GCRODRSolMgr's constructor throws std::logic_error if and
59 // only if its ScalarType (first) template parameter is complex.
60 //
61 // At some point, if GCRODRSolMgr gets fixed so that it works with
62 // complex ScalarType, the second test will no longer pass. This will
63 // be a good thing! The test should still be built in that case, in
64 // order to demonstrate that GCRODRSolMgr compiles for complex
65 // ScalarType. However, in that case, the TEST_THROW macro should be
66 // changed to TEST_NOTHROW, and the macro's second argument should be
67 // removed.
68 //
69 // This test requires that Trilinos was compiled with complex
70 // arithmetic support enabled.
71 //
72 // Updates
73 // 4/15/2014: GCRODR is changed from RealSolverManager to SolverManager,
74 // and the test is changed from THROW to NOTHROW. (phtsuji)
75 
76 TEUCHOS_UNIT_TEST( GCRODR, RealDoesNotThrow )
77 {
78  using Teuchos::RCP;
79  using Teuchos::rcp;
80  typedef double ST;
81  typedef Belos::MultiVec<ST> MV;
82  typedef Belos::Operator<ST> OP;
83  //typedef Belos::MultiVecTraits<ST, MV> MVT;
84  //typedef Belos::OperatorTraits<ST, MV, OP> OPT;
85  typedef Belos::GCRODRSolMgr<ST, MV, OP> sol_mgr_type;
86 
87  RCP<sol_mgr_type> solver;
88  TEST_NOTHROW( solver = rcp (new sol_mgr_type ()) );
89 }
90 
91 TEUCHOS_UNIT_TEST( GCRODR, ComplexDoesNotThrow )
92 {
93  using Teuchos::RCP;
94  using Teuchos::rcp;
95  typedef std::complex<double> ST;
96  typedef Belos::MultiVec<ST> MV;
97  typedef Belos::Operator<ST> OP;
98  // typedef Belos::MultiVecTraits<ST, MV> MVT;
99  // typedef Belos::OperatorTraits<ST, MV, OP> OPT;
100  typedef Belos::GCRODRSolMgr<ST, MV, OP> sol_mgr_type;
101 
102  RCP<sol_mgr_type> solver;
103  TEST_NOTHROW( solver = rcp (new sol_mgr_type ()) );
104 }
#define TEST_NOTHROW(code)
TEUCHOS_UNIT_TEST(Factory, Bug6383)
Definition: Factory.cpp:233
Declaration and definition of Belos::GCRODRSolMgr, which implements the GCRODR (recycling GMRES) solv...
Alternative run-time polymorphic interface for operators.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Class which describes the linear problem to be solved by the iterative solver.
Implementation of the GCRODR (Recycling GMRES) iterative linear solver.
Interface for multivectors used by Belos&#39; linear solvers.
Belos header file which uses auto-configuration information to include necessary C++ headers...