void normaliseP2(double **P, int N) {
  int i, j;
  double sumP = 0.;
  for (i = 0; i < N; i++)
    for (j = 0; j < N; j++)
      sumP += P[i][j];
  for (i = 0; i < N; i++)
    for (j = 0; j < N; j++)
      P[i][j] /= sumP;
  return;
}


x86 assembly (clang, -O0):
<normaliseP2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1179 <normaliseP2+0x79>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1166 <normaliseP2+0x66>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    112d <normaliseP2+0x2d>
jmp    116b <normaliseP2+0x6b>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    111a <normaliseP2+0x1a>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11e3 <normaliseP2+0xe3>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11d0 <normaliseP2+0xd0>
movsd  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1193 <normaliseP2+0x93>
jmp    11d5 <normaliseP2+0xd5>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1180 <normaliseP2+0x80>
pop    %rbp
ret

x86 assembly (clang, -O2):
<normaliseP2>:
push   %r14
push   %rbx
test   %esi,%esi
jle    1286 <normaliseP2+0x186>
mov    %esi,%r14d
lea    -0x1(%r14),%r8
mov    %r14d,%edx
and    $0x7,%edx
mov    %r14d,%eax
and    $0xfffffff8,%eax
xorpd  %xmm0,%xmm0
xor    %r9d,%r9d
jmp    1139 <normaliseP2+0x39>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r9
cmp    %r14,%r9
je     11b0 <normaliseP2+0xb0>
mov    (%rdi,%r9,8),%rcx
xor    %ebx,%ebx
cmp    $0x7,%r8
jb     1188 <normaliseP2+0x88>
cs nopw 0x0(%rax,%rax,1)
nop
addsd  (%rcx,%rbx,8),%xmm0
addsd  0x8(%rcx,%rbx,8),%xmm0
addsd  0x10(%rcx,%rbx,8),%xmm0
addsd  0x18(%rcx,%rbx,8),%xmm0
addsd  0x20(%rcx,%rbx,8),%xmm0
addsd  0x28(%rcx,%rbx,8),%xmm0
addsd  0x30(%rcx,%rbx,8),%xmm0
addsd  0x38(%rcx,%rbx,8),%xmm0
add    $0x8,%rbx
cmp    %rbx,%rax
jne    1150 <normaliseP2+0x50>
test   %rdx,%rdx
je     1130 <normaliseP2+0x30>
lea    (%rcx,%rbx,8),%rcx
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rcx,%rbx,8),%xmm0
add    $0x1,%rbx
cmp    %rbx,%rdx
jne    11a0 <normaliseP2+0xa0>
jmp    1130 <normaliseP2+0x30>
test   %esi,%esi
jle    1286 <normaliseP2+0x186>
mov    %r14d,%r11d
and    $0xfffffffe,%r11d
lea    -0x2(%r11),%r8
mov    %r8,%r10
shr    %r10
add    $0x1,%r10
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mov    %r10,%r9
and    $0xfffffffffffffffe,%r9
xor    %edx,%edx
jmp    11ed <normaliseP2+0xed>
add    $0x1,%rdx
cmp    %r14,%rdx
je     1286 <normaliseP2+0x186>
mov    (%rdi,%rdx,8),%rax
cmp    $0x2,%esi
jae    1200 <normaliseP2+0x100>
xor    %ecx,%ecx
jmp    1260 <normaliseP2+0x160>
nopw   0x0(%rax,%rax,1)
test   %r8,%r8
je     127c <normaliseP2+0x17c>
mov    %r9,%rbx
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movupd (%rax,%rcx,8),%xmm2
movupd 0x10(%rax,%rcx,8),%xmm3
divpd  %xmm1,%xmm2
movupd %xmm2,(%rax,%rcx,8)
divpd  %xmm1,%xmm3
movupd %xmm3,0x10(%rax,%rcx,8)
add    $0x4,%rcx
add    $0xfffffffffffffffe,%rbx
jne    1210 <normaliseP2+0x110>
test   $0x1,%r10b
je     124c <normaliseP2+0x14c>
movupd (%rax,%rcx,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rax,%rcx,8)
mov    %r11,%rcx
cmp    %r14,%r11
je     11e0 <normaliseP2+0xe0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rax,%rcx,8),%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,(%rax,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r14
jne    1260 <normaliseP2+0x160>
jmp    11e0 <normaliseP2+0xe0>
xor    %ecx,%ecx
test   $0x1,%r10b
jne    123e <normaliseP2+0x13e>
jmp    124c <normaliseP2+0x14c>
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @normaliseP2(double** %0, i32 %1) {
%3 = alloca double**
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca double
store double** %0, double*** %3
store i32 %1, i32* %4
store double 0.000000e+00, double* %7
store i32 0, i32* %5
br label %8
8: ; preds = %33, %2
%9 = load i32, i32* %5
%10 = load i32, i32* %4
%11 = icmp slt i32 %9, %10
br i1 %11, label %12, label %36
12: ; preds = %8
store i32 0, i32* %6
br label %13
13: ; preds = %29, %12
%14 = load i32, i32* %6
%15 = load i32, i32* %4
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %32
17: ; preds = %13
%18 = load double**, double*** %3
%19 = load i32, i32* %5
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds double*, double** %18, i64 %20
%22 = load double*, double** %21
%23 = load i32, i32* %6
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %22, i64 %24
%26 = load double, double* %25
%27 = load double, double* %7
%28 = fadd double %27, %26
store double %28, double* %7
br label %29
29: ; preds = %17
%30 = load i32, i32* %6
%31 = add nsw i32 %30, 1
store i32 %31, i32* %6
br label %13
32: ; preds = %13
br label %33
33: ; preds = %32
%34 = load i32, i32* %5
%35 = add nsw i32 %34, 1
store i32 %35, i32* %5
br label %8
36: ; preds = %8
store i32 0, i32* %5
br label %37
37: ; preds = %62, %36
%38 = load i32, i32* %5
%39 = load i32, i32* %4
%40 = icmp slt i32 %38, %39
br i1 %40, label %41, label %65
41: ; preds = %37
store i32 0, i32* %6
br label %42
42: ; preds = %58, %41
%43 = load i32, i32* %6
%44 = load i32, i32* %4
%45 = icmp slt i32 %43, %44
br i1 %45, label %46, label %61
46: ; preds = %42
%47 = load double, double* %7
%48 = load double**, double*** %3
%49 = load i32, i32* %5
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double*, double** %48, i64 %50
%52 = load double*, double** %51
%53 = load i32, i32* %6
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %52, i64 %54
%56 = load double, double* %55
%57 = fdiv double %56, %47
store double %57, double* %55
br label %58
58: ; preds = %46
%59 = load i32, i32* %6
%60 = add nsw i32 %59, 1
store i32 %60, i32* %6
br label %42
61: ; preds = %42
br label %62
62: ; preds = %61
%63 = load i32, i32* %5
%64 = add nsw i32 %63, 1
store i32 %64, i32* %5
br label %37
65: ; preds = %37
ret void
}

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