IFPACK Development
Loading...
Searching...
No Matches
Ifpack_DiagPreconditioner.cpp
1/*@HEADER
2// ***********************************************************************
3//
4// Ifpack: Object-Oriented Algebraic Preconditioner Package
5// Copyright (2002) Sandia Corporation
6//
7// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
8// license for use of this work by or on behalf of the U.S. Government.
9//
10// Redistribution and use in source and binary forms, with or without
11// modification, are permitted provided that the following conditions are
12// met:
13//
14// 1. Redistributions of source code must retain the above copyright
15// notice, this list of conditions and the following disclaimer.
16//
17// 2. Redistributions in binary form must reproduce the above copyright
18// notice, this list of conditions and the following disclaimer in the
19// documentation and/or other materials provided with the distribution.
20//
21// 3. Neither the name of the Corporation nor the names of the
22// contributors may be used to endorse or promote products derived from
23// this software without specific prior written permission.
24//
25// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36//
37// Questions? Contact Michael A. Heroux (maherou@sandia.gov)
38//
39// ***********************************************************************
40//@HEADER
41*/
42
43#include "Ifpack_ConfigDefs.h"
44#include "Ifpack_DiagPreconditioner.h"
45#include "Epetra_BlockMap.h"
46#include "Epetra_Map.h"
47#include "Epetra_MultiVector.h"
48#include "Epetra_Vector.h"
49#include "Epetra_Comm.h"
50
51// ============================================================================
54 const Epetra_Map& RangeMap,
55 const Epetra_Vector& diag) :
56 UseTranspose_(false),
57 DomainMap_(DomainMap),
58 RangeMap_(RangeMap),
59 diag_(diag)
60{ }
61
62// ============================================================================
64{ }
65
66
67// ============================================================================
68int Ifpack_DiagPreconditioner::Apply(const Epetra_MultiVector& /* X */, Epetra_MultiVector& /* Y */) const
69{
70 IFPACK_RETURN(-1); // not defined
71}
72
73// ============================================================================
75{
76 if (X.NumVectors() != Y.NumVectors())
77 IFPACK_CHK_ERR(-1);
78
79 for (int v = 0; v < X.NumVectors(); ++v)
80 for (int i = 0; i < X.MyLength(); ++i)
81 Y[v][i] = diag_[i] * X[v][i];
83
84 return(0);
85}
int NumVectors() const
int MyLength() const
Ifpack_DiagPreconditioner(const Epetra_Map &DomainMap, const Epetra_Map &RangeMap, const Epetra_Vector &diag)
ctor
int ApplyInverse(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const