Sacado_ScalarParameterLibrary.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//                           Sacado Package
//
// Copyright 2006 NTESS and the Sacado contributors.
// SPDX-License-Identifier: LGPL-2.1-or-later
// *****************************************************************************
// @HEADER

#ifndef SACADO_SCALARPARAMETERLIBRARY_HPP
#define SACADO_SCALARPARAMETERLIBRARY_HPP

#include "Sacado_ParameterLibraryBase.hpp"
#include "Sacado_ScalarParameterFamily.hpp"
#include "Sacado_ScalarParameterVector.hpp"

#include "Teuchos_Assert.hpp"

namespace Sacado {

  template <typename EvalTypeTraits = DefaultEvalTypeTraits>
  class ScalarParameterLibrary :
    public ParameterLibraryBase<ScalarParameterFamily<EvalTypeTraits>,
                                ScalarParameterEntry<_,EvalTypeTraits> > {

  public:

    typedef ParameterLibraryBase<ScalarParameterFamily<EvalTypeTraits>,
                                 ScalarParameterEntry<_,EvalTypeTraits> >
    BaseT;

    ScalarParameterLibrary() {}

    virtual ~ScalarParameterLibrary() {}

    void setRealValueForAllTypes(const std::string& name, double value);

    template <class EvalType>
    void
    setRealValue(const std::string& name, double value);

    template <class EvalType>
    void
    setValue(const std::string& name,
             const typename EvalTypeTraits::template apply<EvalType>::type& value);

    template <class EvalType>
    double
    getRealValue(const std::string& name) const;

    template <class EvalType>
    const typename EvalTypeTraits::template apply<EvalType>::type&
    getValue(const std::string& name) const;

    static ScalarParameterLibrary& getInstance() {
      static ScalarParameterLibrary instance;
      return instance;
    }


    template <class EvalType>
    void
    fillVector(const Teuchos::Array<std::string>& names,
               ScalarParameterVector<EvalTypeTraits>& pv);

  private:

    ScalarParameterLibrary(const ScalarParameterLibrary&);

    ScalarParameterLibrary& operator = (const ScalarParameterLibrary&);

  };

}


template <typename EvalTypeTraits>
void
Sacado::ScalarParameterLibrary<EvalTypeTraits>::
setRealValueForAllTypes(const std::string& name, double value)
{
  typename BaseT::FamilyMap::iterator it = this->library.find(name);
  TEUCHOS_TEST_FOR_EXCEPTION(
     it == this->library.end(),
     std::logic_error,
     std::string("Sacado::ScalararameterLibrary::setRealValueForAllTypes():  ")
     + "Invalid parameter family " + name);
  (*it).second->setRealValueForAllTypes(value);
}

template <typename EvalTypeTraits>
template <class EvalType>
void
Sacado::ScalarParameterLibrary<EvalTypeTraits>::
setRealValue(const std::string& name, double value)
{
  typename BaseT::FamilyMap::iterator it = this->library.find(name);
  TEUCHOS_TEST_FOR_EXCEPTION(
     it == this->library.end(),
     std::logic_error,
     std::string("Sacado::ScalarParameterLibrary::setValueAsConstant():  ")
     + "Invalid parameter family " + name);
  (*it).second-> template setRealValue<EvalType>(value);
}

template <typename EvalTypeTraits>
template <class EvalType>
void
Sacado::ScalarParameterLibrary<EvalTypeTraits>::
setValue(
      const std::string& name,
      const typename EvalTypeTraits::template apply<EvalType>::type& value)
{
  typename BaseT::FamilyMap::iterator it = this->library.find(name);
  TEUCHOS_TEST_FOR_EXCEPTION(
      it == this->library.end(),
      std::logic_error,
      std::string("Sacado::ScalarParameterLibrary::setValueAsIndependent():  ")
      + "Invalid parameter family " + name);
  (*it).second-> template setValue<EvalType>(value);
}

template <typename EvalTypeTraits>
template <class EvalType>
double
Sacado::ScalarParameterLibrary<EvalTypeTraits>::
getRealValue(const std::string& name) const
{
  typename BaseT::FamilyMap::const_iterator it = this->library.find(name);
  TEUCHOS_TEST_FOR_EXCEPTION(
                 it == this->library.end(),
                 std::logic_error,
                 std::string("Sacado::ScalarParameterLibrary::getValue():  ")
                 + "Invalid parameter family " + name);
  return (*it).second-> template getRealValue<EvalType>();
}

template <typename EvalTypeTraits>
template <class EvalType>
const typename EvalTypeTraits::template apply<EvalType>::type&
Sacado::ScalarParameterLibrary<EvalTypeTraits>::
getValue(const std::string& name) const
{
  typename BaseT::FamilyMap::const_iterator it = this->library.find(name);
  TEUCHOS_TEST_FOR_EXCEPTION(
                 it == this->library.end(),
                 std::logic_error,
                 std::string("Sacado::ScalarParameterLibrary::getValue():  ")
                 + "Invalid parameter family " + name);
  return (*it).second->template getValue<EvalType>();
}

template <typename EvalTypeTraits>
template <class EvalType>
void
Sacado::ScalarParameterLibrary<EvalTypeTraits>::
fillVector(const Teuchos::Array<std::string>& names,
           Sacado::ScalarParameterVector<EvalTypeTraits>& pv)
{
  typename BaseT::FamilyMap::iterator it;

  // Fill in parameters
  for (unsigned int i=0; i<names.size(); i++) {
    it = this->library.find(names[i]);
    TEUCHOS_TEST_FOR_EXCEPTION(
                   it == this->library.end(),
                   std::logic_error,
                   std::string("Sacado::ParameterLibraryBase::fillVector():  ")
                   + "Invalid parameter family " + names[i]);
    pv.addParam((*it).second, 0.0);
    pv[i].baseValue = (*it).second->template getRealValue<EvalType>();
  }
}


#endif