// File: /workspace/github_c_makefile_star2_O3/cbouilla[P]AES-attacks-finder/SysEqLin.c
// Function: TailleE_Borne
int TailleE_Borne(int * __restrict const X, const int N, SysEqLin E, const int borne) // 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++;
		}
	}
	TriLignesCreusesByNNZ(L,NombreDeLignes);
	
	/* 
	 * On compte le nombre d'éléments non-nuls sur chaque colonne
	 */
	int dim = (1-E->gen)*taille;
	int NombreDeColonnes = 0;
	for (i = 0; i < NombreDeLignes; i++)
	{
		unsigned char pivot_present = 0;
		int j;
		for (j = 0; j < L[i]->nnz; j++)
		{
			if (!nombre_nnz_colonne[L[i]->elements[j].colonne])
			{
				pivot_present = 1;
				NombreDeColonnes++;
			}
			nombre_nnz_colonne[L[i]->elements[j].colonne]++;
		}
		if (pivot_present)
		{
			dim++;
			if (dim>borne) // on est sur que la sortie sera plus grande que borne
			{
				for (j = 0; j <= i ; j++) // on remet le tableau à zero
				{
					int k;
					for (k = 0; k < L[j]->nnz; k++) nombre_nnz_colonne[L[j]->elements[k].colonne]=0;
				}
				//compteur1++;
				//printf("cpt1 = %d, cpt2 = %d, cpt3 = %d, cpt4 = %d\n",compteur1,compteur2,compteur3,compteur4);
				return dim;
			}
		}
	}
	
	dim = (1-E->gen)*taille;
	int ligne = 0;
	while (ligne < NombreDeLignes)
	{
		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++;
				if (dim>borne) goto remise_a_zero;
			}
			NombreDeLignes--;
			L[ligne]=L[NombreDeLignes];
		}
		else ligne++;
	}

	{

		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;
		const int new_borne = borne - dim + 1;
		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]--;
					int colonne;
					colonne_changer[pos]=L[i]->elements[k].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++;
							if (Ncol != new_borne) premier=0;
							else // on est sur que la sortie sera plus grande que borne
							{
								for (i = i; i < NombreDeLignes; i++)
								{
									for (k = 0; k < L[i]->nnz; k++)	nombre_nnz_colonne[L[i]->elements[k].colonne]=0;								
								}
								for (i = 0; i < pos; i++) nouvelle_colonne[colonne_changer[i]]=-1;
								return borne+1;
							}
						}
					}
					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+DimensionBorne_tmp(m,NombreDeLignes,NombreDeColonnes,borne-dim);

	}

	remise_a_zero :
		
	for (ligne = 0; ligne < NombreDeLignes; ligne++)
	{
		int k;
		for (k = 0; k < L[ligne]->nnz; k++) nombre_nnz_colonne[L[ligne]->elements[k].colonne]=0;
	}
	return dim;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x58, %rsp
