Thyra Package Browser (Single Doxygen Collection) Version of the Day
EpetraThyraWrappers_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_EpetraThyraWrappers.hpp"
00045 #include "Thyra_DefaultProductVectorSpace.hpp"
00046 #include "Thyra_VectorStdOps.hpp"
00047 #include "Thyra_MultiVectorStdOps.hpp"
00048 #include "Thyra_TestingTools.hpp"
00049 #include "EpetraThyraAdaptersTestHelpers.hpp"
00050 
00051 #include "Teuchos_UnitTestHarness.hpp"
00052 
00053 
00054 namespace {
00055 
00056 
00057 //
00058 // Unit Tests
00059 //
00060 
00061 
00062 TEUCHOS_UNIT_TEST( get_Epetra_MultiVector, singleBlockProductVector )
00063 {
00064   using Teuchos::Comm;
00065   typedef Teuchos_Ordinal Ordinal;
00066   using Thyra::VectorSpaceBase;
00067   using Thyra::MultiVectorBase;
00068 
00069   const RCP<const Epetra_Comm> epetra_comm = getEpetraComm();
00070   const RCP<const Comm<Ordinal> > comm = Thyra::create_Comm(epetra_comm);
00071   
00072   const RCP<const Epetra_Map> epetra_map = rcp(new Epetra_Map(g_localDim, 0, *epetra_comm));
00073   const RCP<const VectorSpaceBase<double> > vs =  Thyra::create_VectorSpace(epetra_map);
00074 
00075   const RCP<const VectorSpaceBase<double> > pvs = Thyra::productVectorSpace(vs, 1);
00076 
00077   const RCP<MultiVectorBase<double> > pmv = Thyra::createMembers(pvs, 1);
00078 
00079   const double alpha = 3.5;
00080   Thyra::assign<double>( pmv.ptr(), alpha );
00081 
00082   const RCP<Epetra_MultiVector> epetra_mv =
00083     Thyra::get_Epetra_MultiVector(*epetra_map, pmv);
00084 
00085   const RCP<MultiVectorBase<double> > mv2 =
00086     Thyra::create_MultiVector(epetra_mv, pvs);
00087 
00088   Thyra::testRelNormDiffErr<double>(
00089     "*pmv->col(0)", *pmv->col(0),
00090     "*mv2->col(0)", *mv2->col(0),
00091     "max-error", 0.0,
00092     "max-warning", 0.0,
00093     &out
00094     );
00095    
00096 }
00097 
00098 
00099 } // namespace
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines