mirrors/hyperrogue
1
0
Fork 0
mirror of https://github.com/zenorogue/hyperrogue.git synced 2025-10-30 21:42:59 +00:00
Code Issues Releases Wiki Activity

fixed lie_exp to use correct rules in solnih, and also avoid division by h[2] == 0

Browse Source
This commit is contained in:
Zeno Rogue
2022-05-17 09:37:32 +02:00
parent b9fd64e1d9
commit 69a82a8ce1
1 changed files with 17 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
hypgraph.cpp
View File

@@ -2957,25 +2957,31 @@ EX hyperpoint lie_exp(hyperpoint h) {
h[3] = 1; h[3] = 1;
h[2] += h[0] * h[1] / 2; h[2] += h[0] * h[1] / 2;
} }
else if(sol) { else if(sol && !nih) {
h[3] = 1; h[3] = 1;
if(abs(h[2]) > 1e-6) {
h[0] *= (exp(-h[2]) - 1) / -h[2]; h[0] *= (exp(-h[2]) - 1) / -h[2];
h[1] *= (exp(+h[2]) - 1) / h[2]; h[1] *= (exp(+h[2]) - 1) / h[2];
} }
}
else if(sol && nih) { else if(sol && nih) {
h[3] = 1; h[3] = 1;
if(abs(h[2]) > 1e-6) {
ld z = h[2] * log(2); ld z = h[2] * log(2);
h[0] *= (exp(-z) - 1) / -z; h[0] *= (exp(-z) - 1) / -z;
z = h[2] * log(3); z = h[2] * log(3);
h[1] *= (exp(+z) - 1) / z; h[1] *= (exp(+z) - 1) / z;
} }
}
else if(nih) { else if(nih) {
h[3] = 1; h[3] = 1;
if(abs(h[2]) > 1e-6) {
ld z = h[2] * log(2); ld z = h[2] * log(2);
h[0] *= (exp(+z) - 1) / z; h[0] *= (exp(+z) - 1) / z;
z = h[2] * log(3); z = h[2] * log(3);
h[1] *= (exp(+z) - 1) / z; h[1] *= (exp(+z) - 1) / z;
} }
}
else { else {
/* not implemented -- approximate for now */ /* not implemented -- approximate for now */
const int steps = 16; const int steps = 16;