added wolfram link code

2 files changed, 46 insertions, 0 deletions

diff --git a/wolfram_link/Makefile b/wolfram_link/Makefile
new file mode 100644
index 0000000..9d65623
--- /dev/null
+++ b/wolfram_link/Makefile
@@ -0,0 +1,13 @@
+
+CC = clang
+WSPREP = /opt/Mathematica/SystemFiles/Links/WSTP/DeveloperKit/Linux-x86-64/CompilerAdditions/wsprep
+CFLAGS = -g -Os -O3 -Wall -fno-strict-aliasing -Wstrict-overflow -Wno-missing-field-initializers -flto -fopenmp=libiomp5 -I/usr/lib/gcc/x86_64-linux-gnu/4.8/include/ -march=native -I/opt/Mathematica/SystemFiles/Links/WSTP/DeveloperKit/Linux-x86-64/CompilerAdditions/ -I../lib/
+LDFLAGS = -L/opt/Mathematica/SystemFiles/Links/WSTP/DeveloperKit/Linux-x86-64/CompilerAdditions -lm -lpthread -lrt -lstdc++ -ldl -luuid -l WSTP64i4
+
+convex.o: ../lib/convex.c
+	@${CC} -c -o $@ $< ${CFLAGS}
+
+convexminorant: convexminorant.tm convex.o
+	@${WSPREP} $< -o $@.c
+	@${CC} -o $@ convex.o $@.c ${CFLAGS} ${LDFLAGS} 
+
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);
+}
+