00001 // $Id: LOCA_MultiContinuation_ArcLengthGroup.H,v 1.13 2007/06/21 16:22:52 rhoope Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/nox/src-loca/src/LOCA_MultiContinuation_ArcLengthGroup.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_MultiContinuation_ArcLengthGroup.H,v $ 00036 // $Author: rhoope $ 00037 // $Date: 2007/06/21 16:22:52 $ 00038 // $Revision: 1.13 $ 00039 // ************************************************************************ 00040 //@HEADER 00041 00042 #ifndef LOCA_MULTICONTINUATION_ARCLENGTHGROUP_H 00043 #define LOCA_MULTICONTINUATION_ARCLENGTHGROUP_H 00044 00045 #include "LOCA_MultiContinuation_ExtendedGroup.H" // base class 00046 00047 namespace LOCA { 00048 00049 namespace MultiContinuation { 00050 00095 class ArcLengthGroup : 00096 public virtual LOCA::MultiContinuation::ExtendedGroup { 00097 00098 public: 00099 00101 00110 ArcLengthGroup( 00111 const Teuchos::RCP<LOCA::GlobalData>& global_data, 00112 const Teuchos::RCP<LOCA::Parameter::SublistParser>& topParams, 00113 const Teuchos::RCP<Teuchos::ParameterList>& continuationParams, 00114 const Teuchos::RCP<LOCA::MultiContinuation::AbstractGroup>& grp, 00115 const Teuchos::RCP<LOCA::MultiPredictor::AbstractStrategy>& pred, 00116 const vector<int>& paramIDs); 00117 00119 ArcLengthGroup(const ArcLengthGroup& source, 00120 NOX::CopyType type = NOX::DeepCopy); 00121 00123 virtual ~ArcLengthGroup(); 00124 00129 00131 virtual NOX::Abstract::Group& 00132 operator=(const NOX::Abstract::Group& source); 00133 00135 virtual Teuchos::RCP<NOX::Abstract::Group> 00136 clone(NOX::CopyType type = NOX::DeepCopy) const; 00137 00139 00144 00146 virtual void copy(const NOX::Abstract::Group& source); 00147 00149 virtual void scaleTangent(); 00150 00152 virtual double computeScaledDotProduct( 00153 const NOX::Abstract::Vector& x, 00154 const NOX::Abstract::Vector& y) const; 00155 00157 virtual void recalculateScaleFactor(double dpds, double thetaOld, 00158 double& thetaNew); 00159 00161 00162 private: 00163 00165 ArcLengthGroup& operator=(const ArcLengthGroup& source); 00166 00167 protected: 00168 00170 vector<double> theta; 00171 00173 bool doArcLengthScaling; 00174 00176 double gGoal; 00177 00179 double gMax; 00180 00182 double thetaMin; 00183 00185 bool isFirstRescale; 00186 }; 00187 00188 } // namespace MultiContinuation 00189 } // namespace LOCA 00190 00191 #endif