florian/hyperbolic_tilings
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

add documentation

Browse Source
This commit is contained in:
Florian Stecker
2023-01-17 07:15:44 -05:00
parent c62044d637
commit 1b1880fe19
3 changed files with 65 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files

31
hyperbolic.c
View File

@@ -291,24 +291,23 @@ int main(int argc, const char *argv[])
nspecial_hyp = nspecial_rot = 0;
for(int i = 0; i < nspecial; i++) {
LOOP(j) {
int nreal;
int j = 0;
int nreal;
compute_word(ws, special[3*i+j], gen, argv[i+5], j, 0);
nreal = eigenvectors(special[3*i+j], special_eigenvectors[3*i+j], ws);
if(nreal == 3) {
special_attracting[nspecial_hyp] = column(special_eigenvectors[3*i+j], 0);
compute_word(ws, special[3*i+j], gen, argv[i+5], j, 0);
nreal = eigenvectors(special[3*i+j], special_eigenvectors[3*i+j], ws);
if(nreal == 3) {
special_attracting[nspecial_hyp] = column(special_eigenvectors[3*i+j], 0);
// repelling = attracting of inverse
compute_word(ws, special[3*i+j], gen, argv[i+5], j, 1);
eigenvectors(special[3*i+j], special_eigenvectors[3*i+j], ws);
special_repelling[nspecial_hyp] = column(special_eigenvectors[3*i+j], 0);
nspecial_hyp++;
} else {
special_rotation[nspecial_rot] = column(special_eigenvectors[3*i+j], 0);
nspecial_rot++;
}
}
// repelling = attracting of inverse
compute_word(ws, special[3*i+j], gen, argv[i+5], j, 1);
eigenvectors(special[3*i+j], special_eigenvectors[3*i+j], ws);
special_repelling[nspecial_hyp] = column(special_eigenvectors[3*i+j], 0);
nspecial_hyp++;
} else {
special_rotation[nspecial_rot] = column(special_eigenvectors[3*i+j], 0);
nspecial_rot++;
}
}
fprintf(stderr, "%d special elements, %d rotations, %d hyperbolic\n", nspecial, nspecial_rot, nspecial_hyp);