movl %ecx, 0x14(%rsp)
movq %rdx, %r12
movl %esi, %ebx
movq %rdi, %r13
cmpb $0x0, 0x28(%rdx)
jne 0xaa6f
movq %r12, %rdi
callq 0x9eab
movb $0x1, 0x28(%r12)
movl 0x10(%r12), %eax
imull %ebx, %eax
cmpl 0x278c3(%rip), %eax     # 0x32340
jle 0xaaad
movq 0x278c2(%rip), %rdi     # 0x32348
callq 0x2030
movslq 0x10(%r12), %rax
movslq %ebx, %rdi
imulq %rax, %rdi
movl %edi, 0x278a3(%rip)     # 0x32340
shlq $0x3, %rdi
callq 0x21d0
movq %rax, 0x2789b(%rip)     # 0x32348
movslq 0x38(%r12), %rsi
movslq 0x2789f(%rip), %r14     # 0x32358
cmpl %r14d, %esi
jle 0xab23
movq 0x278ab(%rip), %rdi     # 0x32370
shlq $0x2, %rsi
callq 0x2200
movq %rax, 0x2789b(%rip)     # 0x32370
movl 0x38(%r12), %ecx
cmpl %ecx, %r14d
jge 0xaafb
leaq (%rax,%r14,4), %rdi
notl %r14d
addl %r14d, %ecx
leaq 0x4(,%rcx,4), %rdx
movl $0xff, %esi
callq 0x2120
movq 0x2784e(%rip), %rdi     # 0x32350
callq 0x2030
movslq 0x38(%r12), %rdi
movl %edi, 0x27846(%rip)     # 0x32358
movl $0x4, %esi
callq 0x2160
movq %rax, 0x2782d(%rip)     # 0x32350
xorl %r15d, %r15d
movq 0x2781b(%rip), %rax     # 0x32348
movq %rax, 0x8(%rsp)
movl $0x0, %r14d
movl $0x0, %ebp
testl %ebx, %ebx
jle 0xabb6
movl 0x8(%r12), %eax
movl 0x10(%r12), %ecx
movl %ebx, %edx
xorl %esi, %esi
xorl %ebp, %ebp
xorl %r14d, %r14d
movq 0x18(%r12), %r8
movslq (%r13,%rsi,4), %rdi
movl (%r8,%rdi,4), %r8d
imull %ecx, %r8d
cmpl %eax, %r8d
jge 0xabae
movq 0x30(%r12), %r8
movq (%r8,%rdi,8), %r10
cmpl $0x0, (%r10)
jle 0xabac
movslq %r14d, %r8
movq 0x8(%rsp), %r9
leaq (%r9,%r8,8), %r9
xorl %r8d, %r8d
movq 0x8(%r10), %r10
movq (%r10,%r8,8), %r10
movq %r10, (%r9,%r8,8)
incq %r8
movq 0x30(%r12), %r10
movq (%r10,%rdi,8), %r10
movslq (%r10), %r11
cmpq %r11, %r8
jl 0xab89
addl %r8d, %r14d
incl %ebp
incq %rsi
cmpq %rdx, %rsi
jne 0xab54
movl %ebp, 0x28(%rsp)
movq 0x8(%rsp), %rdi
movl %r14d, %esi
callq 0xa4b3
movl 0x10(%r12), %eax
movq 0x2777d(%rip), %rbp     # 0x32350
testl %r14d, %r14d
jle 0xac7b
movl $0x1, %ecx
movl $0x1, %esi
subl %eax, %esi
imull 0x28(%rsp), %esi
movl %r14d, %eax
movq %rax, 0x30(%rsp)
xorl %edi, %edi
movl $0x1, %r8d
xorl %r15d, %r15d
movq 0x8(%rsp), %rax
movq (%rax,%rdi,8), %rax
cmpl $0x0, (%rax)
jle 0xac69
movl %r14d, %edx
movq 0x8(%rax), %r9
xorl %r10d, %r10d
xorl %r11d, %r11d
movslq 0x4(%r9,%r10,8), %rbx
movl (%rbp,%rbx,4), %r13d
xorl %r14d, %r14d
testl %r13d, %r13d
sete %r14b
movzbl %r11b, %r11d
cmovel %ecx, %r11d
addl %r14d, %r15d
incl %r13d
movl %r13d, (%rbp,%rbx,4)
incq %r10
movslq (%rax), %rbx
cmpq %rbx, %r10
jl 0xac1b
testb %r11b, %r11b
je 0xac66
leal 0x1(%rsi), %eax
cmpl 0x14(%rsp), %esi
jge 0xb044
movl %edx, %r14d
movl %eax, %esi
jmp 0xac69
movl %edx, %r14d
incq %rdi
incl %r8d
cmpq 0x30(%rsp), %rdi
jne 0xac00
movl 0x10(%r12), %eax
movl $0x1, %ebx
subl %eax, %ebx
imull 0x28(%rsp), %ebx
testl %r14d, %r14d
jle 0xad59
xorl %ecx, %ecx
movslq %ecx, %rsi
movq 0x8(%rsp), %rax
movq (%rax,%rsi,8), %rax
movl (%rax), %r8d
testl %r8d, %r8d
jle 0xad38
movq 0x8(%rax), %rdi
xorl %r9d, %r9d
xorl %r10d, %r10d
movslq 0x4(%rdi,%r9,8), %r11
movl (%rbp,%r11,4), %r11d
testl %r11d, %r11d
je 0xaccf
cmpl $0x1, %r11d
je 0xacf8
incl %r10d
movl %r9d, %edx
incq %r9
cmpq %r9, %r8
jne 0xacb4
cmpl $0x1, %r10d
jg 0xad4e
jne 0xad38
movq 0x8(%rax), %rax
movslq %edx, %rdi
movslq 0x4(%rax,%rdi,8), %rax
movl $0x0, (%rbp,%rax,4)
decl %r15d
jmp 0xad29
xorl %r9d, %r9d
movslq 0x4(%rdi,%r9,8), %r10
movl (%rbp,%r10,4), %r11d
testl %r11d, %r11d
je 0xad1e
cmpl $0x1, %r11d
jne 0xad13
decl %r15d
decl %r11d
movl %r11d, (%rbp,%r10,4)
movl (%rax), %r8d
incq %r9
movslq %r8d, %r10
cmpq %r10, %r9
jl 0xacfb
leal 0x1(%rbx), %eax
cmpl 0x14(%rsp), %ebx
jge 0xb081
movl %eax, %ebx
movslq %r14d, %rax
decl %r14d
movq 0x8(%rsp), %rdi
movq -0x8(%rdi,%rax,8), %rax
movq %rax, (%rdi,%rsi,8)
jmp 0xad50
incl %ecx
cmpl %r14d, %ecx
jl 0xac92
movq %rbx, 0x38(%rsp)
movl %r14d, 0x10(%rsp)
movl 0x275f2(%rip), %r13d    # 0x3235c
movl 0x275f0(%rip), %ebx     # 0x32360
movl %r15d, %eax
movq %rax, 0x40(%rsp)
cmpl %r13d, %r15d
jle 0xadd7
movl %r15d, 0x275d8(%rip)    # 0x3235c
testl %ebx, %ebx
movq 0x40(%rsp), %r13
jle 0xadb4
movq 0x275d4(%rip), %r14     # 0x32368
xorl %r12d, %r12d
movq (%r14,%r12,8), %rdi
callq 0x2030
movq %r13, %rdi
callq 0x21d0
movq %rax, (%r14,%r12,8)
incq %r12
cmpq %r12, %rbx
jne 0xad97
movq 0x275bd(%rip), %rdi     # 0x32378
callq 0x2030
leaq 0x4(,%r13,4), %rdi
callq 0x21d0
movq %rax, 0x275a4(%rip)     # 0x32378
movl %r15d, %r13d
movl 0x10(%rsp), %r14d
movl %r14d, %ecx
cmpl %ebx, %r14d
movq %rcx, 0x48(%rsp)
jle 0xae30
movslq %ebx, %r12
movq 0x27575(%rip), %rdi     # 0x32368
leaq (,%rcx,8), %rsi
callq 0x2200
movq %rax, %rbx
movq %rax, 0x2755e(%rip)     # 0x32368
movl %r13d, %r13d
movslq %r14d, %r14
movq %r13, %rdi
callq 0x21d0
movq %rax, (%rbx,%r12,8)
incq %r12
cmpq %r12, %r14
jne 0xae10
movl 0x10(%rsp), %r14d
movl %r14d, 0x27530(%rip)    # 0x32360
movl 0x14(%rsp), %eax
movq 0x38(%rsp), %rbx
subl %ebx, %eax
movl %eax, 0x24(%rsp)
testl %r14d, %r14d
jle 0xafc2
leal -0x1(%r15), %r13d
movq 0x27515(%rip), %rax     # 0x32368
movq %rax, 0x50(%rsp)
movq 0x27519(%rip), %rax     # 0x32378
movq %rax, 0x28(%rsp)
movq 0x27505(%rip), %r10     # 0x32370
xorl %r14d, %r14d
xorl %ebx, %ebx
xorl %r12d, %r12d
movq %r10, 0x30(%rsp)
movq 0x50(%rsp), %rax
movq (%rax,%r14,8), %r11
testl %r15d, %r15d
jle 0xaea4
movq %r11, %rdi
xorl %esi, %esi
movq 0x40(%rsp), %rdx
movq %r11, 0x18(%rsp)
callq 0x2120
movq 0x18(%rsp), %r11
movq 0x30(%rsp), %r10
movq 0x8(%rsp), %rax
movq (%rax,%r14,8), %rsi
movl (%rsi), %ecx
testl %ecx, %ecx
jle 0xaf7b
movl $0x1, 0x18(%rsp)
xorl %eax, %eax
movq 0x8(%rsi), %r8
movslq 0x4(%r8,%rax,8), %rdx
movl (%rbp,%rdx,4), %r9d
testl %r9d, %r9d
je 0xaf6c
decl %r9d
movl %r9d, (%rbp,%rdx,4)
movslq 0x4(%r8,%rax,8), %rsi
movslq %r12d, %rdx
movq 0x28(%rsp), %r9
movl %esi, (%r9,%rdx,4)
movl (%r10,%rsi,4), %edi
testl %edi, %edi
js 0xaf07
movl %r12d, %edx
movl %r13d, %r9d
movl %ebx, %r10d
movl %edi, %r13d
jmp 0xaf48
leal 0x1(%r12), %edx
cmpl $0x0, 0x18(%rsp)
jne 0xaf28
movl %r13d, (%r10,%rsi,4)
leal -0x1(%r13), %r9d
movl $0x0, 0x18(%rsp)
movl %ebx, %r10d
jmp 0xaf48
movl %ebx, (%r10,%rsi,4)
cmpl 0x24(%rsp), %ebx
je 0xafdf
leal 0x1(%rbx), %r10d
movl $0x0, 0x18(%rsp)
movl %r13d, %r9d
movl %ebx, %r13d
movb (%r8,%rax,8), %sil
movslq %r13d, %rdi
movb %sil, (%r11,%rdi)
movq 0x8(%rsp), %rsi
movq (%rsi,%r14,8), %rsi
movl (%rsi), %ecx
movl %edx, %r12d
movl %r9d, %r13d
movl %r10d, %ebx
movq 0x30(%rsp), %r10
incq %rax
movslq %ecx, %rdx
cmpq %rdx, %rax
jl 0xaec1
incq %r14
movq 0x48(%rsp), %rcx
cmpq %rcx, %r14
jne 0xae78
testl %r12d, %r12d
movl 0x10(%rsp), %r14d
movq 0x38(%rsp), %rbx
jle 0xafc2
movq 0x273ce(%rip), %rax     # 0x32370
movq 0x273cf(%rip), %rcx     # 0x32378
movl %r12d, %edx
xorl %r8d, %r8d
movslq (%rcx,%r8,4), %rdi
movl $0xffffffff, (%rax,%rdi,4) # imm = 0xFFFFFFFF
incq %r8
cmpq %r8, %rdx
jne 0xafaf
movq 0x2739f(%rip), %rdi     # 0x32368
movl %r14d, %esi
movl %r15d, %edx
movl 0x24(%rsp), %ecx
callq 0x321a
addl %ebx, %eax
jmp 0xb0c1
cmpl %r14d, 0x10(%rsp)
movq 0x48(%rsp), %r8
jle 0xb01f
movq 0x8(%rsp), %rax
movq (%rax,%r14,8), %rax
cmpl $0x0, (%rax)
jle 0xb017
movq 0x8(%rax), %rcx
xorl %esi, %esi
movslq 0x4(%rcx,%rsi,8), %rdi
movl $0x0, (%rbp,%rdi,4)
incq %rsi
movslq (%rax), %rdi
cmpq %rdi, %rsi
jl 0xafff
incq %r14
cmpq %r8, %r14
jne 0xafeb
testl %r12d, %r12d
js 0xb03c
movl %edx, %eax
xorl %ecx, %ecx
movslq (%r9,%rcx,4), %rdx
movl $0xffffffff, (%r10,%rdx,4) # imm = 0xFFFFFFFF
incq %rcx
cmpq %rcx, %rax
jne 0xb028
movl 0x14(%rsp), %eax
incl %eax
jmp 0xb0c1
movl %r8d, %ecx
xorl %edx, %edx
movq 0x8(%rsp), %rsi
movq (%rsi,%rdx,8), %rsi
cmpl $0x0, (%rsi)
jle 0xb077
movq 0x8(%rsi), %rdi
xorl %r8d, %r8d
movslq 0x4(%rdi,%r8,8), %r9
movl $0x0, (%rbp,%r9,4)
incq %r8
movslq (%rsi), %r9
cmpq %r9, %r8
jl 0xb05e
incq %rdx
cmpq %rcx, %rdx
jne 0xb049
jmp 0xb0c1
testl %r14d, %r14d
jle 0xb0c1
movl %r14d, %ecx
xorl %edx, %edx
movq 0x8(%rsp), %rsi
movq (%rsi,%rdx,8), %rsi
cmpl $0x0, (%rsi)
jle 0xb0b9
movq 0x8(%rsi), %rdi
xorl %r8d, %r8d
movslq 0x4(%rdi,%r8,8), %r9
movl $0x0, (%rbp,%r9,4)
incq %r8
movslq (%rsi), %r9
cmpq %r9, %r8
jl 0xb0a0
incq %rdx
cmpq %rcx, %rdx
jne 0xb08b
addq $0x58, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq