wright_omega.h 822 B

1234567891011121314151617181920212223242526
  1. /**
  2. * Santiago Akle
  3. * ICME Stanford University 2014
  4. *
  5. * Computes the value \omega(z) defined as the solution y to
  6. * the equation y+log(y) = z for z real and z>=1.
  7. * Follows the recommendations by
  8. * PIERS W. LAWRENCE, ROBERT M. CORLESS, and DAVID J. JEFFREY.
  9. * Published in:
  10. * Algorithm 917: Complex Double-Precision Evaluation of the Wright \omega Function
  11. * ACM Transactions on Mathematical Software (TOMS) TOMS Homepage table of contents archive
  12. * Volume 38 Issue 3, April 2012
  13. * Article No. 20
  14. * Publication Date 2012-04-01 (yyyy-mm-dd)
  15. * Publisher ACM New York, NY, USA
  16. * ISSN: 0098-3500 EISSN: 1557-7295 doi>10.1145/2168773.2168779
  17. */
  18. #include "glblopts.h"
  19. #if (defined _WIN32 || defined _WIN64 || defined _WINDLL )
  20. #define _USE_MATH_DEFINES
  21. #endif
  22. #include <math.h>
  23. pfloat wrightOmega(pfloat z);