void Mupdate(m, n, alpha, x, y, A, j0) int m, n, j0;
double alpha;
double **A, *x, *y;
{
  register int i, j, n4;
  register double *Aref;
  register double tmp;
  n4 = n / 4;
  n = n % 4;
  for (i = 0; i < m; i++) {
    tmp = alpha * x[i];
    Aref = &(A[i][j0]);
    for (j = 0; j < n4; j++) {
      Aref[4 * j] += tmp * y[4 * j];
      Aref[4 * j + 1] += tmp * y[4 * j + 1];
      Aref[4 * j + 2] += tmp * y[4 * j + 2];
      Aref[4 * j + 3] += tmp * y[4 * j + 3];
    }
    Aref += 4 * n4;
    y += 4 * n4;
    for (j = 0; j < n; j++)
      Aref[j] += tmp * y[j];
  }
}


x86 assembly (clang, -O0):
<Mupdate>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0x8(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,-0x38(%rbp)
mov    -0x8(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %edx,-0x8(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    131b <Mupdate+0x21b>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    128b <Mupdate+0x18b>
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%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    1187 <Mupdate+0x87>
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
mov    -0x40(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1308 <Mupdate+0x208>
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%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    12c2 <Mupdate+0x1c2>
jmp    130d <Mupdate+0x20d>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1142 <Mupdate+0x42>
pop    %rbp
ret

x86 assembly (clang, -O2):
<Mupdate>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
lea    0x3(%rsi),%eax
test   %esi,%esi
cmovns %esi,%eax
test   %edi,%edi
jle    13f5 <Mupdate+0x2f5>
mov    %rcx,%r10
mov    %eax,%ecx
sar    $0x2,%ecx
and    $0xfffffffc,%eax
mov    %esi,%r11d
sub    %eax,%r11d
movslq %r9d,%rbx
lea    0x0(,%rcx,4),%eax
movslq %eax,%r13
mov    %edi,%eax
mov    %rax,-0x8(%rsp)
mov    %rcx,-0x18(%rsp)
mov    %ecx,%ebp
lea    0x0(,%rbp,4),%rax
mov    %rax,-0x28(%rsp)
lea    (%rbx,%rbp,4),%rax
mov    %rax,-0x30(%rsp)
mov    %ebp,%eax
and    $0xfffffffe,%eax
mov    %rax,-0x38(%rsp)
lea    0x0(,%r13,8),%rax
mov    %rax,-0x10(%rsp)
mov    %rbp,-0x40(%rsp)
shl    $0x5,%rbp
and    $0xffffffffffffffc0,%rbp
mov    %rbx,%rcx
lea    0x3(%rbx),%rax
mov    %rax,-0x20(%rsp)
xor    %r14d,%r14d
mov    %r10,%r9
jmp    11af <Mupdate+0xaf>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r14
cmp    -0x8(%rsp),%r14
je     13f5 <Mupdate+0x2f5>
movsd  (%rdx,%r14,8),%xmm9
mulsd  %xmm0,%xmm9
mov    (%r8,%r14,8),%r15
lea    (%r15,%rcx,8),%r12
cmp    $0x4,%esi
jl     1287 <Mupdate+0x187>
cmpl   $0x2,-0x18(%rsp)
jb     11fd <Mupdate+0xfd>
mov    %r14,%rax
imul   %r13,%rax
add    -0x28(%rsp),%rax
lea    (%r10,%rax,8),%rax
cmp    %rax,%r12
jae    12f1 <Mupdate+0x1f1>
mov    -0x30(%rsp),%rax
lea    (%r15,%rax,8),%rax
cmp    %rax,%r9
jae    12f1 <Mupdate+0x1f1>
xor    %edi,%edi
mov    -0x20(%rsp),%rax
lea    (%r15,%rax,8),%rbx
mov    -0x40(%rsp),%rax
sub    %rdi,%rax
shl    $0x5,%rdi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%r9,%rdi,1),%xmm1
mulsd  %xmm9,%xmm1
addsd  -0x18(%rbx,%rdi,1),%xmm1
movsd  %xmm1,-0x18(%rbx,%rdi,1)
movsd  0x8(%r9,%rdi,1),%xmm1
mulsd  %xmm9,%xmm1
addsd  -0x10(%rbx,%rdi,1),%xmm1
movsd  %xmm1,-0x10(%rbx,%rdi,1)
movsd  0x10(%r9,%rdi,1),%xmm1
mulsd  %xmm9,%xmm1
addsd  -0x8(%rbx,%rdi,1),%xmm1
movsd  %xmm1,-0x8(%rbx,%rdi,1)
movsd  0x18(%r9,%rdi,1),%xmm1
mulsd  %xmm9,%xmm1
addsd  (%rbx,%rdi,1),%xmm1
movsd  %xmm1,(%rbx,%rdi,1)
add    $0x20,%rdi
add    $0xffffffffffffffff,%rax
jne    1220 <Mupdate+0x120>
add    -0x10(%rsp),%r9
test   %r11d,%r11d
jle    11a0 <Mupdate+0xa0>
movsd  (%r9),%xmm1
mulsd  %xmm9,%xmm1
addsd  (%r12,%r13,8),%xmm1
movsd  %xmm1,(%r12,%r13,8)
cmp    $0x1,%r11d
je     11a0 <Mupdate+0xa0>
movsd  0x8(%r9),%xmm1
mulsd  %xmm9,%xmm1
addsd  0x8(%r12,%r13,8),%xmm1
movsd  %xmm1,0x8(%r12,%r13,8)
cmp    $0x2,%r11d
je     11a0 <Mupdate+0xa0>
mulsd  0x10(%r9),%xmm9
addsd  0x10(%r12,%r13,8),%xmm9
movsd  %xmm9,0x10(%r12,%r13,8)
jmp    11a0 <Mupdate+0xa0>
mov    %r10,%rax
movapd %xmm9,%xmm8
unpcklpd %xmm9,%xmm8
xor    %r10d,%r10d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movupd (%r9,%r10,1),%xmm3
movupd 0x10(%r9,%r10,1),%xmm4
movupd 0x20(%r9,%r10,1),%xmm6
movupd 0x30(%r9,%r10,1),%xmm7
movapd %xmm3,%xmm5
unpcklpd %xmm6,%xmm5
unpckhpd %xmm6,%xmm3
movapd %xmm4,%xmm6
unpcklpd %xmm7,%xmm6
unpckhpd %xmm7,%xmm4
movupd (%r12,%r10,1),%xmm7
movupd 0x10(%r12,%r10,1),%xmm11
movupd 0x20(%r12,%r10,1),%xmm1
movupd 0x30(%r12,%r10,1),%xmm10
movapd %xmm7,%xmm2
unpcklpd %xmm1,%xmm2
unpckhpd %xmm1,%xmm7
movapd %xmm11,%xmm1
unpcklpd %xmm10,%xmm1
unpckhpd %xmm10,%xmm11
mulpd  %xmm8,%xmm5
addpd  %xmm2,%xmm5
mulpd  %xmm8,%xmm3
addpd  %xmm7,%xmm3
mulpd  %xmm8,%xmm6
addpd  %xmm1,%xmm6
mulpd  %xmm8,%xmm4
addpd  %xmm11,%xmm4
movapd %xmm6,%xmm1
unpcklpd %xmm4,%xmm1
movapd %xmm5,%xmm2
unpcklpd %xmm3,%xmm2
movhlps %xmm6,%xmm4
movhlps %xmm5,%xmm3
movups %xmm3,0x20(%r12,%r10,1)
movups %xmm4,0x30(%r12,%r10,1)
movupd %xmm2,(%r12,%r10,1)
movupd %xmm1,0x10(%r12,%r10,1)
add    $0x40,%r10
cmp    %r10,%rbp
jne    1310 <Mupdate+0x210>
mov    -0x38(%rsp),%rbx
mov    %rbx,%rdi
cmp    -0x40(%rsp),%rbx
mov    %rax,%r10
je     1287 <Mupdate+0x187>
jmp    11ff <Mupdate+0xff>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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