00001 // $Id: LOCA_MultiContinuation_ConstraintInterfaceMVDX.C,v 1.4 2006/08/22 00:01:37 rppawlo Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/nox/src-loca/src/LOCA_MultiContinuation_ConstraintInterfaceMVDX.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_MultiContinuation_ConstraintInterfaceMVDX.C,v $ 00036 // $Author: rppawlo $ 00037 // $Date: 2006/08/22 00:01:37 $ 00038 // $Revision: 1.4 $ 00039 // ************************************************************************ 00040 //@HEADER 00041 00042 #include "LOCA_MultiContinuation_ConstraintInterfaceMVDX.H" 00043 00044 NOX::Abstract::Group::ReturnType 00045 LOCA::MultiContinuation::ConstraintInterfaceMVDX::multiplyDX( 00046 double alpha, 00047 const NOX::Abstract::MultiVector& input_x, 00048 NOX::Abstract::MultiVector::DenseMatrix& result_p) const 00049 { 00050 00051 00052 if (!isDXZero()) { 00053 const NOX::Abstract::MultiVector* dgdx = getDX(); 00054 input_x.multiply(alpha, *dgdx, result_p); 00055 } 00056 else 00057 result_p.putScalar(0.0); 00058 00059 return NOX::Abstract::Group::Ok; 00060 } 00061 00062 NOX::Abstract::Group::ReturnType 00063 LOCA::MultiContinuation::ConstraintInterfaceMVDX::addDX( 00064 Teuchos::ETransp transb, 00065 double alpha, 00066 const NOX::Abstract::MultiVector::DenseMatrix& b, 00067 double beta, 00068 NOX::Abstract::MultiVector& result_x) const 00069 { 00070 if (!isDXZero()) { 00071 const NOX::Abstract::MultiVector* dgdx = getDX(); 00072 result_x.update(transb, alpha, *dgdx, b, beta); 00073 } 00074 else 00075 result_x.scale(beta); 00076 00077 return NOX::Abstract::Group::Ok; 00078 }