From 609fb52b670d8ed74584a988b8c63da82d8d523b Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Tue, 10 Jul 2018 12:36:54 -0400
Subject: added wolfram link code

---
 wolfram_link/convexminorant.tm | 33 +++++++++++++++++++++++++++++++++
 1 file changed, 33 insertions(+)
 create mode 100644 wolfram_link/convexminorant.tm

(limited to 'wolfram_link/convexminorant.tm')

diff --git a/wolfram_link/convexminorant.tm b/wolfram_link/convexminorant.tm
new file mode 100644
index 0000000..48f5650
--- /dev/null
+++ b/wolfram_link/convexminorant.tm
@@ -0,0 +1,33 @@
+
+#include <convex.h>
+#include <wstp.h>
+
+extern int WSMain(int, char **);
+
+void convexminorant(double *, int);
+
+:Begin:
+:Function:      convexminorant
+:Pattern:       GetConvexMinorant[ list:{___Real} ]
+:Arguments:     { list }
+:ArgumentTypes: { Real64List }
+:ReturnType:    Manual
+:End:
+
+:Evaluate:      GetConvexMinorant[ sequence___Float]:= GetConvexMinorant[ {sequence} ]
+
+void convexminorant(double * Gammas, int len) {
+  int i;
+  for (i = 0; i < len; i++) {
+    if (Gammas[i] <= 0) {
+       break;
+    }
+  }
+  double *m = get_convex_minorant(i, Gammas);
+  WSPutReal64List(stdlink, m, i);
+}
+
+int main(int argc, char **argv) {
+  return WSMain(argc, argv);
+}
+
-- 
cgit v1.2.3-70-g09d2