two extra projections in RogueViz

This commit is contained in:
Zeno Rogue 2023-03-16 15:18:02 +01:00
parent eb867ecccc
commit 2370084b11
2 changed files with 59 additions and 1 deletions

View File

@ -211,7 +211,7 @@ bool hypot_zlev(ld zlev, ld& d, ld& df, ld& zf) {
int twopoint_sphere_flips;
bool twopoint_do_flips;
ld find_zlev(hyperpoint& H) {
EX ld find_zlev(hyperpoint& H) {
if(spatial_graphics) {
ld zlev = zlevel(H);

View File

@ -31,6 +31,13 @@ template<class T> void add_complex(const char *name, flagtype flag, const T& f)
});
}
template<class T> void add_extra(const char *name, flagtype flag, const T& f) {
int q = isize(mdinf);
mdinf.emplace_back(modelinfo{name, name, name, mf::euc_boring | mf::broken | flag, 0, 0, 0, 0, 0, nullptr});
while(isize(extra_projections) < q) extra_projections.emplace_back();
extra_projections.emplace_back(f);
}
template<class T> void add_band(const char *name, flagtype flag, const T& f) {
int q = isize(mdinf);
mdinf.emplace_back(modelinfo{name, name, name, mf::euc_boring | mf::broken | flag, 0, 0, 0, 0, 0, nullptr});
@ -181,6 +188,57 @@ void add_extra_projections() {
x = ax;
y = ay;
});
add_extra("double horocyclic", 0, [] (shiftpoint& H_orig, hyperpoint& H, hyperpoint& ret) {
ret = H;
if(hyperbolic) {
find_zlev(H);
models::apply_orientation_yz(H[1], H[2]);
models::apply_orientation(H[0], H[1]);
hyperpoint h1 = deparabolic13(H);
ret[1] = h1[1];
h1 = deparabolic13(cspin90(0, 1) * H);
ret[0] = -h1[1];
if(GDIM == 3) {
h1 = deparabolic13(cspin90(2, 1) * H);
ret[0] = -h1[2];
}
}
ret[LDIM] = 1;
if(hyperbolic) {
models::apply_orientation(ret[1], ret[0]);
models::apply_orientation_yz(ret[2], ret[1]);
}
});
add_extra("azimuthal cylindrical", 0, [] (shiftpoint& H_orig, hyperpoint& H, hyperpoint& ret) {
find_zlev(H);
models::apply_orientation_yz(H[1], H[2]);
models::apply_orientation(H[0], H[1]);
ld x, y;
y = asin_auto(H[1]);
x = asin_auto_clamp(H[0] / cos_auto(y));
if(sphere) {
if(H[LDIM] < 0 && x > 0) x = M_PI - x;
else if(H[LDIM] < 0 && x <= 0) x = -M_PI - x;
}
x += H_orig.shift;
ret[0] = x;
ret[1] = H[1] / H[0] * x;
if(GDIM == 2) ret[2] = H[2] / H[0] * x;
ret[LDIM] = 1;
});
}
void gen_dual() {