CubicSpline.h

Go to the documentation of this file.

// Mantid Repository : https://github.com/mantidproject/mantid
//
// Copyright © 2007 ISIS Rutherford Appleton Laboratory UKRI,
//   NScD Oak Ridge National Laboratory, European Spallation Source,
//   Institut Laue - Langevin & CSNS, Institute of High Energy Physics, CAS
// SPDX - License - Identifier: GPL - 3.0 +
#pragma once
 
//----------------------------------------------------------------------
// Includes
//----------------------------------------------------------------------
#include "MantidCurveFitting/Functions/BackgroundFunction.h"
#include "MantidKernel/GSL_Helpers.h"
 
#include <valarray>
 
namespace Mantid {
namespace CurveFitting {
namespace Functions {
class MANTID_CURVEFITTING_DLL CubicSpline : public BackgroundFunction {
 
public:
  CubicSpline();
 
  std::string name() const override { return "CubicSpline"; }
  void function1D(double *out, const double *xValues, const size_t nData) const override;
  void derivative1D(double *out, const double *xValues, size_t nData, const size_t order) const override;
  void setParameter(size_t i, const double &value, bool explicitlySet = true) override;
  using ParamFunction::setParameter;
 
  void setAttribute(const std::string &attName, const Attribute &) override;
 
  void setXAttribute(const size_t index, double x);
 
private:
  const int m_min_points;
 
  Mantid::Kernel::spline::accel_uptr m_acc;
 
  Mantid::Kernel::spline::spline_uptr m_spline;
 
  mutable bool m_recalculateSpline;
 
  void reallocGSLObjects(const int n);
 
  void setupInput(int const n) const;
 
  void calculateSpline(double *out, const double *xValues, const size_t nData) const;
 
  void calculateDerivative(double *out, const double *xValues, const size_t nData, const size_t order) const;
 
  void checkGSLError(const int status, const int errorType) const;
 
  bool checkXInRange(double x) const;
 
  double splineEval(const double x) const;
};
 
using CubicSpline_sptr = std::shared_ptr<CubicSpline>;
using CubicSpline_const_sptr = const std::shared_ptr<CubicSpline>;
 
} // namespace Functions
} // namespace CurveFitting
} // namespace Mantid