ROL_RandomVector.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//               Rapid Optimization Library (ROL) Package
//
// Copyright 2014 NTESS and the ROL contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef ROL_RANDOMVECTOR_H
#define ROL_RANDOMVECTOR_H
 
#include "ROL_Vector.hpp"
#include "ROL_BoundConstraint.hpp"
#include "ROL_Elementwise_Function.hpp"
 
 
namespace ROL {
 
template<class Real> 
void RandomizeVector( Vector<Real> &x, const Real &lower=0.0, const Real &upper=1.0 ) {
 
  Elementwise::UniformlyRandom<Real> ur(lower,upper);
  x.applyUnary(ur);
}
 
template<class Real>
void RandomizeFeasibleVector( Vector<Real> &x, BoundConstraint<Real> &bnd ) {
  const ROL::Ptr<const Vector<Real> > u = bnd.getUpperBound();
  const ROL::Ptr<const Vector<Real> > l = bnd.getLowerBound();
 
  Elementwise::UniformlyRandomMultiply<Real> urm;
  
  x.set(*u);
  x.axpy(-1.0,*l);
  x.applyUnary(urm);
  x.plus(*l);
}
 
 
} // namespace ROL
 
#endif // ROL_RANDOMVECTOR_H