Thyra_ScaledModelEvaluator_decl.hpp

// @HEADER
// *****************************************************************************
//    Thyra: Interfaces and Support for Abstract Numerical Algorithms
//
// Copyright 2004 NTESS and the Thyra contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER

#ifndef THYRA_SCALED_MODEL_EVALUATOR_HPP
#define THYRA_SCALED_MODEL_EVALUATOR_HPP

#include "Thyra_ModelEvaluatorDelegatorBase.hpp"

namespace Thyra {


template<class Scalar>
class ScaledModelEvaluator : 
    virtual public ModelEvaluatorDelegatorBase<Scalar>
{
public:
  
  ScaledModelEvaluator();
  
  std::string description() const;

  void set_f_scaling(const RCP<const Thyra::VectorBase<Scalar> >& f_scaling);

private:


  void evalModelImpl(
    const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
    const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs
    ) const;

  
private:

  //* Diagonal scaling vector */
  RCP<const Thyra::VectorBase<Scalar> > f_scaling_;

};


template<class Scalar>
RCP<ScaledModelEvaluator<Scalar> >
createNonconstScaledModelEvaluator(const RCP<ModelEvaluator<Scalar > > &model)
{
  RCP<ScaledModelEvaluator<Scalar> > srme(new ScaledModelEvaluator<Scalar>);
  srme->initialize(model);
  return srme;
}


} // namespace Thyra


#endif // THYRA_SCALED_MODEL_EVALUATOR_HPP