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

compute singular values in positive case

Browse Source
This commit is contained in:
Florian Stecker 2020-06-29 18:59:08 -05:00
parent 498634f220
commit 7fb5c1f442
1 changed files with 5 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
hyperbolic.c
View File

@ -191,8 +191,8 @@ int main()
gsl_matrix *coxeter_eigenvectors[3]; gsl_matrix *coxeter_eigenvectors[3];
gsl_matrix *frame; gsl_matrix *frame;
workspace_t *ws; workspace_t *ws;
int elements = 100; int elements = 2000;
int p = 5, q = 5, r = 5; int p = 4, q = 4, r = 4;
group = malloc(elements*sizeof(groupelement_t)); group = malloc(elements*sizeof(groupelement_t));
matrices = malloc(elements*sizeof(gsl_matrix*)); matrices = malloc(elements*sizeof(gsl_matrix*));
@ -257,15 +257,12 @@ int main()
draw_triangle(transformed, frame, "black,fill=black!10,line width=0pt"); draw_triangle(transformed, frame, "black,fill=black!10,line width=0pt");
} }
for(int k = 0; k < 10; k++) { for(int k = 0; k < 2000; k++) {
if(group[k].length % 2) if(group[k].length % 2)
continue; continue;
LOOP(i) transformed[i] = apply(matrices[k], edge_midpoints[(i+2)%3]); LOOP(i) transformed[i] = apply(matrices[k], coxeter_attracting[i]);
LOOP(i) transformed2[i] = apply(matrices[k], coxeter_repelling[i]); // draw_line(transformed[1], transformed[2], frame, "red");
LOOP(i) transformed3[i] = apply(matrices[k], coxeter_repelling[(i+1)%3]);
draw_line(transformed[1], transformed2[1], frame, "red");
draw_line(transformed[1], transformed3[1], frame, "red");
} }
print_tex_footer(); print_tex_footer();