void _inner_ddot_trans(double *Cr, const double *Ar, const double *Br,
                       double alpha, int nVP, int ldB) {
  register int k, iB;
  register double f0, f1, f2, f3, cval0, cval1;
  cval0 = 0.0;
  cval1 = 0.0;
  for (k = 0; k < nVP - 3; k += 4) {
    f0 = Ar[0] * Br[0];
    cval0 += f0;
    f1 = Ar[1] * Br[ldB];
    cval1 += f1;
    f2 = Ar[2] * Br[2 * ldB];
    cval0 += f2;
    f3 = Ar[3] * Br[3 * ldB];
    cval1 += f3;
    Br += 4 * ldB;
    Ar += 4;
  }
  if (k == nVP)
    goto update;
  if (k < nVP - 1) {
    f0 = Ar[0] * Br[0];
    cval0 += f0;
    f1 = Ar[1] * Br[ldB];
    cval1 += f1;
    Br += 2 * ldB;
    Ar += 2;
    k += 2;
  }
  if (k < nVP) {
    f0 = Ar[0] * Br[0];
    cval0 += f0;
    Br += ldB;
    Ar++;
    k++;
  }
update:
  f0 = (cval0 + cval1) * alpha;
  Cr[0] += f0;
}


x86 assembly (clang, -O0):
<_inner_ddot_trans>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8d,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    $0x3,%ecx
cmp    %ecx,%eax
jge    1220 <_inner_ddot_trans+0x120>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x18(%rbp),%rax
imul   $0x3,-0x28(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%ecx
shl    $0x2,%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
add    $0x20,%rax
mov    %rax,-0x10(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x2c(%rbp)
jmp    1133 <_inner_ddot_trans+0x33>
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jne    1231 <_inner_ddot_trans+0x131>
jmp    131c <_inner_ddot_trans+0x21c>
mov    -0x2c(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12bd <_inner_ddot_trans+0x1bd>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
add    $0x10,%rax
mov    %rax,-0x10(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1317 <_inner_ddot_trans+0x217>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x28(%rbp),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x10(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    131c <_inner_ddot_trans+0x21c>
movsd  -0x58(%rbp),%xmm0
addsd  -0x60(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<_inner_ddot_trans>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    $0x4,%ecx
jl     1151 <_inner_ddot_trans+0x51>
movslq %r8d,%r14
lea    (%r14,%r14,1),%eax
movslq %eax,%r9
lea    (%r14,%r14,2),%eax
movslq %eax,%r10
lea    0x0(,%r14,4),%eax
movslq %eax,%r11
lea    -0x4(%rcx),%eax
mov    %eax,%r15d
shr    $0x2,%r15d
add    $0x1,%r15d
cmp    $0x4,%eax
jae    1164 <_inner_ddot_trans+0x64>
xorpd  %xmm1,%xmm1
test   $0x1,%r15b
jne    1246 <_inner_ddot_trans+0x146>
jmp    127d <_inner_ddot_trans+0x17d>
xorpd  %xmm1,%xmm1
xor    %eax,%eax
cmp    %ecx,%eax
jne    1286 <_inner_ddot_trans+0x186>
jmp    12c6 <_inner_ddot_trans+0x1c6>
mov    %r8d,-0x1c(%rsp)
mov    %r15d,%ebp
and    $0xfffffffe,%ebp
lea    0x0(,%r14,8),%r12
mov    %r11,%rbx
shl    $0x4,%rbx
mov    %r9,-0x10(%rsp)
lea    0x0(,%r9,8),%r13
mov    %r10,-0x18(%rsp)
lea    0x0(,%r10,8),%rax
lea    0x0(,%r11,8),%r9
lea    (%rax,%r11,8),%r10
mov    %r11,-0x8(%rsp)
lea    0x0(,%r11,8),%r11
add    %r13,%r11
lea    (%r9,%r14,8),%r8
xorpd  %xmm1,%xmm1
nopl   0x0(%rax)
movsd  (%rdx),%xmm2
movhpd (%rdx,%r12,1),%xmm2
movupd (%rsi),%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm1,%xmm3
movupd 0x10(%rsi),%xmm1
movupd 0x20(%rsi),%xmm2
movsd  (%rdx,%r13,1),%xmm4
movhpd (%rdx,%rax,1),%xmm4
movupd 0x30(%rsi),%xmm5
mulpd  %xmm1,%xmm4
addpd  %xmm3,%xmm4
movsd  (%rdx,%r9,1),%xmm3
movhpd (%rdx,%r8,1),%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm4,%xmm3
movsd  (%rdx,%r11,1),%xmm1
movhpd (%rdx,%r10,1),%xmm1
mulpd  %xmm5,%xmm1
addpd  %xmm3,%xmm1
add    $0x40,%rsi
add    %rbx,%rdx
add    $0xfffffffe,%ebp
jne    11c0 <_inner_ddot_trans+0xc0>
mov    -0x1c(%rsp),%r8d
mov    -0x8(%rsp),%r11
mov    -0x10(%rsp),%r9
mov    -0x18(%rsp),%r10
test   $0x1,%r15b
je     127d <_inner_ddot_trans+0x17d>
movsd  (%rdx),%xmm2
movhpd (%rdx,%r14,8),%xmm2
movupd (%rsi),%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm3,%xmm1
movupd 0x10(%rsi),%xmm2
movsd  (%rdx,%r9,8),%xmm3
movhpd (%rdx,%r10,8),%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm3,%xmm1
lea    (%rdx,%r11,8),%rdx
add    $0x20,%rsi
mov    %ecx,%eax
and    $0xfffffffc,%eax
cmp    %ecx,%eax
je     12c6 <_inner_ddot_trans+0x1c6>
lea    -0x1(%rcx),%ebp
cmp    %ebp,%eax
jge    12b6 <_inner_ddot_trans+0x1b6>
movslq %r8d,%rbp
movsd  (%rdx),%xmm2
movhpd (%rdx,%rbp,8),%xmm2
movupd (%rsi),%xmm3
mulpd  %xmm3,%xmm2
addpd  %xmm2,%xmm1
add    %r8d,%r8d
movslq %r8d,%rbp
lea    (%rdx,%rbp,8),%rdx
add    $0x10,%rsi
or     $0x2,%eax
cmp    %ecx,%eax
jge    12c6 <_inner_ddot_trans+0x1c6>
movsd  (%rsi),%xmm2
mulsd  (%rdx),%xmm2
addsd  %xmm2,%xmm1
movapd %xmm1,%xmm2
unpckhpd %xmm1,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm2
addsd  (%rdi),%xmm2
movsd  %xmm2,(%rdi)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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