Thyra_ScalarProdVectorSpaceBase.hpp

Go to the documentation of this file.
00001 // @HEADER
00002 // ***********************************************************************
00003 // 
00004 //    Thyra: Interfaces and Support for Abstract Numerical Algorithms
00005 //                 Copyright (2004) Sandia Corporation
00006 // 
00007 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00008 // license for use of this work by or on behalf of the U.S. Government.
00009 // 
00010 // This library is free software; you can redistribute it and/or modify
00011 // it under the terms of the GNU Lesser General Public License as
00012 // published by the Free Software Foundation; either version 2.1 of the
00013 // License, or (at your option) any later version.
00014 //  
00015 // This library is distributed in the hope that it will be useful, but
00016 // WITHOUT ANY WARRANTY; without even the implied warranty of
00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00018 // Lesser General Public License for more details.
00019 //  
00020 // You should have received a copy of the GNU Lesser General Public
00021 // License along with this library; if not, write to the Free Software
00022 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
00023 // USA
00024 // Questions? Contact Michael A. Heroux (maherou@sandia.gov) 
00025 // 
00026 // ***********************************************************************
00027 // @HEADER
00028 
00029 #ifndef THYRA_VECTOR_SPACE_STD_BASE_HPP
00030 #define THYRA_VECTOR_SPACE_STD_BASE_HPP
00031 
00032 #include "Thyra_ScalarProdVectorSpaceBaseDecl.hpp"
00033 #include "Thyra_VectorSpaceDefaultBase.hpp"
00034 #include "Thyra_EuclideanScalarProd.hpp"
00035 #include "Thyra_AssertOp.hpp"
00036 
00037 namespace Thyra {
00038 
00039 // Constructors / initializers
00040 
00041 template<class Scalar>
00042 ScalarProdVectorSpaceBase<Scalar>::ScalarProdVectorSpaceBase()
00043   :scalarProd_(Teuchos::rcp(new EuclideanScalarProd<Scalar>()))
00044 {}
00045   
00046 template<class Scalar>
00047 ScalarProdVectorSpaceBase<Scalar>::ScalarProdVectorSpaceBase( const Teuchos::RefCountPtr<const ScalarProdBase<Scalar> > &scalarProd )
00048   :scalarProd_(scalarProd)
00049 {
00050   TEST_FOR_EXCEPT( scalarProd.get()==NULL );
00051 }
00052 
00053 template<class Scalar>
00054 void ScalarProdVectorSpaceBase<Scalar>::setScalarProd( const Teuchos::RefCountPtr<const ScalarProdBase<Scalar> > &scalarProd )
00055 {
00056   TEST_FOR_EXCEPT( scalarProd.get()==NULL );
00057   scalarProd_ = scalarProd;
00058 }
00059 
00060 template<class Scalar>
00061 Teuchos::RefCountPtr<const ScalarProdBase<Scalar> >
00062 ScalarProdVectorSpaceBase<Scalar>::getScalarProd() const
00063 {
00064   return scalarProd_;
00065 }
00066 
00067 // Overridden from VectorSpaceBase
00068 
00069 template<class Scalar>
00070 bool ScalarProdVectorSpaceBase<Scalar>::isEuclidean() const
00071 {
00072   return scalarProd_->isEuclidean();
00073 }
00074 
00075 template<class Scalar>
00076 Scalar ScalarProdVectorSpaceBase<Scalar>::scalarProd( const VectorBase<Scalar>& x, const VectorBase<Scalar>& y ) const
00077 {
00078 #ifdef TEUCHOS_DEBUG
00079   THYRA_ASSERT_VEC_SPACES("ScalarProdVectorSpaceBase<Scalar>::scalarProd(...)",*x.space(),*this);
00080   THYRA_ASSERT_VEC_SPACES("ScalarProdVectorSpaceBase<Scalar>::scalarProd(...)",*y.space(),*this);
00081 #endif
00082   return scalarProd_->scalarProd(x,y);
00083 }
00084 
00085 template<class Scalar>
00086 void ScalarProdVectorSpaceBase<Scalar>::scalarProds( const MultiVectorBase<Scalar>& X, const MultiVectorBase<Scalar>& Y, Scalar scalar_prods[] ) const
00087 {
00088 #ifdef TEUCHOS_DEBUG
00089   THYRA_ASSERT_VEC_SPACES("ScalarProdVectorSpaceBase<Scalar>::scalarProds(...)",*X.range(),*this);
00090   THYRA_ASSERT_VEC_SPACES("ScalarProdVectorSpaceBase<Scalar>::scalarProds(...)",*Y.range(),*this);
00091   THYRA_ASSERT_VEC_SPACES("ScalarProdVectorSpaceBase<Scalar>::scalarProds(...)",*X.domain(),*Y.domain());
00092 #endif
00093   scalarProd_->scalarProds(X,Y,scalar_prods);
00094 }
00095 
00096 } // end namespace Thyra
00097 
00098 #endif  // THYRA_VECTOR_SPACE_STD_BASE_HPP

Generated on Thu Sep 18 12:33:03 2008 for Thyra Package Browser (Single Doxygen Collection) by doxygen 1.3.9.1