EpetraExt
Development
Loading...
Searching...
No Matches
src
transform
EpetraExt_SolverMap_LinearProblem.cpp
Go to the documentation of this file.
1
//@HEADER
2
// ***********************************************************************
3
//
4
// EpetraExt: Epetra Extended - Linear Algebra Services Package
5
// Copyright (2011) Sandia Corporation
6
//
7
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
8
// the U.S. Government retains certain rights in this software.
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
#include <
EpetraExt_SolverMap_LinearProblem.h
>
43
44
#include <Epetra_LinearProblem.h>
45
#include <Epetra_CrsMatrix.h>
46
#include <Epetra_Map.h>
47
48
class
Epetra_MultiVector
;
49
50
namespace
EpetraExt
{
51
52
LinearProblem_SolverMap::
53
~LinearProblem_SolverMap
()
54
{
55
if
(
newObj_
&&
newObj_
!=
origObj_
)
delete
newObj_
;
56
}
57
58
LinearProblem_SolverMap::NewTypeRef
59
LinearProblem_SolverMap::
60
operator()
( OriginalTypeRef orig )
61
{
62
origObj_
= &orig;
63
64
Epetra_CrsMatrix
* OldMatrix =
dynamic_cast<
Epetra_CrsMatrix
*
>
( orig.GetMatrix() );
65
Epetra_MultiVector
* OldRHS = orig.GetRHS();
66
Epetra_MultiVector
* OldLHS = orig.GetLHS();
67
68
Epetra_CrsMatrix
& NewMatrix = SolverMapTrans_( *OldMatrix );
69
70
if
( &NewMatrix == OldMatrix )
//same matrix so use same problem
71
newObj_
=
origObj_
;
72
else
73
newObj_
=
new
Epetra_LinearProblem
( &NewMatrix, OldLHS, OldRHS );
74
75
return
*
newObj_
;
76
}
77
78
}
//namespace EpetraExt
79
EpetraExt_SolverMap_LinearProblem.h
EpetraExt::LinearProblem_SolverMap::~LinearProblem_SolverMap
~LinearProblem_SolverMap()
Destructor.
Definition
EpetraExt_SolverMap_LinearProblem.cpp:53
EpetraExt::LinearProblem_SolverMap::operator()
NewTypeRef operator()(OriginalTypeRef orig)
Constructs "fixed" Epetra_LinearProblem.
Definition
EpetraExt_SolverMap_LinearProblem.cpp:60
EpetraExt::Transform< T, T >::newObj_
NewTypePtr newObj_
Definition
EpetraExt_Transform.h:218
EpetraExt::Transform< T, T >::origObj_
OriginalTypePtr origObj_
Definition
EpetraExt_Transform.h:216
EpetraExt::Transform< T, T >::NewTypeRef
T & NewTypeRef
Definition
EpetraExt_Transform.h:79
Epetra_CrsMatrix
Epetra_LinearProblem
Epetra_MultiVector
EpetraExt
EpetraExt::BlockCrsMatrix: A class for constructing a distributed block matrix.
Definition
EpetraExt_BlockCrsMatrix.cpp:46
Generated by
1.10.0