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/math-lib/cos_Q.c~ | 46 ++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 46 insertions(+)
 create mode 100644 ic-reals-6.3/math-lib/cos_Q.c~

(limited to 'ic-reals-6.3/math-lib/cos_Q.c~')

diff --git a/ic-reals-6.3/math-lib/cos_Q.c~ b/ic-reals-6.3/math-lib/cos_Q.c~
new file mode 100644
index 0000000..b1a4a69
--- /dev/null
+++ b/ic-reals-6.3/math-lib/cos_Q.c~
@@ -0,0 +1,46 @@
+/*
+ * Copyright (C) 2000, Imperial College
+ *
+ * This file is part of the Imperial College Exact Real Arithmetic Library.
+ * See the copyright notice included in the distribution for conditions
+ * of use.
+ */
+
+#include <stdio.h>
+#include "real.h"
+
+Real
+cos_QZ(mpz_t a, mpz_t b)
+{
+	mpz_t x;
+	Real r;
+
+	mpz_init(x);
+	mpz_mul_ui(x, b, 2);
+	r = tan_QZ(a, x);
+	mpz_clear(x);
+	r = tensor_Int(r, r, -1, 1, 0, 0, 0, 0, 1, 1);
+	return r;
+}
+
+Real
+cos_QInt(int a, int b)
+{
+	Real r;
+	mpz_t ap, bp;
+
+	/* check for overflow */
+	if (b > 0x3FFFFFFF || b < 0xC0000000) {
+		mpz_init_set_si(ap, a);
+		mpz_init_set_si(bp, b);
+		mpz_mul_ui(bp, bp, 2);
+		r = tan_QZ(ap, bp);
+		mpz_clear(ap);
+		mpz_clear(bp);
+	}
+	else 
+		r = tan_QInt(a, b * 2);
+
+	r = tensor_Int(r, r, -1, 1, 0, 0, 0, 0, 1, 1);
+	return r;
+}
-- 
cgit v1.2.3