unsigned long add_j_to_k(unsigned long M[][156], unsigned long k,
                         unsigned long j, unsigned long n_blocchi) {
  for (unsigned long I = 0; I < n_blocchi; ++I)
    M[k][I] = M[k][I] ^ M[j][I];
}


x86 assembly (clang, -O0):
<add_j_to_k>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    1180 <add_j_to_k+0x80>
mov    -0x10(%rbp),%rax
imul   $0x4e0,-0x18(%rbp),%rcx
add    %rcx,%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x10(%rbp),%rax
imul   $0x4e0,-0x20(%rbp),%rcx
add    %rcx,%rax
mov    -0x30(%rbp),%rcx
xor    (%rax,%rcx,8),%rdx
mov    -0x10(%rbp),%rax
imul   $0x4e0,-0x18(%rbp),%rcx
add    %rcx,%rax
mov    -0x30(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    111c <add_j_to_k+0x1c>
mov    -0x8(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<add_j_to_k>:
push   %r15
push   %r14
push   %rbx
test   %rcx,%rcx
je     1279 <add_j_to_k+0x179>
cmp    $0x4,%rcx
jb     113c <add_j_to_k+0x3c>
imul   $0x4e0,%rsi,%r15
lea    (%rdi,%r15,1),%r8
imul   $0x4e0,%rdx,%rax
lea    (%rdi,%rax,1),%r9
lea    (%r9,%rcx,8),%rbx
cmp    %rbx,%r8
jae    11b0 <add_j_to_k+0xb0>
lea    (%r8,%rcx,8),%rbx
cmp    %rbx,%r9
jae    11b0 <add_j_to_k+0xb0>
xor    %r10d,%r10d
mov    %r10,%rax
not    %rax
test   $0x1,%cl
je     116a <add_j_to_k+0x6a>
imul   $0x4e0,%rsi,%r8
add    %rdi,%r8
imul   $0x4e0,%rdx,%rbx
add    %rdi,%rbx
mov    (%rbx,%r10,8),%rbx
xor    %rbx,(%r8,%r10,8)
or     $0x1,%r10
add    %rcx,%rax
je     1279 <add_j_to_k+0x179>
imul   $0x4e0,%rdx,%rax
add    %rdi,%rax
add    $0x8,%rax
imul   $0x4e0,%rsi,%rdx
add    %rdi,%rdx
add    $0x8,%rdx
nop
mov    -0x8(%rax,%r10,8),%rsi
xor    %rsi,-0x8(%rdx,%r10,8)
mov    (%rax,%r10,8),%rsi
xor    %rsi,(%rdx,%r10,8)
add    $0x2,%r10
cmp    %r10,%rcx
jne    1190 <add_j_to_k+0x90>
jmp    1279 <add_j_to_k+0x179>
mov    %rcx,%r10
and    $0xfffffffffffffffc,%r10
lea    -0x4(%r10),%rbx
mov    %rbx,%r11
shr    $0x2,%r11
add    $0x1,%r11
test   %rbx,%rbx
je     127f <add_j_to_k+0x17f>
mov    %r11,%r14
and    $0xfffffffffffffffe,%r14
lea    (%rdi,%r15,1),%rbx
add    $0x30,%rbx
lea    (%rdi,%rax,1),%r15
add    $0x30,%r15
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movups -0x30(%r15,%rax,8),%xmm0
movups -0x20(%r15,%rax,8),%xmm1
movups -0x30(%rbx,%rax,8),%xmm2
xorps  %xmm0,%xmm2
movups -0x20(%rbx,%rax,8),%xmm0
xorps  %xmm1,%xmm0
movups -0x10(%rbx,%rax,8),%xmm1
movups (%rbx,%rax,8),%xmm3
movups %xmm2,-0x30(%rbx,%rax,8)
movups %xmm0,-0x20(%rbx,%rax,8)
movups -0x10(%r15,%rax,8),%xmm0
xorps  %xmm1,%xmm0
movups (%r15,%rax,8),%xmm1
xorps  %xmm3,%xmm1
movups %xmm0,-0x10(%rbx,%rax,8)
movups %xmm1,(%rbx,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%r14
jne    11f0 <add_j_to_k+0xf0>
test   $0x1,%r11b
je     1270 <add_j_to_k+0x170>
movups (%r8,%rax,8),%xmm0
movups 0x10(%r8,%rax,8),%xmm1
movups (%r9,%rax,8),%xmm2
xorps  %xmm0,%xmm2
movups 0x10(%r9,%rax,8),%xmm0
xorps  %xmm1,%xmm0
movups %xmm2,(%r8,%rax,8)
movups %xmm0,0x10(%r8,%rax,8)
cmp    %rcx,%r10
jne    113f <add_j_to_k+0x3f>
pop    %rbx
pop    %r14
pop    %r15
ret
xor    %eax,%eax
test   $0x1,%r11b
jne    1249 <add_j_to_k+0x149>
jmp    1270 <add_j_to_k+0x170>

LLVM IR (-O0):
define i64 @add_j_to_k([156 x i64]* %0, i64 %1, i64 %2, i64 %3) {
%5 = alloca i64
%6 = alloca [156 x i64]*
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
store [156 x i64]* %0, [156 x i64]** %6
store i64 %1, i64* %7
store i64 %2, i64* %8
store i64 %3, i64* %9
store i64 0, i64* %10
br label %11
11: ; preds = %34, %4
%12 = load i64, i64* %10
%13 = load i64, i64* %9
%14 = icmp ult i64 %12, %13
br i1 %14, label %15, label %37
15: ; preds = %11
%16 = load [156 x i64]*, [156 x i64]** %6
%17 = load i64, i64* %7
%18 = getelementptr inbounds [156 x i64], [156 x i64]* %16, i64 %17
%19 = load i64, i64* %10
%20 = getelementptr inbounds [156 x i64], [156 x i64]* %18, i64 0, i64 %19
%21 = load i64, i64* %20
%22 = load [156 x i64]*, [156 x i64]** %6
%23 = load i64, i64* %8
%24 = getelementptr inbounds [156 x i64], [156 x i64]* %22, i64 %23
%25 = load i64, i64* %10
%26 = getelementptr inbounds [156 x i64], [156 x i64]* %24, i64 0, i64 %25
%27 = load i64, i64* %26
%28 = xor i64 %21, %27
%29 = load [156 x i64]*, [156 x i64]** %6
%30 = load i64, i64* %7
%31 = getelementptr inbounds [156 x i64], [156 x i64]* %29, i64 %30
%32 = load i64, i64* %10
%33 = getelementptr inbounds [156 x i64], [156 x i64]* %31, i64 0, i64 %32
store i64 %28, i64* %33
br label %34
34: ; preds = %15
%35 = load i64, i64* %10
%36 = add i64 %35, 1
store i64 %36, i64* %10
br label %11
37: ; preds = %11
%38 = load i64, i64* %5
ret i64 %38
}

LLVM IR (-O2):
define i64 @add_j_to_k([156 x i64]* nocapture %0, i64 %1, i64 %2, i64 %3) {
%5 = icmp eq i64 %3, 0
br i1 %5, label %102, label %6
6: ; preds = %4
%7 = icmp ult i64 %3, 4
br i1 %7, label %86, label %8
8: ; preds = %6
%9 = getelementptr [156 x i64], [156 x i64]* %0, i64 %1, i64 0
%10 = getelementptr [156 x i64], [156 x i64]* %0, i64 %1, i64 %3
%11 = getelementptr [156 x i64], [156 x i64]* %0, i64 %2, i64 0
%12 = getelementptr [156 x i64], [156 x i64]* %0, i64 %2, i64 %3
%13 = icmp ult i64* %9, %12
%14 = icmp ult i64* %11, %10
%15 = and i1 %13, %14
br i1 %15, label %86, label %16
16: ; preds = %8
%17 = and i64 %3, -4
%18 = add i64 %17, -4
%19 = lshr exact i64 %18, 2
%20 = add nuw nsw i64 %19, 1
%21 = and i64 %20, 1
%22 = icmp eq i64 %18, 0
br i1 %22, label %64, label %23
23: ; preds = %16
%24 = and i64 %20, 9223372036854775806
br label %25
25: ; preds = %25, %23
%26 = phi i64 [ 0, %23 ], [ %61, %25 ]
%27 = phi i64 [ 0, %23 ], [ %62, %25 ]
%28 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %1, i64 %26
%29 = bitcast i64* %28 to <2 x i64>*
%30 = load <2 x i64>, <2 x i64>* %29
%31 = getelementptr inbounds i64, i64* %28, i64 2
%32 = bitcast i64* %31 to <2 x i64>*
%33 = load <2 x i64>, <2 x i64>* %32
%34 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %2, i64 %26
%35 = bitcast i64* %34 to <2 x i64>*
%36 = load <2 x i64>, <2 x i64>* %35
%37 = getelementptr inbounds i64, i64* %34, i64 2
%38 = bitcast i64* %37 to <2 x i64>*
%39 = load <2 x i64>, <2 x i64>* %38
%40 = xor <2 x i64> %36, %30
%41 = xor <2 x i64> %39, %33
%42 = bitcast i64* %28 to <2 x i64>*
store <2 x i64> %40, <2 x i64>* %42
%43 = bitcast i64* %31 to <2 x i64>*
store <2 x i64> %41, <2 x i64>* %43
%44 = or i64 %26, 4
%45 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %1, i64 %44
%46 = bitcast i64* %45 to <2 x i64>*
%47 = load <2 x i64>, <2 x i64>* %46
%48 = getelementptr inbounds i64, i64* %45, i64 2
%49 = bitcast i64* %48 to <2 x i64>*
%50 = load <2 x i64>, <2 x i64>* %49
%51 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %2, i64 %44
%52 = bitcast i64* %51 to <2 x i64>*
%53 = load <2 x i64>, <2 x i64>* %52
%54 = getelementptr inbounds i64, i64* %51, i64 2
%55 = bitcast i64* %54 to <2 x i64>*
%56 = load <2 x i64>, <2 x i64>* %55
%57 = xor <2 x i64> %53, %47
%58 = xor <2 x i64> %56, %50
%59 = bitcast i64* %45 to <2 x i64>*
store <2 x i64> %57, <2 x i64>* %59
%60 = bitcast i64* %48 to <2 x i64>*
store <2 x i64> %58, <2 x i64>* %60
%61 = add nuw i64 %26, 8
%62 = add i64 %27, 2
%63 = icmp eq i64 %62, %24
br i1 %63, label %64, label %25
64: ; preds = %25, %16
%65 = phi i64 [ 0, %16 ], [ %61, %25 ]
%66 = icmp eq i64 %21, 0
br i1 %66, label %84, label %67
67: ; preds = %64
%68 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %1, i64 %65
%69 = bitcast i64* %68 to <2 x i64>*
%70 = load <2 x i64>, <2 x i64>* %69
%71 = getelementptr inbounds i64, i64* %68, i64 2
%72 = bitcast i64* %71 to <2 x i64>*
%73 = load <2 x i64>, <2 x i64>* %72
%74 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %2, i64 %65
%75 = bitcast i64* %74 to <2 x i64>*
%76 = load <2 x i64>, <2 x i64>* %75
%77 = getelementptr inbounds i64, i64* %74, i64 2
%78 = bitcast i64* %77 to <2 x i64>*
%79 = load <2 x i64>, <2 x i64>* %78
%80 = xor <2 x i64> %76, %70
%81 = xor <2 x i64> %79, %73
%82 = bitcast i64* %68 to <2 x i64>*
store <2 x i64> %80, <2 x i64>* %82
%83 = bitcast i64* %71 to <2 x i64>*
store <2 x i64> %81, <2 x i64>* %83
br label %84
84: ; preds = %64, %67
%85 = icmp eq i64 %17, %3
br i1 %85, label %102, label %86
86: ; preds = %8, %6, %84
%87 = phi i64 [ 0, %8 ], [ 0, %6 ], [ %17, %84 ]
%88 = xor i64 %87, -1
%89 = and i64 %3, 1
%90 = icmp eq i64 %89, 0
br i1 %90, label %98, label %91
91: ; preds = %86
%92 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %1, i64 %87
%93 = load i64, i64* %92
%94 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %2, i64 %87
%95 = load i64, i64* %94
%96 = xor i64 %95, %93
store i64 %96, i64* %92
%97 = or i64 %87, 1
br label %98
98: ; preds = %91, %86
%99 = phi i64 [ %87, %86 ], [ %97, %91 ]
%100 = sub i64 0, %3
%101 = icmp eq i64 %88, %100
br i1 %101, label %102, label %103
102: ; preds = %98, %103, %84, %4
ret i64 undef
103: ; preds = %98, %103
%104 = phi i64 [ %116, %103 ], [ %99, %98 ]
%105 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %1, i64 %104
%106 = load i64, i64* %105
%107 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %2, i64 %104
%108 = load i64, i64* %107
%109 = xor i64 %108, %106
store i64 %109, i64* %105
%110 = add nuw i64 %104, 1
%111 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %1, i64 %110
%112 = load i64, i64* %111
%113 = getelementptr inbounds [156 x i64], [156 x i64]* %0, i64 %2, i64 %110
%114 = load i64, i64* %113
%115 = xor i64 %114, %112
store i64 %115, i64* %111
%116 = add nuw i64 %104, 2
%117 = icmp eq i64 %116, %3
br i1 %117, label %102, label %103
}