summaryrefslogtreecommitdiff
path: root/thirdparty/linux/include/coin/wright_omega.h
diff options
context:
space:
mode:
Diffstat (limited to 'thirdparty/linux/include/coin/wright_omega.h')
-rw-r--r--thirdparty/linux/include/coin/wright_omega.h26
1 files changed, 26 insertions, 0 deletions
diff --git a/thirdparty/linux/include/coin/wright_omega.h b/thirdparty/linux/include/coin/wright_omega.h
new file mode 100644
index 0000000..dd50a77
--- /dev/null
+++ b/thirdparty/linux/include/coin/wright_omega.h
@@ -0,0 +1,26 @@
+/**
+ * Santiago Akle
+ * ICME Stanford University 2014
+ *
+ * Computes the value \omega(z) defined as the solution y to
+ * the equation y+log(y) = z for z real and z>=1.
+ * Follows the recommendations by
+ * PIERS W. LAWRENCE, ROBERT M. CORLESS, and DAVID J. JEFFREY.
+ * Published in:
+ * Algorithm 917: Complex Double-Precision Evaluation of the Wright \omega Function
+ * ACM Transactions on Mathematical Software (TOMS) TOMS Homepage table of contents archive
+ * Volume 38 Issue 3, April 2012
+ * Article No. 20
+ * Publication Date 2012-04-01 (yyyy-mm-dd)
+ * Publisher ACM New York, NY, USA
+ * ISSN: 0098-3500 EISSN: 1557-7295 doi>10.1145/2168773.2168779
+ */
+
+#include "glblopts.h"
+
+#if (defined _WIN32 || defined _WIN64 || defined _WINDLL )
+#define _USE_MATH_DEFINES
+#endif
+#include <math.h>
+
+pfloat wrightOmega(pfloat z);