void add(double *ancMatrixUpdated, double *ancMatrix, int nr, int source,
         int sink) {
  int j, p;
  ancMatrixUpdated[source + nr * sink] = 1;
  for (j = 0; j < nr; j++) {
    if (ancMatrix[sink + j * nr] == 1)
      ancMatrixUpdated[source + j * nr] = 1;
  }
  for (p = 0; p < nr; p++) {
    if (ancMatrix[p + source * nr] == 1) {
      ancMatrixUpdated[p + sink * nr] = 1;
      for (j = 0; j < nr; j++) {
        if (ancMatrix[sink + j * nr] == 1)
          ancMatrixUpdated[p + j * nr] = 1;
      }
    }
  }
}


x86 assembly (clang, -O0):
<add>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xecf(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11ac <add+0xac>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe97(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1199 <add+0x99>
jp     1199 <add+0x99>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
mov    -0x20(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xe6c(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    119e <add+0x9e>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    113d <add+0x3d>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    129d <add+0x19d>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x18(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe21(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    128a <add+0x18a>
jp     128a <add+0x18a>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xdf6(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1285 <add+0x185>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xdbe(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1272 <add+0x172>
jp     1272 <add+0x172>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x20(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xd93(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1277 <add+0x177>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1216 <add+0x116>
jmp    128a <add+0x18a>
jmp    128f <add+0x18f>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11b3 <add+0xb3>
pop    %rbp
ret

x86 assembly (clang, -O2):
<add>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %ecx,%r10d
mov    %r8d,%r9d
imul   %edx,%r9d
lea    (%r9,%r10,1),%eax
cltq
movabs $0x3ff0000000000000,%r15
mov    %rdi,-0x20(%rsp)
mov    %r15,(%rdi,%rax,8)
test   %edx,%edx
jle    131f <add+0x21f>
mov    %edx,%r14d
movslq %r8d,%rbx
movslq %r10d,%rcx
mov    %r14,%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x8(%rsp)
mov    %rbx,-0x18(%rsp)
mov    %rcx,-0x10(%rsp)
je     11c4 <add+0xc4>
mov    %r14d,%r12d
and    $0xfffffffe,%r12d
lea    (%rcx,%r14,1),%rax
mov    -0x20(%rsp),%r8
lea    (%r8,%rax,8),%rdi
mov    %r14,%rax
shl    $0x4,%rax
lea    (%rbx,%r14,1),%rbp
lea    (%rsi,%rbp,8),%rbp
lea    (%r8,%rcx,8),%r13
lea    (%rsi,%rbx,8),%r8
xor    %r11d,%r11d
movsd  0xe74(%rip),%xmm0
xor    %ebx,%ebx
jmp    119c <add+0x9c>
add    $0x2,%rbx
add    %rax,%r11
cmp    %rbx,%r12
je     11c6 <add+0xc6>
movsd  (%r8,%r11,1),%xmm1
ucomisd %xmm0,%xmm1
jne    11af <add+0xaf>
jp     11af <add+0xaf>
mov    %r15,0x0(%r13,%r11,1)
movsd  0x0(%rbp,%r11,1),%xmm1
ucomisd %xmm0,%xmm1
jne    1190 <add+0x90>
jp     1190 <add+0x90>
mov    %r15,(%rdi,%r11,1)
jmp    1190 <add+0x90>
xor    %ebx,%ebx
test   $0x1,%r14b
mov    -0x18(%rsp),%rax
je     11f7 <add+0xf7>
imul   %r14,%rbx
add    %rbx,%rax
movsd  (%rsi,%rax,8),%xmm0
ucomisd 0xe1b(%rip),%xmm0
jne    11f7 <add+0xf7>
jp     11f7 <add+0xf7>
add    -0x10(%rsp),%rbx
mov    -0x20(%rsp),%rax
mov    %r15,(%rax,%rbx,8)
test   %edx,%edx
jle    131f <add+0x21f>
imul   %edx,%r10d
movslq %r10d,%r11
movslq %r9d,%rax
mov    %rax,-0x10(%rsp)
mov    %r14d,%r12d
and    $0xfffffffe,%r12d
mov    -0x20(%rsp),%rbp
lea    0x0(,%r14,8),%rbx
add    %rbp,%rbx
mov    %r14,%r9
shl    $0x4,%r9
mov    -0x18(%rsp),%rcx
lea    (%rcx,%r14,1),%rax
lea    (%rsi,%rax,8),%rax
lea    (%rsi,%rcx,8),%rdx
xor    %r13d,%r13d
movsd  0xdb8(%rip),%xmm0
mov    -0x20(%rsp),%rdi
jmp    1265 <add+0x165>
nop
add    $0x1,%r13
add    $0x8,%rbx
add    $0x8,%rbp
cmp    %r14,%r13
je     131f <add+0x21f>
lea    (%r11,%r13,1),%rcx
movsd  (%rsi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1250 <add+0x150>
jp     1250 <add+0x150>
mov    -0x10(%rsp),%rcx
add    %r13,%rcx
mov    %r15,(%rdi,%rcx,8)
cmpq   $0x0,-0x8(%rsp)
je     12e0 <add+0x1e0>
xor    %r10d,%r10d
xor    %r8d,%r8d
jmp    12ac <add+0x1ac>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x2,%r8
add    %r9,%r10
cmp    %r8,%r12
je     12e3 <add+0x1e3>
movsd  (%rdx,%r10,1),%xmm1
ucomisd %xmm0,%xmm1
jne    12bf <add+0x1bf>
jp     12bf <add+0x1bf>
mov    %r15,0x0(%rbp,%r10,1)
movsd  (%rax,%r10,1),%xmm1
ucomisd %xmm0,%xmm1
jne    12a0 <add+0x1a0>
jp     12a0 <add+0x1a0>
mov    %r15,(%rbx,%r10,1)
jmp    12a0 <add+0x1a0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %r8d,%r8d
test   $0x1,%r14b
je     1250 <add+0x150>
imul   %r14,%r8
mov    -0x18(%rsp),%rcx
add    %r8,%rcx
movsd  (%rsi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1250 <add+0x150>
jp     1250 <add+0x150>
add    %r13,%r8
mov    -0x20(%rsp),%rcx
mov    %r15,(%rcx,%r8,8)
jmp    1250 <add+0x150>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @add(double* %0, double* %1, i32 %2, i32 %3, i32 %4) {
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
store double* %0, double** %6
store double* %1, double** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store i32 %4, i32* %10
%13 = load double*, double** %6
%14 = load i32, i32* %9
%15 = load i32, i32* %8
%16 = load i32, i32* %10
%17 = mul nsw i32 %15, %16
%18 = add nsw i32 %14, %17
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds double, double* %13, i64 %19
store double 1.000000e+00, double* %20
store i32 0, i32* %11
br label %21
21: ; preds = %46, %5
%22 = load i32, i32* %11
%23 = load i32, i32* %8
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %49
25: ; preds = %21
%26 = load double*, double** %7
%27 = load i32, i32* %10
%28 = load i32, i32* %11
%29 = load i32, i32* %8
%30 = mul nsw i32 %28, %29
%31 = add nsw i32 %27, %30
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %26, i64 %32
%34 = load double, double* %33
%35 = fcmp oeq double %34, 1.000000e+00
br i1 %35, label %36, label %45
36: ; preds = %25
%37 = load double*, double** %6
%38 = load i32, i32* %9
%39 = load i32, i32* %11
%40 = load i32, i32* %8
%41 = mul nsw i32 %39, %40
%42 = add nsw i32 %38, %41
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %37, i64 %43
store double 1.000000e+00, double* %44
br label %45
45: ; preds = %36, %25
br label %46
46: ; preds = %45
%47 = load i32, i32* %11
%48 = add nsw i32 %47, 1
store i32 %48, i32* %11
br label %21
49: ; preds = %21
store i32 0, i32* %12
br label %50
50: ; preds = %104, %49
%51 = load i32, i32* %12
%52 = load i32, i32* %8
%53 = icmp slt i32 %51, %52
br i1 %53, label %54, label %107
54: ; preds = %50
%55 = load double*, double** %7
%56 = load i32, i32* %12
%57 = load i32, i32* %9
%58 = load i32, i32* %8
%59 = mul nsw i32 %57, %58
%60 = add nsw i32 %56, %59
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %55, i64 %61
%63 = load double, double* %62
%64 = fcmp oeq double %63, 1.000000e+00
br i1 %64, label %65, label %103
65: ; preds = %54
%66 = load double*, double** %6
%67 = load i32, i32* %12
%68 = load i32, i32* %10
%69 = load i32, i32* %8
%70 = mul nsw i32 %68, %69
%71 = add nsw i32 %67, %70
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %66, i64 %72
store double 1.000000e+00, double* %73
store i32 0, i32* %11
br label %74
74: ; preds = %99, %65
%75 = load i32, i32* %11
%76 = load i32, i32* %8
%77 = icmp slt i32 %75, %76
br i1 %77, label %78, label %102
78: ; preds = %74
%79 = load double*, double** %7
%80 = load i32, i32* %10
%81 = load i32, i32* %11
%82 = load i32, i32* %8
%83 = mul nsw i32 %81, %82
%84 = add nsw i32 %80, %83
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %79, i64 %85
%87 = load double, double* %86
%88 = fcmp oeq double %87, 1.000000e+00
br i1 %88, label %89, label %98
89: ; preds = %78
%90 = load double*, double** %6
%91 = load i32, i32* %12
%92 = load i32, i32* %11
%93 = load i32, i32* %8
%94 = mul nsw i32 %92, %93
%95 = add nsw i32 %91, %94
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds double, double* %90, i64 %96
store double 1.000000e+00, double* %97
br label %98
98: ; preds = %89, %78
br label %99
99: ; preds = %98
%100 = load i32, i32* %11
%101 = add nsw i32 %100, 1
store i32 %101, i32* %11
br label %74
102: ; preds = %74
br label %103
103: ; preds = %102, %54
br label %104
104: ; preds = %103
%105 = load i32, i32* %12
%106 = add nsw i32 %105, 1
store i32 %106, i32* %12
br label %50
107: ; preds = %50
ret void
}

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