00001 /* @HEADER@ */ 00002 /* *********************************************************************** 00003 // 00004 // TSFExtended: Trilinos Solver Framework Extended 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 #include "TSFDenseSerialMatrixFactory.hpp" 00030 #include "TSFDenseSerialMatrix.hpp" 00031 #include "TSFSerialVector.hpp" 00032 #include "TSFVectorSpaceDecl.hpp" 00033 #include "TSFVectorDecl.hpp" 00034 #include "Teuchos_Array.hpp" 00035 #include "Teuchos_MPIComm.hpp" 00036 #include "TSFLinearOperatorDecl.hpp" 00037 00038 #ifndef HAVE_TEUCHOS_EXPLICIT_INSTANTIATION 00039 #include "TSFLinearOperatorImpl.hpp" 00040 #include "TSFVectorImpl.hpp" 00041 #endif 00042 00043 00044 using namespace TSFExtended; 00045 using namespace Teuchos; 00046 using namespace Thyra; 00047 00048 DenseSerialMatrixFactory::DenseSerialMatrixFactory( 00049 const RCP<const SerialVectorSpace>& domain, 00050 const RCP<const SerialVectorSpace>& range) 00051 : range_(range), 00052 domain_(domain) 00053 {} 00054 00055 00056 LinearOperator<double> DenseSerialMatrixFactory::createMatrix() const 00057 { 00058 RCP<LinearOpBase<double> > A 00059 = rcp(new DenseSerialMatrix(domain_, range_)); 00060 return A; 00061 }