// File: /workspace/github_c_makefile_star2_O2/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 $0x48, %rsp
movq %rdi, %r14
movslq 0x8(%rdi), %r15
movslq 0xc(%rdi), %rbp
addq %r15, %rbp
leaq (,%r15,8), %rdi
callq 0x21c0
movq %rax, %rbx
movl %ebp, %eax
imull %r15d, %eax
movslq %eax, %rdi
callq 0x21c0
movq %rbx, %r11
movq %rax, %rdx
xorl %ecx, %ecx
testl %r15d, %r15d
movl $0x0, %eax
movq %r15, 0x10(%rsp)
cmovgl %r15d, %eax
movq %rdx, 0x40(%rsp)
cmpq %rcx, %rax
je 0x91c7
movq %rdx, (%r11,%rcx,8)
incq %rcx
addq %rbp, %rdx
jmp 0x91b6
xorl %ecx, %ecx
cmpq %rax, %rcx
je 0x9206
xorl %edx, %edx
movslq 0xc(%r14), %rsi
cmpq %rsi, %rdx
jge 0x91fc
movq (%r14), %rsi
movq (%rsi,%rcx,8), %rsi
movb (%rsi,%rdx), %sil
movq (%r11,%rcx,8), %rdi
movb %sil, (%rdi,%rdx)
incq %rdx
jmp 0x91d0
movq (%r11,%rcx,8), %rdx
movb $0x0, (%rdx,%rsi)
incq %rsi
cmpq %rbp, %rsi
jl 0x91f1
incq %rcx
jmp 0x91c9
movl 0x8(%r14), %eax
movl 0x10(%r14), %ecx
cltd
idivl %ecx
xorl %edi, %edi
testl %eax, %eax
cmovlel %edi, %eax
cmpl %eax, %edi
je 0x9248
xorl %edx, %edx
cmpl %ecx, %edx
jge 0x9244
imull %edi, %ecx
leal (%rdx,%rcx), %esi
movslq %esi, %rsi
movq (%r11,%rsi,8), %rsi
addl 0xc(%r14), %ecx
addl %edx, %ecx
movslq %ecx, %rcx
movb $0x1, (%rsi,%rcx)
incl %edx
movl 0x10(%r14), %ecx
jmp 0x921e
incl %edi
jmp 0x9218
movl %eax, 0x30(%rsp)
movq %r14, 0x28(%rsp)
movslq 0xc(%r14), %rsi
xorl %eax, %eax
testl %esi, %esi
movl $0x0, %ecx
movq %rsi, 0x38(%rsp)
cmovgl %esi, %ecx
pushq $0x1
popq %rdx
leaq 0xdda0(%rip), %r13      # 0x17010
xorl %r15d, %r15d
movq %r11, 0x18(%rsp)
movslq %r15d, %r10
cmpq %rcx, %rax
movq %r10, (%rsp)
je 0x933b
leaq (%r11,%r10,8), %rsi
xorl %edi, %edi
leaq (%r10,%rdi), %r8
cmpq 0x10(%rsp), %r8
jge 0x9330
movq (%rsi,%rdi,8), %r9
cmpb $0x0, (%r9,%rax)
jne 0x92ad
incq %rdi
jmp 0x928e
testl %edi, %edi
je 0x9322
movq (%rsp), %r10
movq (%r11,%r10,8), %r10
movq (%rsp), %rbx
movq %r9, (%r11,%rbx,8)
movq %r10, (%rsi,%rdi,8)
movq (%rsp), %r10
jmp 0x9322
movq 0x8(%r11,%r8,8), %rdi
cmpb $0x0, (%rdi,%rax)
movq %rsi, %r8
je 0x9322
movq %rdx, %r8
cmpl %r8d, %ebp
jle 0x931f
movq (%r11,%r10,8), %r9
movzbl (%r9,%rax), %r10d
movzbl (%rdi,%r8), %r11d
shll $0x8, %r10d
addq %r13, %r10
movzbl (%rdi,%rax), %ebx
movzbl (%r9,%r8), %r9d
shll $0x8, %ebx
addq %r13, %rbx
movb (%r9,%rbx), %r9b
xorb (%r11,%r10), %r9b
movq (%rsp), %r10
movq 0x18(%rsp), %r11
movb %r9b, (%rdi,%r8)
incq %r8
jmp 0x92dc
movq %rsi, %r8
leaq 0x1(%r8), %rsi
cmpq 0x10(%rsp), %rsi
jl 0x92cb
incl %r15d
incq %rax
incq %rdx
jmp 0x9278
movq 0x38(%rsp), %rbx
leaq 0x1(%rbx), %r10
movl %r15d, 0x34(%rsp)
movl %r15d, %r8d
movslq %r8d, %r12
cmpq %rbp, %rbx
jge 0x9471
leaq (%r11,%r12,8), %rax
xorl %ecx, %ecx
leaq (%r12,%rcx), %rdx
cmpq 0x10(%rsp), %rdx
jge 0x9466
movq %rdx, %r14
movq (%rax,%rcx,8), %rdx
cmpb $0x0, (%rdx,%rbx)
jne 0x937f
incq %rcx
jmp 0x935e
movq %r10, 0x20(%rsp)
movl %r8d, 0xc(%rsp)
testl %ecx, %ecx
je 0x9399
movq (%r11,%r12,8), %rsi
movq %rdx, (%r11,%r12,8)
movq %rsi, (%rax,%rcx,8)
movq (%r11,%r12,8), %r15
movzbl (%r15,%rbx), %edi
callq 0x2433
movzbl %al, %eax
shll $0x8, %eax
movq %rbx, %rcx
movq (%rsp), %rdi
movq %r13, %r9
movq 0x20(%rsp), %r10
cmpq %rbp, %rcx
jge 0x93d6
movzbl (%r15,%rcx), %edx
leaq (%rax,%r9), %rsi
movb (%rdx,%rsi), %dl
movb %dl, (%r15,%rcx)
incq %rcx
jmp 0x93bc
movq %rdi, %rax
movq 0x18(%rsp), %r11
movl 0xc(%rsp), %r8d
cmpq %r12, %rax
jge 0x9458
movq (%r11,%rax,8), %rcx
cmpb $0x0, (%rcx,%rbx)
je 0x941a
movq %r10, %rdx
cmpq %rbp, %rdx
jge 0x9416
movzbl (%rcx,%rbx), %esi
movzbl (%r15,%rdx), %edi
shll $0x8, %esi
addq %r9, %rsi
movb (%rdi,%rsi), %sil
xorb %sil, (%rcx,%rdx)
incq %rdx
jmp 0x93f5
movb $0x0, (%rcx,%rbx)
incq %rax
jmp 0x93e3
movq 0x8(%r11,%r14,8), %rcx
cmpb $0x0, (%rcx,%rbx)
movq %rax, %r14
je 0x9458
movq %r10, %rdx
cmpq %rbp, %rdx
jge 0x9451
movzbl (%rcx,%rbx), %esi
movzbl (%r15,%rdx), %edi
shll $0x8, %esi
addq %r9, %rsi
movb (%rdi,%rsi), %sil
xorb %sil, (%rcx,%rdx)
incq %rdx
jmp 0x9430
movb $0x0, (%rcx,%rbx)
movq %rax, %r14
leaq 0x1(%r14), %rax
cmpq 0x10(%rsp), %rax
jl 0x941f
incl %r8d
incq %rbx
incq %r10
jmp 0x934c
movl %r8d, 0xc(%rsp)
leaq 0x26e3b(%rip), %rax     # 0x302b8
movslq (%rax), %rdi
shlq $0x3, %rdi
callq 0x21c0
movq 0x18(%rsp), %rdi
movq 0x28(%rsp), %r14
movq %rax, 0x30(%r14)
movl 0x30(%rsp), %eax
movq %rax, 0x20(%rsp)
movq (%rsp), %rax
leaq (%rdi,%rax,8), %rbp
xorl %ebx, %ebx
cmpq 0x20(%rsp), %rbx
je 0x960f
movq 0x20(%r14), %rax
movslq (%rax,%rbx,4), %r15
callq 0x9018
movq 0x30(%r14), %rcx
movq %rax, (%rcx,%r15,8)
movslq 0x10(%r14), %r13
movq %r13, %rdi
shlq $0x3, %rdi
callq 0x21c0
movq %r13, %r8
movq (%rsp), %r9
movq 0x18(%rsp), %rdi
movq 0x30(%r14), %rcx
movq (%rcx,%r15,8), %rdx
movq %rax, 0x8(%rdx)
movq (%rcx,%r15,8), %rax
andl $0x0, (%rax)
xorl %esi, %esi
cmpl %r8d, %esi
jge 0x9607
movl %r8d, %eax
imull %ebx, %eax
movq 0x38(%rsp), %rcx
leal (%rsi,%rcx), %r13d
addl %r13d, %eax
cltq
xorl %r14d, %r14d
movq %r9, %rcx
cmpq %r12, %rcx
jge 0x9534
movq (%rdi,%rcx,8), %rdx
cmpb $0x1, (%rdx,%rax)
sbbl $-0x1, %r14d
incq %rcx
jmp 0x951e
testl %r14d, %r14d
je 0x95fb
movq %rsi, 0x10(%rsp)
callq 0x8ef8
movq 0x28(%rsp), %rsi
movq 0x30(%rsi), %rcx
movq (%rcx,%r15,8), %rcx
movq 0x8(%rcx), %rdx
movslq (%rcx), %rcx
movq %rax, (%rdx,%rcx,8)
movl %r14d, %edi
movq %rsi, %r14
shlq $0x3, %rdi
callq 0x21c0
movq (%rsp), %r9
movq 0x30(%r14), %rcx
movq (%rcx,%r15,8), %rdx
movq 0x8(%rdx), %rsi
movslq (%rdx), %rdx
movq (%rsi,%rdx,8), %rdx
movq %rax, 0x8(%rdx)
xorl %edx, %edx
xorl %eax, %eax
leaq (%r9,%rdx), %rsi
cmpq %r12, %rsi
jge 0x95d7
movq (%rbp,%rdx,8), %rsi
movl 0x10(%r14), %edi
imull %ebx, %edi
addl %r13d, %edi
movslq %edi, %rdi
movb (%rsi,%rdi), %sil
testb %sil, %sil
je 0x95d2
movq (%rcx,%r15,8), %rdi
movq 0x8(%rdi), %r8
movslq (%rdi), %rdi
movq (%r8,%rdi,8), %rdi
movq 0x8(%rdi), %rdi
movslq %eax, %r8
movb %sil, (%rdi,%r8,8)
movl %edx, 0x4(%rdi,%r8,8)
incl %eax
incq %rdx
jmp 0x958d
movq (%rcx,%r15,8), %rcx
movq 0x8(%rcx), %rdx
movslq (%rcx), %rsi
movq (%rdx,%rsi,8), %rdx
movl %eax, (%rdx)
leal 0x1(%rsi), %eax
movl %eax, (%rcx)
movl 0x10(%r14), %r8d
movq 0x18(%rsp), %rdi
movq 0x10(%rsp), %rsi
incl %esi
movq 0x28(%rsp), %r14
jmp 0x94fb
incq %rbx
jmp 0x94aa
movl 0xc(%rsp), %eax
subl 0x34(%rsp), %eax
movl %eax, 0x38(%r14)
callq 0x2030
movq 0x40(%rsp), %rdi
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x2030