00001 // $Id: NOX_Belos_PreconditionOperator.C,v 1.3 2006/08/22 00:01:31 rppawlo Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/nox/src-belos/NOX_Belos_PreconditionOperator.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_PreconditionOperator.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_PreconditionOperator.H" 00043 #include "NOX_Belos_MultiVector.H" 00044 #include "NOX_Abstract_MultiVector.H" 00045 #include "NOX_Parameter_List.H" 00046 00047 NOX::Belos::PreconditionOperator::PreconditionOperator( 00048 NOX::Abstract::Group& g, 00049 NOX::Parameter::List& preconditionerParameters) 00050 : grp(g), 00051 precondParams(preconditionerParameters) 00052 { 00053 } 00054 00055 00056 NOX::Belos::PreconditionOperator::~PreconditionOperator() 00057 { 00058 } 00059 00060 ::Belos::ReturnType 00061 NOX::Belos::PreconditionOperator::Apply( 00062 const ::Belos::MultiVec<double>& x, 00063 ::Belos::MultiVec<double>& y, 00064 ::Belos::ETrans trans) const 00065 { 00066 // Cast x and y to NOX::Belos::MultiVec's 00067 const NOX::Belos::MultiVector& nox_belos_x = 00068 dynamic_cast<const NOX::Belos::MultiVector&>(x); 00069 NOX::Belos::MultiVector& nox_belos_y = 00070 dynamic_cast<NOX::Belos::MultiVector&>(y); 00071 00072 // Get underlying NOX::Abstract::MultiVector's 00073 const NOX::Abstract::MultiVector& nox_x = nox_belos_x.getNoxMultiVector(); 00074 NOX::Abstract::MultiVector& nox_y = nox_belos_y.getNoxMultiVector(); 00075 00076 // NOX return type 00077 NOX::Abstract::Group::ReturnType nox_status; 00078 00079 bool useTranspose = false; 00080 if (trans == ::Belos::TRANS) 00081 useTranspose = true; 00082 00083 nox_status = grp.applyRightPreconditioningMultiVector(useTranspose, 00084 precondParams, 00085 nox_x, 00086 nox_y); 00087 00088 return noxReturnTypeToBelos(nox_status); 00089 } 00090 00091 ::Belos::ReturnType 00092 NOX::Belos::PreconditionOperator::ApplyInverse( 00093 const ::Belos::MultiVec<double>& x, 00094 ::Belos::MultiVec<double>& y, 00095 ::Belos::ETrans trans) const 00096 { 00097 return ::Belos::Undefined; 00098 } 00099 00100 ::Belos::ReturnType 00101 NOX::Belos::PreconditionOperator::noxReturnTypeToBelos( 00102 NOX::Abstract::Group::ReturnType noxStatus) const 00103 { 00104 if (noxStatus == NOX::Abstract::Group::Ok || 00105 noxStatus == NOX::Abstract::Group::NotConverged) 00106 return ::Belos::Ok; 00107 else if (noxStatus == NOX::Abstract::Group::NotDefined) 00108 return ::Belos::Undefined; 00109 else 00110 return ::Belos::Error; 00111 } 00112 00113