mirror of
https://github.com/zenorogue/hyperrogue.git
synced 2024-12-23 16:50:27 +00:00
1067 lines
31 KiB
C++
1067 lines
31 KiB
C++
#include "hyper.h"
|
|
|
|
namespace hr {
|
|
|
|
EX namespace intra {
|
|
|
|
EX bool in;
|
|
|
|
#if MAXMDIM >= 4
|
|
#if HDR
|
|
/** information per every space connected with intra-portals */
|
|
struct intra_data {
|
|
gamedata gd;
|
|
geometryinfo gi;
|
|
int wallindex;
|
|
};
|
|
#endif
|
|
|
|
EX vector<intra_data> data;
|
|
|
|
/** index of the space we are currently in */
|
|
EX int current;
|
|
|
|
/** portal debugging flags */
|
|
EX int debug_portal;
|
|
|
|
/** map cells to their intra spaces */
|
|
EX map<cell*, int> intra_id;
|
|
|
|
#if HDR
|
|
/** information about portal (one side) */
|
|
struct portal_data {
|
|
int kind;
|
|
hyperpoint v0;
|
|
ld d;
|
|
transmatrix T;
|
|
transmatrix iT;
|
|
hyperpoint co0;
|
|
hyperpoint co1;
|
|
ld scale;
|
|
/* convert h to portal coordinates ('poco') to usual coordinates */
|
|
hyperpoint to_poco(hyperpoint h) const;
|
|
/* convert h from portal coordinates ('poco') to usual coordinates */
|
|
hyperpoint from_poco(hyperpoint h) const;
|
|
};
|
|
#endif
|
|
|
|
hyperpoint portal_data::to_poco(hyperpoint h) const {
|
|
if(prod && kind == 1) {
|
|
auto dec = product_decompose(h);
|
|
h = dec.second;
|
|
if(bt::in()) {
|
|
h = PIU( deparabolic13(h) );
|
|
h[2] = dec.first - d;
|
|
h[3] = 1;
|
|
}
|
|
else {
|
|
h = T * h;
|
|
h[0] /= h[2];
|
|
h[1] /= h[2];
|
|
h[2] = dec.first - d;
|
|
h[3] = 1;
|
|
}
|
|
if(d<0) h[2] = -h[2], h[0] = -h[0];
|
|
return h;
|
|
}
|
|
else if(prod && kind == 0) {
|
|
h = T * h;
|
|
ld z = product_decompose(h).first;
|
|
h /= exp(z);
|
|
auto h1 = h;
|
|
h[2] = asin_auto_clamp(h1[0]);
|
|
h[1] = z;
|
|
h[0] = asin_auto_clamp(h1[1] / cos_auto(h[2]));
|
|
h[3] = 1;
|
|
return h;
|
|
}
|
|
#if CAP_BT
|
|
else if(hyperbolic && bt::in()) {
|
|
h = deparabolic13(h);
|
|
h[3] = 1;
|
|
tie(h[0], h[1], h[2]) = make_tuple(h[1], h[2], h[0]);
|
|
h = T * h;
|
|
h[2] *= exp(h[1]);
|
|
return h;
|
|
}
|
|
#endif
|
|
else if(sol) {
|
|
h = T * h;
|
|
h[2] *= exp(-h[1]);
|
|
return h;
|
|
}
|
|
else {
|
|
h = T * h;
|
|
h /= h[3];
|
|
if(sphere)
|
|
h[2] /= sqrt(1+h[0]*h[0]+h[1]*h[1]);
|
|
if(hyperbolic)
|
|
h[2] /= sqrt(1-h[0]*h[0]-h[1]*h[1]);
|
|
return h;
|
|
}
|
|
}
|
|
|
|
hyperpoint portal_data::from_poco(hyperpoint h) const {
|
|
if(prod && kind == 1) {
|
|
ld xd = h[2];
|
|
if(d<0) xd = -xd, h[0] = -h[0];
|
|
#if CAP_BT
|
|
if(bt::in()) {
|
|
h[2] = 0;
|
|
return PIU( parabolic13(h) ) * exp(d+xd);
|
|
}
|
|
#endif
|
|
h[2] = 1;
|
|
auto z = product_decompose(h).first;
|
|
return iT * h * exp(d+xd-z);
|
|
}
|
|
else if(prod && kind == 0) {
|
|
auto h0 = h;
|
|
h[0] = sin_auto(h0[2]);
|
|
h[1] = sin_auto(h0[0]) * cos_auto(h0[2]);
|
|
h[2] = cos_auto(h0[0]) * cos_auto(h0[2]);
|
|
h[3] = 1;
|
|
return iT * h * exp(h0[1]);
|
|
}
|
|
#if CAP_BT
|
|
else if(hyperbolic && bt::in()) {
|
|
h[2] *= exp(-h[1]);
|
|
h = iT * h;
|
|
return hr::parabolic13(h[0], h[1]) * xpush0(h[2]);
|
|
}
|
|
else if(sol) {
|
|
h[2] *= exp(h[1]);
|
|
return iT * h;
|
|
}
|
|
#endif
|
|
else {
|
|
h[3] = 1;
|
|
if(sphere)
|
|
h[2] *= sqrt(1+h[0]*h[0]+h[1]*h[1]);
|
|
if(hyperbolic)
|
|
h[2] *= sqrt(1-h[0]*h[0]-h[1]*h[1]);
|
|
return normalize(iT * h);
|
|
}
|
|
}
|
|
|
|
EX portal_data make_portal(cellwalker cw, int spin) {
|
|
if(debug_portal & 289)
|
|
println(hlog, "working in: ", full_geometry_name(), " wall no ", cw.spin, "/", cw.at->type);
|
|
auto& ss = currentmap->get_cellshape(cw.at);
|
|
auto fac = ss.faces_local[cw.spin];
|
|
portal_data id;
|
|
id.scale = 1;
|
|
id.T = Id;
|
|
auto gg = geometry;
|
|
if(prod && cw.spin >= cw.at->type - 2) {
|
|
id.kind = 1;
|
|
id.d = product_decompose(fac[0]).first;
|
|
id.v0 = C0 * exp(id.d);
|
|
#if CAP_BT
|
|
if(bt::in()) {
|
|
for(auto h: fac)
|
|
println(hlog, PIU(deparabolic13(normalize_flat(h))));
|
|
if(cw.spin == cw.at->type - 2)
|
|
fac.pop_back();
|
|
else
|
|
fac.erase(fac.begin() + 1);
|
|
id.scale = log(2)/2;
|
|
}
|
|
#else
|
|
if(false) {}
|
|
#endif
|
|
else {
|
|
hyperpoint ctr = Hypc;
|
|
for(auto p: fac) ctr += product_decompose(p).second;
|
|
ctr = normalize_flat(ctr);
|
|
id.T = gpushxto0(ctr);
|
|
}
|
|
}
|
|
else if(prod) {
|
|
id.kind = 0;
|
|
id.v0 = Hypc;
|
|
id.scale = cgi.plevel;
|
|
for(auto p: fac) id.v0 += p;
|
|
id.v0 = normalize_flat(id.v0);
|
|
hyperpoint h = normalize_flat(fac[0]);
|
|
id.T = cspin(0, 1, -90*degree) * spintox(gpushxto0(id.v0) * h) * gpushxto0(id.v0);
|
|
if((id.T * C0)[0] > 0) id.T = cspin(0, 1, 180*degree) * id.T;
|
|
for(int i=0; i<3; i++) id.T[3][i] = id.T[i][3] = i==3;
|
|
if(debug_portal & 128)
|
|
for(int a=0; a<4; a++) {
|
|
hyperpoint h = fac[a];
|
|
println(hlog, kz(h), " -> ", kz(spintox(id.v0)*h), " -> ", kz(cpush(0, -hdist0(id.v0))) * kz(spintox(id.v0) * h), " -> ", kz(id.to_poco(h)));
|
|
}
|
|
}
|
|
#if CAP_BT
|
|
else if(bt::in()) {
|
|
hyperpoint removed = Hypc;
|
|
|
|
auto facmod = fac;
|
|
if(hyperbolic) for(auto& h: facmod) h = deparabolic13(h);
|
|
|
|
for(int i=0; i<isize(facmod); i++) {
|
|
int i1 = i+1; if(i1 >= isize(facmod)) i1 = 0;
|
|
int i2 = i1+1; if(i2 >= isize(facmod)) i2 = 0;
|
|
if(hypot_d(3, 2*facmod[i1] - facmod[i] - facmod[i2]) < 1e-3) {
|
|
removed = fac[i1];
|
|
facmod.erase(facmod.begin()+i1);
|
|
fac.erase(fac.begin()+i1);
|
|
}
|
|
}
|
|
id.kind = 0;
|
|
id.v0 = Hypc;
|
|
id.T = Id;
|
|
|
|
auto fac1 = fac;
|
|
auto to_coords = [] (hyperpoint& p) {
|
|
if(hyperbolic) {
|
|
p = deparabolic13(p);
|
|
p = hyperpoint(p[1], p[2], p[0], 1);
|
|
}
|
|
};
|
|
for(auto& p: fac1)
|
|
to_coords(p);
|
|
to_coords(removed);
|
|
|
|
for(auto p: fac1) id.v0 += p;
|
|
id.v0 /= isize(fac);
|
|
dynamicval<eGeometry> g(geometry, gCubeTiling);
|
|
id.T = gpushxto0(id.v0);
|
|
for(auto p: fac1) {
|
|
if(abs((id.T * p)[2]) > 1e-3 && abs((id.T * p)[0]) < 1e-3)
|
|
id.T = cspin(2, 0, 90*degree) * id.T;
|
|
if(abs((id.T * p)[2]) > 1e-3 && abs((id.T * p)[1]) < 1e-3)
|
|
id.T = cspin(2, 1, 90*degree) * id.T;
|
|
}
|
|
if((id.T * C03)[2] > 0) id.T = cspin(2, 0, 180*degree) * id.T;
|
|
if(abs((id.T * removed)[0]) > 1e-2) id.T = cspin(0, 1, 90*degree) * id.T;
|
|
if((id.T * removed)[1] < -1e-2) id.T = cspin(0, 1, 180*degree) * id.T;
|
|
vector<hyperpoint> v;
|
|
geometry = gg;
|
|
for(auto f: fac) v.push_back(id.to_poco(f));
|
|
geometry = gCubeTiling;
|
|
ld sca = 1;
|
|
for(int i=0; i<isize(v); i++)
|
|
sca *= sqhypot_d(3, v[i] - v[(1+i) % isize(v)]);
|
|
sca = pow(sca, .5 / isize(v));
|
|
id.scale = sca / 2;
|
|
}
|
|
#endif
|
|
else {
|
|
id.kind = 0;
|
|
id.v0 = project_on_triangle(fac[0], fac[1], fac[2]);
|
|
id.T = cpush(2, -hdist0(id.v0)) * cspin(2, 0, 90*degree) * spintox(id.v0);
|
|
hyperpoint ctr = Hypc;
|
|
for(auto p: fac) ctr += id.T*p;
|
|
ctr = normalize(ctr);
|
|
id.T = gpushxto0(ctr) * id.T;
|
|
}
|
|
if(MDIM == 3) for(int i=0; i<4; i++) id.T[3][i] = id.T[i][3] = i==3;
|
|
id.iT = inverse(id.T);
|
|
if(MDIM == 3) for(int i=0; i<4; i++) id.iT[3][i] = id.iT[i][3] = i==3;
|
|
int first = spin;
|
|
int second = spin + 1;
|
|
first = gmod(first, isize(fac));
|
|
second = gmod(second, isize(fac));
|
|
id.co0 = id.to_poco(fac[first]);
|
|
id.co1 = id.to_poco(fac[second]);
|
|
|
|
if(debug_portal & 32) {
|
|
for(int i=0; i<isize(fac); i++)
|
|
println(hlog, "edge ", i, " length is ", hdist(fac[i], fac[(i+1)%isize(fac)]));
|
|
println(hlog, "chosen edge is ", first, "--", second);
|
|
}
|
|
|
|
if(debug_portal & 256) {
|
|
println(hlog, "portal scale = ", id.scale);
|
|
auto res = [&] (ld x, ld y, ld z) {
|
|
hyperpoint h = hyperpoint(x, y, z, 1);
|
|
return id.from_poco(h);
|
|
};
|
|
for(int x=0; x<5; x++) {
|
|
println(hlog, "horizontal ", x, " = ", hdist(res(x*.1,0,0), res(x*.1+.001,0,0)));
|
|
println(hlog, "vertical ", x, " = ", hdist(res(x*.1,0,0), res(x*.1,0.001,0)));
|
|
println(hlog, "deep ", x, " = ", hdist(res(x*.1,0,0), res(x*.1,0,0.001)));
|
|
}
|
|
hyperpoint a = hyperpoint(.4, .2, .1, 1);
|
|
println(hlog, "a = ", a);
|
|
println(hlog, "b = ", id.from_poco(a));
|
|
println(hlog, "c = ", id.to_poco(id.from_poco(a)));
|
|
}
|
|
|
|
if(debug_portal & 1) {
|
|
for(auto p: fac) {
|
|
auto p2 = id.to_poco(p);
|
|
auto p3 = id.from_poco(p2);
|
|
println(hlog, kz(p), " > ", kz(p2), " > ", kz(p3));
|
|
}
|
|
println(hlog, kz(C0), " > ", kz(id.to_poco(C0)), " > ", kz(id.from_poco(id.to_poco(C0))));
|
|
}
|
|
|
|
return id;
|
|
}
|
|
|
|
#if HDR
|
|
/** information about connection (portal-to-portal) */
|
|
struct connection_data {
|
|
int source_world;
|
|
int target_world;
|
|
cellwalker scw, tcw;
|
|
portal_data id1;
|
|
portal_data id2;
|
|
transmatrix T;
|
|
int spin_value;
|
|
bool mirrored; /* not implemented */
|
|
};
|
|
#endif
|
|
|
|
EX map<cellwalker, connection_data> connections;
|
|
|
|
EX connection_data* find_connection(int a, int b) {
|
|
for(auto& p: connections)
|
|
if(intra_id.at(p.first.at) == a && p.second.target_world == b)
|
|
return &p.second;
|
|
return nullptr;
|
|
}
|
|
|
|
EX void switch_to(int id) {
|
|
if(current == id) return;
|
|
data[current].gd.storegame();
|
|
current = id;
|
|
ginf[gProduct] = data[current].gi;
|
|
data[current].gd.restoregame();
|
|
}
|
|
|
|
void connect_portal_1(cellwalker cw1, cellwalker cw2, int spin) {
|
|
auto& p = connections[cw1];
|
|
p.source_world = intra_id.at(cw1.at);
|
|
p.target_world = intra_id.at(cw2.at);
|
|
p.scw = cw1;
|
|
p.tcw = cw2;
|
|
switch_to(intra_id.at(cw1.at));
|
|
int pspin = 0, nspin = 0;
|
|
if(spin > 0) pspin = spin; else nspin = -spin;
|
|
p.id1 = make_portal(cw1, nspin);
|
|
switch_to(intra_id.at(cw2.at));
|
|
p.id2 = make_portal(cw2, pspin);
|
|
p.spin_value = spin;
|
|
|
|
if(1) {
|
|
dynamicval<eGeometry> g(geometry, gCubeTiling);
|
|
transmatrix T1;
|
|
set_column(T1, 0, p.id1.co0);
|
|
set_column(T1, 1, p.id1.co1);
|
|
set_column(T1, 2, hyperpoint(0,0,p.id1.scale,0));
|
|
set_column(T1, 3, C03);
|
|
transmatrix T2;
|
|
set_column(T2, 0, p.id2.co0);
|
|
set_column(T2, 1, p.id2.co1);
|
|
set_column(T2, 2, hyperpoint(0,0,-p.id2.scale,0));
|
|
set_column(T2, 3, C03);
|
|
if(debug_portal & 2) for(int i=0; i<4; i++)
|
|
println(hlog, "mapping [", p.source_world, "]", get_column(T1, i), " to [", p.target_world, "] ", get_column(T2, i),
|
|
" dists = ", hypot_d(2, get_column(T1,i)), ",", hypot_d(2, get_column(T2,i)));
|
|
p.T = T2 * inverse(T1);
|
|
|
|
if(debug_portal & 2)
|
|
println(hlog, "det = ", det(p.T));
|
|
|
|
if(det(p.T) < 0) {
|
|
set_column(T2, 0, p.id2.co1);
|
|
set_column(T2, 1, p.id2.co0);
|
|
p.T = T2 * inverse(T1);
|
|
}
|
|
}
|
|
|
|
if(debug_portal & 2)
|
|
println(hlog, "got scale = ", tie(p.id1.scale, p.id2.scale));
|
|
|
|
if(debug_portal & 4) for(int i=0; i<5; i++) {
|
|
hyperpoint h = C03;
|
|
if(i == 1) h[0] += .2;
|
|
if(i == 2) h[0] -= .2;
|
|
if(i == 3) h[1] += .2;
|
|
if(i == 4) h[1] -= .2;
|
|
array<hyperpoint, 4> hl;
|
|
ld eps = 1e-5; ld ss = pow(eps, -2);
|
|
hl[0] = h;
|
|
hl[1] = h + point3(eps, 0, 0);
|
|
hl[2] = h + point3(0, eps, 0);
|
|
hl[3] = h + point3(0, 0, eps);
|
|
auto hl1 = hl;
|
|
may_switch_to(cw1.at);
|
|
println(hlog, "checking ", h, " -> L ", kz(p.id1.from_poco(h)), " g: ", full_geometry_name());
|
|
for(auto& hx: hl1) hx = p.id1.from_poco(hx);
|
|
print(hlog, "L side: "); analyze_orthonormal(hl1, ss);
|
|
may_switch_to(cw2.at);
|
|
println(hlog, "checking ", h, " -> R ", kz(p.id2.from_poco(p.T * h)), " g: ", full_geometry_name());
|
|
auto hl2 = hl;
|
|
for(auto& hx: hl2) hx = p.id2.from_poco(p.T * hx);
|
|
print(hlog, "R side: "); analyze_orthonormal(hl2, ss);
|
|
}
|
|
|
|
}
|
|
|
|
EX vector<pair<int, cell*>> full_sample_list;
|
|
|
|
EX void connect_portal(cellwalker cw1, cellwalker cw2, int spin) {
|
|
connect_portal_1(cw1, cw2, spin);
|
|
connect_portal_1(cw2, cw1, -spin);
|
|
}
|
|
|
|
EX void generate_sample_list_for_current() {
|
|
auto v = hybrid::gen_sample_list();
|
|
int gi = 0;
|
|
if(full_sample_list.size()) {
|
|
gi = full_sample_list.back().first;
|
|
full_sample_list.pop_back();
|
|
}
|
|
data[current].wallindex = gi;
|
|
for(auto x: v)
|
|
full_sample_list.emplace_back(x.first + gi, x.second);
|
|
println(hlog, "added ", isize(v)-1, " samples, wallindex = ", data[current].wallindex);
|
|
}
|
|
|
|
EX void regenerate_full_sample_list() {
|
|
resetter ir;
|
|
full_sample_list.clear();
|
|
for(int i=0; i<isize(data); i++) {
|
|
switch_to(i);
|
|
generate_sample_list_for_current();
|
|
}
|
|
println(hlog, isize(full_sample_list), " samples known");
|
|
}
|
|
|
|
/** make currentmap into one of the spaces in intra */
|
|
EX void become() {
|
|
if(intra::in) {
|
|
/* let them add more spaces in this case */
|
|
data[current].gd.storegame();
|
|
intra::in = false;
|
|
return;
|
|
}
|
|
|
|
check_cgi();
|
|
cgi.require_shapes();
|
|
auto& ac = currentmap->allcells();
|
|
current = isize(data);
|
|
for(cell *c: ac)
|
|
intra_id[c] = current;
|
|
for(cell *c: ac)
|
|
currentmap->wall_offset(c);
|
|
for(cell *c: ac) c->item = itNone;
|
|
data.emplace_back();
|
|
data.back().gd.storegame();
|
|
data.back().gi = ginf[gProduct];
|
|
generate_sample_list_for_current();
|
|
sightranges[geometry] = 10;
|
|
}
|
|
|
|
/** after called become() on some spaces, actually start intra */
|
|
EX void start(int id IS(0)) {
|
|
in = true;
|
|
current = id;
|
|
data[current].gd.restoregame();
|
|
ginf[gProduct] = data[current].gi;
|
|
|
|
again:
|
|
int missing = 0;
|
|
for(auto p: intra_id)
|
|
for(int i=0; i<p.first->type; i++) {
|
|
cell *c1 = p.first->move(i);
|
|
if(!c1) continue;
|
|
if(intra_id.count(c1) == 0) {
|
|
intra_id[c1] = p.second;
|
|
missing++;
|
|
}
|
|
}
|
|
if(debug_portal & 64) println(hlog, "missing = ", missing);
|
|
if(missing) goto again;
|
|
}
|
|
|
|
#if HDR
|
|
/** a convenience struct to switch back after a temporary switch_to */
|
|
struct resetter {
|
|
int ic;
|
|
resetter() { ic = current; }
|
|
~resetter() { if(in) switch_to(ic); }
|
|
};
|
|
#endif
|
|
|
|
EX void may_switch_to(cell *c) {
|
|
if(in) switch_to(intra_id.at(c));
|
|
}
|
|
|
|
EX int full_wall_offset(cell *c) {
|
|
int wo = currentmap->wall_offset(c);
|
|
if(in) wo += data[intra_id.at(c)].wallindex;
|
|
return wo;
|
|
}
|
|
|
|
ld dsdet(array<hyperpoint, 4> ds) {
|
|
transmatrix T;
|
|
set_column(T, 0, ds[1]-ds[0]);
|
|
set_column(T, 1, ds[2]-ds[0]);
|
|
set_column(T, 2, ds[3]-ds[0]);
|
|
return det3(T);
|
|
}
|
|
|
|
EX void analyze_orthonormal(array<hyperpoint, 4> ds, ld sca) {
|
|
transmatrix T = gpushxto0(ds[0]);
|
|
vector<ld> orths;
|
|
for(int i: {1,2,3}) {
|
|
ds[i] = T * ds[i];
|
|
if(prod) ds[i][2]--;
|
|
}
|
|
for(int i=0; i<3; i++)
|
|
for(int j=0; j<3; j++)
|
|
orths.push_back(dot_d(3, ds[i+1], ds[j+1]) * sca);
|
|
println(hlog, "orths = ", kz(orths));
|
|
}
|
|
|
|
EX void shift_view_portal(hyperpoint H) {
|
|
shift_view(H);
|
|
if(!through_portal()) return;
|
|
shift_view(-H);
|
|
ld minv = 0, maxv = 1;
|
|
for(int i=0; i<30; i++) {
|
|
ld t = (minv + maxv) / 2;
|
|
shift_view(H * t);
|
|
bool b = through_portal();
|
|
if(b) maxv = t; else minv = t;
|
|
shift_view(H * -t);
|
|
}
|
|
println(hlog, "maxv = ", maxv);
|
|
shift_view(H * maxv);
|
|
check_portal_movement();
|
|
shift_view_portal(H * (1 - maxv));
|
|
}
|
|
|
|
EX const connection_data* through_portal() {
|
|
transmatrix iView = view_inverse(View);
|
|
ld dist = hdist0(iView * C0);
|
|
int nei = -1;
|
|
for(int i=0; i<centerover->type; i++) {
|
|
ld dist1 = hdist0(currentmap->ray_iadj(centerover, i) * iView * C0);
|
|
if(dist1 < dist) nei = i, dist = dist1;
|
|
}
|
|
|
|
auto cw1 = cellwalker(centerover, nei);
|
|
return at_or_null(connections, cw1);
|
|
}
|
|
|
|
EX void check_portal_movement() {
|
|
auto p = through_portal();
|
|
|
|
if(p) {
|
|
ld eps = 1e-5;
|
|
ld ss = pow(eps, -2);
|
|
|
|
array<hyperpoint, 4> ds; /* camera, forward, upward */
|
|
ds[0] = inverse(View) * C0;
|
|
ds[1] = inverse(get_shift_view_of(ctangent(2, -eps), View)) * C0;
|
|
ds[2] = inverse(get_shift_view_of(ctangent(1, +eps), View)) * C0;
|
|
ds[3] = inverse(get_shift_view_of(ctangent(0, +eps), View)) * C0;
|
|
if(debug_portal & 8) {
|
|
#if CAP_BT
|
|
println(hlog, "at = ", ds[0], " det = ", dsdet(ds), " bt = ", bt::minkowski_to_bt(ds[0]));
|
|
#endif
|
|
analyze_orthonormal(ds, ss);
|
|
}
|
|
|
|
for(auto& h: ds) h = p->id1.to_poco(h);
|
|
|
|
if(debug_portal & 8) {
|
|
println(hlog, "poco: at = ", ds[0], " det = ", dsdet(ds));
|
|
if(debug_portal & 16) {
|
|
dynamicval<eGeometry> g(geometry, gCubeTiling);
|
|
analyze_orthonormal(ds, ss);
|
|
}
|
|
}
|
|
|
|
/* reset the original */
|
|
View = Id; NLP = Id;
|
|
|
|
switch_to(p->target_world);
|
|
centerover = p->tcw.at;
|
|
|
|
if(1) {
|
|
dynamicval<eGeometry> g(geometry, gCubeTiling);
|
|
for(auto& h: ds) h = p->T * h;
|
|
}
|
|
|
|
if(debug_portal & 8) {
|
|
println(hlog, "poco2: at = ", ds[0], " det = ", dsdet(ds));
|
|
if(debug_portal & 16) {
|
|
dynamicval<eGeometry> g(geometry, gCubeTiling);
|
|
analyze_orthonormal(ds, ss);
|
|
}
|
|
}
|
|
|
|
for(auto& h: ds) h = p->id2.from_poco(h);
|
|
|
|
if(debug_portal & 8) {
|
|
println(hlog, "goal: at = ", ds[0], " det = ", dsdet(ds));
|
|
analyze_orthonormal(ds, ss);
|
|
}
|
|
|
|
set_view(ds[0], ds[1], ds[2]);
|
|
|
|
if(debug_portal & 8) {
|
|
array<hyperpoint, 4> xds; /* camera, forward, upward */
|
|
xds[0] = inverse(View) * C0;
|
|
xds[1] = inverse(get_shift_view_of(ctangent(2, -eps), View)) * C0;
|
|
xds[2] = inverse(get_shift_view_of(ctangent(1, +eps), View)) * C0;
|
|
xds[3] = inverse(get_shift_view_of(ctangent(0, +eps), View)) * C0;
|
|
#if CAP_BT
|
|
println(hlog, "goal: at = ", xds[0], " det = ", dsdet(xds), " bt = ", bt::minkowski_to_bt(xds[0]));
|
|
#endif
|
|
}
|
|
|
|
ld scale = p->id2.scale / p->id1.scale;
|
|
|
|
camera_speed *= scale;
|
|
anims::cycle_length *= scale;
|
|
#if CAP_VR
|
|
vrhr::absolute_unit_in_meters *= scale;
|
|
#endif
|
|
if(walking::eye_level != -1) walking::eye_level *= scale;
|
|
|
|
walking::floor_dir = -1;
|
|
walking::on_floor_of = nullptr;
|
|
}
|
|
}
|
|
|
|
vector<cellwalker> unconnected;
|
|
|
|
void erase_unconnected(cellwalker cw) {
|
|
for(int i=0; i<isize(unconnected); i++)
|
|
if(unconnected[i] == cw)
|
|
unconnected.erase(unconnected.begin() + i);
|
|
}
|
|
|
|
int edit_spin;
|
|
|
|
EX void show_portals() {
|
|
gamescreen(1);
|
|
|
|
dialog::init(XLAT("manage portals"));
|
|
|
|
cellwalker cw(centerover, point_direction);
|
|
bool valid = point_direction >= 0 && point_direction < centerover->type;
|
|
|
|
dialog::addItem(XLAT("move to the next space"), 'm');
|
|
dialog::add_action([] {
|
|
int ic = (current + 1) % isize(data);
|
|
switch_to(ic);
|
|
});
|
|
|
|
dialog::addSelItem(XLAT("mode"), ray::fixed_map ? "perf" : "edit", 'e');
|
|
dialog::add_action([] { ray::fixed_map = !ray::fixed_map; });
|
|
|
|
if(debug_portal) {
|
|
dialog::addItem(XLAT("debug"), 'd');
|
|
dialog::add_action([] {
|
|
ld eps = 1e-5;
|
|
array<hyperpoint, 4> ds; /* camera, forward, upward */
|
|
ds[0] = inverse(View) * C0;
|
|
ds[1] = inverse(get_shift_view_of(ctangent(2, -eps), View)) * C0;
|
|
ds[2] = inverse(get_shift_view_of(ctangent(1, +eps), View)) * C0;
|
|
ds[3] = inverse(get_shift_view_of(ctangent(0, +eps), View)) * C0;
|
|
set_view(ds[0], ds[1], ds[2]);
|
|
});
|
|
}
|
|
|
|
bool in_list = false; for(cellwalker x: unconnected) if(x == cw) in_list = true;
|
|
|
|
if(!valid) ;
|
|
else if(connections.count(cw)) {
|
|
dialog::addItem(XLAT("disconnect this portal"), 'd');
|
|
dialog::add_action([cw] {
|
|
auto tcw = connections[cw].tcw;
|
|
unconnected.push_back(tcw);
|
|
connections.erase(cw);
|
|
connections.erase(tcw);
|
|
});
|
|
}
|
|
else if(in_list) {
|
|
dialog::addItem(XLAT("remove %1 from the list", lalign(0, cw)), 'r');
|
|
dialog::add_action([cw] {
|
|
erase_unconnected(cw);
|
|
});
|
|
}
|
|
else {
|
|
dialog::addItem(XLAT("add to list"), 'a');
|
|
dialog::add_action([cw] { unconnected.push_back(cw); });
|
|
for(auto p: unconnected) {
|
|
dialog::addItem(XLAT("connect " + lalign(0, p)), '1');
|
|
dialog::add_action([p, cw] {
|
|
connect_portal(cw, p, edit_spin);
|
|
erase_unconnected(p);
|
|
});
|
|
}
|
|
dialog::addSelItem(XLAT("portal orientation"), its(edit_spin), 'o');
|
|
dialog::add_action([] { edit_spin = edit_spin + 1; });
|
|
if(debug_portal) {
|
|
dialog::addItem(XLAT("mirror connection"), 'm');
|
|
dialog::add_action([cw] { connect_portal(cw, cw, edit_spin); });
|
|
}
|
|
if(debug_portal) {
|
|
dialog::addItem(XLAT("test portal here"), 't');
|
|
dialog::add_action([cw] { make_portal(cw, 0); });
|
|
}
|
|
}
|
|
|
|
walking::add_options();
|
|
|
|
dialog::display();
|
|
}
|
|
|
|
#if HDR
|
|
struct portal_to_save {
|
|
cellwalker cw1;
|
|
cellwalker cw2;
|
|
int spin;
|
|
bool mirrored;
|
|
};
|
|
#endif
|
|
|
|
EX vector<portal_to_save> portals_to_save;
|
|
|
|
EX void prepare_to_save() {
|
|
portals_to_save.clear();
|
|
for(auto c: connections) if(c.second.scw < c.second.tcw) {
|
|
portals_to_save.emplace_back(portal_to_save{c.second.scw, c.second.tcw, c.second.spin_value, false});
|
|
}
|
|
}
|
|
|
|
EX void load_saved_portals() {
|
|
for(const auto& p: portals_to_save) connect_portal(p.cw1, p.cw2, p.spin);
|
|
}
|
|
|
|
EX void be_ratio(ld v IS(1)) {
|
|
check_cgi();
|
|
cgi.require_basics();
|
|
PIU( vid.plevel_factor = v * cgi.edgelen / cgi.scalefactor );
|
|
check_cgi();
|
|
cgi.require_basics();
|
|
}
|
|
|
|
EX void be_ratio_edge(int i, ld v IS(1)) {
|
|
start_game();
|
|
ld len = hdist(currentmap->get_corner(cwt.at, i), currentmap->get_corner(cwt.at, (i+1)%cwt.at->type));
|
|
PIU( vid.plevel_factor = v * len / cgi.scalefactor );
|
|
check_cgi();
|
|
cgi.require_basics();
|
|
}
|
|
|
|
/** Remove the space with the given id. Turns off intra */
|
|
EX void kill(int id) {
|
|
if(in) become();
|
|
data.erase(data.begin()+id);
|
|
vector<cellwalker> to_remove;
|
|
for(auto& p: connections) if(intra_id[p.second.scw.at] == id || intra_id[p.second.tcw.at] == id)
|
|
to_remove.push_back(p.first);
|
|
else {
|
|
if(p.second.source_world >= id) p.second.source_world--;
|
|
if(p.second.target_world >= id) p.second.target_world--;
|
|
}
|
|
for(auto r: to_remove) connections.erase(r);
|
|
vector<cell*> to_erase_cell;
|
|
for(auto& p: intra_id)
|
|
if(p.second == id)
|
|
to_erase_cell.push_back(p.first);
|
|
else if(p.second > id)
|
|
p.second--;
|
|
for(auto c: to_erase_cell)
|
|
intra_id.erase(c);
|
|
println(hlog, isize(to_remove), " connections and ", isize(to_erase_cell), " cells erased");
|
|
}
|
|
|
|
EX void erase_all_maps() {
|
|
println(hlog, "erase_all_maps called");
|
|
data[current].gd.storegame();
|
|
in = false;
|
|
for(int i=0; i<isize(data); i++) {
|
|
current = i;
|
|
ginf[gProduct] = data[i].gi;
|
|
data[i].gd.restoregame();
|
|
clearCellMemory();
|
|
}
|
|
intra_id.clear();
|
|
connections.clear();
|
|
data.clear();
|
|
full_sample_list.clear();
|
|
}
|
|
|
|
EX set<cell*> need_to_save;
|
|
|
|
EX void prepare_need_to_save() {
|
|
need_to_save.clear();
|
|
map<cell*, cell*> parent;
|
|
vector<cell*> q;
|
|
cell *s = mapstream::save_start();
|
|
parent[s] = s;
|
|
q = {s};
|
|
for(int i=0; i<isize(q); i++) {
|
|
cell *c = q[i];
|
|
forCellEx(c2, c)
|
|
if(!parent.count(c2)) { parent[c2] = c; q.push_back(c2); }
|
|
}
|
|
for(int i=isize(q)-1; i>=0; i--) {
|
|
cell *c = q[i];
|
|
if(c == cwt.at) need_to_save.insert(c);
|
|
for(auto& p: connections) if(p.first.at == c) need_to_save.insert(c);
|
|
if(need_to_save.count(c)) need_to_save.insert(parent[c]);
|
|
}
|
|
println(hlog, "need to save ", isize(need_to_save), " out of ", isize(q), " cells");
|
|
}
|
|
|
|
auto hooks1 =
|
|
addHook(hooks_o_key, 90, [] (o_funcs& v) {
|
|
if(intra::in) v.push_back(named_dialog(XLAT("manage portals"), show_portals));
|
|
})
|
|
+ arg::add3("-intra-add", [] { start_game(); become(); })
|
|
+ arg::add3("-intra-start", [] { start_game(); become(); start(0); })
|
|
+ arg::add3("-intra-kill", [] { arg::shift(); kill(arg::argi()); start(0); regenerate_full_sample_list(); })
|
|
+ arg::add3("-be-square", [] { be_ratio(); })
|
|
+ arg::add3("-be-square-edge", [] {
|
|
arg::shift(); int i = arg::argi(); be_ratio_edge(i);
|
|
})
|
|
+ arg::add3("-debug-portal", [] { arg::shift(); debug_portal = arg::argi(); });
|
|
#endif
|
|
EX }
|
|
|
|
EX namespace walking {
|
|
|
|
EX bool on;
|
|
|
|
EX bool auto_eyelevel;
|
|
|
|
EX int floor_dir = -1;
|
|
EX cell *on_floor_of = nullptr;
|
|
EX ld eye_level = 0.2174492;
|
|
EX ld eye_angle = 0;
|
|
EX ld eye_angle_scale = 1;
|
|
|
|
int ticks_end, ticks_last;
|
|
|
|
EX set<color_t> colors_of_floors;
|
|
|
|
EX bool isFloor(cell *c) {
|
|
if(!isWall(c)) return false;
|
|
if(colors_of_floors.empty()) return true;
|
|
if(c->wall != waWaxWall) return false;
|
|
return colors_of_floors.count(c->landparam);
|
|
}
|
|
|
|
EX void handle() {
|
|
if(playermoved || !on) return;
|
|
|
|
if(floor_dir == -1 || on_floor_of != centerover) {
|
|
vector<int> choices;
|
|
for(int i=0; i<centerover->type; i++)
|
|
if(isFloor(centerover->cmove(i)))
|
|
choices.push_back(i);
|
|
|
|
if(sol && isize(choices) == 2) choices.pop_back();
|
|
|
|
if(isize(choices) == 1) {
|
|
on_floor_of = centerover;
|
|
floor_dir = choices[0];
|
|
}
|
|
else if(colors_of_floors.empty() && sn::in()) {
|
|
on_floor_of = centerover;
|
|
auto z = inverse(View) * C0;
|
|
switch(geometry) {
|
|
case gSol:
|
|
floor_dir = (z[2] > 0) ? 2 : 6;
|
|
return;
|
|
case gNIH:
|
|
floor_dir = (z[2] > 0) ? 5 : 4;
|
|
return;
|
|
case gSolN:
|
|
floor_dir = (z[2] > 0) ? 4 : 6;
|
|
return;
|
|
default: throw hr_exception("not solnihv");
|
|
}
|
|
}
|
|
#if CAP_BT
|
|
else if(colors_of_floors.empty() && hyperbolic && bt::in()) {
|
|
auto z = bt::minkowski_to_bt(inverse(View) * C0);
|
|
on_floor_of = centerover;
|
|
floor_dir = z[2] > 0 ? bt::updir() : 0;
|
|
println(hlog, "set floor_dir to ", floor_dir);
|
|
}
|
|
#endif
|
|
else {
|
|
println(hlog, "there are ", isize(choices), " choices for floor_dir");
|
|
if(!on_floor_of) return;
|
|
}
|
|
}
|
|
|
|
struct face {
|
|
hyperpoint h0, hx, hy;
|
|
};
|
|
|
|
transmatrix ToOld = currentmap->relative_matrix(on_floor_of, centerover, C0);
|
|
auto& csh = currentmap->get_cellshape(on_floor_of);
|
|
face f;
|
|
f.h0 = ToOld * csh.faces_local[floor_dir][0];
|
|
f.hx = ToOld * csh.faces_local[floor_dir][1];
|
|
f.hy = ToOld * csh.faces_local[floor_dir][2];
|
|
|
|
auto find_nearest = [&] (const face& fac, hyperpoint at) {
|
|
#if CAP_BT
|
|
if(sol) { at[2] = fac.h0[2]; return at; }
|
|
else if(hyperbolic && bt::in()) {
|
|
auto z = bt::minkowski_to_bt(at);
|
|
z[2] = bt::minkowski_to_bt(fac.h0)[2];
|
|
return bt::bt_to_minkowski(z);
|
|
}
|
|
else if(prod && bt::in()) {
|
|
auto dec = product_decompose(at);
|
|
hyperpoint dep = PIU( deparabolic13(dec.second) );
|
|
hyperpoint h = product_decompose(fac.h0).second;
|
|
h = PIU( deparabolic13(h) );
|
|
dep[0] = h[0];
|
|
return zshift(PIU(parabolic13(dep)), dec.first);
|
|
}
|
|
#else
|
|
if(false) {}
|
|
#endif
|
|
else {
|
|
transmatrix M = ray::mirrorize(currentmap->ray_iadj(on_floor_of, floor_dir));
|
|
M = ToOld * M * inverse(ToOld);
|
|
return mid(at, M * at);
|
|
}
|
|
};
|
|
|
|
hyperpoint at = tC0(inverse(View));
|
|
if(invalid_point(at)) {
|
|
println(hlog, "at is invalid!");
|
|
on = false;
|
|
return;
|
|
}
|
|
|
|
auto wallpt = find_nearest(f, at);
|
|
|
|
ld view_eps = 1e-5;
|
|
|
|
transmatrix spin_T;
|
|
bool use_T = false;
|
|
|
|
if(eye_angle) use_T = true, spin_T = cspin(1, 2, -eye_angle * degree);
|
|
#if CAP_VR
|
|
if(vrhr::active() && !vrhr::first && vrhr::hsm != vrhr::eHeadset::none) {
|
|
use_T = true;
|
|
spin_T = vrhr::hmd_ref_at;
|
|
// print(hlog, "HMD seems to be at altitude ", spin_T[1][3], " depth ", spin_T[2][3], " zeros are ", spin_T[3][1], " and ", spin_T[3][2]);
|
|
dynamicval<eGeometry> g(geometry, gCubeTiling);
|
|
spin_T = vrhr::sm * inverse(spin_T);
|
|
eye_level = -spin_T[1][3] / vrhr::absolute_unit_in_meters;
|
|
if(eye_level < .001) eye_level = 0.001;
|
|
vrhr::be_33(spin_T);
|
|
}
|
|
#endif
|
|
|
|
if(use_T) rotate_view(spin_T);
|
|
hyperpoint front = inverse(get_shift_view_of(ctangent(2, -view_eps), View)) * C0;
|
|
hyperpoint up = inverse(get_shift_view_of(ctangent(1, +view_eps), View)) * C0;
|
|
|
|
auto fwallpt = find_nearest(f, front);
|
|
|
|
transmatrix T = nonisotropic ? nisot::translate(wallpt, -1) : gpushxto0(wallpt);
|
|
hyperpoint dx = inverse_exp(shiftless(T * at));
|
|
|
|
transmatrix Tf = nonisotropic ? nisot::translate(fwallpt, -1) : gpushxto0(fwallpt);
|
|
hyperpoint dxf = inverse_exp(shiftless(Tf * front));
|
|
|
|
if(eye_level == -1) eye_level = hypot_d(3, dx);
|
|
|
|
auto smooth = [&] (hyperpoint h1, hyperpoint h2) {
|
|
if(ticks < ticks_end) {
|
|
ld last_t = ilerp(ticks_end-1000, ticks_end, ticks_last);
|
|
ld curr_t = ilerp(ticks_end-1000, ticks_end, ticks);
|
|
last_t = last_t * last_t * (3-2*last_t);
|
|
curr_t = curr_t * curr_t * (3-2*curr_t);
|
|
ld t = ilerp(last_t, 1, curr_t);
|
|
return lerp(h1, h2, t);
|
|
}
|
|
return h2;
|
|
};
|
|
|
|
auto oView = View;
|
|
set_view(
|
|
smooth(at, inverse(T) * direct_exp(dx / hypot_d(3, dx) * eye_level)),
|
|
smooth(front, inverse(Tf) * direct_exp(dxf / hypot_d(3, dxf) * eye_level)),
|
|
smooth(up, inverse(T) * direct_exp(dx / hypot_d(3, dx) * (eye_level + view_eps)))
|
|
);
|
|
if(use_T) rotate_view(inverse(spin_T));
|
|
|
|
auto nat = tC0(inverse(View));
|
|
if(invalid_point(nat)) {
|
|
println(hlog, "at is invalid after fixing!");
|
|
View = oView;
|
|
return;
|
|
}
|
|
|
|
ticks_last = ticks;
|
|
}
|
|
|
|
EX void add_options() {
|
|
dialog::addBoolItem("walking mode", on, 'w');
|
|
dialog::add_action([] {
|
|
on = !on;
|
|
if(on && auto_eyelevel) eye_level = -1;
|
|
floor_dir = -1;
|
|
on_floor_of = nullptr;
|
|
ticks_last = ticks;
|
|
ticks_end = ticks + 1000;
|
|
});
|
|
add_edit(eye_level);
|
|
add_edit(eye_angle);
|
|
if(point_direction >= 0 && point_direction < centerover->type) {
|
|
cell *c = centerover->move(point_direction);
|
|
if(c && c->wall == waWaxWall) {
|
|
color_t col = c->landparam;
|
|
dialog::addBoolItem("we are facing floor (color " + format("%06X", col) + ")", colors_of_floors.count(col), 'n');
|
|
dialog::add_action([col] {
|
|
if(colors_of_floors.count(col)) colors_of_floors.erase(col);
|
|
else colors_of_floors.insert(col);
|
|
});
|
|
}
|
|
}
|
|
}
|
|
|
|
auto a = addHook(hooks_configfile, 100, [] {
|
|
param_b(auto_eyelevel, "auto_eyelevel")
|
|
-> editable("keep eye level when walking enabled", 'L');
|
|
param_f(eye_level, "eye_level")
|
|
-> editable(0, 5, .1, "walking eye level",
|
|
"Distance from the floor to the eye in the walking mode, in absolute units. In VR this is adjusted automatically.",
|
|
'e')
|
|
->set_extra([] { add_edit(auto_eyelevel); });
|
|
param_f(eye_angle, "eye_angle")
|
|
-> editable(-90, 90, 15, "walking eye angle",
|
|
"0 = looking forward, 90 = looking upward. In VR this is adjusted automatically.",
|
|
'k')
|
|
->set_extra([] { add_edit(eye_angle_scale); });
|
|
param_f(eye_angle_scale, "eye_angle_scale")
|
|
-> editable(-2, 0, 2, "eye angle scale",
|
|
"1 = the angle can be changed with keyboard or mouse movements, 0 = the angle is fixed",
|
|
'k');
|
|
})
|
|
+ addHook(hooks_clearmemory, 40, [] { on_floor_of = nullptr; floor_dir = -1; })
|
|
+ arg::add3("-walk-on", [] {
|
|
on = true;
|
|
if(auto_eyelevel) eye_level = -1;
|
|
floor_dir = -1;
|
|
on_floor_of = nullptr;
|
|
ticks_last = ticks_end = ticks;
|
|
});
|
|
|
|
EX }
|
|
}
|