ParDomain.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 +
//----------------------------------------------------------------------
// Includes
//----------------------------------------------------------------------
#include "MantidCurveFitting/ParDomain.h"
#include "MantidKernel/MultiThreaded.h"
 
namespace Mantid::CurveFitting {
 
void ParDomain::getDomainAndValues(size_t i, API::FunctionDomain_sptr &domain, API::FunctionValues_sptr &values) const {
  if (i >= m_creators.size())
    throw std::range_error("Function domain index is out of range.");
  if (!m_domain[i]) {
    m_creators[i]->createDomain(m_domain[i], m_values[i]);
  }
  domain = m_domain[i];
  values = m_values[i];
}
 
void ParDomain::additiveCostFunctionVal(const CostFunctions::CostFuncFitting &costFunction) {
  const int n = static_cast<int>(getNDomains());
  PARALLEL_FOR_NO_WSP_CHECK()
  for (int i = 0; i < n; ++i) {
    API::FunctionDomain_sptr domain;
    API::FunctionValues_sptr values;
    getDomainAndValues(static_cast<size_t>(i), domain, values);
    if (!values) {
      throw std::runtime_error("CostFunction: undefined FunctionValues.");
    }
    costFunction.addVal(domain, values);
  }
}
 
void ParDomain::additiveCostFunctionValDerivHessian(const CostFunctions::CostFuncFitting &costFunction, bool evalDeriv,
                                                    bool evalHessian) {
  const auto n = static_cast<int>(getNDomains());
  PARALLEL_SET_DYNAMIC(0);
  std::vector<API::IFunction_sptr> funs;
  PARALLEL_FOR_NO_WSP_CHECK()
  for (int i = 0; i < n; ++i) {
    API::FunctionDomain_sptr domain;
    API::FunctionValues_sptr values;
    getDomainAndValues(i, domain, values);
    auto simpleValues = std::dynamic_pointer_cast<API::FunctionValues>(values);
    if (!simpleValues) {
      throw std::runtime_error("CostFunction: undefined FunctionValues.");
    }
    std::vector<API::IFunction_sptr>::size_type k = PARALLEL_THREAD_NUMBER;
    PARALLEL_CRITICAL(resize) {
      if (k >= funs.size()) { // cppcheck-suppress unsignedLessThanZero
        funs.resize(k + 1);
      }
      if (!funs[k]) {
        funs[k] = costFunction.getFittingFunction()->clone();
      }
    }
    costFunction.addValDerivHessian(funs[k], domain, simpleValues, evalDeriv, evalHessian);
  }
}
 
} // namespace Mantid::CurveFitting