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