fixed sorting
This commit is contained in:
		
							
								
								
									
										4
									
								
								Makefile
									
									
									
									
									
								
							
							
						
						
									
										4
									
								
								Makefile
									
									
									
									
									
								
							@@ -1,8 +1,8 @@
 | 
				
			|||||||
HEADERS=triangle.h linalg.h queue.h
 | 
					HEADERS=triangle.h linalg.h queue.h
 | 
				
			||||||
 | 
					
 | 
				
			||||||
SPECIAL_OPTIONS=-O0 -g -D_DEBUG -Wall
 | 
					#SPECIAL_OPTIONS=-O0 -g -D_DEBUG
 | 
				
			||||||
#SPECIAL_OPTIONS=-O3 -pg -funroll-loops -fno-inline
 | 
					#SPECIAL_OPTIONS=-O3 -pg -funroll-loops -fno-inline
 | 
				
			||||||
#SPECIAL_OPTIONS=-O3 -flto -funroll-loops -Winline
 | 
					SPECIAL_OPTIONS=-O3 -flto -funroll-loops -Winline
 | 
				
			||||||
#SPECIAL_OPTIONS=
 | 
					#SPECIAL_OPTIONS=
 | 
				
			||||||
 | 
					
 | 
				
			||||||
OPTIONS=-m64 -march=native -mtune=native -std=gnu99 -D_GNU_SOURCE $(SPECIAL_OPTIONS)
 | 
					OPTIONS=-m64 -march=native -mtune=native -std=gnu99 -D_GNU_SOURCE $(SPECIAL_OPTIONS)
 | 
				
			||||||
 
 | 
				
			|||||||
@@ -8,6 +8,7 @@
 | 
				
			|||||||
#define POINCARE 1
 | 
					#define POINCARE 1
 | 
				
			||||||
#define LOOP(i) for(int i = 0; i < 3; i++)
 | 
					#define LOOP(i) for(int i = 0; i < 3; i++)
 | 
				
			||||||
 | 
					
 | 
				
			||||||
 | 
					
 | 
				
			||||||
void cartan_matrix(gsl_matrix *cartan, double a1, double a2, double a3, double s)
 | 
					void cartan_matrix(gsl_matrix *cartan, double a1, double a2, double a3, double s)
 | 
				
			||||||
{
 | 
					{
 | 
				
			||||||
	gsl_matrix_set(cartan, 0, 0, -2);
 | 
						gsl_matrix_set(cartan, 0, 0, -2);
 | 
				
			||||||
@@ -192,7 +193,7 @@ int main()
 | 
				
			|||||||
	gsl_matrix *frame;
 | 
						gsl_matrix *frame;
 | 
				
			||||||
	workspace_t *ws;
 | 
						workspace_t *ws;
 | 
				
			||||||
	int elements = 5000;
 | 
						int elements = 5000;
 | 
				
			||||||
	int p = 3, q = 5, r = 7;
 | 
						int p = 3, q = 3, r = 5;
 | 
				
			||||||
 | 
					
 | 
				
			||||||
	group = malloc(elements*sizeof(groupelement_t));
 | 
						group = malloc(elements*sizeof(groupelement_t));
 | 
				
			||||||
	matrices = malloc(elements*sizeof(gsl_matrix*));
 | 
						matrices = malloc(elements*sizeof(gsl_matrix*));
 | 
				
			||||||
 
 | 
				
			|||||||
		Reference in New Issue
	
	Block a user