// File: /workspace/github_c_makefile_star2_O1/cbouilla[P]AES-attacks-finder/SysEqLin.c
// Function: InitialiseSousMat_tmp
void InitialiseSousMat_tmp(SysEqLin E) // initialise E->sousmat
{

	/*	
	 *	 | * * * * | *          |		 | * * * * | * * * * * * |
	 * 	 | * * * * |   *        |		 |   * * * | * * * * * * |
	 *	 | * * * * |     *      | 		 |     * * | * * * * * * | 
	 *	 | * * * * |       *    |  - >   |       * | * * * * * * |
	 *	 | * * * * |         *  |	     |         | *   * * * * | | -> Sousmat
	 *	 | * * * * |           *|        |         |   * * * * * | | ->
	 * 												 ___.___.___	
	 *									variable:     0   1   2		
	 */

	/* Initialisation Big Matrice */
	const int nombre_ligne_big_matrix = E->a;
	const int nombre_colone_big_matrix = E->a+E->b;
	unsigned char ** __restrict big_matrix=malloc(nombre_ligne_big_matrix*sizeof(unsigned char *));
	unsigned char * __restrict espace_big_matrix=malloc(nombre_colone_big_matrix*nombre_ligne_big_matrix*sizeof(unsigned char));

	int i;
	for (i = 0; i < nombre_ligne_big_matrix; i++) big_matrix[i]=espace_big_matrix + i*nombre_colone_big_matrix;
	for (i = 0; i < nombre_ligne_big_matrix; i++)
	{
		int j;
		for (j = 0; j < E->b; j++) big_matrix[i][j]=E->mat[i][j];
		for (j = E->b; j < nombre_colone_big_matrix; j++) big_matrix[i][j]=0;
	}
	const int nombre_de_variables = E->a/E->gen;
	for (i = 0; i < nombre_de_variables; i++)
	{
		int j;
		for (j = 0; j < E->gen; j++) big_matrix[E->gen*i+j][E->b+E->gen*i+j]=1;
	}
	
	//printf(" Big matrice \n");
	//PrintMat(big_matrix,nombre_ligne_big_matrix,nombre_colone_big_matrix);
	//getchar();

	/* Pivot sur le block gauche de taille E->a x E->b  (E->mat) */
	int dim = 0;
	const int colonne_max_pour_pivot = E->b;
	int pivot_c;
	for (pivot_c = 0; pivot_c < colonne_max_pour_pivot; pivot_c++)
	{
		int new_pivot_l;
		for (new_pivot_l = dim; new_pivot_l < nombre_ligne_big_matrix; new_pivot_l++)
		{
			if (big_matrix[new_pivot_l][pivot_c])
			{
				if (new_pivot_l!=dim)
				{
					unsigned char *tmp = big_matrix[dim];
					big_matrix[dim]=big_matrix[new_pivot_l];
					big_matrix[new_pivot_l]=tmp;
				}
				int ligne;
				for (ligne = new_pivot_l+1; ligne < nombre_ligne_big_matrix; ligne++)
				{
					if (big_matrix[ligne][pivot_c])
					{
						int colonne;
						for (colonne = pivot_c+1; colonne < nombre_colone_big_matrix; colonne++)
						{
							big_matrix[ligne][colonne]=Multiply(big_matrix[dim][pivot_c],big_matrix[ligne][colonne]) ^ Multiply(big_matrix[ligne][pivot_c],big_matrix[dim][colonne]);
						}
						//big_matrix[ligne][pivot_c]=0;
					}
				}
				dim++;
				break;
			}
		}
	}
	//printf(" Big matrice - après premier pivot\n");
	//PrintMat(big_matrix,nombre_ligne_big_matrix,nombre_colone_big_matrix);
	//getchar();
	
	/* Partie sur block en bas à droite de taille (E->a - dim) x E->a */
	const int dim_save = dim;
	for (pivot_c = colonne_max_pour_pivot; pivot_c < nombre_colone_big_matrix; pivot_c++)
	{
		int new_pivot_l;
		for (new_pivot_l = dim; new_pivot_l < nombre_ligne_big_matrix; new_pivot_l++)
		{
			if (big_matrix[new_pivot_l][pivot_c])
			{
				if (new_pivot_l!=dim)
				{
					unsigned char *tmp = big_matrix[dim];
					big_matrix[dim]=big_matrix[new_pivot_l];
					big_matrix[new_pivot_l]=tmp;
				}
				const unsigned char inv = Inverse(big_matrix[dim][pivot_c]);
				int colonne;
				for (colonne = pivot_c; colonne < nombre_colone_big_matrix; colonne++) big_matrix[dim][colonne]=Multiply(inv,big_matrix[dim][colonne]);
				int ligne;
				for (ligne = dim_save; ligne < dim; ligne++)
				{
					if (big_matrix[ligne][pivot_c])
					{
						for (colonne = pivot_c+1; colonne < nombre_colone_big_matrix; colonne++)
						{
							big_matrix[ligne][colonne]^=Multiply(big_matrix[ligne][pivot_c],big_matrix[dim][colonne]);
						}
						big_matrix[ligne][pivot_c]=0;
					}
				}
				for (ligne = new_pivot_l+1; ligne < nombre_ligne_big_matrix; ligne++)
				{
					if (big_matrix[ligne][pivot_c])
					{
						for (colonne = pivot_c+1; colonne < nombre_colone_big_matrix; colonne++)
						{
							big_matrix[ligne][colonne]^=Multiply(big_matrix[ligne][pivot_c],big_matrix[dim][colonne]);
						}
						big_matrix[ligne][pivot_c]=0;
					}
				}
				dim++;
				break;
			}
		}
	}
	//printf(" Big matrice - après second pivot\n");
	//PrintMat(big_matrix,nombre_ligne_big_matrix,nombre_colone_big_matrix);
	//getchar();

	#if 0
	/* On remplit les sous-matrices */
	E->sousmat = malloc(VARIABLES*sizeof(MatriceCreuse));
	for (i = 0; i < nombre_de_variables; i++)
	{
		const int var = E->varligne[i];
		E->sousmat[var] = NouvelleMatriceCreuse();
		/* on compte le nombre d'éléments non nuls */
		E->sousmat[var]->nnz = 0;
		int ligne;
		for (ligne = dim_save; ligne < dim; ligne++)
		{
			int colonne;
			for (colonne = 0; colonne < E->gen ; colonne++) if (big_matrix[ligne][colonne_max_pour_pivot+E->gen*i+colonne]) E->sousmat[var]->nnz++;
		}
		E->sousmat[var]->IA=malloc((E->gen+1)*sizeof(int));
		E->sousmat[var]->A=malloc(E->sousmat[var]->nnz*sizeof(unsigned char));
		E->sousmat[var]->JA=malloc(E->sousmat[var]->nnz*sizeof(int));
		E->sousmat[var]->ligne=0;
		E->sousmat[var]->IA[0]=0;
		int nnz = 0;
		int colonne;
		for (colonne = 0; colonne < E->gen ; colonne++)
		{
			const int nnz_save = nnz;
			for (ligne = dim_save; ligne < dim; ligne++)
			{
				if (big_matrix[ligne][colonne_max_pour_pivot+E->gen*i+colonne])
				{
					E->sousmat[var]->A[nnz]=big_matrix[ligne][colonne_max_pour_pivot+E->gen*i+colonne];
					E->sousmat[var]->JA[nnz]=ligne-dim_save;
					nnz++;
				}
			}
			if (nnz!=nnz_save)
			{
				E->sousmat[var]->ligne++;
				E->sousmat[var]->IA[E->sousmat[var]->ligne] = nnz;
			}
		}
		//printf("Sous-mat %d : \n",i);
		//PrintMat2(E->sousmat[var]);
	}
	#else

	/* On remplit les sous-matrices */
	E->sousmat = malloc(VARIABLES*sizeof(MatriceCreuse));
	for (i = 0; i < nombre_de_variables; i++)
	{
		const int var = E->varligne[i];
		E->sousmat[var] = NouvelleMatriceCreuse();
		E->sousmat[var]->ligne=malloc(E->gen*sizeof(LigneCreuse));
		E->sousmat[var]->nombre_de_lignes = 0;
		int colonne;
		for (colonne = 0; colonne < E->gen ; colonne++)
		{
			int nnz=0;
			/* on compte le nombre d'éléments non nuls */
			int ligne;
			for (ligne = dim_save; ligne < dim; ligne++)
			{
				if (big_matrix[ligne][colonne_max_pour_pivot+E->gen*i+colonne]) nnz++;
			}
			if (nnz)
			{
				E->sousmat[var]->ligne[E->sousmat[var]->nombre_de_lignes]=NouvelleLigneCreuse();
				E->sousmat[var]->ligne[E->sousmat[var]->nombre_de_lignes]->elements=malloc(nnz*sizeof(COL_COEFF));
				nnz=0;
				for (ligne = dim_save; ligne < dim; ligne++)
				{
					if (big_matrix[ligne][colonne_max_pour_pivot+E->gen*i+colonne])
					{
						E->sousmat[var]->ligne[E->sousmat[var]->nombre_de_lignes]->elements[nnz].coeff=big_matrix[ligne][colonne_max_pour_pivot+E->gen*i+colonne];
						E->sousmat[var]->ligne[E->sousmat[var]->nombre_de_lignes]->elements[nnz].colonne=ligne-dim_save;
						nnz++;
					}
				}
				E->sousmat[var]->ligne[E->sousmat[var]->nombre_de_lignes]->nnz=nnz;
				E->sousmat[var]->nombre_de_lignes++;
			}
		}
	}
	#endif
	E->col_max = dim-dim_save;
	//printf("E->col_max = %d",E->col_max);
	//getchar();
	/* Liberation de la memoire */

	free(big_matrix);
	free(espace_big_matrix);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x88, %rsp
movslq 0x8(%rdi), %rbp
movq %rdi, 0x10(%rsp)
movl 0xc(%rdi), %ebx
leal (%rbx,%rbp), %r15d
leaq (,%rbp,8), %rdi
callq 0x21d0
movq %rax, %r14
movl %r15d, %eax
imull %ebp, %eax
movslq %eax, %rdi
callq 0x21d0
movl %r15d, %r8d
movq %rax, %r13
movslq %r15d, %r15
testq %rbp, %rbp
jle 0x9c3a
xorl %eax, %eax
movq %r13, %rcx
movq %rcx, (%r14,%rax,8)
incq %rax
addq %r15, %rcx
cmpq %rax, %rbp
jne 0x9c2b
testl %ebp, %ebp
movq 0x10(%rsp), %rdi
movq %r14, 0x28(%rsp)
movl %r8d, 0x8(%rsp)
jle 0x9ca9
addl %ebp, %ebx
xorl %r12d, %r12d
movl 0xc(%rdi), %edx
testl %edx, %edx
jle 0x9c78
movq (%r14,%r12,8), %rax
xorl %ecx, %ecx
movq (%rdi), %rdx
movq (%rdx,%r12,8), %rdx
movb (%rdx,%rcx), %dl
movb %dl, (%rax,%rcx)
incq %rcx
movslq 0xc(%rdi), %rdx
cmpq %rdx, %rcx
jl 0x9c5f
cmpl %r8d, %edx
jge 0x9ca1
movslq %edx, %rdi
addq (%r14,%r12,8), %rdi
notl %edx
addl %ebx, %edx
incq %rdx
xorl %esi, %esi
callq 0x2120
movl 0x8(%rsp), %r8d
movq 0x28(%rsp), %r14
movq 0x10(%rsp), %rdi
incq %r12
cmpq %rbp, %r12
jne 0x9c52
movl 0x8(%rdi), %eax
cltd
idivl 0x10(%rdi)
testl %eax, %eax
jle 0x9ced
xorl %r9d, %r9d
movl 0x10(%rdi), %edx
testl %edx, %edx
jle 0x9ce5
xorl %ecx, %ecx
imull %r9d, %edx
movslq %edx, %rdx
leaq (%rcx,%rdx), %rsi
movq (%r14,%rsi,8), %rsi
addl 0xc(%rdi), %edx
addl %ecx, %edx
movslq %edx, %rdx
movb $0x1, (%rsi,%rdx)
movl 0x10(%rdi), %edx
incq %rcx
cmpl %edx, %ecx
jl 0x9cc0
incl %r9d
cmpl %eax, %r9d
jne 0x9cb7
movl %eax, 0x4c(%rsp)
movslq 0xc(%rdi), %rcx
movq %rcx, 0x38(%rsp)
testq %rcx, %rcx
movq %rbp, 0x30(%rsp)
movq %r13, 0x58(%rsp)
jle 0x9e05
leaq -0x1(%r15), %rax
xorl %ecx, %ecx
leaq 0xf2f6(%rip), %rdx      # 0x19010
xorl %r11d, %r11d
movl 0x8(%rsp), %r13d
cmpl %ebp, %r11d
jge 0x9df5
movslq %r11d, %rsi
leaq (%r14,%rsi,8), %r8
movl %ebp, %r9d
subl %r11d, %r9d
xorl %edi, %edi
movq (%r8,%rdi,8), %r10
cmpb $0x0, (%r10,%rcx)
jne 0x9d52
incq %rdi
cmpl %edi, %r9d
jne 0x9d3a
jmp 0x9df0
testl %edi, %edi
je 0x9d62
movq (%r14,%rsi,8), %r9
movq %r10, (%r14,%rsi,8)
movq %r9, (%r8,%rdi,8)
movl %r11d, %r8d
incl %r11d
movq %r11, 0x20(%rsp)
addl %edi, %r8d
incl %r8d
cmpl %ebp, %r8d
jge 0x9deb
addq 0x20(%rsp), %rdi
movslq %edi, %rdi
movq (%r14,%rdi,8), %r8
cmpb $0x0, (%r8,%rcx)
je 0x9dd9
cmpq %rax, %rcx
jge 0x9dd9
movq 0x28(%rsp), %r9
movq (%r9,%rsi,8), %r9
movq %rcx, %r10
movzbl (%r9,%rcx), %r11d
movzbl 0x1(%r8,%r10), %ebx
shll $0x8, %r11d
addq %rdx, %r11
movzbl (%r8,%rcx), %r14d
movzbl 0x1(%r9,%r10), %r12d
shll $0x8, %r14d
addq %rdx, %r14
movb (%r12,%r14), %bpl
xorb (%rbx,%r11), %bpl
movb %bpl, 0x1(%r8,%r10)
incq %r10
leal 0x1(%r10), %r11d
cmpl %r11d, %r13d
jg 0x9d9c
incq %rdi
movq 0x30(%rsp), %rbp
cmpq %rbp, %rdi
movq 0x28(%rsp), %r14
jl 0x9d80
movq 0x20(%rsp), %r11
movl 0x8(%rsp), %r8d
incq %rcx
cmpq 0x38(%rsp), %rcx
jne 0x9d22
jmp 0x9e08
xorl %r11d, %r11d
movl %r11d, %edx
cmpl %r8d, 0x38(%rsp)
movq %r11, 0x20(%rsp)
jge 0x9fc6
movslq %r11d, %rax
movq %rax, 0x60(%rsp)
leaq -0x1(%r15), %r9
leaq 0xf1e2(%rip), %r12      # 0x19010
movq 0x38(%rsp), %r10
movq %r9, 0x68(%rsp)
movq %r10, %r13
incq %r10
cmpl %ebp, %r11d
jge 0x9e73
movslq %r10d, %rsi
movslq %r11d, %rdi
leaq (%r14,%rdi,8), %rax
movl %ebp, %ecx
subl %r11d, %ecx
xorl %ebp, %ebp
movq (%rax,%rbp,8), %rdx
cmpb $0x0, (%rdx,%r13)
jne 0x9e7b
incq %rbp
cmpl %ebp, %ecx
jne 0x9e54
movl %r11d, %edx
movq 0x30(%rsp), %rbp
jmp 0x9fba
movl %r11d, %edx
jmp 0x9fba
movq %rsi, 0x78(%rsp)
movq %r10, 0x80(%rsp)
testl %ebp, %ebp
je 0x9e98
movq (%r14,%rdi,8), %rcx
movq %rdx, (%r14,%rdi,8)
movq %rcx, (%rax,%rbp,8)
movq %r11, 0x50(%rsp)
leal 0x1(%r11), %eax
movq %rax, 0x18(%rsp)
addq %rbp, %rax
movq %rax, 0x70(%rsp)
movq %rdi, 0x40(%rsp)
movq (%r14,%rdi,8), %rbx
movzbl (%rbx,%r13), %edi
callq 0x2537
movzbl %al, %eax
shll $0x8, %eax
movq %r13, %rcx
movzbl (%rbx,%rcx), %edx
leaq (%rax,%r12), %rsi
movb (%rdx,%rsi), %dl
movb %dl, (%rbx,%rcx)
incq %rcx
cmpq %r15, %rcx
jl 0x9eca
movq 0x50(%rsp), %rax
cmpl %eax, 0x20(%rsp)
movq 0x28(%rsp), %r14
movl 0x8(%rsp), %r8d
movq 0x18(%rsp), %rdx
movq 0x68(%rsp), %r9
movq 0x80(%rsp), %r10
movq 0x78(%rsp), %r11
movq 0x40(%rsp), %rsi
jge 0x9f5f
movq 0x60(%rsp), %rax
movq (%r14,%rax,8), %rcx
cmpb $0x0, (%rcx,%r13)
je 0x9f57
cmpq %r9, %r13
jge 0x9f48
movq %r11, %rdx
movzbl (%rcx,%r13), %esi
movzbl (%rbx,%rdx), %edi
shll $0x8, %esi
addq %r12, %rsi
movb (%rdi,%rsi), %sil
xorb %sil, (%rcx,%rdx)
incq %rdx
cmpq %r15, %rdx
jl 0x9f29
movb $0x0, (%rcx,%r13)
movq 0x18(%rsp), %rdx
movq 0x40(%rsp), %rsi
incq %rax
cmpq %rsi, %rax
jne 0x9f16
movq 0x50(%rsp), %rax
addl %ebp, %eax
incl %eax
movq 0x30(%rsp), %rbp
cmpl %ebp, %eax
jge 0x9fba
movslq 0x70(%rsp), %rax
movq (%r14,%rax,8), %rcx
cmpb $0x0, (%rcx,%r13)
je 0x9fb2
cmpq %r9, %r13
jge 0x9fa8
movq %r11, %rdx
movzbl (%rcx,%r13), %esi
movzbl (%rbx,%rdx), %edi
shll $0x8, %esi
addq %r12, %rsi
movb (%rdi,%rsi), %sil
xorb %sil, (%rcx,%rdx)
incq %rdx
cmpq %r15, %rdx
jl 0x9f89
movb $0x0, (%rcx,%r13)
movq 0x18(%rsp), %rdx
incq %rax
cmpq %rbp, %rax
jl 0x9f76
movl %edx, %r11d
cmpl %r10d, %r8d
jne 0x9e38
movq %rdx, 0x18(%rsp)
leaq 0x282e6(%rip), %rax     # 0x322b8
movslq (%rax), %rdi
shlq $0x3, %rdi
callq 0x21d0
movq 0x18(%rsp), %rsi
movq 0x20(%rsp), %r11
movq 0x28(%rsp), %rdi
movq 0x10(%rsp), %r10
movq %rax, 0x30(%r10)
movl 0x4c(%rsp), %eax
testl %eax, %eax
jle 0xa1ab
movslq %r11d, %rcx
movl %eax, %eax
movq %rax, 0x40(%rsp)
movslq %esi, %rbp
leaq (%rdi,%rcx,8), %r12
movq %rbp, %r15
movq %rcx, 0x30(%rsp)
subq %rcx, %r15
xorl %ecx, %ecx
movq 0x20(%r10), %rax
movq %rcx, 0x8(%rsp)
movslq (%rax,%rcx,4), %r13
callq 0x9a8a
movq 0x10(%rsp), %rcx
movq 0x30(%rcx), %rcx
movq %rax, (%rcx,%r13,8)
movq 0x10(%rsp), %rax
movslq 0x10(%rax), %rdi
shlq $0x3, %rdi
callq 0x21d0
movq 0x18(%rsp), %rsi
movq 0x20(%rsp), %r11
movq 0x28(%rsp), %rdi
movq 0x10(%rsp), %r10
movq 0x30(%r10), %rcx
movq (%rcx,%r13,8), %rdx
movq %rax, 0x8(%rdx)
movq (%rcx,%r13,8), %rax
movl $0x0, (%rax)
movl 0x10(%r10), %eax
testl %eax, %eax
jle 0xa198
xorl %ebx, %ebx
movl $0x0, %r14d
cmpl %r11d, %esi
jle 0xa0c0
imull 0x8(%rsp), %eax
movq 0x38(%rsp), %rcx
addl %ebx, %ecx
addl %ecx, %eax
cltq
xorl %r14d, %r14d
movq 0x30(%rsp), %rcx
movq (%rdi,%rcx,8), %rdx
cmpb $0x1, (%rdx,%rax)
sbbl $-0x1, %r14d
incq %rcx
cmpq %rcx, %rbp
jne 0xa0ac
testl %r14d, %r14d
je 0xa18a
callq 0x996c
movq 0x10(%rsp), %rcx
movq 0x30(%rcx), %rcx
movq (%rcx,%r13,8), %rcx
movq 0x8(%rcx), %rdx
movslq (%rcx), %rcx
movq %rax, (%rdx,%rcx,8)
movslq %r14d, %rdi
shlq $0x3, %rdi
callq 0x21d0
movq 0x20(%rsp), %r11
movq 0x10(%rsp), %r10
xorl %edx, %edx
movq 0x30(%r10), %rcx
movq (%rcx,%r13,8), %rsi
movq 0x8(%rsi), %rdi
movslq (%rsi), %rsi
movq (%rdi,%rsi,8), %rsi
movq %rax, 0x8(%rsi)
cmpl %r11d, 0x18(%rsp)
jle 0xa16a
movl 0x10(%r10), %eax
imull 0x8(%rsp), %eax
movq 0x38(%rsp), %rdx
addl %ebx, %edx
addl %edx, %eax
cltq
xorl %esi, %esi
xorl %edx, %edx
movq (%r12,%rsi,8), %rdi
movb (%rdi,%rax), %dil
testb %dil, %dil
je 0xa162
movq (%rcx,%r13,8), %r8
movq 0x8(%r8), %r9
movslq (%r8), %r8
movq (%r9,%r8,8), %r8
movq 0x8(%r8), %r8
movslq %edx, %rdx
movb %dil, (%r8,%rdx,8)
movl %esi, 0x4(%r8,%rdx,8)
incl %edx
incq %rsi
cmpq %rsi, %r15
jne 0xa134
movq (%rcx,%r13,8), %rax
movq 0x8(%rax), %rcx
movslq (%rax), %rsi
movq (%rcx,%rsi,8), %rcx
movl %edx, (%rcx)
leal 0x1(%rsi), %ecx
movl %ecx, (%rax)
movq 0x28(%rsp), %rdi
movq 0x18(%rsp), %rsi
incl %ebx
movl 0x10(%r10), %eax
cmpl %eax, %ebx
jl 0xa089
movq 0x8(%rsp), %rcx
incq %rcx
cmpq 0x40(%rsp), %rcx
jne 0xa020
subl %r11d, %esi
movl %esi, 0x38(%r10)
callq 0x2030
movq 0x58(%rsp), %rdi
addq $0x88, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x2030