// File: /workspace/github_c_makefile_star2_O3/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 $0x78, %rsp
movq %rdi, %r15
movslq 0x8(%rdi), %rbx
movl 0xc(%rdi), %ebp
addl %ebx, %ebp
leaq (,%rbx,8), %rdi
callq 0x21d0
movq %rax, %r14
movl %ebp, %eax
imull %ebx, %eax
movslq %eax, %rdi
callq 0x21d0
movl %ebp, %esi
movq %rbx, %r11
movq %rax, 0x50(%rsp)
movslq %ebp, %rbp
movl %r11d, %eax
testq %rbx, %rbx
movq %rbx, (%rsp)
movl %esi, 0x18(%rsp)
jle 0x9f73
xorl %eax, %eax
movq 0x50(%rsp), %rcx
movq %rcx, (%r14,%rax,8)
incq %rax
addq %rbp, %rcx
cmpq %rax, %r11
jne 0x9f06
movl %r11d, %eax
testl %r11d, %r11d
jle 0x9f73
xorl %ebx, %ebx
movl 0xc(%r15), %edx
testl %edx, %edx
jle 0x9f46
movq (%r14,%rbx,8), %rax
xorl %ecx, %ecx
movq (%r15), %rdx
movq (%rdx,%rbx,8), %rdx
movb (%rdx,%rcx), %dl
movb %dl, (%rax,%rcx)
incq %rcx
movslq 0xc(%r15), %rdx
cmpq %rdx, %rcx
jl 0x9f2d
cmpl %esi, %edx
jge 0x9f67
movslq %edx, %rdi
addq (%r14,%rbx,8), %rdi
notl %edx
addl %esi, %edx
incq %rdx
xorl %esi, %esi
callq 0x2120
movl 0x18(%rsp), %esi
movq (%rsp), %r11
incq %rbx
cmpq %r11, %rbx
jne 0x9f1f
movl 0x8(%r15), %eax
movl 0x10(%r15), %ecx
cltd
idivl %ecx
testl %eax, %eax
jle 0x9fb0
xorl %edi, %edi
testl %ecx, %ecx
jle 0x9faa
xorl %edx, %edx
imull %edi, %ecx
leal (%rdx,%rcx), %esi
movslq %esi, %rsi
movq (%r14,%rsi,8), %rsi
addl 0xc(%r15), %ecx
addl %edx, %ecx
movslq %ecx, %rcx
movb $0x1, (%rsi,%rcx)
incl %edx
movl 0x10(%r15), %ecx
cmpl %ecx, %edx
jl 0x9f86
incl %edi
cmpl %eax, %edi
jne 0x9f80
movl %eax, 0x4c(%rsp)
movslq 0xc(%r15), %rbx
testq %rbx, %rbx
movq %r15, 0x38(%rsp)
movq %rbx, 0x40(%rsp)
jle 0xa0d1
leaq -0x1(%rbp), %rax
xorl %ecx, %ecx
leaq 0xf038(%rip), %rdx      # 0x19010
movq $0x0, 0x8(%rsp)
movl %r11d, %r10d
movq 0x8(%rsp), %rsi
subl %esi, %r10d
jle 0xa0c3
movslq %esi, %rdi
incl %esi
leaq (%r14,%rdi,8), %r9
xorl %r8d, %r8d
movq (%r9,%r8,8), %r11
cmpb $0x0, (%r11,%rcx)
jne 0xa01a
incq %r8
cmpl %r8d, %r10d
jne 0x9ffe
movq (%rsp), %r11
jmp 0xa0c3
testl %r8d, %r8d
je 0xa02b
movq (%r14,%rdi,8), %r10
movq %r11, (%r14,%rdi,8)
movq %r10, (%r9,%r8,8)
leal (%rsi,%r8), %r9d
movq (%rsp), %r11
cmpl %r11d, %r9d
jge 0xa0be
movq %rsi, 0x20(%rsp)
movl %esi, %r9d
addq %r8, %r9
movslq %r9d, %r8
movq (%r14,%r8,8), %r9
cmpb $0x0, (%r9,%rcx)
movl 0x18(%rsp), %esi
je 0xa0a1
cmpq %rax, %rcx
jge 0xa0a1
movq (%r14,%rdi,8), %r10
movq %rcx, %r11
movzbl (%r10,%rcx), %ebx
movzbl 0x1(%r9,%r11), %r15d
shll $0x8, %ebx
addq %rdx, %rbx
movzbl (%r9,%rcx), %r12d
movzbl 0x1(%r10,%r11), %r13d
shll $0x8, %r12d
addq %rdx, %r12
movb (%r13,%r12), %r12b
xorb (%r15,%rbx), %r12b
movb %r12b, 0x1(%r9,%r11)
incq %r11
leal 0x1(%r11), %ebx
cmpl %ebx, %esi
jg 0xa065
incq %r8
movq (%rsp), %r11
cmpq %r11, %r8
jl 0xa04a
movq 0x20(%rsp), %rsi
movq %rsi, 0x8(%rsp)
movq 0x40(%rsp), %rbx
jmp 0xa0c3
movq %rsi, 0x8(%rsp)
incq %rcx
cmpq %rbx, %rcx
jne 0x9fe1
jmp 0xa0da
movq $0x0, 0x8(%rsp)
movq 0x8(%rsp), %rax
movl %eax, %ecx
movq %rcx, 0x10(%rsp)
movl 0x18(%rsp), %r8d
cmpl %r8d, %ebx
jge 0xa28a
movslq %eax, %rcx
movq %rcx, 0x60(%rsp)
leaq -0x1(%rbp), %r9
leaq 0xef09(%rip), %r12      # 0x19010
movq 0x40(%rsp), %rsi
movq %rax, 0x10(%rsp)
movq %r9, 0x58(%rsp)
movq %rsi, %r13
incq %rsi
movl %r11d, %eax
movq 0x10(%rsp), %rcx
subl %ecx, %eax
jle 0xa281
movslq %esi, %rdi
movslq %ecx, %rdx
leal 0x1(%rcx), %r10d
movq %rdx, 0x20(%rsp)
leaq (%r14,%rdx,8), %rcx
xorl %r15d, %r15d
movq (%rcx,%r15,8), %rdx
cmpb $0x0, (%rdx,%r13)
jne 0xa15a
incq %r15
cmpl %r15d, %eax
jne 0xa142
jmp 0xa281
movq %rdi, 0x70(%rsp)
movq %rsi, 0x28(%rsp)
testl %r15d, %r15d
movq 0x20(%rsp), %rsi
je 0xa17a
movq (%r14,%rsi,8), %rax
movq %rdx, (%r14,%rsi,8)
movq %rax, (%rcx,%r15,8)
movl %r10d, 0x30(%rsp)
movl %r10d, %eax
addq %r15, %rax
movq %rax, 0x68(%rsp)
movq (%r14,%rsi,8), %rbx
movzbl (%rbx,%r13), %edi
callq 0x2534
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 %rbp, %rcx
jl 0xa1a1
movq 0x60(%rsp), %rax
movq 0x10(%rsp), %rcx
cmpl %ecx, 0x8(%rsp)
movq (%rsp), %r11
movl 0x18(%rsp), %r8d
movq 0x58(%rsp), %r9
movq 0x70(%rsp), %r10
movq 0x20(%rsp), %rdx
jge 0xa223
movq (%r14,%rax,8), %rcx
cmpb $0x0, (%rcx,%r13)
je 0xa21b
cmpq %r9, %r13
jge 0xa211
movq %r10, %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 %rbp, %rdx
jl 0xa1f2
movb $0x0, (%rcx,%r13)
movq 0x20(%rsp), %rdx
incq %rax
cmpq %rdx, %rax
jne 0xa1df
movl 0x30(%rsp), %eax
addl %eax, %r15d
cmpl %r11d, %r15d
jge 0xa277
movslq 0x68(%rsp), %rax
movq (%r14,%rax,8), %rcx
cmpb $0x0, (%rcx,%r13)
je 0xa26b
cmpq %r9, %r13
jge 0xa266
movq %r10, %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 %rbp, %rdx
jl 0xa247
movb $0x0, (%rcx,%r13)
incq %rax
cmpq %r11, %rax
jl 0xa234
movl 0x30(%rsp), %eax
movq %rax, 0x10(%rsp)
movq 0x28(%rsp), %rsi
cmpl %esi, %r8d
jne 0xa116
leaq 0x28027(%rip), %rax     # 0x322b8
movslq (%rax), %rdi
shlq $0x3, %rdi
callq 0x21d0
movq 0x8(%rsp), %r10
movq 0x38(%rsp), %rcx
movq %rax, 0x30(%rcx)
movl 0x4c(%rsp), %eax
testl %eax, %eax
movq 0x10(%rsp), %r11
jle 0xa471
movslq %r10d, %rcx
movl %eax, %eax
movq %rax, 0x28(%rsp)
movslq %r11d, %rbp
leaq (%r14,%rcx,8), %rax
movq %rax, 0x30(%rsp)
movq %rbp, %r15
movq %rcx, 0x20(%rsp)
subq %rcx, %r15
xorl %esi, %esi
movq %rsi, 0x18(%rsp)
movq 0x38(%rsp), %rbx
movq 0x20(%rbx), %rax
movslq (%rax,%rsi,4), %r13
callq 0x9d55
movq 0x30(%rbx), %rcx
movq %rax, (%rcx,%r13,8)
movslq 0x10(%rbx), %r12
leaq (,%r12,8), %rdi
callq 0x21d0
movq 0x18(%rsp), %rsi
movq 0x10(%rsp), %r11
movq 0x8(%rsp), %r10
movq 0x30(%rbx), %rcx
movq (%rcx,%r13,8), %rdx
movq %rax, 0x8(%rdx)
movq %r13, (%rsp)
movq (%rcx,%r13,8), %rax
movl $0x0, (%rax)
testq %r12, %r12
jle 0xa463
xorl %ebx, %ebx
cmpl %r10d, %r11d
jle 0xa458
movl %r12d, %eax
imull %esi, %eax
movq 0x40(%rsp), %rcx
addl %ebx, %ecx
addl %ecx, %eax
cltq
xorl %r13d, %r13d
movq 0x20(%rsp), %rcx
movq (%r14,%rcx,8), %rdx
cmpb $0x1, (%rdx,%rax)
sbbl $-0x1, %r13d
incq %rcx
cmpq %rcx, %rbp
jne 0xa365
testl %r13d, %r13d
je 0xa458
callq 0x9c30
movq 0x38(%rsp), %r12
movq 0x30(%r12), %rcx
movq (%rsp), %rdx
movq (%rcx,%rdx,8), %rcx
movq 0x8(%rcx), %rdx
movslq (%rcx), %rcx
movq %rax, (%rdx,%rcx,8)
movl %r13d, %edi
shlq $0x3, %rdi
callq 0x21d0
movq 0x10(%rsp), %r11
movq 0x8(%rsp), %r10
xorl %edx, %edx
movq 0x30(%r12), %rcx
movq (%rsp), %rsi
movq (%rcx,%rsi,8), %rsi
movq 0x8(%rsi), %rdi
movslq (%rsi), %rsi
movq (%rdi,%rsi,8), %rsi
movq %rax, 0x8(%rsi)
movl 0x10(%r12), %r12d
cmpl %r10d, %r11d
movq 0x30(%rsp), %r13
jle 0xa439
movl %r12d, %eax
imull 0x18(%rsp), %eax
movq 0x40(%rsp), %rdx
addl %ebx, %edx
addl %edx, %eax
cltq
xorl %esi, %esi
xorl %edx, %edx
movq (%r13,%rsi,8), %rdi
movb (%rdi,%rax), %dil
testb %dil, %dil
je 0xa431
movq (%rsp), %r8
movq (%rcx,%r8,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 0xa3fe
movq (%rsp), %rax
movq (%rcx,%rax,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 0x18(%rsp), %rsi
incl %ebx
cmpl %r12d, %ebx
jl 0xa343
incq %rsi
cmpq 0x28(%rsp), %rsi
jne 0xa2df
subl %r10d, %r11d
movq 0x38(%rsp), %rax
movl %r11d, 0x38(%rax)
movq %r14, %rdi
callq 0x2030
movq 0x50(%rsp), %rdi
addq $0x78, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x2030