8 #ifndef IMPSCORE_FUNCTOR_OPEN_CUBIC_SPLINE_H
9 #define IMPSCORE_FUNCTOR_OPEN_CUBIC_SPLINE_H
11 #include <IMP/score_functor/score_functor_config.h>
13 #include "internal/RawOpenCubicSpline.h"
14 IMPSCOREFUNCTOR_BEGIN_NAMESPACE
25 double inverse_spacing_;
26 internal::RawOpenCubicSpline spline_;
42 inverse_spacing_(1.0 / spacing_),
43 spline_(values, spacing_, inverse_spacing_),
45 maxrange_(minrange_ + spacing_ * (values.size() - 1)),
48 "The spacing between values must be positive.");
49 IMP_USAGE_CHECK(values.size() >= 1,
"You must provide at least one value.");
52 template <
unsigned int D>
53 double get_score(
Model *,
55 double distance)
const {
57 if (distance < minrange_ || distance > maxrange_) {
59 if (distance < minrange_)
60 return spline_.get_first();
62 return spline_.get_last();
67 return spline_.evaluate(distance - minrange_, spacing_, inverse_spacing_);
69 template <
unsigned int D>
71 Model *,
const Array<D, ParticleIndex> &,
72 double distance)
const {
74 if (distance < minrange_ || distance > maxrange_) {
76 if (distance < minrange_)
77 return std::make_pair(spline_.get_first(), 0.0);
79 return std::make_pair(spline_.get_last(), 0.0);
81 IMP_THROW(
"Spline out of domain", ModelException);
84 return spline_.evaluate_with_derivative(distance - minrange_, spacing_,
87 template <
unsigned int D>
88 double get_maximum_range(
89 Model *,
const Array<D, ParticleIndex> &)
const {
90 if (!extend_ || spline_.get_last() == 0)
93 return std::numeric_limits<double>::max();
95 template <
unsigned int D>
96 bool get_is_trivially_zero(Model *,
97 const Array<D, ParticleIndex> &,
98 double squared_distance)
const {
99 if (!extend_ || spline_.get_last() == 0) {
100 return squared_distance > algebra::get_squared(maxrange_);
106 IMPSCOREFUNCTOR_END_NAMESPACE
A class to store an fixed array of same-typed values.
Class for storing model, its restraints, constraints, and particles.
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.
Open cubic spline function.
#define IMP_THROW(message, exception_name)
Throw an exception with a message.
double Float
Basic floating-point value (could be float, double...)
#define IMP_USAGE_CHECK(expr, message)
A runtime test for incorrect usage of a class or method.
A functor for computing a distance based score for D particles.
An exception which is thrown when the Model has attributes with invalid values.