From 96f84c77cafdc6e9761e1e17f3bb2a008f489769 Mon Sep 17 00:00:00 2001 From: Zeno Rogue Date: Wed, 23 Dec 2020 23:42:09 +0100 Subject: [PATCH] random_spin3 is now a separate function --- hyperpoint.cpp | 14 ++++++++------ 1 file changed, 8 insertions(+), 6 deletions(-) diff --git a/hyperpoint.cpp b/hyperpoint.cpp index 4afac8d0..5013557f 100644 --- a/hyperpoint.cpp +++ b/hyperpoint.cpp @@ -550,14 +550,16 @@ EX transmatrix cspin(int a, int b, ld alpha) { /** rotate by alpha degrees in the XY plane */ EX transmatrix spin(ld alpha) { return cspin(0, 1, alpha); } +EX transmatrix random_spin3() { + ld alpha2 = asin(randd() * 2 - 1); + ld alpha = randd() * 2 * M_PI; + ld alpha3 = randd() * 2 * M_PI; + return cspin(0, 1, alpha) * cspin(0, 2, alpha2) * cspin(1, 2, alpha3); + } + EX transmatrix random_spin() { if(WDIM == 2) return spin(randd() * 2 * M_PI); - else { - ld alpha2 = asin(randd() * 2 - 1); - ld alpha = randd() * 2 * M_PI; - ld alpha3 = randd() * 2 * M_PI; - return cspin(0, 1, alpha) * cspin(0, 2, alpha2) * cspin(1, 2, alpha3); - } + else return random_spin3(); } EX transmatrix eupush(ld x, ld y) {