// File: /workspace/github_c_makefile_star2_O3/cbouilla[P]AES-attacks-finder/SysEqLin.c
// Function: TailleE
int TailleE(int * __restrict const X, const int N, SysEqLin E) // Renvoi N - dim(E(X)),  N = |X|
{
	//return TailleE(X,N,E);
	/* Beaucoup de variables */
	static int AlloueN = 0;
	static LigneCreuse * __restrict L = NULL;
	static int * __restrict nombre_nnz_colonne = NULL;
	static int AlloueCOL = 0;
	static int AlloueNCOL = 0;
	static int AlloueLIGNE_mat = 0;
	static unsigned char ** __restrict m = NULL;
	static int * __restrict nouvelle_colonne = NULL;
	static int * __restrict colonne_changer = NULL;
	//static int compteur1 = 0, compteur2 = 0, compteur3 = 0, compteur4 = 0;

	
	if (!E->initialise) InitialiseSousMat(E);
		
	if (N*E->gen>AlloueN)
	{
		free(L);
		AlloueN=N*E->gen;
		//printf("AlloueN : %d\n",AlloueN);
		L=malloc(AlloueN*sizeof(LigneCreuse));
	}

	if (E->col_max>AlloueCOL)
	{
		nouvelle_colonne=realloc(nouvelle_colonne,E->col_max*sizeof(int));
		int i;
		for (i = AlloueCOL; i < E->col_max; i++) nouvelle_colonne[i]=-1;
		free(nombre_nnz_colonne);
		AlloueCOL=E->col_max;
		//printf("AlloueCOL : %d\n",AlloueCOL);
		nombre_nnz_colonne=calloc(AlloueCOL,sizeof(int));
	}

	int NombreDeLignes =  0;
	int taille=0,i;
	for (i = 0; i < N; i++)
	{
		if (E->lignevar[X[i]]*E->gen<E->a)
		{
			int j;
			for (j = 0; j < E->sousmat[X[i]]->nombre_de_lignes; j++) L[NombreDeLignes++]=E->sousmat[X[i]]->ligne[j];
			taille++;
		}
	}

	/* 
	 * On compte le nombre d'éléments non-nuls sur chaque colonne
	 */
	int NombreDeColonnes = 0;
	for (i = 0; i < NombreDeLignes; i++)
	{
		int j;
		for (j = 0; j < L[i]->nnz; j++)
		{
			if (!nombre_nnz_colonne[L[i]->elements[j].colonne]) NombreDeColonnes++;
			nombre_nnz_colonne[L[i]->elements[j].colonne]++;
		}
	}

	/*
	 * On va essayer de faire une partie du pivot sur la matrice creuse
	 * (c'est rentable de part la forme des sous-mat)
	 */
	int dim = (1-E->gen)*taille;
	//int cpt=NombreDeLignes;
	int ligne = 0;
	//while (cpt)
	while (ligne<NombreDeLignes)
	{
		//if (ligne==NombreDeLignes) ligne=0;
		int nombre_element_ligne = 0;
		int colonne_a_enlever;
		int j;
		for (j = 0; j < L[ligne]->nnz; j++)
		{
			if (nombre_nnz_colonne[L[ligne]->elements[j].colonne])
			{
				if (nombre_nnz_colonne[L[ligne]->elements[j].colonne]==1) // 1 seul élément non nul sur la colonne -> très bon pivot pour matrice creuse
				{
					int l;
					for (l = 0; l < L[ligne]->nnz; l++)
					{
						if (nombre_nnz_colonne[L[ligne]->elements[l].colonne])
						{
							if (nombre_nnz_colonne[L[ligne]->elements[l].colonne] == 1) NombreDeColonnes--;
							nombre_nnz_colonne[L[ligne]->elements[l].colonne]--; // On met à jour le nombre d'éléments dans les colonnes
						}
					}
					goto suppr_ligne;
				}
				else
				{
					nombre_element_ligne++;
					colonne_a_enlever = j;
				}
			}
		}
		if (nombre_element_ligne<=1)
		{
			if (nombre_element_ligne==1) // 1 seul élément sur la ligne
			{
				nombre_nnz_colonne[L[ligne]->elements[colonne_a_enlever].colonne] = 0; // on supprime la colonne
				NombreDeColonnes--;
				suppr_ligne: // on met à jour la dim et on supprime la ligne
				dim++;
				//cpt=NombreDeLignes;
			}
			NombreDeLignes--;
			L[ligne]=L[NombreDeLignes];
		}
		else ligne++;
		//cpt--;
	}
	
	if (NombreDeColonnes>AlloueNCOL)
	{
		AlloueNCOL=NombreDeColonnes;
		for (i = 0; i < AlloueLIGNE_mat; i++)
		{
			free(m[i]);
			m[i]=malloc(AlloueNCOL*sizeof(unsigned char));
		}
		free(colonne_changer);
		colonne_changer=malloc((AlloueNCOL+1)*sizeof(int));
	}
	if (NombreDeLignes>AlloueLIGNE_mat)
	{
		m = realloc(m,NombreDeLignes*sizeof(unsigned char *));
		for (i = AlloueLIGNE_mat; i < NombreDeLignes; i++) m[i] = malloc(AlloueNCOL*sizeof(unsigned char));
		AlloueLIGNE_mat = NombreDeLignes;
		//printf("AlloueLIGNE_mat : %d\n",AlloueLIGNE_mat);
	}
	
	int pos=0;
	int DerniereColonne = NombreDeColonnes-1;
	int Ncol = 0;
	for (i = 0; i < NombreDeLignes; i++)
	{
		int k;
		unsigned char * __restrict ligne_matrice = m[i];
		for (k = 0; k < NombreDeColonnes; k++) ligne_matrice[k]=0;
		int premier = 1;
		for (k = 0; k < L[i]->nnz; k++)
		{
			if (nombre_nnz_colonne[L[i]->elements[k].colonne])
			{
				nombre_nnz_colonne[L[i]->elements[k].colonne]--;
				colonne_changer[pos]=L[i]->elements[k].colonne;
				int colonne;
				if (nouvelle_colonne[colonne_changer[pos]]>=0) colonne=nouvelle_colonne[colonne_changer[pos]];
				else
				{
					if (!premier)
					{
						nouvelle_colonne[colonne_changer[pos++]]=DerniereColonne;
						colonne=DerniereColonne--;
					}
					else
					{
						nouvelle_colonne[colonne_changer[pos++]]=Ncol;
						colonne=Ncol++;
						premier=0;
					}
				}
				ligne_matrice[colonne]=L[i]->elements[k].coeff;								
			}
		}
	}
			
	for (i = 0; i < pos; i++) nouvelle_colonne[colonne_changer[i]]=-1; // on remet le tableau à zero pour la prochaine execution

	//compteur4++;
	//printf("cpt1 = %d, cpt2 = %d, cpt3 = %d, cpt4 = %d\n",compteur1,compteur2,compteur3,compteur4);
															
	return dim+Dimension_tmp(m,NombreDeLignes,NombreDeColonnes);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x48, %rsp
movq %rdx, %r14
movl %esi, %r12d
movq %rdi, %r15
cmpb $0x0, 0x28(%rdx)
jne 0xa56e
movq %r14, %rdi
callq 0x9eab
movb $0x1, 0x28(%r14)
movl 0x10(%r14), %eax
imull %r12d, %eax
cmpl 0x27d88(%rip), %eax     # 0x32304
jle 0xa5ab
movq 0x27d83(%rip), %rdi     # 0x32308
callq 0x2030
movslq 0x10(%r14), %rax
movslq %r12d, %rdi
imulq %rax, %rdi
movl %edi, 0x27d69(%rip)     # 0x32304
shlq $0x3, %rdi
callq 0x21d0
movq %rax, 0x27d5d(%rip)     # 0x32308
movslq 0x38(%r14), %rsi
movslq 0x27d62(%rip), %rbx     # 0x32318
cmpl %ebx, %esi
jle 0xa61a
movq 0x27d6f(%rip), %rdi     # 0x32330
shlq $0x2, %rsi
callq 0x2200
movq %rax, 0x27d5f(%rip)     # 0x32330
movl 0x38(%r14), %ecx
cmpl %ecx, %ebx
jge 0xa5f3
leaq (%rax,%rbx,4), %rdi
notl %ebx
addl %ebx, %ecx
leaq 0x4(,%rcx,4), %rdx
movl $0xff, %esi
callq 0x2120
movq 0x27d16(%rip), %rdi     # 0x32310
callq 0x2030
movslq 0x38(%r14), %rdi
movl %edi, 0x27d0f(%rip)     # 0x32318
movl $0x4, %esi
callq 0x2160
movq %rax, 0x27cf6(%rip)     # 0x32310
xorl %ebx, %ebx
movl $0x0, %eax
movl $0x0, %r13d
testl %r12d, %r12d
jle 0xa7ce
movl 0x8(%r14), %ecx
movl 0x10(%r14), %edx
movq 0x27cc9(%rip), %rsi     # 0x32308
movl %r12d, %edi
xorl %r8d, %r8d
xorl %eax, %eax
xorl %r13d, %r13d
movq 0x18(%r14), %r10
movslq (%r15,%r8,4), %r9
movl (%r10,%r9,4), %r10d
imull %edx, %r10d
cmpl %ecx, %r10d
jge 0xa6a4
movq 0x30(%r14), %r10
movq (%r10,%r9,8), %r12
cmpl $0x0, (%r12)
jle 0xa6a2
movl %r13d, %ebp
movslq %r13d, %r10
leaq (%rsi,%r10,8), %r11
xorl %r10d, %r10d
movq 0x8(%r12), %r12
movq (%r12,%r10,8), %r12
movq %r12, (%r11,%r10,8)
incq %r10
movq 0x30(%r14), %r12
movq (%r12,%r9,8), %r12
movslq (%r12), %r13
cmpq %r13, %r10
jl 0xa67b
movl %ebp, %r13d
addl %r10d, %r13d
incl %eax
incq %r8
cmpq %rdi, %r8
jne 0xa64a
testl %r13d, %r13d
jle 0xa7ce
movq 0x27c4c(%rip), %rcx     # 0x32308
movq 0x27c4d(%rip), %rdx     # 0x32310
movl %r13d, %esi
xorl %edi, %edi
xorl %ebx, %ebx
movq (%rcx,%rdi,8), %r8
cmpl $0x0, (%r8)
jle 0xa6fb
movq 0x8(%r8), %r9
xorl %r10d, %r10d
movslq 0x4(%r9,%r10,8), %r11
movl (%rdx,%r11,4), %ebp
cmpl $0x1, %ebp
adcl $0x0, %ebx
incl %ebp
movl %ebp, (%rdx,%r11,4)
incq %r10
movslq (%r8), %r11
cmpq %r11, %r10
jl 0xa6db
incq %rdi
cmpq %rsi, %rdi
jne 0xa6ca
movl $0x1, %r15d
subl 0x10(%r14), %r15d
imull %eax, %r15d
testl %r13d, %r13d
jle 0xa7dc
movq 0x27be7(%rip), %rax     # 0x32308
xorl %ecx, %ecx
movq 0x27be6(%rip), %rdx     # 0x32310
movslq %ecx, %rdi
movq (%rax,%rdi,8), %r8
movl (%r8), %r10d
testl %r10d, %r10d
jle 0xa7b0
movq 0x8(%r8), %r9
xorl %r11d, %r11d
xorl %ebp, %ebp
movslq 0x4(%r9,%r11,8), %r14
movl (%rdx,%r14,4), %r14d
testl %r14d, %r14d
je 0xa75b
cmpl $0x1, %r14d
je 0xa782
incl %ebp
movl %r11d, %esi
incq %r11
cmpq %r11, %r10
jne 0xa742
cmpl $0x1, %ebp
jg 0xa7c1
jne 0xa7b0
movq 0x8(%r8), %r8
movslq %esi, %r9
movslq 0x4(%r8,%r9,8), %r8
movl $0x0, (%rdx,%r8,4)
decl %ebx
jmp 0xa7ad
xorl %r11d, %r11d
movslq 0x4(%r9,%r11,8), %r14
movl (%rdx,%r14,4), %ebp
testl %ebp, %ebp
je 0xa7a2
cmpl $0x1, %ebp
jne 0xa799
decl %ebx
decl %ebp
movl %ebp, (%rdx,%r14,4)
movl (%r8), %r10d
incq %r11
movslq %r10d, %r14
cmpq %r14, %r11
jl 0xa785
incl %r15d
movslq %r13d, %r8
decl %r13d
movq -0x8(%rax,%r8,8), %r8
movq %r8, (%rax,%rdi,8)
jmp 0xa7c3
incl %ecx
cmpl %r13d, %ecx
jl 0xa72a
jmp 0xa7dc
movl $0x1, %r15d
subl 0x10(%r14), %r15d
imull %eax, %r15d
movl 0x27b3a(%rip), %ebp     # 0x3231c
movl 0x27b37(%rip), %r14d    # 0x32320
movl %ebx, %eax
movq %rax, 0x18(%rsp)
cmpl %ebp, %ebx
movl %r15d, 0x8(%rsp)
jle 0xa857
movl %ebx, 0x27b1d(%rip)     # 0x3231c
testl %r14d, %r14d
movq 0x18(%rsp), %rbp
jle 0xa830
movq 0x27b18(%rip), %r15     # 0x32328
xorl %r12d, %r12d
movq (%r15,%r12,8), %rdi
callq 0x2030
movq %rbp, %rdi
callq 0x21d0
movq %rax, (%r15,%r12,8)
incq %r12
cmpq %r12, %r14
jne 0xa813
movq 0x27b01(%rip), %rdi     # 0x32338
callq 0x2030
leaq 0x4(,%rbp,4), %rdi
callq 0x21d0
movq %rax, 0x27ae8(%rip)     # 0x32338
movl %ebx, %ebp
movl 0x8(%rsp), %r15d
movl %r13d, %eax
movq %rax, 0x20(%rsp)
cmpl %r14d, %r13d
movl %r13d, 0xc(%rsp)
jle 0xa8ba
movslq %r14d, %r14
movq 0x27ab5(%rip), %rdi     # 0x32328
movq 0x20(%rsp), %rax
leaq (,%rax,8), %rsi
callq 0x2200
movq %rax, %r15
movq %rax, 0x27a99(%rip)     # 0x32328
movl %ebp, %r12d
movslq %r13d, %r13
movq %r12, %rdi
callq 0x21d0
movq %rax, (%r15,%r14,8)
incq %r14
cmpq %r14, %r13
jne 0xa895
movl 0xc(%rsp), %r13d
movl %r13d, 0x27a6b(%rip)    # 0x32320
movl 0x8(%rsp), %r15d
testl %r13d, %r13d
jle 0xaa1d
leal -0x1(%rbx), %r12d
movq 0x27a5a(%rip), %rax     # 0x32328
movq %rax, 0x38(%rsp)
movq 0x27a2e(%rip), %r10     # 0x32308
movq 0x27a2f(%rip), %r14     # 0x32310
movq 0x27a50(%rip), %rax     # 0x32338
movq %rax, 0x40(%rsp)
movq 0x27a3c(%rip), %rax     # 0x32330
movq %rax, 0x10(%rsp)
xorl %ebp, %ebp
xorl %r15d, %r15d
xorl %r13d, %r13d
movq %r10, 0x30(%rsp)
movq 0x38(%rsp), %rax
movq (%rax,%rbp,8), %r11
testl %ebx, %ebx
jle 0xa933
movq %r11, 0x28(%rsp)
movq 0x28(%rsp), %rdi
xorl %esi, %esi
movq 0x18(%rsp), %rdx
callq 0x2120
movq 0x28(%rsp), %r11
movq 0x30(%rsp), %r10
movq (%r10,%rbp,8), %rdx
movl (%rdx), %esi
testl %esi, %esi
jle 0xa9df
movl $0x1, %ecx
xorl %eax, %eax
movq 0x8(%rdx), %rdi
movslq 0x4(%rdi,%rax,8), %r8
movl (%r14,%r8,4), %r9d
testl %r9d, %r9d
je 0xa9d0
decl %r9d
movl %r9d, (%r14,%r8,4)
movslq 0x4(%rdi,%rax,8), %rdx
movslq %r13d, %rsi
movq 0x40(%rsp), %r8
movl %edx, (%r8,%rsi,4)
movq 0x10(%rsp), %rsi
movl (%rsi,%rdx,4), %esi
testl %esi, %esi
js 0xa989
movl %r12d, %r8d
movl %r15d, %r9d
movl %esi, %r12d
jmp 0xa9ba
incl %r13d
testl %ecx, %ecx
jne 0xa9a5
movq 0x10(%rsp), %rcx
movl %r12d, (%rcx,%rdx,4)
leal -0x1(%r12), %r8d
xorl %ecx, %ecx
movl %r15d, %r9d
jmp 0xa9ba
movq 0x10(%rsp), %rcx
movl %r15d, (%rcx,%rdx,4)
leal 0x1(%r15), %r9d
xorl %ecx, %ecx
movl %r12d, %r8d
movl %r15d, %r12d
movb (%rdi,%rax,8), %dl
movslq %r12d, %rsi
movb %dl, (%r11,%rsi)
movq (%r10,%rbp,8), %rdx
movl (%rdx), %esi
movl %r8d, %r12d
movl %r9d, %r15d
incq %rax
movslq %esi, %rdi
cmpq %rdi, %rax
jl 0xa948
incq %rbp
cmpq 0x20(%rsp), %rbp
jne 0xa906
testl %r13d, %r13d
movl 0x8(%rsp), %r15d
jle 0xaa1d
movq 0x27932(%rip), %rax     # 0x32330
movq 0x27933(%rip), %rcx     # 0x32338
movl %r13d, %edx
xorl %esi, %esi
movslq (%rcx,%rsi,4), %rdi
movl $0xffffffff, (%rax,%rdi,4) # imm = 0xFFFFFFFF
incq %rsi
cmpq %rsi, %rdx
jne 0xaa0a
movq 0x27904(%rip), %rdi     # 0x32328
movl 0xc(%rsp), %esi
movl %ebx, %edx
callq 0x3374
addl %r15d, %eax
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq