HcalCubicInterpolator Class Reference

#include <HcalCubicInterpolator.h>

Inheritance diagram for HcalCubicInterpolator:

Public Types
typedef std::tuple< double, double, double >	Triple

Public Member Functions
HcalCubicInterpolator	approximateInverse () const
	HcalCubicInterpolator ()
	HcalCubicInterpolator (const std::vector< Triple > &points)
double	operator() (double x) const override
double	xmax () const override
double	xmin () const override
	~HcalCubicInterpolator () override
Public Member Functions inherited from AbsHcalFunctor
bool	operator!= (const AbsHcalFunctor &r) const
bool	operator== (const AbsHcalFunctor &r) const
virtual	~AbsHcalFunctor ()

Protected Member Functions
bool	isEqual (const AbsHcalFunctor &other) const override

Private Member Functions
template<class Archive >
void	serialize (Archive &ar, unsigned)

Private Attributes
std::vector< double >	abscissae_
std::vector< double >	derivatives_
std::vector< double >	values_

Friends
class	boost::serialization::access

Additional Inherited Members
Static Protected Member Functions inherited from AbsHcalFunctor
template<class Iter >
static bool	isStrictlyDecreasing (Iter begin, Iter const end)
template<class Iter >
static bool	isStrictlyIncreasing (Iter begin, Iter const end)

Detailed Description

Definition at line 17 of file HcalCubicInterpolator.h.

Member Typedef Documentation

typedef std::tuple<double, double, double> HcalCubicInterpolator::Triple

Definition at line 20 of file HcalCubicInterpolator.h.

Constructor & Destructor Documentation

HcalCubicInterpolator::HcalCubicInterpolator

(

)

Definition at line 6 of file HcalCubicInterpolator.cc.

Referenced by approximateInverse().

{}

HcalCubicInterpolator::HcalCubicInterpolator ( const std::vector< Triple > &  points )

explicit

Definition at line 8 of file HcalCubicInterpolator.cc.

References abscissae_, derivatives_, Exception, mps_fire::i, submitPVValidationJobs::t, createJobs::tmp, and values_.

                                                                            {
  const std::size_t sz = points.size();
  if (sz) {
    std::vector<Triple> tmp(points);
    std::sort(tmp.begin(), tmp.end());
    abscissae_.reserve(sz);
    values_.reserve(sz);
    derivatives_.reserve(sz);
    for (std::size_t i = 0; i < sz; ++i) {
      const Triple& t(tmp[i]);
      abscissae_.push_back(std::get<0>(t));
      values_.push_back(std::get<1>(t));
      derivatives_.push_back(std::get<2>(t));
    }
    const std::size_t szm1 = sz - 1;
    for (std::size_t i = 0; i < szm1; ++i)
      if (abscissae_[i] == abscissae_[i + 1])
        throw cms::Exception(
            "In HcalCubicInterpolator constructor:"
            " abscissae must not coincide");
  }
}

HcalCubicInterpolator::~HcalCubicInterpolator ( )

inlineoverride

Definition at line 30 of file HcalCubicInterpolator.h.

{}

Member Function Documentation

HcalCubicInterpolator HcalCubicInterpolator::approximateInverse

(

)

const

Definition at line 74 of file HcalCubicInterpolator.cc.

References abscissae_, derivatives_, Exception, HcalCubicInterpolator(), mps_fire::i, AbsHcalFunctor::isStrictlyDecreasing(), AbsHcalFunctor::isStrictlyIncreasing(), and values_.

                                                                      {
  const bool monotonous =
      isStrictlyIncreasing(values_.begin(), values_.end()) || isStrictlyDecreasing(values_.begin(), values_.end());
  if (!monotonous)
    throw cms::Exception(
        "In HcalCubicInterpolator::inverse:"
        " can't invert non-monotonous functor");
  const std::size_t sz = abscissae_.size();
  std::vector<Triple> points;
  points.reserve(sz);
  for (std::size_t i = 0; i < sz; ++i) {
    const double dydx = derivatives_[i];
    if (dydx == 0.0)
      throw cms::Exception(
          "In HcalCubicInterpolator::inverse:"
          " can't invert functor with derivatives of 0");
    points.push_back(Triple(values_[i], abscissae_[i], 1.0 / dydx));
  }
  return HcalCubicInterpolator(points);
}

