DiffSphere.cpp

Go to the documentation of this file.

// Mantid Repository : https://github.com/mantidproject/mantid
//
// Copyright © 2018 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 +
#include "MantidCurveFitting/Functions/DiffSphere.h"
#include "MantidAPI/FunctionFactory.h"
#include "MantidAPI/IFunction.h"
#include "MantidAPI/ImmutableCompositeFunction.h"
#include "MantidCurveFitting/Functions/ElasticDiffSphere.h"
#include "MantidCurveFitting/Functions/InelasticDiffSphere.h"
 
#include <boost/math/special_functions/bessel.hpp>
 
#include <cmath>
#include <limits>
 
namespace {
Mantid::Kernel::Logger g_log("DiffSphere");
}
 
namespace Mantid::CurveFitting::Functions {
 
DECLARE_FUNCTION(DiffSphere)
 
 
void DiffSphere::trickleDownAttribute(const std::string &name) {
  for (size_t iFun = 0; iFun < this->nFunctions(); iFun++) {
    API::IFunction_sptr fun = this->getFunction(iFun);
    if (fun->hasAttribute(name)) {
      fun->setAttribute(name, this->getAttribute(name));
    }
  }
}
 
void DiffSphere::declareAttribute(const std::string &name, const API::IFunction::Attribute &defaultValue) {
  API::ImmutableCompositeFunction::declareAttribute(name, defaultValue);
  this->trickleDownAttribute(name);
}
 
void DiffSphere::setAttribute(const std::string &name, const API::IFunction::Attribute &att) {
  API::ImmutableCompositeFunction::setAttribute(name, att);
  this->trickleDownAttribute(name);
}
 
void DiffSphere::init() {
  m_elastic = std::dynamic_pointer_cast<ElasticDiffSphere>(
      API::FunctionFactory::Instance().createFunction("ElasticDiffSphere"));
  this->addFunction(m_elastic);
  m_inelastic = std::dynamic_pointer_cast<InelasticDiffSphere>(
      API::FunctionFactory::Instance().createFunction("InelasticDiffSphere"));
  this->addFunction(m_inelastic);
 
  this->setAttributeValue("NumDeriv", true);
  this->declareAttribute("Q", API::IFunction::Attribute(1.0));
 
  // Set the ties between Elastic and Inelastic parameters
  this->addDefaultTies("f0.Height=f1.Intensity,f0.Radius=f1.Radius,f0.Centre=0");
  this->applyTies();
}
 
} // namespace Mantid::CurveFitting::Functions