simplified ggpushxto0

2020-09-22 15:51:37 +02:00 · 2020-09-22 15:51:37 +02:00 · 70d92bb157
parent f4a596ed87
commit 70d92bb157
1 changed files with 1 additions and 1 deletions
--- a/hyperpoint.cpp
+++ b/hyperpoint.cpp
@ -803,7 +803,7 @@ EX transmatrix ggpushxto0(const hyperpoint& H, ld co) {
    }
  transmatrix res = Id;
  if(sqhypot_d(GDIM, H) < 1e-12) return res;
-  ld fac = (H[LDIM]-1) / sqhypot_d(GDIM, H);
+  ld fac = -curvature()/(H[LDIM]+1);
  for(int i=0; i<GDIM; i++)
  for(int j=0; j<GDIM; j++)
    res[i][j] += H[i] * H[j] * fac;