void diffusionmarine(double pyZ[], int pyBord[], int pyDep[], int pyNgbs[][20],
                     double pyEdge[][20], double pyDist[][20], int pyGIDs[],
                     double pyDiff[], int pylocalNb, int pyglobalNb) {
  int i;
  for (i = 0; i < pyglobalNb; i++) {
    pyDiff[i] = 0.;
  }
  int k;
  for (k = 0; k < pylocalNb; k++) {
    int gid = pyGIDs[k];
    int p;
    if (pyBord[gid] > 0.) {
      for (p = 0; p < 20; p++) {
        int ngbid = pyNgbs[gid][p];
        if (ngbid < 0) {
          break;
        }
        if (pyBord[ngbid] > 0) {
          if (pyDep[gid] > 0 && pyZ[gid] > pyZ[ngbid] && pyDist[gid][p] > 0.) {
            pyDiff[gid] +=
                pyEdge[gid][p] * (pyZ[ngbid] - pyZ[gid]) / pyDist[gid][p];
          }
          if (pyDep[ngbid] > 0 && pyZ[gid] < pyZ[ngbid] &&
              pyDist[gid][p] > 0.) {
            pyDiff[gid] +=
                pyEdge[gid][p] * (pyZ[ngbid] - pyZ[gid]) / pyDist[gid][p];
          }
        }
        if (pyBord[ngbid] < 1) {
          if (pyZ[ngbid] < pyZ[gid] && pyDep[gid] > 0 && pyDist[gid][p] > 0.) {
            pyDiff[gid] +=
                pyEdge[gid][p] * (pyZ[ngbid] - pyZ[gid]) / pyDist[gid][p];
          }
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<diffusionmarine>:
push   %rbp
mov    %rsp,%rbp
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x28(%rbp),%eax
jge    115b <diffusionmarine+0x5b>
mov    0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1131 <diffusionmarine+0x31>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    1470 <diffusionmarine+0x370>
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    145d <diffusionmarine+0x35d>
movl   $0x0,-0x40(%rbp)
cmpl   $0x14,-0x40(%rbp)
jge    1458 <diffusionmarine+0x358>
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0x50,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x44(%rbp)
cmpl   $0x0,-0x44(%rbp)
jge    11cf <diffusionmarine+0xcf>
jmp    1458 <diffusionmarine+0x358>
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jle    136a <diffusionmarine+0x26a>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jle    12a1 <diffusionmarine+0x1a1>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    12a1 <diffusionmarine+0x1a1>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    12a1 <diffusionmarine+0x1a1>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jle    1365 <diffusionmarine+0x265>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1365 <diffusionmarine+0x265>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1365 <diffusionmarine+0x265>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    136a <diffusionmarine+0x26a>
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
cmpl   $0x1,(%rax,%rcx,4)
jge    1445 <diffusionmarine+0x345>
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1440 <diffusionmarine+0x340>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jle    1440 <diffusionmarine+0x340>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1440 <diffusionmarine+0x340>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1445 <diffusionmarine+0x345>
jmp    144a <diffusionmarine+0x34a>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    119d <diffusionmarine+0x9d>
jmp    145d <diffusionmarine+0x35d>
jmp    1462 <diffusionmarine+0x362>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1162 <diffusionmarine+0x62>
pop    %rbp
ret

x86 assembly (clang, -O2):
<diffusionmarine>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r9,(%rsp)
mov    %r8,%r15
mov    %rcx,%r12
mov    %rdx,%rbx
mov    %rsi,%r13
mov    %rdi,%rbp
mov    0x58(%rsp),%eax
mov    0x50(%rsp),%ecx
mov    0x48(%rsp),%r14
test   %eax,%eax
jle    1153 <diffusionmarine+0x43>
mov    %eax,%edx
shl    $0x3,%rdx
mov    %r14,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
mov    0x50(%rsp),%ecx
test   %ecx,%ecx
jle    127e <diffusionmarine+0x16e>
mov    0x40(%rsp),%r8
mov    %ecx,%r9d
xor    %r11d,%r11d
xorpd  %xmm0,%xmm0
jmp    117d <diffusionmarine+0x6d>
nopl   0x0(%rax)
add    $0x1,%r11
cmp    %r9,%r11
je     127e <diffusionmarine+0x16e>
movslq (%r8,%r11,4),%rsi
cmpl   $0x0,0x0(%r13,%rsi,4)
jle    1170 <diffusionmarine+0x60>
mov    %rsi,%rax
shl    $0x5,%rax
lea    (%rax,%rax,4),%rdi
lea    (%r15,%rdi,1),%r10
add    (%rsp),%rdi
lea    (%rsi,%rsi,4),%rcx
shl    $0x4,%rcx
add    %r12,%rcx
xor    %eax,%eax
jmp    11ba <diffusionmarine+0xaa>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rax
cmp    $0x14,%rax
je     1170 <diffusionmarine+0x60>
mov    (%rcx,%rax,4),%edx
test   %edx,%edx
js     1170 <diffusionmarine+0x60>
cmpl   $0x0,0x0(%r13,%rdx,4)
jle    1230 <diffusionmarine+0x120>
cmpl   $0x0,(%rbx,%rsi,4)
jle    1210 <diffusionmarine+0x100>
movsd  0x0(%rbp,%rsi,8),%xmm2
movsd  0x0(%rbp,%rdx,8),%xmm1
ucomisd %xmm1,%xmm2
jbe    1210 <diffusionmarine+0x100>
movsd  (%rdi,%rax,8),%xmm3
ucomisd %xmm0,%xmm3
jbe    1210 <diffusionmarine+0x100>
subsd  %xmm2,%xmm1
mulsd  (%r10,%rax,8),%xmm1
divsd  %xmm3,%xmm1
addsd  (%r14,%rsi,8),%xmm1
movsd  %xmm1,(%r14,%rsi,8)
cs nopw 0x0(%rax,%rax,1)
cmpl   $0x0,(%rbx,%rdx,4)
jle    11b0 <diffusionmarine+0xa0>
movsd  0x0(%rbp,%rsi,8),%xmm2
movsd  0x0(%rbp,%rdx,8),%xmm1
ucomisd %xmm2,%xmm1
jbe    11b0 <diffusionmarine+0xa0>
jmp    1250 <diffusionmarine+0x140>
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rdx,8),%xmm1
movsd  0x0(%rbp,%rsi,8),%xmm2
ucomisd %xmm1,%xmm2
jbe    11b0 <diffusionmarine+0xa0>
cmpl   $0x0,(%rbx,%rsi,4)
jle    11b0 <diffusionmarine+0xa0>
movsd  (%rdi,%rax,8),%xmm3
ucomisd %xmm0,%xmm3
jbe    11b0 <diffusionmarine+0xa0>
subsd  %xmm2,%xmm1
mulsd  (%r10,%rax,8),%xmm1
divsd  %xmm3,%xmm1
addsd  (%r14,%rsi,8),%xmm1
movsd  %xmm1,(%r14,%rsi,8)
jmp    11b0 <diffusionmarine+0xa0>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @diffusionmarine(double* %0, i32* %1, i32* %2, [20 x i32]* %3, [20 x double]* %4, [20 x double]* %5, i32* %6, double* %7, i32 %8, i32 %9) {
%11 = alloca double*
%12 = alloca i32*
%13 = alloca i32*
%14 = alloca [20 x i32]*
%15 = alloca [20 x double]*
%16 = alloca [20 x double]*
%17 = alloca i32*
%18 = alloca double*
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
store double* %0, double** %11
store i32* %1, i32** %12
store i32* %2, i32** %13
store [20 x i32]* %3, [20 x i32]** %14
store [20 x double]* %4, [20 x double]** %15
store [20 x double]* %5, [20 x double]** %16
store i32* %6, i32** %17
store double* %7, double** %18
store i32 %8, i32* %19
store i32 %9, i32* %20
store i32 0, i32* %21
br label %26
26: ; preds = %35, %10
%27 = load i32, i32* %21
%28 = load i32, i32* %20
%29 = icmp slt i32 %27, %28
br i1 %29, label %30, label %38
30: ; preds = %26
%31 = load double*, double** %18
%32 = load i32, i32* %21
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %31, i64 %33
store double 0.000000e+00, double* %34
br label %35
35: ; preds = %30
%36 = load i32, i32* %21
%37 = add nsw i32 %36, 1
store i32 %37, i32* %21
br label %26
38: ; preds = %26
store i32 0, i32* %22
br label %39
39: ; preds = %289, %38
%40 = load i32, i32* %22
%41 = load i32, i32* %19
%42 = icmp slt i32 %40, %41
br i1 %42, label %43, label %292
43: ; preds = %39
%44 = load i32*, i32** %17
%45 = load i32, i32* %22
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds i32, i32* %44, i64 %46
%48 = load i32, i32* %47
store i32 %48, i32* %23
%49 = load i32*, i32** %12
%50 = load i32, i32* %23
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds i32, i32* %49, i64 %51
%53 = load i32, i32* %52
%54 = sitofp i32 %53 to double
%55 = fcmp ogt double %54, 0.000000e+00
br i1 %55, label %56, label %288
56: ; preds = %43
store i32 0, i32* %24
br label %57
57: ; preds = %284, %56
%58 = load i32, i32* %24
%59 = icmp slt i32 %58, 20
br i1 %59, label %60, label %287
60: ; preds = %57
%61 = load [20 x i32]*, [20 x i32]** %14
%62 = load i32, i32* %23
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds [20 x i32], [20 x i32]* %61, i64 %63
%65 = load i32, i32* %24
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds [20 x i32], [20 x i32]* %64, i64 0, i64 %66
%68 = load i32, i32* %67
store i32 %68, i32* %25
%69 = load i32, i32* %25
%70 = icmp slt i32 %69, 0
br i1 %70, label %71, label %72
71: ; preds = %60
br label %287
72: ; preds = %60
%73 = load i32*, i32** %12
%74 = load i32, i32* %25
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds i32, i32* %73, i64 %75
%77 = load i32, i32* %76
%78 = icmp sgt i32 %77, 0
br i1 %78, label %79, label %210
79: ; preds = %72
%80 = load i32*, i32** %13
%81 = load i32, i32* %23
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds i32, i32* %80, i64 %82
%84 = load i32, i32* %83
%85 = icmp sgt i32 %84, 0
br i1 %85, label %86, label %144
86: ; preds = %79
%87 = load double*, double** %11
%88 = load i32, i32* %23
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %87, i64 %89
%91 = load double, double* %90
%92 = load double*, double** %11
%93 = load i32, i32* %25
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds double, double* %92, i64 %94
%96 = load double, double* %95
%97 = fcmp ogt double %91, %96
br i1 %97, label %98, label %144
98: ; preds = %86
%99 = load [20 x double]*, [20 x double]** %16
%100 = load i32, i32* %23
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds [20 x double], [20 x double]* %99, i64 %101
%103 = load i32, i32* %24
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds [20 x double], [20 x double]* %102, i64 0, i64 %104
%106 = load double, double* %105
%107 = fcmp ogt double %106, 0.000000e+00
br i1 %107, label %108, label %144
108: ; preds = %98
%109 = load [20 x double]*, [20 x double]** %15
%110 = load i32, i32* %23
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds [20 x double], [20 x double]* %109, i64 %111
%113 = load i32, i32* %24
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds [20 x double], [20 x double]* %112, i64 0, i64 %114
%116 = load double, double* %115
%117 = load double*, double** %11
%118 = load i32, i32* %25
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds double, double* %117, i64 %119
%121 = load double, double* %120
%122 = load double*, double** %11
%123 = load i32, i32* %23
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds double, double* %122, i64 %124
%126 = load double, double* %125
%127 = fsub double %121, %126
%128 = fmul double %116, %127
%129 = load [20 x double]*, [20 x double]** %16
%130 = load i32, i32* %23
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds [20 x double], [20 x double]* %129, i64 %131
%133 = load i32, i32* %24
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds [20 x double], [20 x double]* %132, i64 0, i64 %134
%136 = load double, double* %135
%137 = fdiv double %128, %136
%138 = load double*, double** %18
%139 = load i32, i32* %23
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds double, double* %138, i64 %140
%142 = load double, double* %141
%143 = fadd double %142, %137
store double %143, double* %141
br label %144
144: ; preds = %108, %98, %86, %79
%145 = load i32*, i32** %13
%146 = load i32, i32* %25
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds i32, i32* %145, i64 %147
%149 = load i32, i32* %148
%150 = icmp sgt i32 %149, 0
br i1 %150, label %151, label %209
151: ; preds = %144
%152 = load double*, double** %11
%153 = load i32, i32* %23
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds double, double* %152, i64 %154
%156 = load double, double* %155
%157 = load double*, double** %11
%158 = load i32, i32* %25
%159 = sext i32 %158 to i64
%160 = getelementptr inbounds double, double* %157, i64 %159
%161 = load double, double* %160
%162 = fcmp olt double %156, %161
br i1 %162, label %163, label %209
163: ; preds = %151
%164 = load [20 x double]*, [20 x double]** %16
%165 = load i32, i32* %23
%166 = sext i32 %165 to i64
%167 = getelementptr inbounds [20 x double], [20 x double]* %164, i64 %166
%168 = load i32, i32* %24
%169 = sext i32 %168 to i64
%170 = getelementptr inbounds [20 x double], [20 x double]* %167, i64 0, i64 %169
%171 = load double, double* %170
%172 = fcmp ogt double %171, 0.000000e+00
br i1 %172, label %173, label %209
173: ; preds = %163
%174 = load [20 x double]*, [20 x double]** %15
%175 = load i32, i32* %23
%176 = sext i32 %175 to i64
%177 = getelementptr inbounds [20 x double], [20 x double]* %174, i64 %176
%178 = load i32, i32* %24
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds [20 x double], [20 x double]* %177, i64 0, i64 %179
%181 = load double, double* %180
%182 = load double*, double** %11
%183 = load i32, i32* %25
%184 = sext i32 %183 to i64
%185 = getelementptr inbounds double, double* %182, i64 %184
%186 = load double, double* %185
%187 = load double*, double** %11
%188 = load i32, i32* %23
%189 = sext i32 %188 to i64
%190 = getelementptr inbounds double, double* %187, i64 %189
%191 = load double, double* %190
%192 = fsub double %186, %191
%193 = fmul double %181, %192
%194 = load [20 x double]*, [20 x double]** %16
%195 = load i32, i32* %23
%196 = sext i32 %195 to i64
%197 = getelementptr inbounds [20 x double], [20 x double]* %194, i64 %196
%198 = load i32, i32* %24
%199 = sext i32 %198 to i64
%200 = getelementptr inbounds [20 x double], [20 x double]* %197, i64 0, i64 %199
%201 = load double, double* %200
%202 = fdiv double %193, %201
%203 = load double*, double** %18
%204 = load i32, i32* %23
%205 = sext i32 %204 to i64
%206 = getelementptr inbounds double, double* %203, i64 %205
%207 = load double, double* %206
%208 = fadd double %207, %202
store double %208, double* %206
br label %209
209: ; preds = %173, %163, %151, %144
br label %210
210: ; preds = %209, %72
%211 = load i32*, i32** %12
%212 = load i32, i32* %25
%213 = sext i32 %212 to i64
%214 = getelementptr inbounds i32, i32* %211, i64 %213
%215 = load i32, i32* %214
%216 = icmp slt i32 %215, 1
br i1 %216, label %217, label %283
217: ; preds = %210
%218 = load double*, double** %11
%219 = load i32, i32* %25
%220 = sext i32 %219 to i64
%221 = getelementptr inbounds double, double* %218, i64 %220
%222 = load double, double* %221
%223 = load double*, double** %11
%224 = load i32, i32* %23
%225 = sext i32 %224 to i64
%226 = getelementptr inbounds double, double* %223, i64 %225
%227 = load double, double* %226
%228 = fcmp olt double %222, %227
br i1 %228, label %229, label %282
229: ; preds = %217
%230 = load i32*, i32** %13
%231 = load i32, i32* %23
%232 = sext i32 %231 to i64
%233 = getelementptr inbounds i32, i32* %230, i64 %232
%234 = load i32, i32* %233
%235 = icmp sgt i32 %234, 0
br i1 %235, label %236, label %282
236: ; preds = %229
%237 = load [20 x double]*, [20 x double]** %16
%238 = load i32, i32* %23
%239 = sext i32 %238 to i64
%240 = getelementptr inbounds [20 x double], [20 x double]* %237, i64 %239
%241 = load i32, i32* %24
%242 = sext i32 %241 to i64
%243 = getelementptr inbounds [20 x double], [20 x double]* %240, i64 0, i64 %242
%244 = load double, double* %243
%245 = fcmp ogt double %244, 0.000000e+00
br i1 %245, label %246, label %282
246: ; preds = %236
%247 = load [20 x double]*, [20 x double]** %15
%248 = load i32, i32* %23
%249 = sext i32 %248 to i64
%250 = getelementptr inbounds [20 x double], [20 x double]* %247, i64 %249
%251 = load i32, i32* %24
%252 = sext i32 %251 to i64
%253 = getelementptr inbounds [20 x double], [20 x double]* %250, i64 0, i64 %252
%254 = load double, double* %253
%255 = load double*, double** %11
%256 = load i32, i32* %25
%257 = sext i32 %256 to i64
%258 = getelementptr inbounds double, double* %255, i64 %257
%259 = load double, double* %258
%260 = load double*, double** %11
%261 = load i32, i32* %23
%262 = sext i32 %261 to i64
%263 = getelementptr inbounds double, double* %260, i64 %262
%264 = load double, double* %263
%265 = fsub double %259, %264
%266 = fmul double %254, %265
%267 = load [20 x double]*, [20 x double]** %16
%268 = load i32, i32* %23
%269 = sext i32 %268 to i64
%270 = getelementptr inbounds [20 x double], [20 x double]* %267, i64 %269
%271 = load i32, i32* %24
%272 = sext i32 %271 to i64
%273 = getelementptr inbounds [20 x double], [20 x double]* %270, i64 0, i64 %272
%274 = load double, double* %273
%275 = fdiv double %266, %274
%276 = load double*, double** %18
%277 = load i32, i32* %23
%278 = sext i32 %277 to i64
%279 = getelementptr inbounds double, double* %276, i64 %278
%280 = load double, double* %279
%281 = fadd double %280, %275
store double %281, double* %279
br label %282
282: ; preds = %246, %236, %229, %217
br label %283
283: ; preds = %282, %210
br label %284
284: ; preds = %283
%285 = load i32, i32* %24
%286 = add nsw i32 %285, 1
store i32 %286, i32* %24
br label %57
287: ; preds = %71, %57
br label %288
288: ; preds = %287, %43
br label %289
289: ; preds = %288
%290 = load i32, i32* %22
%291 = add nsw i32 %290, 1
store i32 %291, i32* %22
br label %39
292: ; preds = %39
ret void
}

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