mirrors/hyperrogue
1
0
Fork 0
mirror of https://github.com/zenorogue/hyperrogue.git synced 2026-01-25 12:23:33 +00:00
Code Issues Releases Wiki Activity

moved precise_hdist to main

Browse Source
This commit is contained in:
Zeno Rogue
2025-12-05 01:39:39 +01:00
parent 1e30a539de
commit 7338a510d3
2 changed files with 21 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
hyperpoint.cpp
View File

@@ -1395,6 +1395,27 @@ EX ld hdist(const shiftpoint& h1, const shiftpoint& h2) {
return hdist(h1.h, unshift(h2, h1.shift));
}
EX ld precise_hdist(hyperpoint vi, hyperpoint vj) {
int n = MDIM-1;
hassert(n == 2 || n == 3);
ld da = acosh(vi[n]);
ld db = acosh(vj[n]);
ld rs = sqhypot_d(n, vi) * sqhypot_d(n, vj);
if(!rs) return da + db;
ld cosphi = 0;
for(int i=0; i<n; i++) cosphi += vi[i] * vj[i];
cosphi /= sqrt(rs);
ld co = sinh(da) * sinh(db) * (1 - cosphi);
ld v = cosh(da - db) + co;
if(v < 1) return 0;
return acosh(v);
}
/** like orthogonal_move but fol may be factor (in 2D graphics) or level (elsewhere) */
EX hyperpoint orthogonal_move_fol(const hyperpoint& h, double fol) {
if(GDIM == 2) return scale_point(h, fol);

16
rogueviz/dhrg/loglik.cpp
View File

@@ -195,20 +195,4 @@ template<class T> auto parallelize(long long N, T action) -> decltype(action(0,0
#endif
}
ld precise_hdist(hyperpoint vi, hyperpoint vj) {
ld da = acosh(vi[2]);
ld db = acosh(vj[2]);
ld phia = atan2(vi[0], vi[1]);
ld phib = atan2(vj[0], vj[1]);
ld co = sinh(da) * sinh(db) * (1 - cos(phia-phib));
// - (vi[0]*vj[0] + vi[1]*vj[1]);
ld v = cosh(da - db) + co;
if(v < 1) return 0;
return acosh(v);
}
}