// File: /workspace/github_c_makefile_star2_O1/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 0xa2aa
movq %r14, %rdi
callq 0x9bd9
movb $0x1, 0x28(%r14)
movl 0x10(%r14), %eax
imull %r12d, %eax
cmpl 0x2804c(%rip), %eax     # 0x32304
jle 0xa2e7
movq 0x28047(%rip), %rdi     # 0x32308
callq 0x2030
movslq 0x10(%r14), %rax
movslq %r12d, %rdi
imulq %rax, %rdi
movl %edi, 0x2802d(%rip)     # 0x32304
shlq $0x3, %rdi
callq 0x21d0
movq %rax, 0x28021(%rip)     # 0x32308
movslq 0x38(%r14), %rsi
movslq 0x28026(%rip), %rbx     # 0x32318
cmpl %ebx, %esi
jle 0xa356
movq 0x28033(%rip), %rdi     # 0x32330
shlq $0x2, %rsi
callq 0x2200
movq %rax, 0x28023(%rip)     # 0x32330
movl 0x38(%r14), %ecx
cmpl %ecx, %ebx
jge 0xa32f
leaq (%rax,%rbx,4), %rdi
notl %ebx
addl %ebx, %ecx
leaq 0x4(,%rcx,4), %rdx
movl $0xff, %esi
callq 0x2120
movq 0x27fda(%rip), %rdi     # 0x32310
callq 0x2030
movslq 0x38(%r14), %rdi
movl %edi, 0x27fd3(%rip)     # 0x32318
movl $0x4, %esi
callq 0x2160
movq %rax, 0x27fba(%rip)     # 0x32310
xorl %ebx, %ebx
movl $0x0, %r13d
movl $0x0, %eax
testl %r12d, %r12d
jle 0xa3e4
movl 0x8(%r14), %ecx
movl 0x10(%r14), %edx
movq 0x27f91(%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 0xa3dc
movq 0x30(%r14), %r10
movq (%r10,%r9,8), %r12
cmpl $0x0, (%r12)
jle 0xa3da
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 0xa3b3
movl %ebp, %r13d
addl %r10d, %r13d
incl %eax
incq %r8
cmpq %rdi, %r8
jne 0xa382
testl %r13d, %r13d
jle 0xa437
movq 0x27f18(%rip), %rcx     # 0x32308
movq 0x27f19(%rip), %rdx     # 0x32310
movl %r13d, %esi
xorl %edi, %edi
xorl %ebx, %ebx
movq (%rcx,%rdi,8), %r8
cmpl $0x0, (%r8)
jle 0xa42f
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 0xa40f
incq %rdi
cmpq %rsi, %rdi
jne 0xa3fe
movl $0x1, %r12d
subl 0x10(%r14), %r12d
imull %eax, %r12d
testl %r13d, %r13d
jle 0xa50e
movq 0x27eb3(%rip), %rax     # 0x32308
xorl %ecx, %ecx
movq 0x27eb2(%rip), %rdx     # 0x32310
movslq %ecx, %r8
movq (%rax,%r8,8), %rdi
movslq (%rdi), %r10
testq %r10, %r10
jle 0xa499
movq 0x8(%rdi), %r9
xorl %r14d, %r14d
xorl %r11d, %r11d
movslq 0x4(%r9,%r14,8), %r15
movl (%rdx,%r15,4), %ebp
testl %ebp, %ebp
je 0xa48f
cmpl $0x1, %ebp
je 0xa4c3
incl %r11d
movl %r14d, %esi
incq %r14
cmpq %r14, %r10
jne 0xa477
jmp 0xa49c
xorl %r11d, %r11d
cmpl $0x1, %r11d
jg 0xa4bf
leaq (%rax,%r8,8), %r8
jne 0xa4f7
movq 0x8(%rdi), %rdi
movslq %esi, %r9
movslq 0x4(%rdi,%r9,8), %rdi
movl $0x0, (%rdx,%rdi,4)
decl %ebx
jmp 0xa4f4
incl %ecx
jmp 0xa505
leaq (%rax,%r8,8), %r8
cmpl $0x0, (%rdi)
jle 0xa4f4
xorl %r10d, %r10d
movslq 0x4(%r9,%r10,8), %r11
movl (%rdx,%r11,4), %ebp
testl %ebp, %ebp
je 0xa4e9
cmpl $0x1, %ebp
jne 0xa4e3
decl %ebx
decl %ebp
movl %ebp, (%rdx,%r11,4)
incq %r10
movslq (%rdi), %r11
cmpq %r11, %r10
jl 0xa4cf
incl %r12d
movslq %r13d, %rdi
decl %r13d
movq -0x8(%rax,%rdi,8), %rdi
movq %rdi, (%r8)
cmpl %r13d, %ecx
jl 0xa45e
movl %r12d, 0xc(%rsp)
movl %ebx, %ebp
cmpl 0x27e01(%rip), %ebx     # 0x3231c
jle 0xa576
movl %ebx, 0x27df9(%rip)     # 0x3231c
movslq 0x27df6(%rip), %r14     # 0x32320
testq %r14, %r14
jle 0xa556
movq 0x27df2(%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 0xa539
movq 0x27ddb(%rip), %rdi     # 0x32338
callq 0x2030
leaq 0x4(,%rbp,4), %rdi
callq 0x21d0
movq %rax, 0x27dc2(%rip)     # 0x32338
movslq 0x27da3(%rip), %r14     # 0x32320
movl %r13d, %eax
movq %rax, 0x10(%rsp)
cmpl %r14d, %r13d
movl %r13d, 0x8(%rsp)
jle 0xa5dc
movq 0x27d92(%rip), %rdi     # 0x32328
movq 0x10(%rsp), %rax
leaq (,%rax,8), %rsi
callq 0x2200
movq %rax, %r15
movq %rax, 0x27d76(%rip)     # 0x32328
movl 0x27d63(%rip), %r12d    # 0x3231c
movslq %r13d, %r13
movq %r12, %rdi
callq 0x21d0
movq %rax, (%r15,%r14,8)
incq %r14
cmpq %r14, %r13
jne 0xa5bc
movl 0x8(%rsp), %r13d
movl %r13d, 0x27d44(%rip)    # 0x32320
testl %r13d, %r13d
jle 0xa703
movq %rbp, 0x30(%rsp)
leal -0x1(%rbx), %r15d
movq 0x27d33(%rip), %rax     # 0x32328
movq %rax, 0x38(%rsp)
movq 0x27d07(%rip), %r9      # 0x32308
movq 0x27d08(%rip), %r12     # 0x32310
movq 0x27d29(%rip), %r10     # 0x32338
movq 0x27d1a(%rip), %r11     # 0x32330
xorl %ebp, %ebp
xorl %r14d, %r14d
xorl %r13d, %r13d
movq %r9, 0x28(%rsp)
movq %r10, 0x20(%rsp)
movq %r11, 0x18(%rsp)
movq 0x38(%rsp), %rax
movq (%rax,%rbp,8), %rdi
testl %ebx, %ebx
movq %rdi, 0x40(%rsp)
jle 0xa65a
xorl %esi, %esi
movq 0x30(%rsp), %rdx
callq 0x2120
movq 0x18(%rsp), %r11
movq 0x20(%rsp), %r10
movq 0x28(%rsp), %r9
movq (%r9,%rbp,8), %rdx
cmpl $0x0, (%rdx)
jle 0xa6f3
movl $0x1, %ecx
xorl %eax, %eax
movq 0x8(%rdx), %rdx
movslq 0x4(%rdx,%rax,8), %rsi
movl (%r12,%rsi,4), %edi
testl %edi, %edi
je 0xa6e0
decl %edi
movl %edi, (%r12,%rsi,4)
movslq 0x4(%rdx,%rax,8), %rsi
movslq %r13d, %rdi
movl %esi, (%r10,%rdi,4)
movl (%r11,%rsi,4), %r8d
testl %r8d, %r8d
js 0xa6a5
movl %r15d, %esi
movl %r14d, %edi
movl %r8d, %r15d
jmp 0xa6cb
incl %r13d
testl %ecx, %ecx
jne 0xa6bb
movl %r15d, (%r11,%rsi,4)
leal -0x1(%r15), %esi
xorl %ecx, %ecx
movl %r14d, %edi
jmp 0xa6cb
movl %r14d, (%r11,%rsi,4)
leal 0x1(%r14), %edi
xorl %ecx, %ecx
movl %r15d, %esi
movl %r14d, %r15d
movb (%rdx,%rax,8), %dl
movslq %r15d, %r8
movq 0x40(%rsp), %r14
movb %dl, (%r14,%r8)
movl %esi, %r15d
movl %edi, %r14d
incq %rax
movq (%r9,%rbp,8), %rdx
movslq (%rdx), %rsi
cmpq %rsi, %rax
jl 0xa66e
incq %rbp
cmpq 0x10(%rsp), %rbp
jne 0xa62d
jmp 0xa706
xorl %r13d, %r13d
testl %r13d, %r13d
jle 0xa731
movq 0x27c1e(%rip), %rax     # 0x32330
movq 0x27c1f(%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 0xa71e
movq 0x27bf0(%rip), %rdi     # 0x32328
movl 0x8(%rsp), %esi
movl %ebx, %edx
callq 0x3048
addl 0xc(%rsp), %eax
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq