43#ifndef RTOPPACK_TOP_RANDOMIZE_HPP
44#define RTOPPACK_TOP_RANDOMIZE_HPP
46#include "RTOpPack_RTOpTHelpers.hpp"
75 const Scalar& l = -ScalarTraits<Scalar>::one(),
76 const Scalar& u = +ScalarTraits<Scalar>::one()
97 const Ptr<ReductTarget> &reduct_obj_inout
100 typedef typename Teuchos::ArrayRCP<Scalar>::iterator iter_t;
103 validate_apply_op<Scalar>(*
this, 0, 1,
false,
104 sub_vecs, targ_sub_vecs, reduct_obj_inout.getConst());
107 (void)reduct_obj_inout;
110 const index_type subDim = targ_sub_vecs[0].subDim();
111 const index_type globalOffset = targ_sub_vecs[0].globalOffset();
113 iter_t z0_val = targ_sub_vecs[0].values().begin();
114 const ptrdiff_t z0_s = targ_sub_vecs[0].stride();
117 const Scalar a = Scalar(0.5)*(
u_ -
l_);
118 const Scalar b = Scalar(0.5)*(
u_ +
l_);
119 for(
index_type i = 0; i < subDim; ++i, z0_val += z0_s )
121 Teuchos::ScalarTraits<Scalar>::seedrandom(
seed_+globalOffset+i);
122 *z0_val = a * Teuchos::ScalarTraits<Scalar>::random() + b;
135template<
class Scalar>
Class for a non-changeable sub-vector.
Templated interface to vector reduction/transformation operators {abstract}.
void apply_op(const ArrayView< const ConstSubVectorView< Scalar > > &sub_vecs, const ArrayView< const SubVectorView< Scalar > > &targ_sub_vecs, const Ptr< ReductTarget > &reduct_obj) const
Apply the reduction/transformation operator to a set of sub-vectors.
void setOpNameBase(const std::string &op_name_base)
Just set the operator name.
Class for a changeable sub-vector.
Generate a random vector in the range [l,u]: z0[i] = 0.5*((u-l)*Teuchos::ScalarTraits<Scalar>::random...
unsigned int get_seed() const
void set_bounds(const Scalar &l, const Scalar &u)
TOpRandomize(const Scalar &l=-ScalarTraits< Scalar >::one(), const Scalar &u=+ScalarTraits< Scalar >::one())
void set_seed(const unsigned int seed)
static unsigned int static_seed_
void apply_op_impl(const ArrayView< const ConstSubVectorView< Scalar > > &sub_vecs, const ArrayView< const SubVectorView< Scalar > > &targ_sub_vecs, const Ptr< ReductTarget > &reduct_obj_inout) const
static void set_static_seed(const unsigned int static_seed)
static unsigned int get_static_seed()
Teuchos_Ordinal index_type