mirrors/hyperrogue
1
0
Fork 0
mirror of https://github.com/zenorogue/hyperrogue.git synced 2024-11-20 11:54:48 +00:00
Code Issues Releases Wiki Activity

fixed a bug in inverse caused by precision issues

Browse Source
This commit is contained in:
? 2019-02-25 04:05:56 +01:00 committed by Zeno Rogue
parent bbffc21a95
commit 3b030471b0
1 changed files with 13 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
hyperpoint.cpp
View File

@ -515,13 +515,18 @@ transmatrix inverse(const transmatrix& T) {
transmatrix T2 = Id; transmatrix T2 = Id;
for(int a=0; a<MDIM; a++) { for(int a=0; a<MDIM; a++) {
for(int b=a; b<=DIM; b++) int best = a;
if(T1[b][a]) {
if(b != a) for(int b=a+1; b<MDIM; b++)
for(int c=0; c<MDIM; c++) if(abs(T1[b][a]) > abs(T1[best][a]))
swap(T1[b][c], T1[a][c]), swap(T2[b][c], T2[a][c]); best = b;
break;
} int b = best;
if(b != a)
for(int c=0; c<MDIM; c++)
swap(T1[b][c], T1[a][c]), swap(T2[b][c], T2[a][c]);
if(!T1[a][a]) { inverse_error(T); return Id; } if(!T1[a][a]) { inverse_error(T); return Id; }
for(int b=a+1; b<=DIM; b++) { for(int b=a+1; b<=DIM; b++) {
ld co = -T1[b][a] / T1[a][a]; ld co = -T1[b][a] / T1[a][a];
@ -575,7 +580,7 @@ ld circlelength(ld r) {
// distance between two points // distance between two points
double hdist(const hyperpoint& h1, const hyperpoint& h2) { double hdist(const hyperpoint& h1, const hyperpoint& h2) {
return hdist0(gpushxto0(h1) * h2); // return hdist0(gpushxto0(h1) * h2);
ld iv = intval(h1, h2); ld iv = intval(h1, h2);
switch(cgclass) { switch(cgclass) {
case gcEuclid: case gcEuclid: