From ff07260b38d521c771346b147c1b43601ebbdf02 Mon Sep 17 00:00:00 2001
From: = <computergeek3.142@gmail.com>
Date: Sun, 23 Aug 2020 16:17:00 -0700
Subject: [PATCH] Adding convergents based search

---
 Makefile         |  6 ++++
 cf_convergents.c | 81 ++++++++++++++++++++++++++++++++++++++++++++++++
 2 files changed, 87 insertions(+)
 create mode 100644 Makefile
 create mode 100644 cf_convergents.c

diff --git a/Makefile b/Makefile
new file mode 100644
index 0000000..9af3437
--- /dev/null
+++ b/Makefile
@@ -0,0 +1,6 @@
+CC=gcc
+L=-lm -lgmp
+
+default: cf_convergents.c
+	$(CC) cf_convergents.c $(L) -o convergents
+
diff --git a/cf_convergents.c b/cf_convergents.c
new file mode 100644
index 0000000..f9bc41c
--- /dev/null
+++ b/cf_convergents.c
@@ -0,0 +1,81 @@
+#include <stdlib.h>
+#include <stdio.h>
+#include <math.h>
+#include <gmp.h>
+
+mpz_t temp_h;
+mpz_t temp_k;
+mpf_t temp_f;
+mpf_t temp_f2;
+
+void get_cf(mpf_t r, mpz_t h, mpz_t k, mpz_t h_old, mpz_t k_old, mpz_t a){
+	mpz_set_f(a, r);
+	mpf_set_z(temp_f, a);
+	mpf_sub(r, r, temp_f);
+	mpf_ui_div(r, 1, r);
+
+	mpz_mul(temp_h, a, h);
+	mpz_add(temp_h, temp_h, h_old);
+	mpz_mul(temp_k, a, k);
+	mpz_add(temp_k, temp_k, k_old);
+
+	mpz_set(h_old, h);
+	mpz_set(k_old, k);
+	mpz_set(h, temp_h);
+	mpz_set(k, temp_k);
+}
+
+int main(int argc, char **argv){
+	mpf_t r;
+	mpz_t h;
+	mpz_t k;
+	mpz_t h_old;
+	mpz_t k_old;
+	mpz_t a;
+	int i;
+	double d;
+
+	if(argc > 1){
+		//Allow the user to specify a default precision using a command-line argument
+		mpf_set_default_prec(atoi(argv[1]));
+	}
+
+	mpf_init(r);
+	mpf_init(temp_f);
+	mpf_init(temp_f2);
+
+	mpz_init(temp_h);
+	mpz_init(temp_k);
+	mpz_init(h);
+	mpz_init(k);
+	mpz_init(h_old);
+	mpz_init(k_old);
+	mpz_init(a);
+
+	mpf_set_d(r, M_PI/2);
+
+	mpz_set_si(h, 1);
+	mpz_set_si(k_old, 1);
+
+	for(i = 0; i < 10; i++){
+		get_cf(r, h, k, h_old, k_old, a);
+		if(mpz_odd_p(k)){
+			printf("candidate: ");
+			mpz_out_str(stdout, 0, h);
+			printf("\n");
+		}
+	}
+
+	mpf_clear(r);
+	mpf_clear(temp_f);
+	mpf_clear(temp_f2);
+
+	mpz_clear(temp_h);
+	mpz_clear(temp_k);
+	mpz_clear(h);
+	mpz_clear(k);
+	mpz_clear(h_old);
+	mpz_clear(k_old);
+	mpz_clear(a);
+}
+