void ScatterUpdateFO(long dimi, long *structi, long dimj, long *structj,
                     long destdim, double *oldupdate, double *newupdate) {
  long i, j;
  double *srccol, *destcol;
  for (j = 0; j < dimj; j++) {
    if (structj)
      destcol = &newupdate[structj[j] * destdim];
    else
      destcol = &newupdate[j * destdim];
    srccol = &oldupdate[j * dimi];
    if (structi)
      for (i = 0; i < dimi; i++) {
        destcol[structi[i]] += srccol[i];
        srccol[i] = 0.0;
      }
    else
      for (i = 0; i < dimi; i++) {
        destcol[i] += srccol[i];
        srccol[i] = 0.0;
      }
  }
}


x86 assembly (clang, -O0):
<ScatterUpdateFO>:
push   %rbp
mov    %rsp,%rbp
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)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    1275 <ScatterUpdateFO+0x175>
cmpq   $0x0,-0x20(%rbp)
je     1166 <ScatterUpdateFO+0x66>
mov    0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x40(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
imul   -0x28(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
jmp    117e <ScatterUpdateFO+0x7e>
mov    0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
imul   -0x28(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x40(%rbp),%rcx
imul   -0x8(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
cmpq   $0x0,-0x10(%rbp)
je     1204 <ScatterUpdateFO+0x104>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
jge    11ff <ScatterUpdateFO+0xff>
mov    -0x48(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x38(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
mov    -0x38(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11a9 <ScatterUpdateFO+0xa9>
jmp    125f <ScatterUpdateFO+0x15f>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
jge    125a <ScatterUpdateFO+0x15a>
mov    -0x48(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x38(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
mov    -0x38(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    120c <ScatterUpdateFO+0x10c>
jmp    125f <ScatterUpdateFO+0x15f>
jmp    1264 <ScatterUpdateFO+0x164>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1128 <ScatterUpdateFO+0x28>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ScatterUpdateFO>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8,-0x8(%rsp)
test   %rdx,%rdx
jle    146c <ScatterUpdateFO+0x36c>
mov    %r9,%r14
mov    %rdx,%r8
mov    0x38(%rsp),%rax
lea    (%rax,%rdi,8),%rdx
mov    %rdx,-0x28(%rsp)
lea    -0x4(%rdi),%rdx
mov    %rdx,-0x30(%rsp)
mov    %rdx,%rbx
shr    $0x2,%rbx
add    $0x1,%rbx
mov    %rdi,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %rdi,%rbp
and    $0xfffffffffffffffc,%rbp
mov    %rbp,-0x48(%rsp)
mov    %rbx,%rbp
mov    %rbx,-0x38(%rsp)
and    $0xfffffffffffffffe,%rbx
mov    %rbx,-0x50(%rsp)
mov    %rdi,%rbp
neg    %rbp
mov    %rbp,-0x18(%rsp)
lea    0x8(%r9),%rbp
lea    0x0(,%rdi,8),%rbx
mov    %rbx,-0x10(%rsp)
lea    0x30(%r9),%r11
lea    0x30(%rax),%rbx
mov    %rbx,-0x58(%rsp)
add    $0x8,%rax
mov    %rax,-0x20(%rsp)
xor    %r12d,%r12d
xorps  %xmm0,%xmm0
mov    %r9,%r10
mov    %r8,-0x40(%rsp)
jmp    11cb <ScatterUpdateFO+0xcb>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r12
mov    -0x10(%rsp),%rax
add    %rax,%rbp
add    %rax,%r11
add    %rax,%r10
cmp    %r8,%r12
je     146c <ScatterUpdateFO+0x36c>
mov    %r12,%rbx
imul   %rdi,%rbx
mov    %r12,%r9
test   %rcx,%rcx
je     11de <ScatterUpdateFO+0xde>
mov    (%rcx,%r12,8),%r9
lea    (%r14,%rbx,8),%rax
mov    %rax,-0x60(%rsp)
imul   -0x8(%rsp),%r9
mov    0x38(%rsp),%rax
lea    (%rax,%r9,8),%r13
test   %rsi,%rsi
je     1290 <ScatterUpdateFO+0x190>
test   %rdi,%rdi
jle    11b0 <ScatterUpdateFO+0xb0>
xor    %eax,%eax
cmp    $0x1,%rdi
je     125c <ScatterUpdateFO+0x15c>
nopl   0x0(%rax)
movsd  -0x8(%rbp,%rax,8),%xmm1
mov    (%rsi,%rax,8),%rbx
addsd  0x0(%r13,%rbx,8),%xmm1
movsd  %xmm1,0x0(%r13,%rbx,8)
movq   $0x0,-0x8(%rbp,%rax,8)
movsd  0x0(%rbp,%rax,8),%xmm1
mov    0x8(%rsi,%rax,8),%rbx
addsd  0x0(%r13,%rbx,8),%xmm1
movsd  %xmm1,0x0(%r13,%rbx,8)
movq   $0x0,0x0(%rbp,%rax,8)
add    $0x2,%rax
cmp    %rax,%rdx
jne    1210 <ScatterUpdateFO+0x110>
test   $0x1,%dil
je     11b0 <ScatterUpdateFO+0xb0>
mov    -0x60(%rsp),%r9
movsd  (%r9,%rax,8),%xmm1
mov    (%rsi,%rax,8),%rbx
addsd  0x0(%r13,%rbx,8),%xmm1
movsd  %xmm1,0x0(%r13,%rbx,8)
movq   $0x0,(%r9,%rax,8)
jmp    11b0 <ScatterUpdateFO+0xb0>
test   %rdi,%rdi
jle    11b0 <ScatterUpdateFO+0xb0>
cmp    $0x4,%rdi
jb     12c3 <ScatterUpdateFO+0x1c3>
add    %rdi,%rbx
lea    (%r14,%rbx,8),%rbx
cmp    %rbx,%r13
jae    1352 <ScatterUpdateFO+0x252>
mov    -0x28(%rsp),%rax
lea    (%rax,%r9,8),%rbx
cmp    %rbx,-0x60(%rsp)
jae    1352 <ScatterUpdateFO+0x252>
xor    %r15d,%r15d
mov    %r15,%rbx
test   $0x1,%dil
je     12f7 <ScatterUpdateFO+0x1f7>
mov    -0x60(%rsp),%rax
movsd  (%rax,%r15,8),%xmm1
addsd  0x0(%r13,%r15,8),%xmm1
movsd  %xmm1,0x0(%r13,%r15,8)
movq   $0x0,(%rax,%r15,8)
mov    %r15,%rbx
or     $0x1,%rbx
not    %r15
cmp    -0x18(%rsp),%r15
je     11b0 <ScatterUpdateFO+0xb0>
mov    -0x20(%rsp),%rax
lea    (%rax,%r9,8),%rax
xchg   %ax,%ax
movsd  (%r10,%rbx,8),%xmm1
addsd  -0x8(%rax,%rbx,8),%xmm1
movsd  %xmm1,-0x8(%rax,%rbx,8)
movq   $0x0,(%r10,%rbx,8)
movsd  0x8(%r10,%rbx,8),%xmm1
addsd  (%rax,%rbx,8),%xmm1
movsd  %xmm1,(%rax,%rbx,8)
movq   $0x0,0x8(%r10,%rbx,8)
add    $0x2,%rbx
cmp    %rbx,%rdi
jne    1310 <ScatterUpdateFO+0x210>
jmp    11b0 <ScatterUpdateFO+0xb0>
mov    %r14,%rax
cmpq   $0x4,-0x30(%rsp)
jae    1365 <ScatterUpdateFO+0x265>
xor    %r15d,%r15d
jmp    1406 <ScatterUpdateFO+0x306>
mov    -0x58(%rsp),%rbx
lea    (%rbx,%r9,8),%r14
mov    -0x50(%rsp),%r8
xor    %r15d,%r15d
cs nopw 0x0(%rax,%rax,1)
movupd -0x30(%r11,%r15,8),%xmm1
movupd -0x20(%r11,%r15,8),%xmm2
movupd -0x30(%r14,%r15,8),%xmm3
addpd  %xmm1,%xmm3
movupd -0x20(%r14,%r15,8),%xmm1
addpd  %xmm2,%xmm1
movupd %xmm3,-0x30(%r14,%r15,8)
movupd %xmm1,-0x20(%r14,%r15,8)
movups %xmm0,-0x30(%r11,%r15,8)
movups %xmm0,-0x20(%r11,%r15,8)
movupd -0x10(%r11,%r15,8),%xmm1
movupd (%r11,%r15,8),%xmm2
movupd -0x10(%r14,%r15,8),%xmm3
addpd  %xmm1,%xmm3
movupd (%r14,%r15,8),%xmm1
addpd  %xmm2,%xmm1
movupd %xmm3,-0x10(%r14,%r15,8)
movupd %xmm1,(%r14,%r15,8)
movups %xmm0,-0x10(%r11,%r15,8)
movups %xmm0,(%r11,%r15,8)
add    $0x8,%r15
add    $0xfffffffffffffffe,%r8
jne    1380 <ScatterUpdateFO+0x280>
testb  $0x1,-0x38(%rsp)
je     144e <ScatterUpdateFO+0x34e>
mov    -0x60(%rsp),%rbx
movupd (%rbx,%r15,8),%xmm1
movupd 0x10(%rbx,%r15,8),%xmm2
movupd 0x0(%r13,%r15,8),%xmm3
addpd  %xmm1,%xmm3
movupd 0x10(%r13,%r15,8),%xmm1
addpd  %xmm2,%xmm1
movupd %xmm3,0x0(%r13,%r15,8)
movupd %xmm1,0x10(%r13,%r15,8)
movups %xmm0,(%rbx,%r15,8)
movups %xmm0,0x10(%rbx,%r15,8)
mov    -0x48(%rsp),%rbx
mov    %rbx,%r15
cmp    %rdi,%rbx
mov    %rax,%r14
mov    -0x40(%rsp),%r8
je     11b0 <ScatterUpdateFO+0xb0>
jmp    12c6 <ScatterUpdateFO+0x1c6>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @ScatterUpdateFO(i64 %0, i64* %1, i64 %2, i64* %3, i64 %4, double* nocapture %5, double* nocapture %6) {
%8 = icmp sgt i64 %2, 0
br i1 %8, label %9, label %183
9: ; preds = %7
%10 = icmp eq i64* %3, null
%11 = icmp eq i64* %1, null
%12 = icmp sgt i64 %0, 0
%13 = getelementptr double, double* %6, i64 %0
%14 = add i64 %0, -4
%15 = lshr i64 %14, 2
%16 = add nuw nsw i64 %15, 1
%17 = and i64 %0, 1
%18 = icmp eq i64 %0, 1
%19 = and i64 %0, -2
%20 = icmp eq i64 %17, 0
%21 = icmp ult i64 %0, 4
%22 = and i64 %0, -4
%23 = and i64 %16, 1
%24 = icmp ult i64 %14, 4
%25 = and i64 %16, 9223372036854775806
%26 = icmp eq i64 %23, 0
%27 = icmp eq i64 %22, %0
%28 = and i64 %0, 1
%29 = icmp eq i64 %28, 0
%30 = sub i64 0, %0
br label %31
31: ; preds = %9, %180
%32 = phi i64 [ 0, %9 ], [ %181, %180 ]
%33 = mul i64 %32, %0
%34 = getelementptr double, double* %5, i64 %33
%35 = add i64 %33, %0
%36 = getelementptr double, double* %5, i64 %35
br i1 %10, label %40, label %37
37: ; preds = %31
%38 = getelementptr inbounds i64, i64* %3, i64 %32
%39 = load i64, i64* %38
br label %40
40: ; preds = %31, %37
%41 = phi i64 [ %39, %37 ], [ %32, %31 ]
%42 = mul i64 %41, %4
%43 = getelementptr double, double* %6, i64 %42
%44 = mul nsw i64 %32, %0
%45 = getelementptr inbounds double, double* %5, i64 %44
br i1 %11, label %48, label %46
46: ; preds = %40
br i1 %12, label %47, label %180
47: ; preds = %46
br i1 %18, label %170, label %134
48: ; preds = %40
br i1 %12, label %49, label %180
49: ; preds = %48
br i1 %21, label %121, label %50
50: ; preds = %49
%51 = getelementptr double, double* %13, i64 %42
%52 = icmp ult double* %43, %36
%53 = icmp ult double* %34, %51
%54 = and i1 %52, %53
br i1 %54, label %121, label %55
55: ; preds = %50
br i1 %24, label %99, label %56
56: ; preds = %55, %56
%57 = phi i64 [ %96, %56 ], [ 0, %55 ]
%58 = phi i64 [ %97, %56 ], [ 0, %55 ]
%59 = getelementptr inbounds double, double* %45, i64 %57
%60 = bitcast double* %59 to <2 x double>*
%61 = load <2 x double>, <2 x double>* %60
%62 = getelementptr inbounds double, double* %59, i64 2
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = getelementptr inbounds double, double* %43, i64 %57
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = getelementptr inbounds double, double* %65, i64 2
%69 = bitcast double* %68 to <2 x double>*
%70 = load <2 x double>, <2 x double>* %69
%71 = fadd <2 x double> %61, %67
%72 = fadd <2 x double> %64, %70
%73 = bitcast double* %65 to <2 x double>*
store <2 x double> %71, <2 x double>* %73
%74 = bitcast double* %68 to <2 x double>*
store <2 x double> %72, <2 x double>* %74
%75 = bitcast double* %59 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %75
%76 = bitcast double* %62 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %76
%77 = or i64 %57, 4
%78 = getelementptr inbounds double, double* %45, i64 %77
%79 = bitcast double* %78 to <2 x double>*
%80 = load <2 x double>, <2 x double>* %79
%81 = getelementptr inbounds double, double* %78, i64 2
%82 = bitcast double* %81 to <2 x double>*
%83 = load <2 x double>, <2 x double>* %82
%84 = getelementptr inbounds double, double* %43, i64 %77
%85 = bitcast double* %84 to <2 x double>*
%86 = load <2 x double>, <2 x double>* %85
%87 = getelementptr inbounds double, double* %84, i64 2
%88 = bitcast double* %87 to <2 x double>*
%89 = load <2 x double>, <2 x double>* %88
%90 = fadd <2 x double> %80, %86
%91 = fadd <2 x double> %83, %89
%92 = bitcast double* %84 to <2 x double>*
store <2 x double> %90, <2 x double>* %92
%93 = bitcast double* %87 to <2 x double>*
store <2 x double> %91, <2 x double>* %93
%94 = bitcast double* %78 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %94
%95 = bitcast double* %81 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %95
%96 = add nuw i64 %57, 8
%97 = add i64 %58, 2
%98 = icmp eq i64 %97, %25
br i1 %98, label %99, label %56
99: ; preds = %56, %55
%100 = phi i64 [ 0, %55 ], [ %96, %56 ]
br i1 %26, label %120, label %101
101: ; preds = %99
%102 = getelementptr inbounds double, double* %45, i64 %100
%103 = bitcast double* %102 to <2 x double>*
%104 = load <2 x double>, <2 x double>* %103
%105 = getelementptr inbounds double, double* %102, i64 2
%106 = bitcast double* %105 to <2 x double>*
%107 = load <2 x double>, <2 x double>* %106
%108 = getelementptr inbounds double, double* %43, i64 %100
%109 = bitcast double* %108 to <2 x double>*
%110 = load <2 x double>, <2 x double>* %109
%111 = getelementptr inbounds double, double* %108, i64 2
%112 = bitcast double* %111 to <2 x double>*
%113 = load <2 x double>, <2 x double>* %112
%114 = fadd <2 x double> %104, %110
%115 = fadd <2 x double> %107, %113
%116 = bitcast double* %108 to <2 x double>*
store <2 x double> %114, <2 x double>* %116
%117 = bitcast double* %111 to <2 x double>*
store <2 x double> %115, <2 x double>* %117
%118 = bitcast double* %102 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %118
%119 = bitcast double* %105 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %119
br label %120
120: ; preds = %99, %101
br i1 %27, label %180, label %121
121: ; preds = %50, %49, %120
%122 = phi i64 [ 0, %50 ], [ 0, %49 ], [ %22, %120 ]
%123 = xor i64 %122, -1
br i1 %29, label %131, label %124
124: ; preds = %121
%125 = getelementptr inbounds double, double* %45, i64 %122
%126 = load double, double* %125
%127 = getelementptr inbounds double, double* %43, i64 %122
%128 = load double, double* %127
%129 = fadd double %126, %128
store double %129, double* %127
store double 0.000000e+00, double* %125
%130 = or i64 %122, 1
br label %131
131: ; preds = %124, %121
%132 = phi i64 [ %122, %121 ], [ %130, %124 ]
%133 = icmp eq i64 %123, %30
br i1 %133, label %180, label %155
134: ; preds = %47, %134
%135 = phi i64 [ %152, %134 ], [ 0, %47 ]
%136 = phi i64 [ %153, %134 ], [ 0, %47 ]
%137 = getelementptr inbounds double, double* %45, i64 %135
%138 = load double, double* %137
%139 = getelementptr inbounds i64, i64* %1, i64 %135
%140 = load i64, i64* %139
%141 = getelementptr inbounds double, double* %43, i64 %140
%142 = load double, double* %141
%143 = fadd double %138, %142
store double %143, double* %141
store double 0.000000e+00, double* %137
%144 = or i64 %135, 1
%145 = getelementptr inbounds double, double* %45, i64 %144
%146 = load double, double* %145
%147 = getelementptr inbounds i64, i64* %1, i64 %144
%148 = load i64, i64* %147
%149 = getelementptr inbounds double, double* %43, i64 %148
%150 = load double, double* %149
%151 = fadd double %146, %150
store double %151, double* %149
store double 0.000000e+00, double* %145
%152 = add nuw nsw i64 %135, 2
%153 = add i64 %136, 2
%154 = icmp eq i64 %153, %19
br i1 %154, label %170, label %134
155: ; preds = %131, %155
%156 = phi i64 [ %168, %155 ], [ %132, %131 ]
%157 = getelementptr inbounds double, double* %45, i64 %156
%158 = load double, double* %157
%159 = getelementptr inbounds double, double* %43, i64 %156
%160 = load double, double* %159
%161 = fadd double %158, %160
store double %161, double* %159
store double 0.000000e+00, double* %157
%162 = add nuw nsw i64 %156, 1
%163 = getelementptr inbounds double, double* %45, i64 %162
%164 = load double, double* %163
%165 = getelementptr inbounds double, double* %43, i64 %162
%166 = load double, double* %165
%167 = fadd double %164, %166
store double %167, double* %165
store double 0.000000e+00, double* %163
%168 = add nuw nsw i64 %156, 2
%169 = icmp eq i64 %168, %0
br i1 %169, label %180, label %155
170: ; preds = %134, %47
%171 = phi i64 [ 0, %47 ], [ %152, %134 ]
br i1 %20, label %180, label %172
172: ; preds = %170
%173 = getelementptr inbounds double, double* %45, i64 %171
%174 = load double, double* %173
%175 = getelementptr inbounds i64, i64* %1, i64 %171
%176 = load i64, i64* %175
%177 = getelementptr inbounds double, double* %43, i64 %176
%178 = load double, double* %177
%179 = fadd double %174, %178
store double %179, double* %177
store double 0.000000e+00, double* %173
br label %180
180: ; preds = %172, %170, %131, %155, %120, %46, %48
%181 = add nuw nsw i64 %32, 1
%182 = icmp eq i64 %181, %2
br i1 %182, label %183, label %31
183: ; preds = %180, %7
ret void
}