42 #ifndef STOKHOS_IFPACK2_UQ_PCE_HPP
43 #define STOKHOS_IFPACK2_UQ_PCE_HPP
57 template <
typename XV,
class SizeType>
61 template <
typename S,
typename ... P,
class SizeType>
63 Kokkos::
View< Sacado::UQ::PCE<S>*,P... >,
66 typedef Kokkos::View< Sacado::UQ::PCE<S>*,P... >
XVector;
71 typedef typename array_type::non_const_value_type
value_type;
72 typedef Kokkos::Details::ArithTraits<value_type>
KAT;
82 const typename XVector::non_const_value_type& min_val) :
85 m_min_val_mag(
KAT::
abs(m_min_val)),
89 KOKKOS_INLINE_FUNCTION
101 template<
class XV,
class SizeType>
104 template<
class S,
class ...P,
class SizeType>
106 Kokkos::
View< Sacado::UQ::PCE<S>*,P... >, SizeType > {
107 typedef Kokkos::View< Sacado::UQ::PCE<S>*,P... >
XV;
111 const typename XV::non_const_value_type& minVal)
115 "LocalReciprocalThreshold not implemented for non-constant minVal");
126 Kokkos::parallel_for( X.extent(0), op );
134 #endif // STOKHOS_IFPACK2_UQ_PCE_HPP
array_type::execution_space execution_space
array_type::non_const_value_type value_type
const value_type m_min_val_mag
Kokkos::DefaultExecutionSpace execution_space
Kokkos::Details::ArithTraits< value_type > KAT
Kokkos::View< Sacado::UQ::PCE< S > *, P... > XVector
KOKKOS_INLINE_FUNCTION void raise_error(const char *msg)
KOKKOS_INLINE_FUNCTION constexpr std::enable_if< is_view_uq_pce< View< T, P...> >::value, unsigned >::type dimension_scalar(const View< T, P...> &view)
Kokkos::View< Sacado::UQ::PCE< S > *, P... > XV
V_ReciprocalThresholdSelfFunctor(const XVector &x, const typename XVector::non_const_value_type &min_val)
const value_type m_min_val
XVector::array_type array_type
KOKKOS_INLINE_FUNCTION bool is_constant(const T &x)
KOKKOS_INLINE_FUNCTION PCE< Storage > abs(const PCE< Storage > &a)
expr1 expr1 expr1 expr2 expr1 expr1 c expr2 expr1 c fastAccessCoeff(j)-expr2.val(j)
static void compute(const XV &X, const typename XV::non_const_value_type &minVal)
KOKKOS_INLINE_FUNCTION void operator()(const size_type i) const