IMP  2.2.1
The Integrative Modeling Platform
score_functor/OpenCubicSpline.h
Go to the documentation of this file.
1 /**
2  * \file IMP/score_functor/OpenCubicSpline.h
3  * \brief A Score on the distance between a pair of particles.
4  *
5  * Copyright 2007-2014 IMP Inventors. All rights reserved.
6  */
7 
8 #ifndef IMPSCORE_FUNCTOR_OPEN_CUBIC_SPLINE_H
9 #define IMPSCORE_FUNCTOR_OPEN_CUBIC_SPLINE_H
10 
11 #include <IMP/score_functor/score_functor_config.h>
12 #include "Score.h"
13 #include "internal/RawOpenCubicSpline.h"
14 IMPSCOREFUNCTOR_BEGIN_NAMESPACE
15 
16 //! Open cubic spline function.
17 /** This function interpolates a set of evenly spaced values using a spline.
18  The the second derivative at the termini is set to zero. See
19  core::ClosedCubicSpline for a periodic version.
20 
21  \see core::ClosedCubicSpline
22  */
23 class OpenCubicSpline : public Score {
24  double spacing_;
25  double inverse_spacing_;
26  internal::RawOpenCubicSpline spline_;
27  double minrange_;
28  double maxrange_;
29  bool extend_;
30 
31  public:
32  //! Constructor.
33  /** \param[in] values Score value at each spline point
34  \param[in] minrange Feature value at first spline point
35  \param[in] spacing Distance (in feature space) between points
36  \param[in] extend If extend is true, the nearest value is used
37  as the value for any feature outside of the range.
38  */
39  OpenCubicSpline(const Floats &values, Float minrange, Float spacing,
40  bool extend = false)
41  : spacing_(spacing),
42  inverse_spacing_(1.0 / spacing_),
43  spline_(values, spacing_, inverse_spacing_),
44  minrange_(minrange),
45  maxrange_(minrange_ + spacing_ * (values.size() - 1)),
46  extend_(extend) {
47  IMP_USAGE_CHECK(spacing > 0,
48  "The spacing between values must be positive.");
49  IMP_USAGE_CHECK(values.size() >= 1, "You must provide at least one value.");
50  }
51  // depend on get_is_trivially_zero
52  template <unsigned int D>
53  double get_score(kernel::Model *,
55  double distance) const {
56  // check for feature in range
57  if (distance < minrange_ || distance > maxrange_) {
58  if (extend_) {
59  if (distance < minrange_)
60  return spline_.get_first();
61  else
62  return spline_.get_last();
63  } else {
64  IMP_THROW("Spline out of domain", ModelException);
65  }
66  }
67  return spline_.evaluate(distance - minrange_, spacing_, inverse_spacing_);
68  }
69  template <unsigned int D>
70  DerivativePair get_score_and_derivative(
71  kernel::Model *, const base::Array<D, kernel::ParticleIndex> &,
72  double distance) const {
73  // check for distance in range
74  if (distance < minrange_ || distance > maxrange_) {
75  if (extend_) {
76  if (distance < minrange_)
77  return std::make_pair(spline_.get_first(), 0.0);
78  else
79  return std::make_pair(spline_.get_last(), 0.0);
80  } else {
81  IMP_THROW("Spline out of domain", ModelException);
82  }
83  }
84  return spline_.evaluate_with_derivative(distance - minrange_, spacing_,
85  inverse_spacing_);
86  }
87  template <unsigned int D>
88  double get_maximum_range(
89  kernel::Model *, const base::Array<D, kernel::ParticleIndex> &) const {
90  if (!extend_ || spline_.get_last() == 0)
91  return maxrange_;
92  else
93  return std::numeric_limits<double>::max();
94  }
95  template <unsigned int D>
96  bool get_is_trivially_zero(kernel::Model *,
97  const base::Array<D, kernel::ParticleIndex> &,
98  double squared_distance) const {
99  if (!extend_ || spline_.get_last() == 0) {
100  return squared_distance > algebra::get_squared(maxrange_);
101  } else
102  return false;
103  }
104 };
105 
106 IMPSCOREFUNCTOR_END_NAMESPACE
107 
108 #endif /* IMPSCORE_FUNCTOR_OPEN_CUBIC_SPLINE_H */
#define IMP_USAGE_CHECK(expr, message)
A runtime test for incorrect usage of a class or method.
A class to store an fixed array of same-typed values.
Definition: base/Array.h:33
IMP::kernel::Model Model
OpenCubicSpline(const Floats &values, Float minrange, Float spacing, bool extend=false)
Constructor.
A Score on the distance between a pair of particles.
std::pair< double, double > DerivativePair
A pair representing a function value with its first derivative.
Definition: base/types.h:23
#define IMP_THROW(message, exception_name)
Throw an exception with a message.
double Float
Basic floating-point value (could be float, double...)
Definition: base/types.h:20
Class for storing model, its restraints, constraints, and particles.
Definition: kernel/Model.h:72