Stokhos Package Browser (Single Doxygen Collection) Version of the Day
Loading...
Searching...
No Matches
Tpetra_Details_fill_UQ_PCE.hpp
Go to the documentation of this file.
1// @HEADER
2// ***********************************************************************
3//
4// Tpetra: Templated Linear Algebra Services Package
5// Copyright (2008) 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#ifndef TPETRA_DETAILS_FILL_UQ_PCE_HPP
43#define TPETRA_DETAILS_FILL_UQ_PCE_HPP
44
45#include "Tpetra_Details_fill.hpp"
47
48namespace Tpetra {
49namespace Details {
50namespace Blas {
51
52template<class DT, class ... DP,
53 class ValueType,
54 class IndexType,
55 class ExecutionSpace>
56typename std::enable_if<
57 Kokkos::is_view_uq_pce< Kokkos::View<DT,DP...> >::value >::type
58fill (const ExecutionSpace& execSpace,
59 const Kokkos::View<DT,DP...>& X,
60 const ValueType& alpha,
61 const IndexType numRows,
62 const IndexType numCols)
63{
64 static_assert (std::is_integral<IndexType>::value,
65 "IndexType must be a built-in integer type.");
66 Kokkos::deep_copy(execSpace, X, alpha);
67}
68
69template<class DT, class ... DP,
70 class ValueType,
71 class IndexType,
72 class ExecutionSpace>
73typename std::enable_if<
74 Kokkos::is_view_uq_pce< Kokkos::View<DT,DP...> >::value >::type
75fill (const ExecutionSpace& execSpace,
76 const Kokkos::View<DT,DP...>& X,
77 const ValueType& alpha,
78 const IndexType numRows,
79 const IndexType numCols,
80 const size_t whichVectors[])
81{
82 typedef Kokkos::View<DT,DP...> ViewType;
83 static_assert (ViewType::Rank == 2, "ViewType must be a rank-2 "
84 "Kokkos::View in order to call the \"whichVectors\" "
85 "specialization of fill.");
86 static_assert (std::is_integral<IndexType>::value,
87 "IndexType must be a built-in integer type.");
88 for (IndexType k = 0; k < numCols; ++k) {
89 const IndexType j = whichVectors[k];
90 auto X_j = Kokkos::subview (X, Kokkos::ALL (), j);
91 Kokkos::deep_copy(execSpace, X_j, alpha);
92 }
93}
94
95} // namespace Blas
96} // namespace Details
97} // namespace Tpetra
98
99#endif // TPETRA_DETAILS_FILL_UQ_PCE_HPP
std::enable_if< Kokkos::is_view_uq_pce< Kokkos::View< DT, DP... > >::value >::type fill(const ExecutionSpace &execSpace, const Kokkos::View< DT, DP... > &X, const ValueType &alpha, const IndexType numRows, const IndexType numCols)