summaryrefslogtreecommitdiff
path: root/math/goldenSectionSearch.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'math/goldenSectionSearch.cpp')
-rw-r--r--math/goldenSectionSearch.cpp14
1 files changed, 14 insertions, 0 deletions
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;
+}