void svd_daxpy(long n, double da, double *dx, long incx, double *dy,
               long incy) {
  long i;
  if (n <= 0 || incx == 0 || incy == 0 || da == 0.0)
    return;
  if (incx == 1 && incy == 1)
    for (i = 0; i < n; i++) {
      *dy += da * (*dx++);
      dy++;
    }
  else {
    if (incx < 0)
      dx += (-n + 1) * incx;
    if (incy < 0)
      dy += (-n + 1) * incy;
    for (i = 0; i < n; i++) {
      *dy += da * (*dx);
      dx += incx;
      dy += incy;
    }
  }
  return;
}


x86 assembly (clang, -O0):
<svd_daxpy>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
cmpq   $0x0,-0x8(%rbp)
jle    1156 <svd_daxpy+0x56>
cmpq   $0x0,-0x20(%rbp)
je     1156 <svd_daxpy+0x56>
cmpq   $0x0,-0x30(%rbp)
je     1156 <svd_daxpy+0x56>
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    115b <svd_daxpy+0x5b>
jp     115b <svd_daxpy+0x5b>
jmp    1293 <svd_daxpy+0x193>
cmpq   $0x1,-0x20(%rbp)
jne    11d5 <svd_daxpy+0xd5>
cmpq   $0x1,-0x30(%rbp)
jne    11d5 <svd_daxpy+0xd5>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
jge    11d0 <svd_daxpy+0xd0>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x18(%rbp)
movsd  (%rax),%xmm2
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x28(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1179 <svd_daxpy+0x79>
jmp    128e <svd_daxpy+0x18e>
cmpq   $0x0,-0x20(%rbp)
jge    11fb <svd_daxpy+0xfb>
xor    %eax,%eax
sub    -0x8(%rbp),%rax
add    $0x1,%rax
imul   -0x20(%rbp),%rax
shl    $0x3,%rax
add    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
cmpq   $0x0,-0x30(%rbp)
jge    1221 <svd_daxpy+0x121>
xor    %eax,%eax
sub    -0x8(%rbp),%rax
add    $0x1,%rax
imul   -0x30(%rbp),%rax
shl    $0x3,%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x28(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
jge    1289 <svd_daxpy+0x189>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
shl    $0x3,%rax
add    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
mov    -0x30(%rbp),%rax
shl    $0x3,%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1229 <svd_daxpy+0x129>
jmp    128e <svd_daxpy+0x18e>
jmp    1293 <svd_daxpy+0x193>
pop    %rbp
ret

x86 assembly (clang, -O2):
<svd_daxpy>:
push   %r15
push   %r14
push   %r12
push   %rbx
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1117 <svd_daxpy+0x17>
jnp    133c <svd_daxpy+0x23c>
test   %rdi,%rdi
jle    133c <svd_daxpy+0x23c>
test   %rdx,%rdx
je     133c <svd_daxpy+0x23c>
test   %r8,%r8
je     133c <svd_daxpy+0x23c>
cmp    $0x1,%rdx
jne    1223 <svd_daxpy+0x123>
cmp    $0x1,%r8
jne    1223 <svd_daxpy+0x123>
cmp    $0x4,%rdi
jb     1166 <svd_daxpy+0x66>
lea    (%rsi,%rdi,8),%rax
cmp    %rcx,%rax
jbe    1344 <svd_daxpy+0x244>
lea    (%rcx,%rdi,8),%rax
cmp    %rsi,%rax
jbe    1344 <svd_daxpy+0x244>
xor    %r9d,%r9d
mov    %r9,%rdx
not    %rdx
add    %rdi,%rdx
mov    %rdi,%rbx
and    $0x3,%rbx
je     11a4 <svd_daxpy+0xa4>
neg    %rbx
xor    %eax,%eax
movsd  (%rsi),%xmm1
mulsd  %xmm0,%xmm1
add    $0x8,%rsi
addsd  (%rcx),%xmm1
movsd  %xmm1,(%rcx)
add    $0x8,%rcx
add    $0xffffffffffffffff,%rax
cmp    %rax,%rbx
jne    1180 <svd_daxpy+0x80>
sub    %rax,%r9
cmp    $0x3,%rdx
jb     133c <svd_daxpy+0x23c>
sub    %r9,%rdi
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rsi,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rcx,%rax,8),%xmm1
movsd  %xmm1,(%rcx,%rax,8)
movsd  0x8(%rsi,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x8(%rcx,%rax,8),%xmm1
movsd  %xmm1,0x8(%rcx,%rax,8)
movsd  0x10(%rsi,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x10(%rcx,%rax,8),%xmm1
movsd  %xmm1,0x10(%rcx,%rax,8)
movsd  0x18(%rsi,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x18(%rcx,%rax,8),%xmm1
movsd  %xmm1,0x18(%rcx,%rax,8)
add    $0x4,%rax
cmp    %rax,%rdi
jne    11c0 <svd_daxpy+0xc0>
jmp    133c <svd_daxpy+0x23c>
mov    $0x1,%ebx
sub    %rdi,%rbx
mov    %rbx,%r9
imul   %r8,%r9
mov    %r8,%rax
sar    $0x3f,%rax
and    %r9,%rax
lea    (%rcx,%rax,8),%rax
imul   %rdx,%rbx
mov    %rdx,%rcx
sar    $0x3f,%rcx
and    %rbx,%rcx
lea    (%rsi,%rcx,8),%rcx
lea    -0x1(%rdi),%rsi
mov    %edi,%r11d
and    $0x3,%r11d
cmp    $0x3,%rsi
jb     1310 <svd_daxpy+0x210>
and    $0xfffffffffffffffc,%rdi
mov    %rdx,%r9
shl    $0x5,%r9
mov    %r8,%r10
shl    $0x5,%r10
lea    0x0(,%rdx,8),%r15
lea    0x0(,%r8,8),%rsi
nopl   0x0(%rax)
mov    %rcx,%r14
mov    %rax,%rbx
movsd  (%rcx),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rax),%xmm1
movsd  %xmm1,(%rax)
add    %r15,%rcx
add    %rsi,%rax
movsd  (%r14,%r15,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rbx,%rsi,1),%xmm1
movsd  %xmm1,(%rbx,%rsi,1)
lea    (%rax,%rsi,1),%r12
movsd  (%r15,%rcx,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rsi,%rax,1),%xmm1
add    %r15,%rcx
movsd  %xmm1,(%rsi,%rax,1)
lea    (%r12,%rsi,1),%rax
movsd  (%r15,%rcx,1),%xmm1
mulsd  %xmm0,%xmm1
add    %r15,%rcx
addsd  (%rsi,%r12,1),%xmm1
movsd  %xmm1,(%rsi,%r12,1)
add    %r15,%rcx
add    %rsi,%rax
add    $0xfffffffffffffffc,%rdi
jne    1290 <svd_daxpy+0x190>
add    %r9,%r14
add    %r10,%rbx
mov    %r14,%rcx
mov    %rbx,%rax
test   %r11,%r11
je     133c <svd_daxpy+0x23c>
shl    $0x3,%rdx
shl    $0x3,%r8
nopl   (%rax)
movsd  (%rcx),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rax),%xmm1
movsd  %xmm1,(%rax)
add    %rdx,%rcx
add    %r8,%rax
add    $0xffffffffffffffff,%r11
jne    1320 <svd_daxpy+0x220>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
mov    %rdi,%r9
and    $0xfffffffffffffffc,%r9
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%r9),%rax
mov    %rax,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rax,%rax
je     143c <svd_daxpy+0x33c>
mov    %r8,%rbx
and    $0xfffffffffffffffe,%rbx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd (%rsi,%rdx,8),%xmm2
movupd 0x10(%rsi,%rdx,8),%xmm3
mulpd  %xmm1,%xmm2
movupd (%rcx,%rdx,8),%xmm4
addpd  %xmm2,%xmm4
movupd 0x10(%rcx,%rdx,8),%xmm2
movupd 0x20(%rcx,%rdx,8),%xmm5
movupd 0x30(%rcx,%rdx,8),%xmm6
mulpd  %xmm1,%xmm3
addpd  %xmm2,%xmm3
movupd %xmm4,(%rcx,%rdx,8)
movupd %xmm3,0x10(%rcx,%rdx,8)
movupd 0x20(%rsi,%rdx,8),%xmm2
movupd 0x30(%rsi,%rdx,8),%xmm3
mulpd  %xmm1,%xmm2
addpd  %xmm5,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm6,%xmm3
movupd %xmm2,0x20(%rcx,%rdx,8)
movupd %xmm3,0x30(%rcx,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rbx
jne    1380 <svd_daxpy+0x280>
test   $0x1,%r8b
je     1426 <svd_daxpy+0x326>
movupd (%rsi,%rdx,8),%xmm2
movupd 0x10(%rsi,%rdx,8),%xmm3
movupd (%rcx,%rdx,8),%xmm4
movupd 0x10(%rcx,%rdx,8),%xmm5
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm2,(%rcx,%rdx,8)
movupd %xmm3,0x10(%rcx,%rdx,8)
cmp    %rdi,%r9
je     133c <svd_daxpy+0x23c>
lea    (%rsi,%r9,8),%rsi
lea    (%rcx,%r9,8),%rcx
jmp    1169 <svd_daxpy+0x69>
xor    %edx,%edx
test   $0x1,%r8b
jne    13f5 <svd_daxpy+0x2f5>
jmp    1426 <svd_daxpy+0x326>

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

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