#ifndef _utl_Moyal_h_
#define _utl_Moyal_h_

// $Id$

#include <cmath>


namespace utl {

  /** \brief Moyal function
    \image html utl_Moyal.png

    Moyal function is defined as
    \f[
      {\rm M}(x) = \exp(-{\textstyle\frac12}(x+{\rm e}^{-x}))
    \f]
    \ingroup math
  */

  inline double Moyal(const double x)
  { return std::exp(-0.5*(x + std::exp(-x))); }


  /** \brief Inverse of the Moyal function
    \image html utl_InverseMoyal.png
    \author Darko Veberic
    \ingroup math
  */

  template<int branch>
  double InverseMoyal(const double y, const double eps = 1e-6);

}


#endif

// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// End: