00001 // $Id: LOCA_MultiPredictor_AbstractStrategy.C,v 1.3 2006/08/22 00:01:37 rppawlo Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/nox/src-loca/src/LOCA_MultiPredictor_AbstractStrategy.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-loca/src/LOCA_MultiPredictor_AbstractStrategy.C,v $ 00036 // $Author: rppawlo $ 00037 // $Date: 2006/08/22 00:01:37 $ 00038 // $Revision: 1.3 $ 00039 // ************************************************************************ 00040 //@HEADER 00041 00042 #include "LOCA_MultiPredictor_AbstractStrategy.H" 00043 #include "LOCA_MultiContinuation_ExtendedVector.H" 00044 #include "LOCA_MultiContinuation_ExtendedMultiVector.H" 00045 #include "LOCA_MultiContinuation_ExtendedGroup.H" 00046 00047 void 00048 LOCA::MultiPredictor::AbstractStrategy::setPredictorOrientation( 00049 bool baseOnSecant, 00050 const vector<double>& stepSize, 00051 const LOCA::MultiContinuation::ExtendedGroup& grp, 00052 const LOCA::MultiContinuation::ExtendedVector& prevXVec, 00053 const LOCA::MultiContinuation::ExtendedVector& xVec, 00054 LOCA::MultiContinuation::ExtendedVector& secant, 00055 LOCA::MultiContinuation::ExtendedMultiVector& tangent) 00056 { 00057 // If orientation is not based on a secant vector (i.e., first or last 00058 // steps in a continuation run) make parameter component of predictor 00059 // positive 00060 if (!baseOnSecant) { 00061 for (int i=0; i<tangent.numVectors(); i++) 00062 if (tangent.getScalar(i,i) < 0.0) 00063 tangent[i].scale(-1.0); 00064 return; 00065 } 00066 00067 // compute secant 00068 secant.update(1.0, xVec, -1.0, prevXVec, 0.0); 00069 00070 for (int i=0; i<tangent.numVectors(); i++) 00071 if (grp.computeScaledDotProduct(secant, tangent[i])*stepSize[i] < 0.0) 00072 tangent[i].scale(-1.0); 00073 }