EpetraExt Package Browser (Single Doxygen Collection)
Development
Loading...
Searching...
No Matches
src
transform
EpetraExt_Overlap_CrsGraph.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_Overlap_CrsGraph.h
>
43
44
#include <
Epetra_Import.h
>
45
#include <
Epetra_CrsGraph.h
>
46
#include <
Epetra_Map.h
>
47
48
namespace
EpetraExt
{
49
50
CrsGraph_Overlap::
51
~CrsGraph_Overlap
()
52
{
53
if
(
newObj_
)
delete
newObj_
;
54
55
if
(
OverlapMap_
)
delete
OverlapMap_
;
56
}
57
58
CrsGraph_Overlap::NewTypeRef
59
CrsGraph_Overlap::
60
operator()
( OriginalTypeRef orig )
61
{
62
origObj_
= &orig;
63
64
//check that this is a distributed graph and overlap level is not zero
65
if
( orig.DomainMap().DistributedGlobal() &&
levelOverlap_
)
66
{
67
Epetra_CrsGraph
* OverlapGraph =
new
Epetra_CrsGraph
( orig );
68
OverlapMap_
=
new
Epetra_BlockMap
( orig.RowMap() );
69
70
Epetra_BlockMap
* DomainMap = &(
const_cast<
Epetra_BlockMap
&
>
(orig.DomainMap()));
71
Epetra_BlockMap
* RangeMap = &(
const_cast<
Epetra_BlockMap
&
>
(orig.RangeMap()));
72
73
for
(
int
level = 0; level <
levelOverlap_
; ++level )
74
{
75
Epetra_BlockMap
* OldRowMap =
OverlapMap_
;
76
Epetra_CrsGraph
* OldGraph = OverlapGraph;
77
78
Epetra_Import
& OverlapImporter = *(
const_cast<
Epetra_Import
*
>
( OldGraph->Importer() ));
79
OverlapMap_
=
new
Epetra_BlockMap
( OverlapImporter.TargetMap() );
80
81
//filter to local square block on last level if required
82
if
(
squareLocalBlock_
&& level==(
levelOverlap_
-1) )
83
OverlapGraph =
new
Epetra_CrsGraph
(
Copy
, *
OverlapMap_
, *
OverlapMap_
, 0 );
84
else
85
OverlapGraph =
new
Epetra_CrsGraph
(
Copy
, *
OverlapMap_
, 0 );
86
87
OverlapGraph->Import( *OldGraph, OverlapImporter,
Insert
);
88
OverlapGraph->FillComplete( *DomainMap, *RangeMap );
89
90
delete
OldGraph;
91
delete
OldRowMap;
92
}
93
94
newObj_
= OverlapGraph;
95
}
96
else
//just create a copy since this is not a InPlaceTransform
97
newObj_
=
new
Epetra_CrsGraph
( orig );
98
99
return
*
newObj_
;
100
}
101
102
}
// namespace EpetraExt
EpetraExt_Overlap_CrsGraph.h
Insert
Insert
Epetra_CrsGraph.h
Copy
Copy
Epetra_Import.h
Epetra_Map.h
EpetraExt::CrsGraph_Overlap::squareLocalBlock_
const bool squareLocalBlock_
Definition
EpetraExt_Overlap_CrsGraph.h:60
EpetraExt::CrsGraph_Overlap::operator()
NewTypeRef operator()(OriginalTypeRef orig)
Constructs "overlapped" Epetra_CrsGraph from original.
Definition
EpetraExt_Overlap_CrsGraph.cpp:60
EpetraExt::CrsGraph_Overlap::~CrsGraph_Overlap
~CrsGraph_Overlap()
Destructor.
Definition
EpetraExt_Overlap_CrsGraph.cpp:51
EpetraExt::CrsGraph_Overlap::OverlapMap_
Epetra_BlockMap * OverlapMap_
Definition
EpetraExt_Overlap_CrsGraph.h:62
EpetraExt::CrsGraph_Overlap::levelOverlap_
const int levelOverlap_
Definition
EpetraExt_Overlap_CrsGraph.h:58
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_BlockMap
Epetra_CrsGraph
Epetra_Import
EpetraExt
EpetraExt::BlockCrsMatrix: A class for constructing a distributed block matrix.
Definition
EpetraExt_BlockCrsMatrix.cpp:46
Generated by
1.10.0