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 "TSFDenseLUSolver.hpp" 00030 #include "TSFDenseSerialMatrix.hpp" 00031 #include "TSFLinearOperatorDecl.hpp" 00032 00033 #ifndef HAVE_TEUCHOS_EXPLICIT_INSTANTIATION 00034 #include "TSFLinearOperatorImpl.hpp" 00035 #include "TSFVectorImpl.hpp" 00036 #endif 00037 00038 using namespace TSFExtended; 00039 using namespace Teuchos; 00040 using namespace Thyra; 00041 using std::setw; 00042 00043 00044 DenseLUSolver::DenseLUSolver() 00045 : LinearSolverBase<double>(ParameterList()) 00046 { 00047 } 00048 00049 SolverState<double> DenseLUSolver::solve(const LinearOperator<double>& op, 00050 const Vector<double>& rhs, 00051 Vector<double>& soln) const 00052 { 00053 return denseSolve(op, rhs, soln); 00054 }