From 11da511c784eca003deb90c23570f0873954e0de Mon Sep 17 00:00:00 2001
From: Duncan Wilkie <antigravityd@gmail.com>
Date: Sat, 18 Nov 2023 06:11:09 -0600
Subject: Initial commit.

---
 ic-reals-6.3/tests/t2.c | 42 ++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 42 insertions(+)
 create mode 100644 ic-reals-6.3/tests/t2.c

(limited to 'ic-reals-6.3/tests/t2.c')

diff --git a/ic-reals-6.3/tests/t2.c b/ic-reals-6.3/tests/t2.c
new file mode 100644
index 0000000..a41eed0
--- /dev/null
+++ b/ic-reals-6.3/tests/t2.c
@@ -0,0 +1,42 @@
+#include <stdio.h>
+#include "real.h"
+#include <math.h>
+
+/*
+ * computes log(exp(x))
+ */
+main(int argc, char *argv[])
+{
+	Real x, y, z;
+	double f;
+	Real makeRealSignCNQInt(Sign, char *, int, int, int);
+
+	MyName = argv[0];
+
+	if (argc != 7) {
+		fprintf(stderr, "%s <sign> <c> <n> <a> <b> <ndigits>\n", MyName);
+		exit(1);
+	}
+
+	initReals();
+
+	y = makeRealSignCNQInt(
+				atoi(argv[1]),	/* sign */
+				argv[2],		/* c */
+				atoi(argv[3]), 	/* n */
+				atoi(argv[4]),	/* a */
+				atoi(argv[5]));	/* b */
+
+	print_R_Dec(y, atoi(argv[6]));
+	printf("\n");
+
+	f = realToDouble(y);
+	printf("x=%f\n",f);
+
+	x = exp_R(y);
+	x = log_R(x);
+
+	print_R_Dec(x, atoi(argv[6]));
+	printf("\n");
+	printf("t2(x)=%f\n", log(exp(f)));
+}
-- 
cgit v1.2.3