void babl_matrix_mul_vector_buf4(const double *mat, const double *v_in,
                                 double *v_out, int samples) {
  int i;
  const double m_0_0 = mat[0 * 3 + 0];
  const double m_0_1 = mat[0 * 3 + 1];
  const double m_0_2 = mat[0 * 3 + 2];
  const double m_1_0 = mat[1 * 3 + 0];
  const double m_1_1 = mat[1 * 3 + 1];
  const double m_1_2 = mat[1 * 3 + 2];
  const double m_2_0 = mat[2 * 3 + 0];
  const double m_2_1 = mat[2 * 3 + 1];
  const double m_2_2 = mat[2 * 3 + 2];
  for (i = 0; i < samples; i++) {
    const double a = v_in[0], b = v_in[1], c = v_in[2];
    v_out[0] = m_0_0 * a + m_0_1 * b + m_0_2 * c;
    v_out[1] = m_1_0 * a + m_1_1 * b + m_1_2 * c;
    v_out[2] = m_2_0 * a + m_2_1 * b + m_2_2 * c;
    v_out[3] = v_in[3];
    v_in += 4;
    v_out += 4;
  }
}


x86 assembly (clang, -O0):
<babl_matrix_mul_vector_buf4>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x40(%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12a8 <babl_matrix_mul_vector_buf4+0x1a8>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x38(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x40(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x58(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm2
movsd  -0x60(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x10(%rax)
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x18(%rax)
mov    -0x10(%rbp),%rax
add    $0x20,%rax
mov    %rax,-0x10(%rbp)
mov    -0x18(%rbp),%rax
add    $0x20,%rax
mov    %rax,-0x18(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1197 <babl_matrix_mul_vector_buf4+0x97>
pop    %rbp
ret

x86 assembly (clang, -O2):
<babl_matrix_mul_vector_buf4>:
sub    $0x38,%rsp
test   %ecx,%ecx
jle    120b <babl_matrix_mul_vector_buf4+0x10b>
movsd  (%rdi),%xmm0
movsd  0x8(%rdi),%xmm8
movsd  0x10(%rdi),%xmm9
movsd  0x18(%rdi),%xmm3
movsd  0x20(%rdi),%xmm4
movsd  0x28(%rdi),%xmm11
movsd  0x30(%rdi),%xmm10
movsd  0x38(%rdi),%xmm1
movsd  0x40(%rdi),%xmm2
xor    %r8d,%r8d
mov    %ecx,%r9d
add    $0xffffffff,%r9d
je     1170 <babl_matrix_mul_vector_buf4+0x70>
mov    %r9,%rax
shl    $0x5,%rax
lea    (%rsi,%rax,1),%rdi
add    $0x20,%rdi
cmp    %rdx,%rdi
jbe    1210 <babl_matrix_mul_vector_buf4+0x110>
add    %rdx,%rax
add    $0x20,%rax
cmp    %rsi,%rax
jbe    1210 <babl_matrix_mul_vector_buf4+0x110>
mov    %rsi,%r11
mov    %rdx,%r10
unpcklpd %xmm4,%xmm8
unpcklpd %xmm3,%xmm0
unpcklpd %xmm11,%xmm9
sub    %r8d,%ecx
xor    %eax,%eax
nopl   0x0(%rax)
movsd  (%r11,%rax,1),%xmm3
movsd  0x8(%r11,%rax,1),%xmm6
movsd  0x10(%r11,%rax,1),%xmm7
movapd %xmm1,%xmm4
mulsd  %xmm6,%xmm4
unpcklpd %xmm6,%xmm6
mulpd  %xmm8,%xmm6
movapd %xmm10,%xmm5
mulsd  %xmm3,%xmm5
unpcklpd %xmm3,%xmm3
mulpd  %xmm0,%xmm3
addpd  %xmm6,%xmm3
movapd %xmm2,%xmm6
mulsd  %xmm7,%xmm6
unpcklpd %xmm7,%xmm7
mulpd  %xmm9,%xmm7
addpd  %xmm3,%xmm7
movupd %xmm7,(%r10,%rax,1)
addsd  %xmm4,%xmm5
addsd  %xmm6,%xmm5
movsd  %xmm5,0x10(%r10,%rax,1)
movsd  0x18(%r11,%rax,1),%xmm3
movsd  %xmm3,0x18(%r10,%rax,1)
add    $0x20,%rax
add    $0xffffffff,%ecx
jne    1190 <babl_matrix_mul_vector_buf4+0x90>
add    $0x38,%rsp
ret
add    $0x1,%r9
mov    %r9,%r8
and    $0xfffffffffffffffe,%r8
mov    %r8,%r10
shl    $0x5,%r10
lea    (%rsi,%r10,1),%r11
add    %rdx,%r10
movapd %xmm8,%xmm5
unpcklpd %xmm8,%xmm5
movapd %xmm5,0x20(%rsp)
movapd %xmm0,-0x30(%rsp)
movapd %xmm0,%xmm7
unpcklpd %xmm0,%xmm7
movapd %xmm9,%xmm0
unpcklpd %xmm9,%xmm0
movapd %xmm0,0x10(%rsp)
movapd %xmm4,-0x60(%rsp)
unpcklpd %xmm4,%xmm4
movapd %xmm4,(%rsp)
movapd %xmm3,-0x40(%rsp)
unpcklpd %xmm3,%xmm3
movapd %xmm3,-0x10(%rsp)
movapd %xmm11,-0x50(%rsp)
movapd %xmm11,%xmm0
unpcklpd %xmm11,%xmm0
movapd %xmm1,-0x20(%rsp)
unpcklpd %xmm1,%xmm1
movapd %xmm10,-0x80(%rsp)
movapd %xmm10,%xmm5
unpcklpd %xmm10,%xmm5
movapd %xmm2,-0x70(%rsp)
movapd %xmm2,%xmm4
movapd %xmm7,%xmm2
unpcklpd %xmm4,%xmm4
xor    %edi,%edi
mov    %r8,%rax
nopl   0x0(%rax)
movupd (%rsi,%rdi,1),%xmm14
movupd 0x10(%rsi,%rdi,1),%xmm12
movupd 0x20(%rsi,%rdi,1),%xmm6
movupd 0x30(%rsi,%rdi,1),%xmm11
movapd %xmm14,%xmm13
unpcklpd %xmm6,%xmm13
unpckhpd %xmm6,%xmm14
movapd %xmm12,%xmm10
unpcklpd %xmm11,%xmm10
movapd 0x20(%rsp),%xmm6
mulpd  %xmm14,%xmm6
movapd %xmm2,%xmm15
mulpd  %xmm13,%xmm15
addpd  %xmm6,%xmm15
movapd 0x10(%rsp),%xmm6
mulpd  %xmm10,%xmm6
addpd  %xmm15,%xmm6
movapd (%rsp),%xmm7
mulpd  %xmm14,%xmm7
movapd -0x10(%rsp),%xmm3
mulpd  %xmm13,%xmm3
addpd  %xmm7,%xmm3
movapd %xmm0,%xmm7
mulpd  %xmm10,%xmm7
addpd  %xmm3,%xmm7
mulpd  %xmm1,%xmm14
mulpd  %xmm5,%xmm13
addpd  %xmm14,%xmm13
mulpd  %xmm4,%xmm10
addpd  %xmm13,%xmm10
movapd %xmm6,%xmm3
unpcklpd %xmm7,%xmm3
movhlps %xmm6,%xmm7
movsd  %xmm10,%xmm12
unpckhpd %xmm11,%xmm10
movups %xmm7,0x20(%rdx,%rdi,1)
movupd %xmm10,0x30(%rdx,%rdi,1)
movupd %xmm3,(%rdx,%rdi,1)
movupd %xmm12,0x10(%rdx,%rdi,1)
add    $0x40,%rdi
add    $0xfffffffffffffffe,%rax
jne    12c0 <babl_matrix_mul_vector_buf4+0x1c0>
cmp    %r8,%r9
movapd -0x80(%rsp),%xmm10
movapd -0x20(%rsp),%xmm1
movapd -0x70(%rsp),%xmm2
movapd -0x30(%rsp),%xmm0
movapd -0x40(%rsp),%xmm3
movapd -0x50(%rsp),%xmm11
movapd -0x60(%rsp),%xmm4
jne    1176 <babl_matrix_mul_vector_buf4+0x76>
jmp    120b <babl_matrix_mul_vector_buf4+0x10b>

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

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