EpetraExt Package Browser (Single Doxygen Collection)
Development
Loading...
Searching...
No Matches
src
transform
EpetraExt_Reindex_MultiVector.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_Reindex_MultiVector.h
>
43
44
#include <vector>
45
46
#include <
Epetra_Map.h
>
47
#include <
Epetra_MultiVector.h
>
48
49
namespace
EpetraExt
{
50
51
MultiVector_Reindex::
52
~MultiVector_Reindex
()
53
{
54
if
(
newObj_
)
delete
newObj_
;
55
}
56
57
MultiVector_Reindex::NewTypeRef
58
MultiVector_Reindex::
59
operator()
( OriginalTypeRef orig )
60
{
61
origObj_
= &orig;
62
63
//test std::map, must have same number of local and global elements as original row std::map
64
assert( orig.Map().NumMyElements() ==
NewRowMap_
.NumMyElements() );
65
66
std::vector<double*> MyValues(1);
67
int
MyLDA;
68
int
NumVectors = orig.NumVectors();
69
orig.ExtractView( &MyValues[0], &MyLDA );
70
71
Epetra_MultiVector
* NewMV =
new
Epetra_MultiVector
(
View
,
NewRowMap_
, MyValues[0], MyLDA, NumVectors );
72
73
newObj_
= NewMV;
74
75
return
*NewMV;
76
}
77
78
MultiVector_Reindex::NewTypeRCP
79
MultiVector_Reindex::transform
(OriginalTypeRef orig)
80
{
81
//test std::map, must have same number of local and global elements as original row std::map
82
assert( orig.Map().NumMyElements() ==
NewRowMap_
.NumMyElements() );
83
84
std::vector<double*> MyValues(1);
85
int
MyLDA;
86
int
NumVectors = orig.NumVectors();
87
orig.ExtractView( &MyValues[0], &MyLDA );
88
89
return
Teuchos::rcp(
new
Epetra_MultiVector
(
View
,
NewRowMap_
, MyValues[0], MyLDA, NumVectors ));
90
}
91
92
}
// namespace EpetraExt
93
EpetraExt_Reindex_MultiVector.h
View
View
Epetra_Map.h
Epetra_MultiVector.h
EpetraExt::MultiVector_Reindex::operator()
NewTypeRef operator()(OriginalTypeRef orig)
Constructs a "reindexed" view of the original using the given NewRowMap.
Definition
EpetraExt_Reindex_MultiVector.cpp:59
EpetraExt::MultiVector_Reindex::transform
NewTypeRCP transform(OriginalTypeRef orig)
Definition
EpetraExt_Reindex_MultiVector.cpp:79
EpetraExt::MultiVector_Reindex::~MultiVector_Reindex
~MultiVector_Reindex()
Destructor.
Definition
EpetraExt_Reindex_MultiVector.cpp:52
EpetraExt::MultiVector_Reindex::NewRowMap_
const Epetra_Map & NewRowMap_
Definition
EpetraExt_Reindex_MultiVector.h:57
EpetraExt::Transform< T, T >::NewTypeRCP
Teuchos::RCP< T > NewTypeRCP
Definition
EpetraExt_Transform.h:78
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_MultiVector
EpetraExt
EpetraExt::BlockCrsMatrix: A class for constructing a distributed block matrix.
Definition
EpetraExt_BlockCrsMatrix.cpp:46
Generated by
1.10.0