bool HcalCubicInterpolator::isEqual ( const AbsHcalFunctor &  other ) const

inlineoverrideprotectedvirtual

Implements AbsHcalFunctor.

Definition at line 41 of file HcalCubicInterpolator.h.

References abscissae_, derivatives_, alignCSCRings::r, and values_.

                                                                  {
    const HcalCubicInterpolator& r = static_cast<const HcalCubicInterpolator&>(other);
    return abscissae_ == r.abscissae_ && values_ == r.values_ && derivatives_ == r.derivatives_;
  }

double HcalCubicInterpolator::operator() ( double  x ) const

overridevirtual

Implements AbsHcalFunctor.

Definition at line 31 of file HcalCubicInterpolator.cc.

References abscissae_, derivatives_, PVValHelper::dx, mps_fire::result, submitPVValidationJobs::t, cuda_std::upper_bound(), values_, and x.

                                                             {
  double result = 0.0;
  const std::size_t sz = abscissae_.size();
  if (sz) {
    if (sz > 1) {
      const std::size_t szm1 = sz - 1;
      if (x >= abscissae_[szm1])
        result = values_[szm1] + derivatives_[szm1] * (x - abscissae_[szm1]);
      else if (x <= abscissae_[0])
        result = values_[0] + derivatives_[0] * (x - abscissae_[0]);
      else {
        const std::size_t cell = std::upper_bound(abscissae_.begin(), abscissae_.end(), x) - abscissae_.begin() - 1;
        const std::size_t cellp1 = cell + 1;
        const double dx = abscissae_[cellp1] - abscissae_[cell];
        const double t = (x - abscissae_[cell]) / dx;
        const double onemt = 1.0 - t;
        const double h00 = onemt * onemt * (1.0 + 2.0 * t);
        const double h10 = onemt * onemt * t;
        const double h01 = t * t * (3.0 - 2.0 * t);
        const double h11 = t * t * onemt;
        result = h00 * values_[cell] + h10 * dx * derivatives_[cell] + h01 * values_[cellp1] -
                 h11 * dx * derivatives_[cellp1];
      }
    } else
      result = values_[0] + derivatives_[0] * (x - abscissae_[0]);
  }
  return result;
}

template<class Archive >

void HcalCubicInterpolator::serialize ( Archive &  ar, unsigned    )

inlineprivate

Definition at line 54 of file HcalCubicInterpolator.h.

References abscissae_, derivatives_, and values_.

                                                 {
    boost::serialization::base_object<AbsHcalFunctor>(*this);
    ar& abscissae_& values_& derivatives_;
  }

double HcalCubicInterpolator::xmax ( ) const

overridevirtual

Reimplemented from AbsHcalFunctor.

Definition at line 67 of file HcalCubicInterpolator.cc.

References abscissae_, and mps_fire::result.

Referenced by svgfig.XAxis::__repr__().

                                         {
  double result = 0.0;
  if (!abscissae_.empty())
    result = abscissae_.back();
  return result;
}

double HcalCubicInterpolator::xmin ( ) const

overridevirtual

Reimplemented from AbsHcalFunctor.

Definition at line 60 of file HcalCubicInterpolator.cc.

References abscissae_, and mps_fire::result.

Referenced by svgfig.XAxis::__repr__(), svgfig.Axes::__repr__(), svgfig.HGrid::__repr__(), svgfig.Grid::__repr__(), and svgfig.Axes::SVG().

                                         {
  double result = 0.0;
  if (!abscissae_.empty())
    result = abscissae_[0];
  return result;
}

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 51 of file HcalCubicInterpolator.h.

Member Data Documentation

std::vector<double> HcalCubicInterpolator::abscissae_

private

Definition at line 47 of file HcalCubicInterpolator.h.

Referenced by approximateInverse(), HcalCubicInterpolator(), isEqual(), operator()(), serialize(), xmax(), and xmin().

std::vector<double> HcalCubicInterpolator::derivatives_

private

Definition at line 49 of file HcalCubicInterpolator.h.

Referenced by approximateInverse(), HcalCubicInterpolator(), isEqual(), operator()(), and serialize().

std::vector<double> HcalCubicInterpolator::values_

private

Definition at line 48 of file HcalCubicInterpolator.h.

Referenced by approximateInverse(), HcalCubicInterpolator(), isEqual(), operator()(), and serialize().