void addMul2Tri(double base[], const double b[], const double a[],
                const int n_ra, const int n_ca, const double c[],
                const int n_rc, const int n_cc, double bc[]) {
  int ci, ri, si;
  double d;
  for (ri = 0; ri < n_ra; ri++) {
    for (ci = 0; ci < n_cc; ci++) {
      bc[((ri) + (ci) * (n_ra))] = 0.0;
      for (si = 0; si < n_rc; si++) {
        bc[((ri) + (ci) * (n_ra))] +=
            b[((ri > si) ? (((ri) * ((ri) + 1)) / 2 + (si))
                         : (((si) * ((si) + 1)) / 2 + (ri)))] *
            c[si + ci * n_rc];
      }
    }
  }
  for (ri = 0; ri < n_ca; ri++) {
    for (ci = 0; ci < n_cc; ci++) {
      d = 0.0;
      for (si = 0; si < n_ra; si++) {
        d += a[((si) + (ri) * (n_ra))] * bc[((si) + (ci) * (n_ra))];
      }
      base[((ri) + (ci) * (n_ca))] += d;
    }
  }
}


x86 assembly (clang, -O0):
<addMul2Tri>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1246 <addMul2Tri+0x146>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    1233 <addMul2Tri+0x133>
mov    0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    1220 <addMul2Tri+0x120>
mov    -0x10(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x34(%rbp),%eax
jle    11ac <addMul2Tri+0xac>
mov    -0x30(%rbp),%eax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x34(%rbp),%eax
mov    %eax,-0x4c(%rbp)
jmp    11c6 <addMul2Tri+0xc6>
mov    -0x34(%rbp),%eax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x30(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x48(%rbp),%rax
mov    -0x4c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    116d <addMul2Tri+0x6d>
jmp    1225 <addMul2Tri+0x125>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    113f <addMul2Tri+0x3f>
jmp    1238 <addMul2Tri+0x138>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    112c <addMul2Tri+0x2c>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    131a <addMul2Tri+0x21a>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    1307 <addMul2Tri+0x207>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12d7 <addMul2Tri+0x1d7>
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    127b <addMul2Tri+0x17b>
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1260 <addMul2Tri+0x160>
jmp    130c <addMul2Tri+0x20c>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    124d <addMul2Tri+0x14d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<addMul2Tri>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r9,-0x20(%rsp)
mov    %r8d,-0x14(%rsp)
mov    %rdx,-0x10(%rsp)
mov    %rdi,-0x8(%rsp)
mov    0x48(%rsp),%r14
mov    %ecx,%eax
mov    %ecx,-0x18(%rsp)
test   %ecx,%ecx
mov    %rax,%rcx
jle    11fb <addMul2Tri+0xfb>
mov    0x38(%rsp),%eax
mov    %eax,%edx
movslq %eax,%r11
mov    %r11d,%ebp
shl    $0x3,%r11
mov    $0x2,%edi
xor    %eax,%eax
mov    %rax,-0x28(%rsp)
xor    %r10d,%r10d
jmp    1176 <addMul2Tri+0x76>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r10
add    %rdi,-0x28(%rsp)
add    $0x2,%rdi
cmp    %rcx,%r10
je     11fb <addMul2Tri+0xfb>
cmpl   $0x0,0x40(%rsp)
jle    1160 <addMul2Tri+0x60>
mov    -0x28(%rsp),%r13
shr    %r13d
mov    -0x20(%rsp),%r9
xor    %eax,%eax
jmp    11a0 <addMul2Tri+0xa0>
xchg   %ax,%ax
add    $0x1,%rax
add    %r11,%r9
mov    0x40(%rsp),%ebx
cmp    %rbx,%rax
je     1160 <addMul2Tri+0x60>
mov    %rax,%r8
imul   %rcx,%r8
add    %r10,%r8
movq   $0x0,(%r14,%r8,8)
test   %edx,%edx
jle    1190 <addMul2Tri+0x90>
xorpd  %xmm0,%xmm0
xor    %r12d,%r12d
jmp    11f0 <addMul2Tri+0xf0>
nop
lea    0x1(%r12),%r15d
imul   %r12d,%r15d
shr    %r15d
add    %r10d,%r15d
movslq %r15d,%rbx
movsd  (%rsi,%rbx,8),%xmm1
mulsd  (%r9,%r12,8),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r14,%r8,8)
add    $0x1,%r12
cmp    %r12,%rbp
je     1190 <addMul2Tri+0x90>
cmp    %r12,%r10
jbe    11c0 <addMul2Tri+0xc0>
lea    (%r12,%r13,1),%r15d
jmp    11cf <addMul2Tri+0xcf>
mov    %rcx,%r11
mov    -0x14(%rsp),%eax
test   %eax,%eax
mov    -0x18(%rsp),%edi
mov    -0x8(%rsp),%rcx
jle    132c <addMul2Tri+0x22c>
movslq %edi,%rdx
mov    %eax,%r8d
mov    %r11d,%esi
and    $0xfffffffe,%esi
mov    -0x10(%rsp),%rax
lea    0x8(%rax),%rbx
mov    %rdx,-0x28(%rsp)
lea    0x0(,%rdx,8),%r15
lea    0x8(%r14),%rax
mov    %rax,-0x20(%rsp)
xor    %r13d,%r13d
mov    -0x10(%rsp),%r9
jmp    1260 <addMul2Tri+0x160>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r13
add    %r15,%rbx
cmp    %r8,%r13
je     132c <addMul2Tri+0x22c>
cmpl   $0x0,0x40(%rsp)
jle    1250 <addMul2Tri+0x150>
mov    %r13,%r10
imul   -0x28(%rsp),%r10
mov    -0x20(%rsp),%rbp
xor    %r12d,%r12d
jmp    12a8 <addMul2Tri+0x1a8>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
mov    %r12,%rax
imul   %r8,%rax
add    %r13,%rax
addsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%rax,8)
add    $0x1,%r12
add    %r15,%rbp
mov    0x40(%rsp),%eax
cmp    %rax,%r12
je     1250 <addMul2Tri+0x150>
test   %edi,%edi
jle    1280 <addMul2Tri+0x180>
cmp    $0x1,%edi
jne    12c0 <addMul2Tri+0x1c0>
xorpd  %xmm0,%xmm0
xor    %edx,%edx
jmp    12fc <addMul2Tri+0x1fc>
nopl   0x0(%rax)
xorpd  %xmm1,%xmm1
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movsd  -0x8(%rbx,%rdx,8),%xmm2
movsd  (%rbx,%rdx,8),%xmm0
mulsd  -0x8(%rbp,%rdx,8),%xmm2
mulsd  0x0(%rbp,%rdx,8),%xmm0
addsd  %xmm1,%xmm2
addsd  %xmm2,%xmm0
add    $0x2,%rdx
movapd %xmm0,%xmm1
cmp    %rdx,%rsi
jne    12d0 <addMul2Tri+0x1d0>
test   $0x1,%r11b
je     1284 <addMul2Tri+0x184>
mov    %r12,%rcx
imul   -0x28(%rsp),%rcx
lea    (%rdx,%r10,1),%rax
movsd  (%r9,%rax,8),%xmm1
add    %rdx,%rcx
mulsd  (%r14,%rcx,8),%xmm1
mov    -0x8(%rsp),%rcx
addsd  %xmm1,%xmm0
jmp    1284 <addMul2Tri+0x184>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @addMul2Tri(double* %0, double* %1, double* %2, i32 %3, i32 %4, double* %5, i32 %6, i32 %7, double* %8) {
%10 = alloca double*
%11 = alloca double*
%12 = alloca double*
%13 = alloca i32
%14 = alloca i32
%15 = alloca double*
%16 = alloca i32
%17 = alloca i32
%18 = alloca double*
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca double
store double* %0, double** %10
store double* %1, double** %11
store double* %2, double** %12
store i32 %3, i32* %13
store i32 %4, i32* %14
store double* %5, double** %15
store i32 %6, i32* %16
store i32 %7, i32* %17
store double* %8, double** %18
store i32 0, i32* %20
br label %23
23: ; preds = %98, %9
%24 = load i32, i32* %20
%25 = load i32, i32* %13
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %101
27: ; preds = %23
store i32 0, i32* %19
br label %28
28: ; preds = %94, %27
%29 = load i32, i32* %19
%30 = load i32, i32* %17
%31 = icmp slt i32 %29, %30
br i1 %31, label %32, label %97
32: ; preds = %28
%33 = load double*, double** %18
%34 = load i32, i32* %20
%35 = load i32, i32* %19
%36 = load i32, i32* %13
%37 = mul nsw i32 %35, %36
%38 = add nsw i32 %34, %37
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %33, i64 %39
store double 0.000000e+00, double* %40
store i32 0, i32* %21
br label %41
41: ; preds = %90, %32
%42 = load i32, i32* %21
%43 = load i32, i32* %16
%44 = icmp slt i32 %42, %43
br i1 %44, label %45, label %93
45: ; preds = %41
%46 = load double*, double** %11
%47 = load i32, i32* %20
%48 = load i32, i32* %21
%49 = icmp sgt i32 %47, %48
br i1 %49, label %50, label %58
50: ; preds = %45
%51 = load i32, i32* %20
%52 = load i32, i32* %20
%53 = add nsw i32 %52, 1
%54 = mul nsw i32 %51, %53
%55 = sdiv i32 %54, 2
%56 = load i32, i32* %21
%57 = add nsw i32 %55, %56
br label %66
58: ; preds = %45
%59 = load i32, i32* %21
%60 = load i32, i32* %21
%61 = add nsw i32 %60, 1
%62 = mul nsw i32 %59, %61
%63 = sdiv i32 %62, 2
%64 = load i32, i32* %20
%65 = add nsw i32 %63, %64
br label %66
66: ; preds = %58, %50
%67 = phi i32 [ %57, %50 ], [ %65, %58 ]
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %46, i64 %68
%70 = load double, double* %69
%71 = load double*, double** %15
%72 = load i32, i32* %21
%73 = load i32, i32* %19
%74 = load i32, i32* %16
%75 = mul nsw i32 %73, %74
%76 = add nsw i32 %72, %75
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %71, i64 %77
%79 = load double, double* %78
%80 = load double*, double** %18
%81 = load i32, i32* %20
%82 = load i32, i32* %19
%83 = load i32, i32* %13
%84 = mul nsw i32 %82, %83
%85 = add nsw i32 %81, %84
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %80, i64 %86
%88 = load double, double* %87
%89 = call double @llvm.fmuladd.f64(double %70, double %79, double %88)
store double %89, double* %87
br label %90
90: ; preds = %66
%91 = load i32, i32* %21
%92 = add nsw i32 %91, 1
store i32 %92, i32* %21
br label %41
93: ; preds = %41
br label %94
94: ; preds = %93
%95 = load i32, i32* %19
%96 = add nsw i32 %95, 1
store i32 %96, i32* %19
br label %28
97: ; preds = %28
br label %98
98: ; preds = %97
%99 = load i32, i32* %20
%100 = add nsw i32 %99, 1
store i32 %100, i32* %20
br label %23
101: ; preds = %23
store i32 0, i32* %20
br label %102
102: ; preds = %156, %101
%103 = load i32, i32* %20
%104 = load i32, i32* %14
%105 = icmp slt i32 %103, %104
br i1 %105, label %106, label %159
106: ; preds = %102
store i32 0, i32* %19
br label %107
107: ; preds = %152, %106
%108 = load i32, i32* %19
%109 = load i32, i32* %17
%110 = icmp slt i32 %108, %109
br i1 %110, label %111, label %155
111: ; preds = %107
store double 0.000000e+00, double* %22
store i32 0, i32* %21
br label %112
112: ; preds = %137, %111
%113 = load i32, i32* %21
%114 = load i32, i32* %13
%115 = icmp slt i32 %113, %114
br i1 %115, label %116, label %140
116: ; preds = %112
%117 = load double*, double** %12
%118 = load i32, i32* %21
%119 = load i32, i32* %20
%120 = load i32, i32* %13
%121 = mul nsw i32 %119, %120
%122 = add nsw i32 %118, %121
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds double, double* %117, i64 %123
%125 = load double, double* %124
%126 = load double*, double** %18
%127 = load i32, i32* %21
%128 = load i32, i32* %19
%129 = load i32, i32* %13
%130 = mul nsw i32 %128, %129
%131 = add nsw i32 %127, %130
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %126, i64 %132
%134 = load double, double* %133
%135 = load double, double* %22
%136 = call double @llvm.fmuladd.f64(double %125, double %134, double %135)
store double %136, double* %22
br label %137
137: ; preds = %116
%138 = load i32, i32* %21
%139 = add nsw i32 %138, 1
store i32 %139, i32* %21
br label %112
140: ; preds = %112
%141 = load double, double* %22
%142 = load double*, double** %10
%143 = load i32, i32* %20
%144 = load i32, i32* %19
%145 = load i32, i32* %14
%146 = mul nsw i32 %144, %145
%147 = add nsw i32 %143, %146
%148 = sext i32 %147 to i64
%149 = getelementptr inbounds double, double* %142, i64 %148
%150 = load double, double* %149
%151 = fadd double %150, %141
store double %151, double* %149
br label %152
152: ; preds = %140
%153 = load i32, i32* %19
%154 = add nsw i32 %153, 1
store i32 %154, i32* %19
br label %107
155: ; preds = %107
br label %156
156: ; preds = %155
%157 = load i32, i32* %20
%158 = add nsw i32 %157, 1
store i32 %158, i32* %20
br label %102
159: ; preds = %102
ret void
}

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