mirrors/hyperrogue
1
0
Fork 0
mirror of https://github.com/zenorogue/hyperrogue.git synced 2025-10-31 05:52:59 +00:00
Code Issues Releases Wiki Activity
Files
56ee09fe1ce143311e87a90783a7e0d404c9c310
hyperrogue/geom-exp.cpp

1411 lines
43 KiB
C++
Raw Blame History

// Hyperbolic Rogue -- the 'experiments with geometry' menu
// Copyright (C) 2011-2019 Zeno Rogue, see 'hyper.cpp' for details
/** \file geom-exp.cpp
* \brief The 'experiments with geometry' menu
*
* Implementation of this menu, and computation of the statistics shown there
*/
#include "hyper.h"
namespace hr {
int eupage = 0;
int euperpage = 9;
string euchelp =
"If you want to know how much the gameplay is affected by the "
"hyperbolic geometry in HyperRogue, this mode is for you!\n\n"
"You can try many different geometries here. We start by gluing "
"n-gons in such a way that k of them meet in every vertex. "
"Depending on n and k, this either folds into a sphere, unfolds into a plane, "
"or requires a hyperbolic space. The result may be then 'bitruncated' by "
"replacing each vertex by a 2k-gon. Furthermore, you can play "
"with quotient geometries. For example, the elliptic geometry is "
"obtained from the sphere by making the antipodes be the same point, "
"so you return to the same spot (but as a mirror image) after going there. "
"Have fun experimenting! "
"Achievements and leaderboards do not work in geometry experiments, "
"except some specific ones.\n\n"
"In standard geometry (bitruncated or not), you can play the full game, but in other geometries "
"you select a particular land. Lands are unlocked by visiting them in this "
"session, or permanently by collecting 25 treasure. Try Crossroads in Euclidean "
"or chaos mode in non-standard non-quotient hyperbolic to visit many lands. "
"Highlights:\n"
"* Crystal World and Warped Coast can be understood as extra geometries.\n"
"* Halloween is specially designed for spherical geometry.\n"
"* To see the difference, try Hunting Grounds in Euclidean -- it is impossible.\n";
#if CAP_FIELD
void showQuotientConfig() {
using namespace fieldpattern;
cmode = sm::SIDE | sm::MAYDARK;
gamescreen();
dialog::init(XLAT("field quotient"));
fgeomextra& gxcur = fgeomextras[current_extra];
for(int i=0; i<isize(fgeomextras); i++) {
auto& g = fgeomextras[i];
dialog::addBoolItem(ginf[g.base].tiling_name, g.base == gxcur.base, 'a'+i);
dialog::add_action([i] { current_extra = i; });
}
dialog::addBreak(100);
nextPrimes(gxcur);
string stars[3] = {"", "*", "**"};
for(int i=0; i<isize(gxcur.primes); i++) {
auto& p = gxcur.primes[i];
int sides = ginf[gxcur.base].sides;
bool sa = shapename_available(sides);
string s = sa ? "order %1%2 (%3 %4s)" : "order %1%2 (%3 %4-gons)";
s = XLAT(s, its(p.p), p.squared ? "²" : "", its(p.cells), sa ? shapename(sides) : its(sides));
s += stars[gxcur.dualval[i]];
dialog::addBoolItem(s, i == gxcur.current_prime_id, 'A'+i);
dialog::add_action([&gxcur, i] { gxcur.current_prime_id = i; });
}
if(isize(gxcur.primes) <= 6) {
dialog::addBreak(100);
dialog::addHelp(
"This geometry is obtained by applying the same 'generators' which "
"lead to creating the given basic hyperbolic geometry, "
"but using a fixed finite field instead of the field of reals. "
"It can be also interpreted as a quotient of the given basic geometry. "
"Warning: field patterns based on large primes might generate for a long time."
);
dialog::addBreak(100);
}
dialog::addItem("find the next prime", 'p');
dialog::add_action([&gxcur] { nextPrime(gxcur); });
dialog::addItem("activate", 'x');
dialog::add_action_confirmed([&gxcur] {
set_geometry(gxcur.base);
enableFieldChange();
set_geometry(gFieldQuotient);
start_game();
});
dialog::addItem("find alternate manifolds", 'y');
dialog::add_action_confirmed([&gxcur] {
set_geometry(gxcur.base);
triplet_id = 0;
enableFieldChange();
set_geometry(gFieldQuotient);
start_game();
auto& cfp = currfp;
auto triplets = cfp.find_triplets();
pushScreen([triplets] {
cmode = sm::SIDE | sm::MAYDARK;
gamescreen();
dialog::init(XLAT("alternate manifolds"));
int id = 0;
for(auto t: triplets) {
dialog::addItem(XLAT("generators (%1,%2), size %3", its(t.i), its(t.j), its(t.size/S7)), 'a'+id);
dialog::add_action([id] {
stop_game();
triplet_id = id;
fieldpattern::enableFieldChange();
start_game();
});
id++;
}
dialog::addBreak(100);
dialog::addHelp(XLAT(
"This option finds alternate solutions. For example, there are three {7,3} manifolds with 156 heptagons each (\"first Hurwitz triplet\").")
);
dialog::display();
});
});
dialog::addItem("default", 'z');
dialog::add_action_confirmed([] {
set_geometry(gEuclid);
fieldpattern::quotient_field_changed = false;
set_geometry(gFieldQuotient);
start_game();
});
dialog::addBack();
dialog::display();
}
#endif
EX string bitruncnames[5] = {" (b)", " (n)", " (g)", " (i)", " (d)"};
void validity_info() {
int vccolors[4] = {0xFF0000, 0xFF8000, 0xFFFF00, 0x00FF00};
auto lv = land_validity(specialland);
if(lv.flags & lv::display_error_message)
dialog::addInfo(XLAT(lv.msg), vccolors[lv.quality_level]);
else
dialog::addBreak(100);
}
EX bool showquotients;
string validclasses[4] = {" (X)", " (½)", "", " (!)"};
EX void gen_landvisited() {
if(cheater) for(int i=0; i<landtypes; i++) landvisited[i] = true;
for(int i=0; i<landtypes; i++)
if(hiitemsMax(treasureType(eLand(i))) >= 25) landvisited[i] = true;
landvisited[laCrossroads] = true;
landvisited[laCrossroads4] = true;
landvisited[laIce] = true;
landvisited[laHunting] = true;
landvisited[laMirrorOld] = true;
landvisited[laPrincessQuest] = cheater || princess::everSaved;
landvisited[laWildWest] = true;
landvisited[laHalloween] = true;
landvisited[laWarpCoast] = true;
landvisited[laGraveyard] = true;
landvisited[laDual] = true;
landvisited[laDocks] |= landvisited[laWarpCoast];
landvisited[laSnakeNest] |= landvisited[laRedRock];
landvisited[laCamelot] |= hiitemsMax(treasureType(laCamelot)) >= 1;
landvisited[laCA] = true;
landvisited[laAsteroids] = true;
}
EX void ge_land_selection() {
cmode = sm::SIDE | sm::MAYDARK;
gamescreen();
gen_landvisited();
dialog::init(XLAT("select the starting land"));
if(dialog::infix != "") mouseovers = dialog::infix;
generateLandList([] (eLand l) {
if(dialog::infix != "" && !dialog::hasInfix(linf[l].name)) return false;
return !!(land_validity(l).flags & lv::appears_in_geom_exp);
});
stable_sort(landlist.begin(), landlist.end(), [] (eLand l1, eLand l2) { return land_validity(l1).quality_level > land_validity(l2).quality_level; });
dialog::start_list(900, 900, '1');
for(auto& l: landlist) {
string s = XLAT1(linf[l].name);
if(landvisited[l] || unlock_all) {
dialog::addBoolItem(s, l == specialland, dialog::list_fake_key++);
}
else {
dialog::addSelItem(s, XLAT("(locked)"), dialog::list_fake_key++);
}
dialog::lastItem().color = linf[l].color;
dialog::lastItem().value += validclasses[land_validity(l).quality_level];
dialog::add_action([l] {
if(landvisited[l]) dialog::do_if_confirmed(dual::mayboth([l] {
stop_game_and_switch_mode(tactic::on ? rg::tactic : rg::nothing);
firstland = specialland = l;
if(l == laCanvas || l == laAsteroids || (land_validity(l).flags & lv::switch_to_single))
land_structure = lsSingle;
else if(among(l, laCrossroads, laCrossroads2))
land_structure = lsNiceWalls;
else if(among(l, laCrossroads4))
land_structure = lsNoWalls;
start_game();
popScreen();
}));
});
}
dialog::end_list();
dialog::addBreak(25);
validity_info();
dialog::addBreak(25);
dialog::addInfo(XLAT("press letters to search"));
dual::add_choice();
dialog::addBack();
dialog::display();
keyhandler = [] (int sym, int uni) {
dialog::handleNavigation(sym, uni);
if(dialog::editInfix(uni)) dialog::list_skip = 0;
else if(doexiton(sym, uni)) popScreen();
};
}
EX void activate_ge_land_selection() {
dialog::infix = "";
eupage = 0;
pushScreen(ge_land_selection);
}
#if HDR
struct geometry_filter {
string name;
int hotkey;
/** test if the current geometry matches the filter */
function<bool()> test;
};
#endif
EX geometry_filter *current_filter;
bool forced_quotient() { return quotient && !(cgflags & qOPTQ); }
EX geometry_filter gf_hyperbolic = {"hyperbolic", 'h', [] { return (arcm::in() || arb::in() || hyperbolic) && !forced_quotient(); }};
EX geometry_filter gf_spherical = {"spherical", 's', [] { return (arcm::in() || arb::in() || sphere) && !forced_quotient(); }};
EX geometry_filter gf_euclidean = {"Euclidean", 'e', [] { return (arcm::in() || arb::in() || euclid) && !forced_quotient(); }};
EX geometry_filter gf_other = {"non-isotropic", 'n', [] { return mproduct || mtwisted || nonisotropic; }};
EX geometry_filter gf_regular_2d = {"regular 2D tesselations", 'r', [] {
return standard_tiling() && WDIM == 2 && !forced_quotient();
}};
EX geometry_filter gf_regular_3d = {"regular 3D honeycombs", '3', [] {
if(euclid) return geometry == gCubeTiling;
return !bt::in() && !kite::in() && WDIM == 3 && !forced_quotient() && !nonisotropic && !mproduct;
}};
EX geometry_filter gf_quotient = {"interesting quotient spaces", 'q', [] {
return forced_quotient() && !elliptic;
}};
EX geometry_filter gf_tes_file = {"load from file", 'f', [] { return arb::in(); }};
EX geometry_filter gf_no_filters = {"no filters", '-', [] { return true; }};
EX vector<geometry_filter*> available_filters = { &gf_hyperbolic, &gf_spherical, &gf_euclidean, &gf_other, &gf_regular_2d, &gf_regular_3d, &gf_quotient, &gf_tes_file, &gf_no_filters };
void ge_select_filter() {
cmode = sm::SIDE | sm::MAYDARK;
gamescreen();
dialog::init(XLAT("geometries"));
for(auto f: available_filters) {
if(current_filter)
dialog::addBoolItem(XLAT(f->name), f == current_filter, f->hotkey);
else
dialog::addItem(XLAT(f->name), f->hotkey);
dialog::add_action([f] { current_filter = f; popScreen(); });
}
dialog::addBack();
dialog::display();
}
void set_default_filter() {
current_filter = hyperbolic ? &gf_hyperbolic : euclid ? &gf_euclidean : sphere ? &gf_spherical : &gf_other;
for(auto f: available_filters) if(f->test()) current_filter = f;
}
void set_or_configure_geometry(eGeometry g) {
if(0) ;
#if CAP_CRYSTAL
else if(g == gCrystal)
pushScreen(crystal::show);
#endif
#if CAP_ARCM
else if(g == gArchimedean)
pushScreen(arcm::show);
#endif
else if(g == gArbitrary)
arb::choose();
else {
bool quo = false;
if(among(g, gProduct, gTwistedProduct)) {
if(WDIM == 3) {
addMessage(XLAT("Only works with 2D geometries"));
return;
}
if(g == gTwistedProduct) {
if(nonorientable) {
addMessage(XLAT("Only works in orientable spaces"));
return;
}
bool ok = true;
quo = sphere || quotient;
if(arcm::in()) ok = PURE;
else if(bt::in() || aperiodic) ok = false;
else if(GOLDBERG && S3 == 4 && gp::param == gp::loc{1,1}) ok = true;
else if(UNRECTIFIED && gp::param == gp::loc{1,1}) ok = true;
else ok = PURE || BITRUNCATED;
if(!ok) {
addMessage(XLAT("Only works with (semi-)regular tilings"));
return;
}
#if CAP_ARCM
if(arcm::in()) {
int steps, single_step;
if(!arcm::current.get_step_values(steps, single_step)) {
addMessage(XLAT("That would have %1/%2 levels", its(steps), its(single_step)));
return;
}
}
#endif
}
}
dual::may_split_or_do([g, quo] {
set_geometry(g);
if(quo) hybrid::fixup_csteps();
});
start_game();
}
}
/** is g2 the same tiling as the current geometry (geometry)? */
bool same_tiling(eGeometry g2) {
/* no quotients for fractals */
if(cgflags & qFRACTAL) return g2 == geometry;
if(g2 == gCrystal)
return S3 == 4;
if(g2 == gFieldQuotient && (hyperbolic || (geometry == gCubeTiling && reg3::cubes_reg3)) && standard_tiling())
return true;
if(g2 == gFieldQuotient && geometry != gFieldQuotient) {
int ce = 0;
for(auto& ge: fieldpattern::fgeomextras) {
if(ginf[ge.base].tiling_name == ginf[geometry].tiling_name) {
fieldpattern::current_extra = ce;
return true;
}
ce++;
}
}
return ginf[g2].tiling_name == ginf[geometry].tiling_name;
}
void ge_select_tiling() {
if(current_filter == &gf_tes_file) { popScreen(); set_or_configure_geometry(gArbitrary); }
cmode = sm::SIDE | sm::MAYDARK;
gamescreen();
if(!current_filter) { popScreen(); return; }
dialog::init();
dialog::addItem(XLAT(current_filter->name), 'x');
dialog::add_action_push(ge_select_filter);
vector<eGeometry> geometries;
dialog::addBreak(100);
dialog::start_list(500, 1500, 'a');
for(int i=0; i<isize(ginf); i++) {
eGeometry g = eGeometry(i);
if(among(g, gProduct, gTwistedProduct)) hybrid::configure(g);
bool orig_el = elliptic;
bool on = geometry == g;
bool in_2d = WDIM == 2;
dynamicval<eGeometry> cg(geometry, g);
if(cgflags & qDEPRECATED) continue;
if(arcm::in() && !CAP_ARCM) continue;
if(cryst && !CAP_CRYSTAL) continue;
if(sol && !CAP_SOLV) continue;
if(arb::in() && (ISMOBILE || ISWEB)) continue;
if(WDIM == 3 && MAXMDIM == 3) continue;
if(geometry == gFieldQuotient && !CAP_FIELD) continue;
if(geometry == gFake) continue;
if(!current_filter->test()) continue;
if(orig_el) {
for(int j=0; j<isize(ginf); j++)
if(ginf[j].tiling_name == ginf[i].tiling_name)
geometry = g = eGeometry(j);
}
bool is_product = (geometry == gProduct && in_2d);
bool is_twisted = (geometry == gTwistedProduct && in_2d);
dialog::addBoolItem(
is_product ? XLAT("current geometry x E") :
is_twisted ? XLAT("twisted current x E") :
XLAT(ginf[g].menu_displayed_name), on, dialog::list_fake_key++);
dialog::lastItem().value += validclasses[land_validity(specialland).quality_level];
dialog::add_action([g] { set_or_configure_geometry(g); });
}
dialog::end_list();
dialog::addBreak(100);
dual::add_choice();
dialog::addBack();
dialog::display();
}
EX string current_proj_name() {
bool h = hyperbolic || sn::in();
if(vpconf.model == mdPanini && vpconf.alpha == 1)
return XLAT("stereographic Panini");
else if(vpconf.model != mdDisk)
return models::get_model_name(vpconf.model);
else if(h && vpconf.alpha == 1)
return XLAT("Poincaré model");
else if(h && vpconf.alpha == 0)
return XLAT("Klein-Beltrami model");
else if(h && vpconf.alpha == -1)
return XLAT("inverted Poincaré model");
else if(sphere && vpconf.alpha == 1)
return XLAT("stereographic projection");
else if(sphere && vpconf.alpha == 0)
return XLAT("gnomonic projection");
else if(sphere && vpconf.alpha >= 999)
return XLAT("orthographic projection");
else if(h && vpconf.alpha >= 999)
return XLAT("Gans model");
else
return XLAT("general perspective");
}
EX string dim_name() {
return " (" + its(WDIM) + "D)";
}
#if CAP_THREAD && MAXMDIM >= 4
EX void showQuotientConfig3() {
using namespace fieldpattern;
cmode = sm::SIDE | sm::MAYDARK;
gamescreen();
dialog::init(XLAT("field quotient"));
auto& ds = discoveries[cginf.tiling_name];
if(!ds.discoverer) {
dialog::addItem("start discovery", 's');
dialog::add_action([&ds] { ds.activate(); });
}
else if(ds.is_suspended) {
dialog::addItem("resume discovery", 's');
dialog::add_action([&ds] { ds.activate(); });
}
else {
dialog::addItem("suspend discovery", 's');
dialog::add_action([&ds] { ds.suspend(); });
}
auto& e = ds.experiment;
if(!e.Prime)
dialog::addBreak(100);
else {
string s = its(e.Prime);
if(e.wsquare) s += "²";
dialog::addInfo(s);
}
dialog::addBreak(100);
if(1) {
std::unique_lock<std::mutex> lk(ds.lock);
auto&l = ds.hashes_found;
for(auto& v: l) {
char x = 'a';
string s = XLAT("#%1, cells: %2, p=%3", itsh(v.first), its(get<5>(v.second)), its(get<0>(v.second)) + (get<1>(v.second) ? "²" : ""));
dialog::addItem(s, x++);
dialog::add_action([&v] {
stop_game();
int tmp;
tie(currfp.Prime, currfp.wsquare, currfp.R, currfp.P, currfp.X, tmp) = v.second;
currfp.Field = currfp.wsquare ? currfp.Prime * currfp.Prime : currfp.Prime;
currfp.generate_all3();
currfp.analyze();
start_game();
});
}
}
dialog::addBreak(100);
dialog::addBack();
dialog::display();
}
#endif
EX string geometry_name(eGeometryClass gc) {
switch(gc) {
case gcHyperbolic:
return XLAT("hyperbolic");
case gcEuclid:
if(cgflags & qAFFINE)
return XLAT("affine");
return XLAT("flat");
case gcSphere:
return XLAT("spherical");
case gcSol:
return XLAT("Sol");
case gcNIH:
return XLAT("hyperbolic (3:2)");
case gcSolN:
return XLAT("Sol (3:2)");
case gcNil:
return XLAT("Nil");
case gcSL2:
return XLAT("~SL(2,R)~");
case gcProduct:
return XLAT("%1 x E", PIU(geometry_name()));
}
return "?";
}
EX string geometry_name() {
if(cgi.emb && cgi.emb->is_same_in_same())
return geometry_name(geom3::mgclass());
else if(embedded_plane && gproduct)
return geometry_name(geom3::mgclass()) + " (x E)";
else if(embedded_plane)
return geometry_name(geom3::mgclass()) + " @ " + geometry_name(geom3::ggclass());
else if(among(ginf[geometry].cclass, gcHyperbolic, gcEuclid, gcSphere))
return geometry_name(ginf[geometry].cclass) + dim_name();
else
return geometry_name(ginf[geometry].cclass);
}
EX void select_quotient_screen() {
cmode = sm::SIDE | sm::MAYDARK;
gamescreen();
dialog::init(XLAT("quotient spaces in ") + ginf[geometry].tiling_name);
char key = 'a';
for(int i=0; i<isize(ginf); i++) {
auto g = eGeometry(i);
if(ginf[g].flags & qDEPRECATED) continue;
if(same_tiling(g)) {
dialog::addBoolItem(
(ginf[g].flags & qANYQ) ?
XLAT(ginf[g].menu_displayed_name) :
"no quotient",
g == geometry, key++);
dialog::add_action([g] {
if(g == gFieldQuotient && WDIM == 3) {
if(geometry != gFieldQuotient) {
stop_game();
fieldpattern::field_from_current();
set_geometry(gFieldQuotient);
for(int p=2;; p++) { currfp.Prime = p; currfp.force_hash = 0; if(!currfp.solve()) break; }
println(hlog, "set prime = ", currfp.Prime);
start_game();
}
#if CAP_THREAD && MAXMDIM >= 4
pushScreen(showQuotientConfig3);
#endif
}
else if(g == gFieldQuotient)
pushScreen(showQuotientConfig);
#if CAP_CRYSTAL
else if(g == gCrystal)
pushScreen(crystal::show);
#endif
else {
dual::may_split_or_do([g] { set_geometry(g); });
start_game();
}
});
}
}
dialog::addBack();
dialog::display();
}
EX void select_quotient() {
if(meuclid && !aperiodic && !arcm::in() && !reg3::cubes_reg3 && !(cgflags & qFRACTAL)) {
euc::prepare_torus3();
pushScreen(euc::show_torus3);
}
else if(nil && !mhybrid) {
nilv::prepare_niltorus3(),
pushScreen(nilv::show_niltorus3);
}
#if CAP_SOLV
else if(asonov::in()) {
asonov::prepare_config();
pushScreen(asonov::show_config);
}
#endif
else if(mproduct)
pushScreen(product::show_config);
else if(mtwisted)
hybrid::configure_period();
else {
vector<eGeometry> choices;
for(int i=0; i<isize(ginf); i++) if(same_tiling(eGeometry(i))) choices.push_back(eGeometry(i));
println(hlog, "choices = ", choices);
if(isize(choices) > 1)
pushScreen(select_quotient_screen);
/* else if(isize(choices) > 1) {
set_geometry(choices[choices[0] == geometry ? 1 : 0]);
start_game();
} */
else
addMessage("No quotient spaces available in the current tiling.");
}
}
EX string full_geometry_name() {
string qstring = ginf[geometry].quotient_name;
bool variable =
!(mproduct || mhybrid || bt::in() || (WDIM == 3 && !reg3::in()) || aperiodic || arb::in());
string fgname = XLAT(ginf[geometry].tiling_name);
if(qstring != "none") fgname += " " + XLAT(qstring);
#if CAP_ARCM
if(arcm::in()) fgname = arcm::current.symbol;
#endif
if(variable) fgname = gp::operation_name() + " " + fgname;
return fgname;
}
void action_change_variation() {
if(0) ;
#if CAP_ARCM
else if(arcm::in()) arcm::next_variation();
#endif
#if MAXMDIM >= 4
else if(reg3::in() || geometry == gCubeTiling) reg3::configure_variation();
#endif
else if(euc::in(2,4) || !CAP_GP) dialog::do_if_confirmed([] {
set_variation(PURE ? eVariation::bitruncated : eVariation::pure);
start_game();
});
#if CAP_GP
else // if(S3 == 3)
gp::configure();
#endif
}
EX void menuitem_change_variation(key_type key) {
dialog::addSelItem(XLAT("variations"), gp::operation_name(), key);
dialog::add_action(action_change_variation);
}
EX void menuitem_change_geometry(key_type key) {
dialog::addSelItem(XLAT("geometry/topology/tiling"), full_geometry_name(), key);
dialog::add_action_push(current_filter ? ge_select_tiling : ge_select_filter);
}
EX void menuitem_projection(key_type key) {
dialog::addSelItem(XLAT("projection"), current_proj_name(), key);
dialog::add_action_push(models::model_menu);
}
EX void menuitem_binary_width(key_type key) {
dialog::addSelItem(XLAT("binary tiling width"), fts(vid.binary_width), key);
dialog::add_action([] {
dialog::editNumber(vid.binary_width, 0, 2, 0.1, 1, XLAT("binary tiling width"), "");
dialog::get_ne().reaction = [] () {
#if CAP_TEXTURE
texture::config.remap();
#endif
#if CAP_SOLV
if(asonov::in()) asonov::prepare();
#endif
};
});
}
EX void menuitem_nilwidth(key_type key) {
dialog::addSelItem(XLAT("Nil width"), fts(nilv::nilwidth), key);
dialog::add_action([] {
dialog::editNumber(nilv::nilwidth, 0.01, 2, 0.1, 1, XLAT("Nil width"), "");
dialog::get_ne().reaction = [] {
ray::reset_raycaster();
twist::clear_twisted_matrices();
};
dialog::bound_low(0.01);
});
}
EX void edit_stretch() {
stretch::mstretch = false;
ray::reset_raycaster();
dialog::editNumber(stretch::factor, -1, 9, 0.1, 0, XLAT("stretched geometry"),
XLAT(
"Stretch the metric along the fibers. This can currently be done in rotation spaces and in 8-cell, 24-cell and 120-cell. "
"Value of 0 means not stretched, -1 means S2xE or H2xE (works only in the limit). (Must be > -1)"
)
);
dialog::get_ne().reaction = [] { if(abs(stretch::factor+1) < 1e-3) stretch::factor = -.9; ray::reset_raycaster(); };
}
#if HDR
struct geometry_data {
int ts, tv, nom, denom, euler, demigenus, worldsize;
int area;
string spf;
string size_str;
};
#endif
EX geometry_data compute_geometry_data() {
geometry_data gd;
auto& ts = gd.ts;
auto& tv = gd.tv;
ts = ginf[geometry].sides;
tv = ginf[geometry].vertex;
gd.nom = (BITRUNCATED ? tv+ts : tv) * 4;
gd.denom = (2*ts + 2*tv - ts * tv);
#if CAP_GP
if(GOLDBERG || INVERSE) {
gd.area = PIU(cgi.gpdata->area);
if(GOLDBERG || WARPED) {
gd.nom = 2 * (2*tv + (S3-2) * ts * (gd.area-1));
}
else if(UNRECTIFIED) {
if((gp::param.first + gp::param.second) % 2 == 0)
gd.nom = ts * 2 * gd.area;
else
gd.nom = (2*tv + (S3-2) * ts * (gd.area-1));
}
else if(UNTRUNCATED) {
if((gp::param.first - gp::param.second) % 3 == 0) {
gd.nom = ts * 4 * gd.area;
gd.denom *= 3;
}
else {
gd.nom = 2 * (2*tv + (S3-2) * ts * (gd.area-1));
gd.denom *= 3;
}
}
}
else
#endif
gd.area = PURE ? 1 : 3;
gd.euler = 0;
if(meuclid) gd.euler = 0;
else if(msphere && nonorientable) gd.euler = 1;
else if(msphere) gd.euler = 2;
else if(!closed_manifold) gd.euler = -2;
else if(WDIM == 3) gd.euler = 0;
else switch(geometry) {
case gFieldQuotient:
gd.worldsize = isize(currentmap->allcells());
gd.euler = 2 * gd.worldsize * gd.denom / gd.nom;
break;
case gMinimal:
gd.euler = -1;
break;
case gZebraQuotient:
case gBolza:
gd.euler = -2;
break;
case gKleinQuartic:
case gSchmutzM2:
case gBolza2:
gd.euler = -4;
break;
case gSchmutzM3:
case gBring:
gd.euler = -6;
break;
case gMacbeath:
gd.euler = -12;
break;
default:
gd.worldsize = isize(currentmap->allcells());
println(hlog, "warning: Euler characteristics unknown, worldsize = ", gd.worldsize);
gd.euler = 2 * gd.worldsize * gd.denom / gd.nom;
break;
}
gd.nom *= gd.euler;
gd.denom *= 2;
if(mhybrid) gd.nom *= hybrid::csteps, gd.denom *= cgi.single_step;
int g = gcd(gd.nom, gd.denom);
if(g) {
gd.nom /= g;
gd.denom /= g;
}
if(euclid && closed_manifold) {
gd.worldsize = euc::eu.det;
if(BITRUNCATED) gd.worldsize *= (a4 ? 2 : 3);
if(GOLDBERG) gd.worldsize *= cgi.gpdata->area;
#if CAP_IRR
if(IRREGULAR) gd.worldsize *= isize(irr::cells) / isize(irr::cells_of_heptagon);
#endif
}
else
gd.worldsize = gd.denom ? gd.nom / gd.denom : 0;
if(gd.euler < 0 && !closed_manifold)
gd.worldsize = -gd.worldsize;
string spf = its(ts);
if(0) ;
#if CAP_ARCM
else if(arcm::in()) {
spf = "";
for(int i: arcm::current.faces) {
if(spf != "") spf += ",";
spf += its(i);
}
if(BITRUNCATED) spf = "[" + spf + "]," + its(arcm::current.N * 2) + "," + its(arcm::current.N * 2);
if(DUAL) spf = its(arcm::current.N) + "^[" + spf + "]";
}
#endif
#if CAP_BT
else if(bt::in()) switch(geometry) {
case gBinaryTiling: spf = "6,[6,7],7"; break;
case gBinary4: spf = "5,5,5[,5]"; break;
case gTernary: spf = "6,6,6[,6]"; break;
default: spf = "?"; break;
}
#endif
else if(BITRUNCATED)
spf = spf + "," + its(S6) + "," + its(S6);
#if CAP_IRR
else if(IRREGULAR && irr::bitruncations_performed)
spf = "[4..8],6,6";
else if(IRREGULAR)
spf = "[4..8]^3";
#endif
#if CAP_GP
else if(UNRECTIFIED || UNTRUNCATED) {
if(UNRECTIFIED && (gp::param.first-gp::param.second) % 2 != 0)
spf = "(?)";
else if(UNTRUNCATED && (gp::param.first-gp::param.second) % 3 != 0)
spf = "(?)";
else {
string each = UNRECTIFIED ? "4" : "3";
int first = UNRECTIFIED ? 4 : 6;
int second = S7;
if(gp::param == gp::loc(1, 1))
first = second;
else if(second < first)
swap(first, second);
spf = each;
for(int z=1; z<first; z++) spf += "," + each;
if(first != second) {
spf += "[";
for(int z=first; z<second; z++) spf += "," + each;
spf += "]";
}
}
}
else if(GOLDBERG && S3 == 4 && gp::param == gp::loc(1, 1))
spf = spf + ",4," + spf + ",4";
else if(GOLDBERG && S3 == 4 && gp::param == gp::loc(2, 0))
spf = spf + ",4,4,4";
else if(GOLDBERG && S3 == 4)
spf = "[" + spf + ",4],4,4,4";
else if(WARPED && S3 == 3 && gp::param == gp::loc(1,1))
spf = spf + ",3,3";
else if(WARPED && S3 == 3)
spf = "[" + spf + ",6],3,3";
else if(GOLDBERG && S3 == 3)
spf = "[" + spf + ",6],6,6";
#endif
else {
string spf0 = spf;
for(int z=1; z<S3; z++) spf = spf + "," + spf0;
}
gd.size_str =
disksize ? its(isize(currentmap->allcells())) :
#if CAP_BT
bt::in() ? fts(1440._deg * sqrt(2) * log(2) / pow(vid.binary_width, WDIM-1), 4) + " exp(∞)" :
#endif
#if CAP_ARCM
arcm::in() && (WDIM == 2) ? arcm::current.world_size() :
(arcm::in() && sphere) ? its(isize(currentmap->allcells())) :
#endif
#if CAP_CRYSTAL
cryst ? "∞^" + its(ts/2) :
#endif
WDIM == 3 && closed_manifold ? its(isize(currentmap->allcells())) :
WDIM == 3 && euclid ? "∞" :
gd.worldsize < 0 ? (gd.nom%gd.denom ? its(gd.nom)+"/"+its(gd.denom) : its(-gd.worldsize)) + " exp(∞)":
(euclid && quotient && !closed_manifold) ? "∞" :
gd.worldsize == 0 ? "∞²" :
its(gd.worldsize);
#if CAP_IRR
if(hyperbolic && IRREGULAR) {
gd.nom = isize(irr::cells);
// both Klein Quartic and Bolza2 are double the Zebra quotiennt
gd.denom = -2;
if(!quotient) gd.worldsize = gd.nom / gd.denom;
}
#endif
if(WDIM == 3) gd.euler = 0;
gd.demigenus = 2 - gd.euler;
return gd;
}
EX void add_size_action() {
dialog::add_action([] {
if(!viewdists) { enable_viewdists(); pushScreen(viewdist_configure_dialog); }
else if(viewdists) viewdists = false;
});
}
EX void showEuclideanMenu() {
// for(int i=2; i<lt; i++) landvisited[i] = true;
cmode = sm::SIDE | sm::MAYDARK;
gamescreen();
dialog::init(XLAT("experiment with geometry"));
dialog::addSelItem(XLAT("geometry"), geometry_name(), 'd');
dialog::add_action([] { pushScreen(ge_select_tiling); pushScreen(ge_select_filter); });
dialog::addSelItem(XLAT("basic tiling"), XLAT(ginf[geometry].tiling_name), 't');
dialog::add_action([] {
if(!current_filter || !current_filter->test()) set_default_filter();
pushScreen(ge_select_tiling);
});
string qstring = ginf[geometry].quotient_name;
if(meuclid && quotient) qstring = ONOFF(true);
if(qstring == "none")
dialog::addBoolItem(XLAT("quotient space"), false, 'q');
else
dialog::addSelItem(XLAT("quotient space"), XLAT(qstring), 'q');
dialog::add_action(select_quotient);
#if CAP_ARCM
if(arcm::in()) {
dialog::addItem(XLAT("advanced parameters"), '4');
dialog::add_action_push(arcm::show);
}
#endif
#if CAP_CRYSTAL
if(cryst) {
dialog::addItem(XLAT("advanced parameters"), '4');
dialog::add_action_push(crystal::show);
}
#endif
if(fake::available()) {
dialog::addItem(XLAT("fake curvature"), '4');
dialog::add_action([] {
if(fake::in()) fake::configure();
else dialog::cheat_if_confirmed(
fake::configure
);
});
}
if(arb::in() && !(arb::current.sliders.empty() && arb::current.intsliders.empty())) {
dialog::addItem(XLAT("tessellation sliders"), '4');
dialog::add_action_push(arb::set_sliders);
}
if(arb::is_apeirogonal()) add_edit(arb::apeirogon_simplified_display);
if(arb::in()) add_edit(debug_tiles);
#if MAXMDIM >= 4
if(cgflags & qULTRA) {
dialog::addBoolItem(XLAT("truncate ultra-vertices with mirrors"), reg3::ultra_mirror_on, 'Z');
dialog::add_action([] {
reg3::ultra_mirror_on = !reg3::ultra_mirror_on;
ray::reset_raycaster();
});
}
#endif
if(mproduct) {
dialog::addSelItem(XLAT("Z-level height factor"), fts(vid.plevel_factor), 'Z');
dialog::add_action([] {
dialog::editNumber(vid.plevel_factor, 0, 2, 0.1, 0.7, XLAT("Z-level height factor"), "");
dialog::get_ne().reaction = ray::reset_raycaster;
});
}
else if(mtwisted) {
dialog::addSelItem(XLAT("number of levels"), its(hybrid::csteps / cgi.single_step), 'L');
dialog::add_action(hybrid::configure_period);
if(nil) menuitem_nilwidth('v');
}
else if(bt::in()) {
menuitem_binary_width('v');
add_edit_wall_quality('W');
}
else if(hat::in()) {
add_edit(hat::hat_param);
add_edit(hat::hat_param_imag);
}
else if(nil) {
menuitem_nilwidth('v');
}
else if((WDIM == 3 || aperiodic || arb::in()) && !reg3::in() && geometry != gCubeTiling) dialog::addBreak(100);
else if(cgclass & qFRACTAL)
dialog::addBreak(100);
else
menuitem_change_variation('v');
if(in_s2xe()) {
dialog::addSelItem(XLAT("precision of S2xE rings"), its(s2xe::qrings), '5');
dialog::add_action([] {
dialog::editNumber(s2xe::qrings, 1, 256, 4, 32, XLAT("precision of S2xE rings"),
XLAT(
"In S2xE, objects at spherical distances which are multiples of π will look like "
"rings, and objects close to these will look like crescents. "
"This setting controls the quality of rendering these rings and crescents.")
);
dialog::bound_low(1);
dialog::bound_up(256);
});
}
#if MAXMDIM >= 4
if(mhybrid) {
auto r = hybrid::underlying_scale;
dialog::addSelItem(XLAT("view the underlying geometry"), r > 0 ? fts(r)+"x" : ONOFF(false), '6');
dialog::add_action([] {
dialog::editNumber(hybrid::underlying_scale, 0, 1, 0.05, 0.25, XLAT("view the underlying geometry"),
geometry == gTwistedProduct ?
XLAT("The space you are currently in a twisted product space. ")
: XLAT("You are currently in a product space.") +
XLAT(
"This option lets you see the underlying space. Lands and some walls (e.g. in the Graveyard) are based on "
"the respective features in the underlying world, but details such as monsters or items are ignored."
)
);
dialog::bound_low(0);
dialog::bound_up(1);
dialog::get_di().extra_options = [] () { hybrid::draw_underlying(true); };
});
}
#endif
if(stretch::applicable()) {
dialog::addSelItem(XLAT("stretched geometry"), fts(stretch::factor), 'S');
dialog::add_action(edit_stretch);
}
dialog::addBreak(100);
menuitem_land_structure('l');
add_edit(req_disksize);
if(specialland == laMinefield && closed_or_bounded) {
dialog::addSelItem(XLAT("number of mines"), its(bounded_mine_quantity), 'm');
dialog::add_action([] {
dialog::editNumber(bounded_mine_quantity, 0, bounded_mine_max, 1, (bounded_mine_max+5)/10,
XLAT("number of mines"), "");
dialog::get_ne().reaction = [] {
if(bounded_mine_quantity < 0) bounded_mine_quantity = 0;
if(bounded_mine_quantity > bounded_mine_max) bounded_mine_quantity = bounded_mine_max;
};
dialog::get_ne().reaction_final = [] {
bounded_mine_percentage = bounded_mine_quantity * 1. / bounded_mine_max;
stop_game();
start_game();
};
});
}
if(geometry_has_alt_mine_rule()) {
dialog::addSelItem(XLAT("adjacency rule"), mine_adjacency_rule ? XLAT("vertex") : WDIM == 3 ? XLAT("face") : XLAT("edge"), 'M');
dialog::add_action([] {
stop_game();
mine_adjacency_rule = !mine_adjacency_rule;
start_game();
addMessage(XLAT("Note: adjacency rule affects environmental effects, but not movement."));
});
}
dialog::addBoolItem(XLAT("pattern"), specialland == laCanvas, 'p');
if(specialland == laCanvas) dialog::lastItem().value = ccolor::which->name;
dialog::add_action_push(patterns::showPrePattern);
validity_info();
if(WDIM == 3) {
dialog::addItem(XLAT("3D configuration"), '9');
dialog::add_action_push(show3D);
}
menuitem_projection('1');
if(lie_movement_available()) {
add_edit(nisot::geodesic_movement);
}
#if CAP_CRYSTAL && MAXMDIM >= 4
crystal::add_crystal_transform('x');
#endif
dialog::addBreak(50);
#if CAP_SHOT
dialog::addItem(XLAT("take screenshot"), 's');
dialog::add_action_push(shot::menu);
#endif
dialog::addHelp();
dialog::addBack();
dialog::addBreak(150);
dialog::addTitle(XLAT("info about: %1", full_geometry_name()), 0xFFFFFF, 150);
auto gd = compute_geometry_data();
if(WDIM == 2 && !arb::in() && !aperiodic) dialog::addSelItem(XLAT("faces per vertex"), gd.spf, 0);
if(arb::in() && arb::current.comment != "") {
dialog::addBreak(100);
dialog::addHelp(arb::current.comment);
}
if(WDIM == 2 && quotient && closed_manifold) {
dialog::addItem(XLAT("fundamental domain"), 'F');
dialog::add_action_push(fundamental::showMenu);
}
dialog::addSelItem(XLAT("size of the world"), gd.size_str, '3');
add_size_action();
if(closed_manifold) {
dialog::addSelItem(XLAT("Euler characteristics"), its(gd.euler), 0);
if(WDIM == 3) ;
else if(nonorientable)
dialog::addSelItem(XLAT("demigenus"), its(gd.demigenus), 0);
else
dialog::addSelItem(XLAT("genus"), its(gd.demigenus/2), 0);
}
else dialog::addBreak(200);
dialog::display();
}
EX void runGeometryExperiments() {
if(!geometry && specialland == laIce)
specialland = getLandForList(cwt.at);
pushScreen(showEuclideanMenu);
}
#if CAP_COMMANDLINE
EX eGeometry readGeo(const string& ss) {
for(int i=0; i<isize(ginf); i++) if(ginf[i].shortname == ss) return eGeometry(i);
bool numeric = true;
for(char c: ss) if(c < '0' || c > '9') numeric = false;
if(numeric) return eGeometry(atoi(ss.c_str()));
for(int i=0; i<isize(ginf); i++) if(appears(ginf[i].menu_displayed_name, ss)) {
return eGeometry(i);
break;
}
return gNormal;
}
EX void field_quotient_3d(int p, unsigned hash) {
check_cgi();
cgi.require_basics();
stop_game_and_switch_mode(rg::nothing);
fieldpattern::field_from_current();
set_geometry(gFieldQuotient);
for(;; p++) {
println(hlog, "trying p = ", p);
currfp.Prime = p; currfp.force_hash = hash;
if(!currfp.solve()) break;
}
println(hlog, "set prime = ", currfp.Prime);
}
EX void field_quotient_3d(string code) {
check_cgi();
cgi.require_basics();
stop_game_and_switch_mode(rg::nothing);
fieldpattern::field_from_current();
set_geometry(gFieldQuotient);
shstream ins(code);
hread_fpattern(ins, currfp);
}
EX void field_quotient_2d(int group, int id, int triplet) {
using namespace fieldpattern;
current_extra = group;
auto& gxcur = fgeomextras[current_extra];
while(id >= isize(gxcur.primes)) nextPrime(gxcur);
fgeomextras[current_extra].current_prime_id = id;
if(triplet != -1)
triplet_id = triplet;
enableFieldChange();
set_geometry(gFieldQuotient);
}
int read_geom_args() {
using namespace arg;
if(0) ;
#if CAP_FIELD
else if(argis("-qpar")) {
int p;
shift(); sscanf(argcs(), "%d,%d,%d",
&p, &quotientspace::rvadd, &quotientspace::rvdir
);
autocheat = true;
currfp.init(p);
}
else if(argis("-qpar2")) {
stop_game_and_switch_mode(rg::nothing);
int a, b;
shift(); sscanf(argcs(), "%d,%d", &a, &b);
field_quotient_2d(a, b, -1);
}
else if(argis("-triplet")) {
stop_game();
shift(); fieldpattern::triplet_id = argi();
fieldpattern::enableFieldChange();
}
else if(argis("-to-fq")) {
int p = 2;
shift();
if(args() == "p") {
shift(); p = argi();
shift();
}
unsigned hash = arghex();
field_quotient_3d(p, hash);
}
else if(argis("-cs")) {
shift(); cheat();
fieldpattern::matrix M = currfp.strtomatrix(args());
fieldpattern::subpathid = currfp.matcode[M];
fieldpattern::subpathorder = currfp.order(M);
}
else if(argis("-fwrite")) {
shstream hs;
hwrite_fpattern(hs, currfp);
println(hlog, "current fieldpattern hash: ", currfp.hashv);
println(hlog, "current fieldpattern: ", as_cstring(hs.s));
}
else if(argis("-csp")) {
cheat();
currfp.findsubpath();
}
#endif
else if(argis("-mineadj")) {
shift(); mine_adjacency_rule = argi();
}
TOGGLE('7', PURE, set_variation(PURE ? eVariation::bitruncated : eVariation::pure))
else if(argis("-geo")) {
PHASEFROM(2);
shift();
set_geometry(readGeo(args()));
}
#if CAP_GP
else if(argis("-gp")) {
PHASEFROM(2);
shift(); gp::param.first = argi();
shift(); gp::param.second = argi();
set_variation(eVariation::goldberg);
}
else if(argis("-gph")) {
PHASEFROM(2);
shift(); int x = argi();
shift(); int y = argi();
gp::whirl_set({x, y});
}
else if(argis("-unrectified")) {
PHASEFROM(2);
gp::param = gp::univ_param();
set_variation(eVariation::unrectified);
}
else if(argis("-untruncated")) {
PHASEFROM(2);
gp::param = gp::univ_param();
set_variation(eVariation::untruncated);
}
else if(argis("-warped")) {
PHASEFROM(2);
gp::param = gp::univ_param();
set_variation(eVariation::warped);
}
#if MAXMDIM >= 4
else if(argis("-subcubes")) {
PHASEFROM(2);
stop_game();
set_variation(eVariation::subcubes);
shift(); reg3::subcube_count = argi();
}
else if(argis("-dual-subcubes")) {
PHASEFROM(2);
stop_game();
set_variation(eVariation::dual_subcubes);
shift(); reg3::subcube_count = argi();
}
else if(argis("-bch-subcubes")) {
PHASEFROM(2);
stop_game();
set_variation(eVariation::bch);
shift(); reg3::subcube_count = argi();
}
else if(argis("-bch-oct")) {
PHASEFROM(2);
stop_game();
set_variation(eVariation::bch_oct);
shift(); reg3::subcube_count = argi();
}
else if(argis("-coxeter")) {
PHASEFROM(2);
stop_game();
set_variation(eVariation::coxeter);
shift(); reg3::coxeter_param = argi();
}
#endif
#endif
#if CAP_FIELD
else if(argis("-fi")) {
fieldpattern::info();
exit(0);
}
else if(argis("-fi-at")) {
geometry = gNormal;
shift(); dynamicval<int> s7(S7, argi());
shift(); dynamicval<int> s3(S3, argi());
fieldpattern::info();
exit(0);
}
else if(argis("-fi-geo")) {
fieldpattern::info();
exit(0);
}
else if(argis("-qs")) {
cheat();
shift(); currfp.qpaths.push_back(args());
}
#if MAXMDIM >= 4
else if(argis("-truncate-ultra")) {
shift(); reg3::ultra_mirror_on = argi();
}
#endif
else if(argis("-d:quotient"))
launch_dialog(WDIM == 2 ? showQuotientConfig : showQuotientConfig3);
else if(argis("-uqf"))
fieldpattern::use_quotient_fp = true;
#endif
else if(argis("-d:geom"))
launch_dialog(showEuclideanMenu);
else return 1;
return 0;
}
auto ah_geom = addHook(hooks_args, 0, read_geom_args);
#endif
}
Reference in New Issue View Git Blame Copy Permalink