SigmoidConverter.cpp

Go to the documentation of this file.

#include <cmath>
#include <iostream>
 
#include <ElementsKernel/Logging.h>
 
#include "ModelFitting/Parameters/SigmoidConverter.h"
 
namespace ModelFitting {
 
using namespace std;
 
static Elements::Logging logger = Elements::Logging::getLogger("ModelFitting");
 
SigmoidConverter::~SigmoidConverter() = default;
 
double SigmoidConverter::worldToEngine(const double world_value) const {
  if (world_value < m_min_value || world_value > m_max_value) {
    logger.warn() << "WorldToEngine SigmoidConverter: world values outside of possible range";
  }
 
  double num = world_value - m_min_value;
  double den = m_max_value - world_value;
  return (num > 1e-50 ? (den > 1e-50 ? log(num/den) : 50.0) : -50.0);
}
 
double SigmoidConverter::engineToWorld(const double engine_value) const {
  auto clamped_value = std::max(-50.0, std::min(50.0, engine_value));
  return m_min_value + (m_max_value - m_min_value) / (1 + exp(-clamped_value));
}
 
double SigmoidConverter::getEngineToWorldDerivative(const double value) const {
  return (value - m_min_value) * (m_max_value - value) / (m_max_value - m_min_value);
}
 
 
}// namespace ModelFitting