IFPACK Development
Ifpack_DiagPreconditioner.cpp
00001 /*@HEADER
00002 // ***********************************************************************
00003 //
00004 //       Ifpack: Object-Oriented Algebraic Preconditioner Package
00005 //                 Copyright (2002) Sandia Corporation
00006 //
00007 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00008 // license for use of this work by or on behalf of the U.S. Government.
00009 //
00010 // This library is free software; you can redistribute it and/or modify
00011 // it under the terms of the GNU Lesser General Public License as
00012 // published by the Free Software Foundation; either version 2.1 of the
00013 // License, or (at your option) any later version.
00014 //
00015 // This library is distributed in the hope that it will be useful, but
00016 // WITHOUT ANY WARRANTY; without even the implied warranty of
00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00018 // Lesser General Public License for more details.
00019 //
00020 // You should have received a copy of the GNU Lesser General Public
00021 // License along with this library; if not, write to the Free Software
00022 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
00023 // USA
00024 // Questions? Contact Michael A. Heroux (maherou@sandia.gov)
00025 //
00026 // ***********************************************************************
00027 //@HEADER
00028 */
00029 
00030 #include "Ifpack_ConfigDefs.h"
00031 #include "Ifpack_DiagPreconditioner.h"
00032 #include "Epetra_BlockMap.h"
00033 #include "Epetra_Map.h"
00034 #include "Epetra_MultiVector.h"
00035 #include "Epetra_Vector.h"
00036 #include "Epetra_Comm.h"
00037 
00038 // ============================================================================ 
00039 Ifpack_DiagPreconditioner::
00040 Ifpack_DiagPreconditioner(const Epetra_Map& DomainMap,
00041                           const Epetra_Map& RangeMap,
00042                           const Epetra_Vector& diag) : 
00043   UseTranspose_(false),
00044   DomainMap_(DomainMap),
00045   RangeMap_(RangeMap),
00046   diag_(diag)
00047 { }
00048 
00049 // ============================================================================ 
00050 Ifpack_DiagPreconditioner::~Ifpack_DiagPreconditioner()
00051 { }
00052 
00053 
00054 // ============================================================================ 
00055 int Ifpack_DiagPreconditioner::Apply(const Epetra_MultiVector& X, Epetra_MultiVector& Y) const
00056 {
00057   IFPACK_RETURN(-1); // not defined
00058 }
00059 
00060 // ============================================================================ 
00061  int Ifpack_DiagPreconditioner::ApplyInverse(const Epetra_MultiVector& X, Epetra_MultiVector& Y) const
00062 {
00063   if (X.NumVectors() != Y.NumVectors())
00064     IFPACK_CHK_ERR(-1); 
00065 
00066   for (int v = 0; v < X.NumVectors(); ++v)
00067     for (int i = 0; i < X.MyLength(); ++i)
00068       Y[v][i] = diag_[i] * X[v][i];
00070 
00071   return(0);
00072 }
 All Classes Files Functions Variables Enumerations Friends