mirrors/hyperrogue
1
0
Fork 0
mirror of https://github.com/zenorogue/hyperrogue.git synced 2025-01-12 10:20:32 +00:00
Code Issues Releases Wiki Activity

rogueviz::sag:: automatically set RT to average values

Browse Source
This commit is contained in:
Zeno Rogue 2022-10-23 15:49:59 +02:00
parent c8a1c8b5d7
commit a1695a82d0
1 changed files with 20 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
rogueviz/sag.cpp
View File

@ -478,6 +478,21 @@ namespace sag {
} }
} }
void compute_auto_rt() {
ld sum0 = 0, sum1 = 0, sum2 = 0;
for(auto& tab: sagdist) for(auto i: tab) {
sum0 ++;
sum1 += i;
sum2 += i*i;
}
lgsag.R = sum1 / sum0;
lgsag.T = sqrt((sum2 - sum1*sum1/sum0) / sum0);
println(hlog, "automatically set R = ", lgsag.R, " and ", lgsag.T, " max_sag_dist = ", max_sag_dist);
if(logistic_cost) compute_loglik_tab();
}
void optimize_sag_loglik() { void optimize_sag_loglik() {
vector<int> indist(max_sag_dist, 0); vector<int> indist(max_sag_dist, 0);
@ -1057,6 +1072,11 @@ int readArgs() {
shift(); sag::lgsag.T = argf(); shift(); sag::lgsag.T = argf();
if(sag::logistic_cost) compute_loglik_tab(); if(sag::logistic_cost) compute_loglik_tab();
} }
else if(argis("-sagrt-auto")) {
compute_auto_rt();
}
else if(argis("-sag_use_loglik")) { else if(argis("-sag_use_loglik")) {
shift(); sag::logistic_cost = argi(); shift(); sag::logistic_cost = argi();
if(sag::logistic_cost) compute_loglik_tab(); if(sag::logistic_cost) compute_loglik_tab();