From 7fb5c1f442fd74f4a208ad13515b52b22bc9895c Mon Sep 17 00:00:00 2001 From: Florian Stecker Date: Mon, 29 Jun 2020 18:59:08 -0500 Subject: [PATCH] compute singular values in positive case --- hyperbolic.c | 13 +++++-------- 1 file changed, 5 insertions(+), 8 deletions(-) diff --git a/hyperbolic.c b/hyperbolic.c index 1c0ee19..bb13fc1 100644 --- a/hyperbolic.c +++ b/hyperbolic.c @@ -191,8 +191,8 @@ int main() gsl_matrix *coxeter_eigenvectors[3]; gsl_matrix *frame; workspace_t *ws; - int elements = 100; - int p = 5, q = 5, r = 5; + int elements = 2000; + int p = 4, q = 4, r = 4; group = malloc(elements*sizeof(groupelement_t)); matrices = malloc(elements*sizeof(gsl_matrix*)); @@ -257,15 +257,12 @@ int main() 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) continue; - LOOP(i) transformed[i] = apply(matrices[k], edge_midpoints[(i+2)%3]); - LOOP(i) transformed2[i] = apply(matrices[k], coxeter_repelling[i]); - 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"); + LOOP(i) transformed[i] = apply(matrices[k], coxeter_attracting[i]); +// draw_line(transformed[1], transformed[2], frame, "red"); } print_tex_footer();