// File: /workspace/github_c_makefile_star2_O2/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 $0x48, %rsp
movl %ecx, 0x14(%rsp)
movq %rdx, %r13
movl %esi, %r12d
movq %rdi, %rbp
cmpb $0x0, 0x28(%rdx)
jne 0x9b8e
movq %r13, %rdi
callq 0x9638
movl 0x10(%r13), %eax
imull %r12d, %eax
cmpl 0x267a4(%rip), %eax     # 0x30340
jle 0x9bcb
movq 0x267a3(%rip), %rdi     # 0x30348
callq 0x2030
movslq 0x10(%r13), %rax
movslq %r12d, %rdi
imulq %rax, %rdi
movl %edi, 0x26785(%rip)     # 0x30340
shlq $0x3, %rdi
callq 0x21c0
movq %rax, 0x2677d(%rip)     # 0x30348
movslq 0x38(%r13), %rsi
movslq 0x26782(%rip), %rbx     # 0x30358
cmpl %ebx, %esi
jle 0x9c28
movq 0x2678f(%rip), %rdi     # 0x30370
shlq $0x2, %rsi
callq 0x21f0
movq %rax, 0x2677f(%rip)     # 0x30370
movslq 0x38(%r13), %rcx
cmpq %rcx, %rbx
jge 0x9c03
orl $-0x1, (%rax,%rbx,4)
incq %rbx
jmp 0x9bf5
movq 0x26746(%rip), %rdi     # 0x30350
callq 0x2030
movslq 0x38(%r13), %rdi
movl %edi, 0x2673f(%rip)     # 0x30358
pushq $0x4
popq %rsi
callq 0x2150
movq %rax, 0x26728(%rip)     # 0x30350
movq 0x26719(%rip), %rax     # 0x30348
movq %rax, 0x8(%rsp)
xorl %eax, %eax
testl %r12d, %r12d
cmovlel %eax, %r12d
xorl %r15d, %r15d
xorl %ebx, %ebx
cmpq %r12, %rax
je 0x9c97
movq 0x18(%r13), %rdx
movslq (%rbp,%rax,4), %rcx
movl 0x10(%r13), %esi
imull (%rdx,%rcx,4), %esi
cmpl 0x8(%r13), %esi
jge 0x9c92
movslq %r15d, %rdx
movq 0x8(%rsp), %rsi
leaq (%rsi,%rdx,8), %rsi
xorl %edx, %edx
movq 0x30(%r13), %rdi
movq (%rdi,%rcx,8), %rdi
movslq (%rdi), %r8
cmpq %r8, %rdx
jge 0x9c8d
movq 0x8(%rdi), %rdi
movq (%rdi,%rdx,8), %rdi
movq %rdi, (%rsi,%rdx,8)
incq %rdx
jmp 0x9c6c
addl %edx, %r15d
incl %ebx
incq %rax
jmp 0x9c42
movq 0x8(%rsp), %rdi
movq %r15, %r14
movl %r14d, %esi
callq 0x964e
pushq $0x1
popq %rcx
movl %ecx, %edx
subl 0x10(%r13), %edx
movl %ebx, 0x18(%rsp)
imull %ebx, %edx
movq 0x26692(%rip), %rbx     # 0x30350
xorl %esi, %esi
testl %r14d, %r14d
movl $0x0, %edi
cmovgl %r14d, %edi
movl %ecx, %r8d
xorl %r12d, %r12d
cmpq %rdi, %rsi
je 0x9d36
movq 0x8(%rsp), %rax
movq (%rax,%rsi,8), %rax
xorl %r9d, %r9d
xorl %r10d, %r10d
movslq (%rax), %r11
cmpq %r11, %r9
jge 0x9d1a
movq 0x8(%rax), %r11
movslq 0x4(%r11,%r9,8), %r11
movl (%rbx,%r11,4), %ebp
xorl %r14d, %r14d
testl %ebp, %ebp
sete %r14b
movzbl %r10b, %r10d
cmovel %ecx, %r10d
addl %r14d, %r12d
incl %ebp
movl %ebp, (%rbx,%r11,4)
incq %r9
jmp 0x9ce6
testb %r10b, %r10b
je 0x9d2e
leal 0x1(%rdx), %eax
cmpl 0x14(%rsp), %edx
jge 0x9e9c
movl %eax, %edx
incq %rsi
incl %r8d
jmp 0x9cd2
pushq $0x1
popq %rax
subl 0x10(%r13), %eax
imull 0x18(%rsp), %eax
movq %rax, 0x38(%rsp)
xorl %ecx, %ecx
movslq %r15d, %r13
cmpl %r15d, %ecx
jge 0x9e28
movslq %ecx, %rsi
movq 0x8(%rsp), %rax
movq (%rax,%rsi,8), %rax
movl (%rax), %edi
xorl %r9d, %r9d
testl %edi, %edi
movl $0x0, %r10d
cmovgl %edi, %r10d
xorl %r11d, %r11d
cmpq %r9, %r10
je 0x9d9b
movq 0x8(%rax), %r8
movslq 0x4(%r8,%r9,8), %r14
movl (%rbx,%r14,4), %ebp
testl %ebp, %ebp
je 0x9d96
cmpl $0x1, %ebp
je 0x9da5
incl %r11d
movl %r9d, %edx
incq %r9
jmp 0x9d75
cmpl $0x1, %r11d
jle 0x9dde
incl %ecx
jmp 0x9d49
movq 0x8(%rsp), %r9
leaq (%r9,%rsi,8), %rsi
xorl %r9d, %r9d
movslq %edi, %r10
cmpq %r10, %r9
jge 0x9dfc
movslq 0x4(%r8,%r9,8), %r10
movl (%rbx,%r10,4), %r11d
testl %r11d, %r11d
je 0x9dd9
cmpl $0x1, %r11d
jne 0x9dd0
decl %r12d
decl %r11d
movl %r11d, (%rbx,%r10,4)
movl (%rax), %edi
incq %r9
jmp 0x9db1
movq 0x8(%rsp), %rdi
leaq (%rdi,%rsi,8), %rsi
jne 0x9e13
movq 0x8(%rax), %rax
movslq %edx, %rdi
movslq 0x4(%rax,%rdi,8), %rax
andl $0x0, (%rbx,%rax,4)
decl %r12d
movq 0x38(%rsp), %rdi
leal 0x1(%rdi), %eax
cmpl %edi, 0x14(%rsp)
jle 0xa0af
movq %rax, 0x38(%rsp)
decl %r15d
movq 0x8(%rsp), %rax
movq -0x8(%rax,%r13,8), %rax
movq %rax, (%rsi)
jmp 0x9d49
movq %r15, 0x30(%rsp)
movl 0x26529(%rip), %ebp     # 0x3035c
movslq 0x26526(%rip), %r15     # 0x30360
movq 0x26527(%rip), %rax     # 0x30368
movq %rax, 0x20(%rsp)
cmpl %ebp, %r12d
jle 0x9efd
movl %r12d, 0x26506(%rip)    # 0x3035c
movl %r12d, %eax
movq %rax, 0x28(%rsp)
xorl %r14d, %r14d
testl %r15d, %r15d
movl $0x0, %eax
cmovgl %r15d, %eax
movq %rax, 0x18(%rsp)
movq 0x20(%rsp), %rbp
cmpq %r14, 0x18(%rsp)
je 0x9ed5
movq (%rbp,%r14,8), %rdi
callq 0x2030
movq 0x28(%rsp), %rdi
callq 0x21c0
movq %rax, (%rbp,%r14,8)
incq %r14
jmp 0x9e77
movl %r8d, %ecx
xorl %edx, %edx
cmpq %rcx, %rdx
je 0xa11f
movq 0x8(%rsp), %rsi
movq (%rsi,%rdx,8), %rsi
xorl %edi, %edi
movslq (%rsi), %r8
cmpq %r8, %rdi
jge 0x9ed0
movq 0x8(%rsi), %r8
movslq 0x4(%r8,%rdi,8), %r8
andl $0x0, (%rbx,%r8,4)
incq %rdi
jmp 0x9eb5
incq %rdx
jmp 0x9ea1
movq 0x2649c(%rip), %rdi     # 0x30378
callq 0x2030
movq 0x28(%rsp), %rax
leaq 0x4(,%rax,4), %rdi
callq 0x21c0
movq %rax, 0x2647e(%rip)     # 0x30378
movl %r12d, %ebp
movq 0x30(%rsp), %rax
cmpl %r15d, %eax
jle 0x9f4b
movl %eax, %esi
shlq $0x3, %rsi
movq 0x20(%rsp), %rdi
callq 0x21f0
movq %rax, 0x20(%rsp)
movq %rax, 0x26445(%rip)     # 0x30368
movl %ebp, %ebp
cmpq %r13, %r15
jge 0x9f40
movq %rbp, %rdi
callq 0x21c0
movq 0x20(%rsp), %rcx
movq %rax, (%rcx,%r15,8)
incq %r15
jmp 0x9f25
movq 0x30(%rsp), %rax
movl %eax, 0x26415(%rip)     # 0x30360
leal -0x1(%r12), %r9d
movl 0x14(%rsp), %ecx
subl 0x38(%rsp), %ecx
movl %ecx, 0x28(%rsp)
movq 0x26415(%rip), %rcx     # 0x30378
movq %rcx, 0x18(%rsp)
movq 0x26401(%rip), %rdx     # 0x30370
xorl %esi, %esi
testl %r12d, %r12d
movl $0x0, %r8d
cmovgl %r12d, %r8d
testl %eax, %eax
movl $0x0, %ecx
cmovgl %eax, %ecx
movq %rcx, 0x40(%rsp)
xorl %edi, %edi
xorl %r15d, %r15d
cmpq 0x40(%rsp), %rsi
je 0xa0e7
movq 0x20(%rsp), %rax
movq (%rax,%rsi,8), %r11
xorl %eax, %eax
cmpq %rax, %r8
je 0x9fb7
movb $0x0, (%r11,%rax)
incq %rax
jmp 0x9fa8
xorl %r13d, %r13d
pushq $0x1
popq %rbp
movq 0x8(%rsp), %rax
movq (%rax,%rsi,8), %rax
movslq (%rax), %rcx
cmpq %rcx, %r13
jge 0xa04b
movq 0x8(%rax), %r14
movslq 0x4(%r14,%r13,8), %rax
movl (%rbx,%rax,4), %ecx
testl %ecx, %ecx
je 0xa043
decl %ecx
movl %ecx, (%rbx,%rax,4)
movslq 0x4(%r14,%r13,8), %rax
movslq %edi, %r10
movq 0x18(%rsp), %rcx
movl %eax, (%rcx,%r10,4)
movl (%rdx,%rax,4), %ecx
testl %ecx, %ecx
js 0xa006
movl %r9d, %r10d
movl %r15d, %eax
movl %ecx, %r9d
jmp 0xa032
incl %edi
testl %ebp, %ebp
jne 0xa01b
movl %r9d, (%rdx,%rax,4)
leal -0x1(%r9), %r10d
xorl %ebp, %ebp
movl %r15d, %eax
jmp 0xa032
movl %r15d, (%rdx,%rax,4)
cmpl 0x28(%rsp), %r15d
je 0xa053
leal 0x1(%r15), %eax
xorl %ebp, %ebp
movl %r9d, %r10d
movl %r15d, %r9d
movb (%r14,%r13,8), %cl
movslq %r9d, %r9
movb %cl, (%r11,%r9)
movl %r10d, %r9d
movl %eax, %r15d
incq %r13
jmp 0x9fbd
incq %rsi
jmp 0x9f92
movq 0x30(%rsp), %r8
cmpl %esi, %r8d
cmovlel %esi, %r8d
cmpq %r8, %rsi
je 0xa08e
movq 0x8(%rsp), %rax
movq (%rax,%rsi,8), %rax
xorl %edi, %edi
movslq (%rax), %rcx
cmpq %rcx, %rdi
jge 0xa089
movq 0x8(%rax), %rcx
movslq 0x4(%rcx,%rdi,8), %rcx
andl $0x0, (%rbx,%rcx,4)
incq %rdi
jmp 0xa06f
incq %rsi
jmp 0xa05f
xorl %eax, %eax
cmpq %r10, %rax
jg 0xa0a7
movq 0x18(%rsp), %rcx
movslq (%rcx,%rax,4), %rcx
orl $-0x1, (%rdx,%rcx,4)
incq %rax
jmp 0xa090
movl 0x14(%rsp), %eax
incl %eax
jmp 0xa11f
xorl %ecx, %ecx
testl %r15d, %r15d
cmovlel %ecx, %r15d
cmpq %r15, %rcx
je 0xa11f
movq 0x8(%rsp), %rdx
movq (%rdx,%rcx,8), %rdx
xorl %esi, %esi
movslq (%rdx), %rdi
cmpq %rdi, %rsi
jge 0xa0e2
movq 0x8(%rdx), %rdi
movslq 0x4(%rdi,%rsi,8), %rdi
andl $0x0, (%rbx,%rdi,4)
incq %rsi
jmp 0xa0c8
incq %rcx
jmp 0xa0b8
xorl %eax, %eax
testl %edi, %edi
cmovlel %eax, %edi
movq 0x30(%rsp), %rsi
cmpq %rax, %rdi
je 0xa10a
movq 0x18(%rsp), %rcx
movslq (%rcx,%rax,4), %rcx
orl $-0x1, (%rdx,%rcx,4)
incq %rax
jmp 0xa0f3
movq 0x20(%rsp), %rdi
movl %r12d, %edx
movl 0x28(%rsp), %ecx
callq 0x2c07
addl 0x38(%rsp), %eax
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq