00001 // $Id: LOCA_MultiPredictor_Tangent.H,v 1.5 2007/06/21 16:22:52 rhoope Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/nox/src-loca/src/LOCA_MultiPredictor_Tangent.H,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-loca/src/LOCA_MultiPredictor_Tangent.H,v $ 00036 // $Author: rhoope $ 00037 // $Date: 2007/06/21 16:22:52 $ 00038 // $Revision: 1.5 $ 00039 // ************************************************************************ 00040 //@HEADER 00041 00042 #ifndef LOCA_MULTIPREDICTOR_TANGENT_H 00043 #define LOCA_MULTIPREDICTOR_TANGENT_H 00044 00045 #include "LOCA_MultiPredictor_AbstractStrategy.H" // base class 00046 00047 // forward declarations 00048 namespace LOCA{ 00049 class GlobalData; 00050 } 00051 00052 namespace LOCA { 00053 00054 namespace MultiPredictor { 00055 00057 00067 class Tangent : public LOCA::MultiPredictor::AbstractStrategy { 00068 00069 public: 00070 00072 00079 Tangent(const Teuchos::RCP<LOCA::GlobalData>& global_data, 00080 const Teuchos::RCP<Teuchos::ParameterList>& predParams, 00081 const Teuchos::RCP<Teuchos::ParameterList>& solverParams); 00082 00084 virtual ~Tangent(); 00085 00087 Tangent(const Tangent& source, NOX::CopyType type = NOX::DeepCopy); 00088 00090 virtual LOCA::MultiPredictor::AbstractStrategy& 00091 operator=(const LOCA::MultiPredictor::AbstractStrategy& source); 00092 00094 virtual Teuchos::RCP<LOCA::MultiPredictor::AbstractStrategy> 00095 clone(NOX::CopyType type = NOX::DeepCopy) const; 00096 00106 virtual NOX::Abstract::Group::ReturnType 00107 compute(bool baseOnSecant, const vector<double>& stepSize, 00108 LOCA::MultiContinuation::ExtendedGroup& grp, 00109 const LOCA::MultiContinuation::ExtendedVector& prevXVec, 00110 const LOCA::MultiContinuation::ExtendedVector& xVec); 00111 00113 00118 virtual NOX::Abstract::Group::ReturnType 00119 evaluate(const vector<double>& stepSize, 00120 const LOCA::MultiContinuation::ExtendedVector& xVec, 00121 LOCA::MultiContinuation::ExtendedMultiVector& result) const; 00122 00124 virtual NOX::Abstract::Group::ReturnType 00125 computeTangent(LOCA::MultiContinuation::ExtendedMultiVector& tangent); 00126 00128 00131 virtual bool isTangentScalable() const; 00132 00133 protected: 00134 00136 Teuchos::RCP<LOCA::GlobalData> globalData; 00137 00139 Teuchos::RCP<Teuchos::ParameterList> linSolverParams; 00140 00142 Teuchos::RCP<NOX::Abstract::MultiVector> fdfdp; 00143 00145 Teuchos::RCP<LOCA::MultiContinuation::ExtendedMultiVector> tangent; 00146 00148 Teuchos::RCP<LOCA::MultiContinuation::ExtendedVector> secant; 00149 00151 bool initialized; 00152 00153 }; 00154 } // namespace Predictor 00155 } // namespace LOCA 00156 00157 #endif