math/goldenSectionSearch.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

ld gss(ld l, ld r, function<ld(ld)> f) {
	ld inv = (sqrt(5.0l) - 1) / 2;
	ld x1 = r - inv*(r-l), x2 = l + inv*(r-l);
	ld f1 = f(x1), f2 = f(x2);
	for (int i = 0; i < 200; i++) {
		if (f1 < f2) { //change to > to find maximum
			u = x2; x2 = x1; f2 = f1;
			x1 = r - inv*(r-l); f1 = f(x1);
		} else {
			l = x1; x1 = x2; f1 = f2;
			x2 = l + inv*(r-l); f2 = f(x2);
		}
	}
	return l;
}