mfexpcon.cpp

Go to the documentation of this file.

/*
\file Robust exponential functions for constant scalar objects.

Author: David Fournier
Copyright (c) 2008-2015 Regents of the University of California
*/
#include <fvar.hpp>

double mfexp(double x, double b)
{
  if (x > b)
  {
    return exp(b)*(1.+2.*(x-b))/(1.+x-b);
  }
  else if (x < -b)
  {
    return exp(-b)*(1.-x-b)/(1.+2.*(-x-b));
  }
  return exp(x);
}
double mfexp(double x)
{
  double b = 60;
  return mfexp(x, b);
}