Xpetra
Version of the Day
Loading...
Searching...
No Matches
Xpetra_MultiVector_def.hpp
Go to the documentation of this file.
1
// @HEADER
2
//
3
// ***********************************************************************
4
//
5
// Xpetra: A linear algebra interface package
6
// Copyright 2012 Sandia Corporation
7
//
8
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
9
// the U.S. Government retains certain rights in this software.
10
//
11
// Redistribution and use in source and binary forms, with or without
12
// modification, are permitted provided that the following conditions are
13
// met:
14
//
15
// 1. Redistributions of source code must retain the above copyright
16
// notice, this list of conditions and the following disclaimer.
17
//
18
// 2. Redistributions in binary form must reproduce the above copyright
19
// notice, this list of conditions and the following disclaimer in the
20
// documentation and/or other materials provided with the distribution.
21
//
22
// 3. Neither the name of the Corporation nor the names of the
23
// contributors may be used to endorse or promote products derived from
24
// this software without specific prior written permission.
25
//
26
// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27
// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29
// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30
// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31
// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33
// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34
// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35
// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36
// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37
//
38
// Questions? Contact
39
// Jonathan Hu (jhu@sandia.gov)
40
// Andrey Prokopenko (aprokop@sandia.gov)
41
// Ray Tuminaro (rstumin@sandia.gov)
42
//
43
// ***********************************************************************
44
//
45
// @HEADER
46
#ifndef XPETRA_MULTIVECTOR_DEF_HPP
47
#define XPETRA_MULTIVECTOR_DEF_HPP
48
49
#include "
Xpetra_MultiVector_decl.hpp
"
50
51
52
namespace
Xpetra
{
53
54
55
template
<
class
Scalar,
class
LocalOrdinal,
class
GlobalOrdinal,
class
Node>
56
MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
57
~MultiVector
()
58
{
59
}
60
61
62
template
<
class
Scalar,
class
LocalOrdinal,
class
GlobalOrdinal,
class
Node>
63
MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>
&
64
MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
65
operator=
(
const
MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>
& rhs)
66
{
67
assign(rhs);
// dispatch to protected virtual method
68
return
*
this
;
69
}
70
71
72
template
<
class
Scalar,
class
LocalOrdinal,
class
GlobalOrdinal,
class
Node>
73
void
74
MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
75
Xpetra_randomize
()
76
{
77
typedef
Teuchos::ScalarTraits<Scalar>
SCT;
78
79
const
size_t
numVectors = getNumVectors();
80
for
(
size_t
i = 0; i < numVectors; i++)
81
{
82
Teuchos::ArrayRCP<Scalar>
datai = getDataNonConst(i);
83
84
const
size_t
myLength = getLocalLength();
85
for
(
size_t
j = 0; j < myLength; j++)
86
{
87
datai[ j ] = SCT::random();
88
}
89
}
90
}
91
92
template
<
class
Scalar,
class
LocalOrdinal,
class
GlobalOrdinal,
class
Node>
93
void
94
MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
95
Xpetra_randomize
(
const
Scalar& minVal,
const
Scalar& maxVal)
96
{
97
typedef
Teuchos::ScalarTraits<Scalar>
SCT;
98
Scalar point5 = SCT::one()/ (SCT::one() + SCT::one());
99
100
const
size_t
numVectors = getNumVectors();
101
for
(
size_t
i = 0; i < numVectors; i++)
102
{
103
Teuchos::ArrayRCP<Scalar>
datai = getDataNonConst(i);
104
105
const
size_t
myLength = getLocalLength();
106
for
(
size_t
j = 0; j < myLength; j++)
107
{
108
datai[ j ] = point5*(maxVal-minVal)*SCT::random()+point5*(maxVal+minVal);
109
}
110
}
111
}
112
113
114
}
// namespace Xpetra
115
116
#endif
// XPETRA_MULTIVECTOR_DEF_HPP
117
Xpetra_MultiVector_decl.hpp
Teuchos::ArrayRCP
Xpetra::MultiVector
Definition
Xpetra_MultiVector_decl.hpp:80
Xpetra::MultiVector::operator=
MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > & operator=(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &rhs)
Assignment operator: Does a deep copy.
Definition
Xpetra_MultiVector_def.hpp:65
Xpetra::MultiVector::Xpetra_randomize
virtual void Xpetra_randomize()
Set multi-vector values to random numbers. XPetra implementation.
Definition
Xpetra_MultiVector_def.hpp:75
Xpetra::MultiVector::~MultiVector
virtual ~MultiVector()
Destructor.
Definition
Xpetra_MultiVector_def.hpp:57
Xpetra
Xpetra namespace
Definition
Xpetra_BlockedCrsMatrix.hpp:89
Teuchos::ScalarTraits
src
MultiVector
Xpetra_MultiVector_def.hpp
Generated by
1.10.0