void rtrend(float *y, int n) {
  int i;
  double a, b, a11, a12, a22, y1, y2;
  y1 = y2 = 0.;
  for (i = 0; i < n; i++) {
    y1 += i * y[i];
    y2 += y[i];
  }
  a12 = 0.5 * n * (n - 1);
  a11 = a12 * (2 * n - 1) / 3.;
  a22 = n;
  b = a11 * a22 - a12 * a12;
  a = (a22 * y1 - a12 * y2) / b;
  b = (a11 * y2 - a12 * y1) / b;
  for (i = 0; i < n; i++) {
    y[i] = y[i] - a * i - b;
  }
}


x86 assembly (clang, -O0):
<rtrend>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1177 <rtrend+0x77>
cvtsi2ssl -0x10(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1122 <rtrend+0x22>
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  0xe84(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe47(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  -0x30(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  -0x48(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12ef <rtrend+0x1ef>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  -0x18(%rbp),%xmm0
cvtsi2sdl -0x10(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
subsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1285 <rtrend+0x185>
pop    %rbp
ret

x86 assembly (clang, -O2):
<rtrend>:
mov    %esi,%r9d
test   %esi,%esi
jle    1114 <rtrend+0x14>
cmp    $0x1,%esi
jne    111a <rtrend+0x1a>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    1171 <rtrend+0x71>
xorpd  %xmm0,%xmm0
jmp    1191 <rtrend+0x91>
mov    %r9d,%edx
and    $0xfffffffe,%edx
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
movss  (%rdi,%rcx,4),%xmm2
movss  0x4(%rdi,%rcx,4),%xmm3
mulss  %xmm2,%xmm1
unpcklps %xmm2,%xmm1
cvtps2pd %xmm1,%xmm1
addpd  %xmm0,%xmm1
lea    0x1(%rcx),%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
mulss  %xmm3,%xmm0
unpcklps %xmm3,%xmm0
cvtps2pd %xmm0,%xmm0
addpd  %xmm1,%xmm0
add    $0x2,%rcx
cmp    %rdx,%rcx
jne    1130 <rtrend+0x30>
test   $0x1,%r9b
je     1191 <rtrend+0x91>
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
movss  (%rdi,%rcx,4),%xmm2
mulss  %xmm2,%xmm1
unpcklps %xmm2,%xmm1
cvtps2pd %xmm1,%xmm1
addpd  %xmm1,%xmm0
test   %esi,%esi
jle    13a1 <rtrend+0x2a1>
lea    -0x1(%rsi),%eax
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
movsd  0xe55(%rip),%xmm1
mulsd  %xmm2,%xmm1
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
mulsd  %xmm1,%xmm3
lea    -0x1(%rsi,%rsi,1),%eax
cvtsi2sd %eax,%xmm4
mulsd  %xmm3,%xmm4
divsd  0xe3a(%rip),%xmm4
movapd 0xe4a(%rip),%xmm5
xorpd  %xmm3,%xmm5
movapd %xmm0,%xmm6
mulsd  %xmm3,%xmm6
mulsd  %xmm3,%xmm3
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
mulsd  %xmm1,%xmm5
mulsd  %xmm4,%xmm1
mulsd  %xmm2,%xmm4
subsd  %xmm3,%xmm4
subsd  %xmm6,%xmm1
divsd  %xmm4,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm5,%xmm0
xorpd  0xe0a(%rip),%xmm0
divsd  %xmm4,%xmm0
cmp    $0x4,%esi
jae    1226 <rtrend+0x126>
xor    %ecx,%ecx
jmp    1370 <rtrend+0x270>
mov    %r9d,%ecx
and    $0xfffffffc,%ecx
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
movapd %xmm1,%xmm10
unpcklpd %xmm1,%xmm10
lea    -0x4(%rcx),%rax
mov    %rax,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rax,%rax
je     13a2 <rtrend+0x2a2>
mov    %r8,%rsi
and    $0xfffffffffffffffe,%rsi
movdqa 0xddb(%rip),%xmm4
xor    %edx,%edx
movdqa 0xde0(%rip),%xmm8
movdqa 0xde7(%rip),%xmm9
nopl   0x0(%rax)
cvtps2pd 0x8(%rdi,%rdx,4),%xmm7
cvtps2pd (%rdi,%rdx,4),%xmm5
pshufd $0xee,%xmm4,%xmm6
cvtdq2pd %xmm6,%xmm6
cvtdq2pd %xmm4,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm5,%xmm2
mulpd  %xmm3,%xmm6
addpd  %xmm7,%xmm6
subpd  %xmm10,%xmm6
subpd  %xmm10,%xmm2
cvtpd2ps %xmm2,%xmm2
cvtpd2ps %xmm6,%xmm5
unpcklpd %xmm5,%xmm2
movupd %xmm2,(%rdi,%rdx,4)
movdqa %xmm4,%xmm2
cvtps2pd 0x18(%rdi,%rdx,4),%xmm5
paddd  %xmm8,%xmm2
cvtps2pd 0x10(%rdi,%rdx,4),%xmm6
pshufd $0xee,%xmm2,%xmm7
cvtdq2pd %xmm7,%xmm7
cvtdq2pd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm6,%xmm2
mulpd  %xmm3,%xmm7
addpd  %xmm5,%xmm7
subpd  %xmm10,%xmm2
cvtpd2ps %xmm2,%xmm2
subpd  %xmm10,%xmm7
cvtpd2ps %xmm7,%xmm5
unpcklpd %xmm5,%xmm2
movupd %xmm2,0x10(%rdi,%rdx,4)
add    $0x8,%rdx
paddd  %xmm9,%xmm4
add    $0xfffffffffffffffe,%rsi
jne    1280 <rtrend+0x180>
cvtdq2pd %xmm4,%xmm5
pshufd $0xee,%xmm4,%xmm2
cvtdq2pd %xmm2,%xmm4
test   $0x1,%r8b
je     1367 <rtrend+0x267>
cvtps2pd 0x8(%rdi,%rdx,4),%xmm2
cvtps2pd (%rdi,%rdx,4),%xmm6
mulpd  %xmm3,%xmm5
addpd  %xmm6,%xmm5
mulpd  %xmm3,%xmm4
addpd  %xmm2,%xmm4
subpd  %xmm10,%xmm4
subpd  %xmm10,%xmm5
cvtpd2ps %xmm5,%xmm2
cvtpd2ps %xmm4,%xmm3
unpcklpd %xmm3,%xmm2
movupd %xmm2,(%rdi,%rdx,4)
cmp    %r9,%rcx
je     13a1 <rtrend+0x2a1>
nopl   0x0(%rax)
movss  (%rdi,%rcx,4),%xmm2
cvtss2sd %xmm2,%xmm2
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
subsd  %xmm1,%xmm3
xorps  %xmm2,%xmm2
cvtsd2ss %xmm3,%xmm2
movss  %xmm2,(%rdi,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%r9
jne    1370 <rtrend+0x270>
ret
xorpd  %xmm5,%xmm5
movhpd 0xc62(%rip),%xmm5
movapd 0xc7a(%rip),%xmm4
xor    %edx,%edx
test   $0x1,%r8b
jne    1333 <rtrend+0x233>
jmp    1367 <rtrend+0x267>

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

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