Thyra Package Browser (Single Doxygen Collection) Version of the Day
EpetraLinearOpParallel_UnitTests.cpp
Go to the documentation of this file.
00001 /*
00002 // @HEADER
00003 // ***********************************************************************
00004 // 
00005 //    Thyra: Interfaces and Support for Abstract Numerical Algorithms
00006 //                 Copyright (2004) Sandia Corporation
00007 // 
00008 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00009 // license for use of this work by or on behalf of the U.S. Government.
00010 // 
00011 // Redistribution and use in source and binary forms, with or without
00012 // modification, are permitted provided that the following conditions are
00013 // met:
00014 //
00015 // 1. Redistributions of source code must retain the above copyright
00016 // notice, this list of conditions and the following disclaimer.
00017 //
00018 // 2. Redistributions in binary form must reproduce the above copyright
00019 // notice, this list of conditions and the following disclaimer in the
00020 // documentation and/or other materials provided with the distribution.
00021 //
00022 // 3. Neither the name of the Corporation nor the names of the
00023 // contributors may be used to endorse or promote products derived from
00024 // this software without specific prior written permission.
00025 //
00026 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
00027 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
00028 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
00029 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
00030 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
00031 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
00032 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
00033 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
00034 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
00035 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
00036 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00037 //
00038 // Questions? Contact Roscoe A. Bartlett (bartlettra@ornl.gov) 
00039 // 
00040 // ***********************************************************************
00041 // @HEADER
00042 */
00043 
00044 #include "Thyra_EpetraLinearOp.hpp"
00045 #include "Thyra_LinearOpTester.hpp"
00046 #include "Thyra_DefaultBlockedLinearOp.hpp"
00047 #include "Thyra_MultiVectorBase.hpp"
00048 #include "Thyra_MultiVectorStdOps.hpp"
00049 #include "EpetraThyraAdaptersTestHelpers.hpp"
00050 
00051 #include "Teuchos_UnitTestHarness.hpp"
00052 
00053 
00054 namespace Thyra {
00055 
00056 
00057 //
00058 // Unit Tests
00059 //
00060 
00061 
00062 TEUCHOS_UNIT_TEST( EpetraLinearOpParallel, blocked_op )
00063 {
00064 
00065   using Teuchos::describe;
00066 
00067   const int scale=30;
00068   
00069   // build sub operators
00070   RCP<const LinearOpBase<double> > A00 = epetraLinearOp(getEpetraMatrix(scale*4,scale*4,0));
00071   RCP<const LinearOpBase<double> > A01 = epetraLinearOp(getEpetraMatrix(scale*4,scale*3,1));
00072   RCP<const LinearOpBase<double> > A02 = epetraLinearOp(getEpetraMatrix(scale*4,scale*2,2));
00073   RCP<const LinearOpBase<double> > A10 = epetraLinearOp(getEpetraMatrix(scale*3,scale*4,3));
00074   RCP<const LinearOpBase<double> > A11 = epetraLinearOp(getEpetraMatrix(scale*3,scale*3,4));
00075   RCP<const LinearOpBase<double> > A12 = epetraLinearOp(getEpetraMatrix(scale*3,scale*2,5));
00076   RCP<const LinearOpBase<double> > A20 = epetraLinearOp(getEpetraMatrix(scale*2,scale*4,6));
00077   RCP<const LinearOpBase<double> > A21 = epetraLinearOp(getEpetraMatrix(scale*2,scale*3,8));
00078   RCP<const LinearOpBase<double> > A22 = epetraLinearOp(getEpetraMatrix(scale*2,scale*2,8));
00079   
00080   out << "Sub operators built" << std::endl;
00081 
00082   typedef double S;
00083   const RCP<PhysicallyBlockedLinearOpBase<S> > M00 = defaultBlockedLinearOp<S>();
00084   M00->beginBlockFill(2,2);
00085   M00->setBlock(0,0,A00);
00086   M00->setBlock(0,1,A02);
00087   M00->setBlock(1,0,A20);
00088   M00->setBlock(1,1,A22);
00089   M00->endBlockFill();
00090 
00091   out << "Built M00" << std::endl;
00092 
00093   const RCP<PhysicallyBlockedLinearOpBase<S> > M10 = defaultBlockedLinearOp<S>();
00094   M10->beginBlockFill(2,1);
00095   M10->setBlock(0,0,A01);
00096   M10->setBlock(1,0,A21);
00097   M10->endBlockFill();
00098 
00099   out << "Built M10" << std::endl;
00100 
00101   const RCP<PhysicallyBlockedLinearOpBase<S> > M01 = defaultBlockedLinearOp<S>();
00102   M01->beginBlockFill(1,2);
00103   M01->setBlock(0,0,A10);
00104   M01->setBlock(0,1,A12);
00105   M01->endBlockFill();
00106 
00107   out << "Built M01" << std::endl;
00108 
00109   const RCP<PhysicallyBlockedLinearOpBase<S> > M = defaultBlockedLinearOp<S>();
00110   M->beginBlockFill(2,2);
00111   M->setBlock(0,0,A11);
00112   M->setBlock(0,1,M01);
00113   M->setBlock(1,0,M10);
00114   M->setBlock(1,1,M00);
00115   M->endBlockFill();
00116 
00117   out << "Built M" << std::endl;
00118 
00119   out << "M = " << describe(*M, Teuchos::VERB_MEDIUM);
00120   out << "M->range() = " << describe(*M->range(), Teuchos::VERB_MEDIUM);
00121   out << "M->domain() = " << describe(*M->range(), Teuchos::VERB_MEDIUM);
00122 
00123   out << "Test complete" << std::endl;
00124 
00125 }
00126 
00127 
00128 } // namespace Thyra
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines