add hoc

1 files changed, 75 insertions, 0 deletions

diff --git a/src/cmd/hoc/math.c b/src/cmd/hoc/math.c
new file mode 100644
index 00000000..a4545838
--- /dev/null
+++ b/src/cmd/hoc/math.c
@@ -0,0 +1,75 @@
+#include <u.h>
+#include <libc.h>
+
+#include "hoc.h"
+
+double	errcheck(double, char*);
+
+double
+Log(double x)
+{
+	return errcheck(log(x), "log");
+}
+double
+Log10(double x)
+{
+	return errcheck(log10(x), "log10");
+}
+
+double
+Sqrt(double x)
+{
+	return errcheck(sqrt(x), "sqrt");
+}
+
+double
+Exp(double x)
+{
+	return errcheck(exp(x), "exp");
+}
+
+double
+Asin(double x)
+{
+	return errcheck(asin(x), "asin");
+}
+
+double
+Acos(double x)
+{
+	return errcheck(acos(x), "acos");
+}
+
+double
+Sinh(double x)
+{
+	return errcheck(sinh(x), "sinh");
+}
+double
+Cosh(double x)
+{
+	return errcheck(cosh(x), "cosh");
+}
+double
+Pow(double x, double y)
+{
+	return errcheck(pow(x,y), "exponentiation");
+}
+
+double
+integer(double x)
+{
+	if(x<-2147483648.0 || x>2147483647.0)
+		execerror("argument out of domain", 0);
+	return (double)(long)x;
+}
+
+double
+errcheck(double d, char* s)	/* check result of library call */
+{
+	if(isNaN(d))
+		execerror(s, "argument out of domain");
+	if(isInf(d, 0))
+		execerror(s, "result out of range");
+	return d;
+}