Thyra Package Browser (Single Doxygen Collection) Version of the Day
Loading...
Searching...
No Matches
EpetraLinearOpParallel_UnitTests.cpp
Go to the documentation of this file.
1/*
2// @HEADER
3// ***********************************************************************
4//
5// Thyra: Interfaces and Support for Abstract Numerical Algorithms
6// Copyright (2004) Sandia Corporation
7//
8// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
9// license for use of this work by or on behalf of the U.S. Government.
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 Roscoe A. Bartlett (bartlettra@ornl.gov)
39//
40// ***********************************************************************
41// @HEADER
42*/
43
45#include "Thyra_LinearOpTester.hpp"
46#include "Thyra_DefaultBlockedLinearOp.hpp"
47#include "Thyra_MultiVectorBase.hpp"
48#include "Thyra_MultiVectorStdOps.hpp"
50
52
53
54namespace Thyra {
55
56
57//
58// Unit Tests
59//
60
61
63{
64
65 using Teuchos::describe;
66
67 const int scale=30;
68
69 // build sub operators
70 RCP<const LinearOpBase<double> > A00 = epetraLinearOp(getEpetraMatrix(scale*4,scale*4,0));
71 RCP<const LinearOpBase<double> > A01 = epetraLinearOp(getEpetraMatrix(scale*4,scale*3,1));
72 RCP<const LinearOpBase<double> > A02 = epetraLinearOp(getEpetraMatrix(scale*4,scale*2,2));
73 RCP<const LinearOpBase<double> > A10 = epetraLinearOp(getEpetraMatrix(scale*3,scale*4,3));
74 RCP<const LinearOpBase<double> > A11 = epetraLinearOp(getEpetraMatrix(scale*3,scale*3,4));
75 RCP<const LinearOpBase<double> > A12 = epetraLinearOp(getEpetraMatrix(scale*3,scale*2,5));
76 RCP<const LinearOpBase<double> > A20 = epetraLinearOp(getEpetraMatrix(scale*2,scale*4,6));
77 RCP<const LinearOpBase<double> > A21 = epetraLinearOp(getEpetraMatrix(scale*2,scale*3,8));
78 RCP<const LinearOpBase<double> > A22 = epetraLinearOp(getEpetraMatrix(scale*2,scale*2,8));
79
80 out << "Sub operators built" << std::endl;
81
82 typedef double S;
84 M00->beginBlockFill(2,2);
85 M00->setBlock(0,0,A00);
86 M00->setBlock(0,1,A02);
87 M00->setBlock(1,0,A20);
88 M00->setBlock(1,1,A22);
89 M00->endBlockFill();
90
91 out << "Built M00" << std::endl;
92
94 M10->beginBlockFill(2,1);
95 M10->setBlock(0,0,A01);
96 M10->setBlock(1,0,A21);
97 M10->endBlockFill();
98
99 out << "Built M10" << std::endl;
100
102 M01->beginBlockFill(1,2);
103 M01->setBlock(0,0,A10);
104 M01->setBlock(0,1,A12);
105 M01->endBlockFill();
106
107 out << "Built M01" << std::endl;
108
110 M->beginBlockFill(2,2);
111 M->setBlock(0,0,A11);
112 M->setBlock(0,1,M01);
113 M->setBlock(1,0,M10);
114 M->setBlock(1,1,M00);
115 M->endBlockFill();
116
117 out << "Built M" << std::endl;
118
119 out << "M = " << describe(*M, Teuchos::VERB_MEDIUM);
120 out << "M->range() = " << describe(*M->range(), Teuchos::VERB_MEDIUM);
121 out << "M->domain() = " << describe(*M->range(), Teuchos::VERB_MEDIUM);
122
123 out << "Test complete" << std::endl;
124
125}
126
127
128} // namespace Thyra
TEUCHOS_UNIT_TEST(Comm, Issue1029)