00001 // $Id: NOX_Belos_JacobianOperator.C,v 1.3 2006/08/22 00:01:31 rppawlo Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/nox/src-belos/NOX_Belos_JacobianOperator.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-belos/NOX_Belos_JacobianOperator.C,v $ 00036 // $Author: rppawlo $ 00037 // $Date: 2006/08/22 00:01:31 $ 00038 // $Revision: 1.3 $ 00039 // ************************************************************************ 00040 //@HEADER 00041 00042 #include "NOX_Belos_JacobianOperator.H" 00043 #include "NOX_Belos_MultiVector.H" 00044 #include "NOX_Abstract_MultiVector.H" 00045 00046 NOX::Belos::JacobianOperator::JacobianOperator(NOX::Abstract::Group& g) 00047 : grp(g) 00048 { 00049 } 00050 00051 00052 NOX::Belos::JacobianOperator::~JacobianOperator() 00053 { 00054 } 00055 00056 ::Belos::ReturnType 00057 NOX::Belos::JacobianOperator::Apply( 00058 const ::Belos::MultiVec<double>& x, 00059 ::Belos::MultiVec<double>& y, 00060 ::Belos::ETrans trans) const 00061 { 00062 // Cast x and y to NOX::Belos::MultiVec's 00063 const NOX::Belos::MultiVector& nox_belos_x = 00064 dynamic_cast<const NOX::Belos::MultiVector&>(x); 00065 NOX::Belos::MultiVector& nox_belos_y = 00066 dynamic_cast<NOX::Belos::MultiVector&>(y); 00067 00068 // Get underlying NOX::Abstract::MultiVector's 00069 const NOX::Abstract::MultiVector& nox_x = nox_belos_x.getNoxMultiVector(); 00070 NOX::Abstract::MultiVector& nox_y = nox_belos_y.getNoxMultiVector(); 00071 00072 // NOX return type 00073 NOX::Abstract::Group::ReturnType nox_status; 00074 00075 if (trans == ::Belos::NOTRANS) 00076 nox_status = grp.applyJacobianMultiVector(nox_x,nox_y); 00077 else 00078 nox_status = grp.applyJacobianTransposeMultiVector(nox_x,nox_y); 00079 00080 return noxReturnTypeToBelos(nox_status); 00081 } 00082 00083 ::Belos::ReturnType 00084 NOX::Belos::JacobianOperator::ApplyInverse( 00085 const ::Belos::MultiVec<double>& x, 00086 ::Belos::MultiVec<double>& y, 00087 ::Belos::ETrans trans) const 00088 { 00089 return ::Belos::Undefined; 00090 } 00091 00092 ::Belos::ReturnType 00093 NOX::Belos::JacobianOperator::noxReturnTypeToBelos( 00094 NOX::Abstract::Group::ReturnType noxStatus) const 00095 { 00096 if (noxStatus == NOX::Abstract::Group::Ok || 00097 noxStatus == NOX::Abstract::Group::NotConverged) 00098 return ::Belos::Ok; 00099 else if (noxStatus == NOX::Abstract::Group::NotDefined) 00100 return ::Belos::Undefined; 00101 else 00102 return ::Belos::Error; 00103 } 00104 00105