00001 // $Id: NOX_Solver_Factory.C,v 1.1 2007/08/01 20:50:40 rppawlo Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/nox/src/NOX_Solver_Factory.C,v $ 00003 00004 //@HEADER 00005 // ************************************************************************ 00006 // 00007 // NOX: An Object-Oriented Nonlinear Solver Package 00008 // Copyright (2002) Sandia Corporation 00009 // 00010 // LOCA: Library of Continuation Algorithms Package 00011 // Copyright (2005) Sandia Corporation 00012 // 00013 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive 00014 // license for use of this work by or on behalf of the U.S. Government. 00015 // 00016 // This library is free software; you can redistribute it and/or modify 00017 // it under the terms of the GNU Lesser General Public License as 00018 // published by the Free Software Foundation; either version 2.1 of the 00019 // License, or (at your option) any later version. 00020 // 00021 // This library is distributed in the hope that it will be useful, but 00022 // WITHOUT ANY WARRANTY; without even the implied warranty of 00023 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00024 // Lesser General Public License for more details. 00025 // 00026 // You should have received a copy of the GNU Lesser General Public 00027 // License along with this library; if not, write to the Free Software 00028 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 00029 // USA 00030 // 00031 // Questions? Contact Roger Pawlowski (rppawlo@sandia.gov) or 00032 // Eric Phipps (etphipp@sandia.gov), Sandia National Laboratories. 00033 // ************************************************************************ 00034 // CVS Information 00035 // $Source: /space/CVS/Trilinos/packages/nox/src/NOX_Solver_Factory.C,v $ 00036 // $Author: rppawlo $ 00037 // $Date: 2007/08/01 20:50:40 $ 00038 // $Revision: 1.1 $ 00039 // ************************************************************************ 00040 //@HEADER 00041 00042 #include "NOX_Common.H" 00043 #include "NOX_Solver_Factory.H" // class definition 00044 #include "NOX_Solver_Generic.H" 00045 00046 // Header files for different solvers 00047 #include "NOX_Solver_LineSearchBased.H" // LineSearch method 00048 #include "NOX_Solver_TrustRegionBased.H" // Trust region method 00049 #include "NOX_Solver_InexactTrustRegionBased.H" // Inexact Trust region method 00050 #include "NOX_Solver_TensorBased.H" // Tensor method 00051 #ifdef WITH_PRERELEASE 00052 #include "NOX_Solver_TensorBasedTest.H" // Tensor-Krylov method 00053 #endif 00054 00055 // ************************************************************************ 00056 // ************************************************************************ 00057 NOX::Solver::Factory::Factory() 00058 { } 00059 00060 // ************************************************************************ 00061 // ************************************************************************ 00062 NOX::Solver::Factory::~Factory() 00063 { } 00064 00065 // ************************************************************************ 00066 // ************************************************************************ 00067 Teuchos::RCP<NOX::Solver::Generic> NOX::Solver::Factory:: 00068 buildSolver(const Teuchos::RCP<NOX::Abstract::Group>& grp, 00069 const Teuchos::RCP<NOX::StatusTest::Generic>& tests, 00070 const Teuchos::RCP<Teuchos::ParameterList>& params) 00071 { 00072 using Teuchos::RCP; 00073 using Teuchos::rcp; 00074 RCP<NOX::Solver::Generic> solver; 00075 00076 string method = params->get("Nonlinear Solver", "Line Search Based"); 00077 00078 if ((method == "Newton") || (method == "Line Search Based")) 00079 solver = rcp(new LineSearchBased(grp, tests, params)); 00080 else if (method == "Trust Region Based") 00081 solver = rcp(new TrustRegionBased(grp, tests, params)); 00082 else if (method == "Inexact Trust Region Based") 00083 solver = rcp(new InexactTrustRegionBased(grp, tests, params)); 00084 else if (method == "Tensor Based") 00085 solver = rcp(new TensorBased(grp, tests, params)); 00086 #ifdef WITH_PRERELEASE 00087 else if (method == "Tensor-Krylov Based") 00088 solver = rcp(new TensorBasedTest(grp, tests, params)); 00089 #endif 00090 else { 00091 std::ostringstream msg; 00092 msg << "Error - NOX::Solver::Manager::buildSolver() - The \"Nonlinear Solver\" parameter \"" << method << "\" is not a valid solver option. Please fix your parameter list!"; 00093 TEST_FOR_EXCEPTION(true, std::logic_error, msg.str()); 00094 } 00095 00096 return solver; 00097 } 00098 00099 // ************************************************************************ 00100 // ************************************************************************ 00101 // Nonmember function 00102 Teuchos::RCP<NOX::Solver::Generic> 00103 NOX::Solver::buildSolver(const Teuchos::RCP<NOX::Abstract::Group>& grp, 00104 const Teuchos::RCP<NOX::StatusTest::Generic>& tests, 00105 const Teuchos::RCP<Teuchos::ParameterList>& params) 00106 { 00107 NOX::Solver::Factory factory; 00108 return factory.buildSolver(grp, tests, params); 00109 } 00110 00111 // ************************************************************************ 00112 // ************************************************************************ 00113