summaryrefslogtreecommitdiff
path: root/test/math/goldenSectionSearch.cpp
diff options
context:
space:
mode:
authormzuenni <mzuenni@users.noreply.github.com>2024-07-28 22:54:40 +0200
committerGitHub <noreply@github.com>2024-07-28 22:54:40 +0200
commit8d11c6c8213f46f0fa19826917c255edd5d43cb1 (patch)
tree96d75baff33d5a04b5a60f1a41f514a26c716874 /test/math/goldenSectionSearch.cpp
parent8c33b4e0d3030cfed17fc64b4fe41133339f6d87 (diff)
Test (#4)
* update * moved content in subdir * rename file * add test setup * add test setup * add github action * automaticly test all cpp files * timeout after 10s * setulimit and dont zero memory * test build pdf * install latexmk * update * update * ngerman * fonts * removed old code * add first test * added tests * test in sorted order * more tests * simplified test * more tests * fix suffix tree * fixes and improvements * done ust lst directly * fix swap * add links to pdf * fix constants * add primorial * add comment * various improvements * more tests * added missing stuf * more tests * fix tests * more tests * more tests * more tests * fix recursion? * test trie * more tests * only use python temporarily for listings * only use python temporarily for listings * more tests * fix longestCommonSubstring * more tests * more tests * made code more similiar * fix? * more tests * more tests * more tests * add ahoCorasick test + limit 4GB stack size * more tests * fix test * add additional test * more tests * more tests * fix? * better fix * fix virtual tree * more tests * more tests * recursive closest pair * more tests * decrease limit * new tests * more tests * fix name * more tests * add test * new test * more tests * more tests * more tests * more tests * new test and content * new code * new code * larger tests * fix and test * new test * new test * update pdf * remove comments * new test * more tests * more testcases * more tests * increased limit * more tests * more tests * more tests * new tests * more tests * shortened code * new test * add basic tests for bigint * more tests * removed old files * new test * ignore some files * more auto more ccw * fix test * more tests * fix * new tests * more tests * more tests * stronger test * actually verify delaunay... * more tests * fix header * more tests * run tests parallel? * test parralel? * add --missing * separate workflows * test * is the pdf checked? * separate workflows * fix workflow * more workflows --------- Co-authored-by: Yidi <noob999noob999@gmail.com>
Diffstat (limited to 'test/math/goldenSectionSearch.cpp')
-rw-r--r--test/math/goldenSectionSearch.cpp74
1 files changed, 74 insertions, 0 deletions
diff --git a/test/math/goldenSectionSearch.cpp b/test/math/goldenSectionSearch.cpp
new file mode 100644
index 0000000..565a21c
--- /dev/null
+++ b/test/math/goldenSectionSearch.cpp
@@ -0,0 +1,74 @@
+#include "../util.h"
+#include <math/goldenSectionSearch.cpp>
+
+struct RandomFunction {
+ ld min;
+ vector<pair<ld, int>> polys;
+ RandomFunction(ld l, ld r) : min(Random::real<ld>(l, r)) {
+ do {
+ polys.emplace_back(
+ Random::real<ld>(0, 1e9),
+ 2 * Random::integer<int>(1, 5)
+ );
+ } while(false && Random::integer<int>(4) != 0);
+ }
+
+ ld operator()(ld x){
+ ld res = 0;
+ for (auto [m, p] : polys) {
+ res += m * pow(x - min, p);
+ }
+ return res;
+ }
+
+ friend ostream& operator<<(ostream& os, const RandomFunction& f) {
+ string plus = "";
+ for (auto [m, p] : f.polys) {
+ os << setprecision(15) << plus << m << "*(x-" << f.min << ")**" << p;
+ plus = "+";
+ }
+ return os;
+ }
+};
+
+void stress_test() {
+ int queries = 0;
+ for (int i = 0; i < 50'000; i++) {
+ ld l = Random::real<double>(-200, 200);
+ ld r = Random::real<double>(-200, 200);
+ if (l > r) swap(l, r);
+
+ RandomFunction f(l, r);
+
+ ld got = gss(l, r, f);
+ ld expected = f.min;
+ if (float_error(got, expected) > 1e-6) {
+ cerr << f << endl;
+ cerr << "got: " << got << ", expected: " << expected << FAIL;
+ }
+ queries++;
+ }
+ cerr << "tested random queries: " << queries << endl;
+}
+
+constexpr int N = 10'000;
+void performance_test() {
+ timer t;
+ RandomFunction f(-200, 200);
+ f.polys.resize(1);
+
+ hash_t hash = 0;
+ for (int i = 0; i < N; i++) {
+ t.start();
+ hash += gss(-200, 200, f);
+ t.stop();
+ }
+ if (t.time > 1000) cerr << "too slow: " << t.time << FAIL;
+ cerr << "tested performance: " << t.time << "ms (hash: " << hash << ")" << endl;
+}
+
+int main() {
+ stress_test();
+ performance_test();
+}
+