00001 /* @HEADER@ */ 00002 // ************************************************************************ 00003 // 00004 // Playa: Programmable Linear Algebra 00005 // Copyright 2012 Sandia Corporation 00006 // 00007 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation, 00008 // the U.S. Government retains certain rights in this software. 00009 // 00010 // Redistribution and use in source and binary forms, with or without 00011 // modification, are permitted provided that the following conditions are 00012 // met: 00013 // 00014 // 1. Redistributions of source code must retain the above copyright 00015 // notice, this list of conditions and the following disclaimer. 00016 // 00017 // 2. Redistributions in binary form must reproduce the above copyright 00018 // notice, this list of conditions and the following disclaimer in the 00019 // documentation and/or other materials provided with the distribution. 00020 // 00021 // 3. Neither the name of the Corporation nor the names of the 00022 // contributors may be used to endorse or promote products derived from 00023 // this software without specific prior written permission. 00024 // 00025 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY 00026 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 00027 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR 00028 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE 00029 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, 00030 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, 00031 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 00032 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF 00033 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING 00034 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS 00035 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 00036 // 00037 // Questions? Contact Kevin Long (kevin.long@ttu.edu) 00038 // 00039 00040 /* @HEADER@ */ 00041 00042 #ifndef PLAYA_EPETRAVECTORSPACE_HPP 00043 #define PLAYA_EPETRAVECTORSPACE_HPP 00044 00045 #include "PlayaDefs.hpp" 00046 #include "Epetra_Map.h" 00047 #include "Teuchos_RefCountPtr.hpp" 00048 #include "PlayaMPIComm.hpp" 00049 #include "PlayaVectorSpaceBaseDecl.hpp" 00050 00051 namespace Playa 00052 { 00053 using namespace Teuchos; 00054 00055 00056 /** 00057 * Adaptor wrapping Epetra map in the Thyra vector space system. 00058 */ 00059 class EpetraVectorSpace : public VectorSpaceBase<double> 00060 { 00061 public: 00062 /** */ 00063 EpetraVectorSpace(const RCP<const Epetra_Map>& map); 00064 00065 /** */ 00066 RCP<VectorBase<double> > 00067 createMember(const VectorSpace<double>& self) const ; 00068 00069 /** @name Overridden form Teuchos::Describable */ 00070 //@{ 00071 /** \brief . */ 00072 std::string description() const; 00073 //@} 00074 00075 /** @name Public overridden from VectorSpace */ 00076 //@{ 00077 00078 /** */ 00079 int dim() const {return globalDim_;} 00080 00081 /** */ 00082 int numLocalElements() const {return numLocalElements_;} 00083 00084 /** */ 00085 int baseGlobalNaturalIndex() const {return baseGlobalNaturalIndex_;} 00086 00087 /** */ 00088 bool isCompatible(const VectorSpaceBase<double>* other) const ; 00089 00090 int numBlocks() const {return 1;} 00091 00092 //@} 00093 00094 /** */ 00095 const RCP<const Epetra_Map>& epetraMap() const 00096 {return epetraMap_;} 00097 00098 00099 /** */ 00100 virtual const MPIComm& comm() const {return comm_;} 00101 00102 protected: 00103 00104 MPIComm epetraCommToTeuchosMPIComm(const Epetra_Comm& epComm) ; 00105 00106 private: 00107 RCP<const Epetra_Map> epetraMap_; 00108 00109 MPIComm comm_; 00110 00111 int globalDim_; 00112 00113 int baseGlobalNaturalIndex_; 00114 00115 int numLocalElements_; 00116 }; 00117 00118 } 00119 00120 #endif