From 5ab8a5088b729a9953b8dff1b2a985dc8fb2098b Mon Sep 17 00:00:00 2001
From: mzuenni <michi.zuendorf@gmail.com>
Date: Mon, 27 Jun 2022 17:19:28 +0200
Subject: updated tcr

---
 math/goldenSectionSearch.cpp | 14 ++++++++++++++
 1 file changed, 14 insertions(+)
 create mode 100644 math/goldenSectionSearch.cpp

(limited to 'math/goldenSectionSearch.cpp')

diff --git a/math/goldenSectionSearch.cpp b/math/goldenSectionSearch.cpp
new file mode 100644
index 0000000..2be8c2d
--- /dev/null
+++ b/math/goldenSectionSearch.cpp
@@ -0,0 +1,14 @@
+ld gss(ld l, ld r, function<ld(ld)> f) {
+	ld r = (sqrtl(5.0l)-1)/2;
+	ld x1 = b - r*(b-a), x2 = a + r*(b-a);
+	ld f1 = f(x1), f2 = f(x2);
+	for (int i = 0; i < 200; i++) {
+		if (f1 < f2) { //change to > to find maximum
+			b = x2; x2 = x1; f2 = f1;
+			x1 = b - r*(b-a); f1 = f(x1);
+		} else {
+			a = x1; x1 = x2; f1 = f2;
+			x2 = a + r*(b-a); f2 = f(x2);
+	}}
+	return a;
+}
-- 
cgit v1.2.3