int parfit(float *y, int ip, int *n) {
  int i, imax, n2, i1, i2;
  float b, c, sumx2, sumx4, sumy, sumxy, sumx2y;
  n2 = (*n) / 2;
  *n = 2 * n2 + 1;
  i1 = ip - n2;
  i2 = ip + n2;
  sumx2 = 0.0;
  sumx4 = 0.0;
  sumy = 0.0;
  sumxy = 0.0;
  sumx2y = 0.0;
  for (i = i1; i <= i2; i++) {
    int ii, ix;
    float yi;
    yi = y[i];
    ix = i - ip;
    ii = ix * ix;
    sumx2 += ii;
    sumx4 += ii * ii;
    sumy += yi;
    sumxy += ix * yi;
    sumx2y += ii * yi;
  }
  b = sumxy * ((*n) * sumx4 - sumx2 * sumx2);
  c = sumx2 * ((*n) * sumx2y - sumx2 * sumy);
  imax = ip - 0.5 * b / c + 0.5;
  if (c > 0)
    imax = -imax;
  return imax;
}


x86 assembly (clang, -O0):
<parfit>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0xc(%rbp),%eax
sub    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0xc(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x28(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1219 <parfit+0x119>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x54(%rbp)
mov    -0x1c(%rbp),%eax
sub    -0xc(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
imul   -0x50(%rbp),%eax
mov    %eax,-0x4c(%rbp)
cvtsi2ssl -0x4c(%rbp),%xmm0
addss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x4c(%rbp),%eax
imul   -0x4c(%rbp),%eax
cvtsi2ss %eax,%xmm0
addss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x54(%rbp),%xmm0
addss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
cvtsi2ssl -0x50(%rbp),%xmm0
movss  -0x54(%rbp),%xmm2
movss  -0x44(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
cvtsi2ssl -0x4c(%rbp),%xmm0
movss  -0x54(%rbp),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    116f <parfit+0x6f>
movss  -0x44(%rbp),%xmm0
mov    -0x18(%rbp),%rax
cvtsi2ssl (%rax),%xmm1
movss  -0x3c(%rbp),%xmm3
movss  -0x38(%rbp),%xmm2
mulss  -0x38(%rbp),%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
cvtsi2ssl (%rax),%xmm1
movss  -0x48(%rbp),%xmm3
movss  -0x38(%rbp),%xmm2
mulss  -0x40(%rbp),%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
movss  -0x30(%rbp),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  0xd5d(%rip),%xmm1
mulsd  %xmm2,%xmm1
movss  -0x34(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  0xd40(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x20(%rbp)
movss  -0x34(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    12e4 <parfit+0x1e4>
xor    %eax,%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<parfit>:
mov    (%rdx),%eax
mov    %eax,%r8d
shr    $0x1f,%r8d
add    %eax,%r8d
mov    %r8d,%ecx
sar    %ecx
or     $0x1,%r8d
mov    %r8d,(%rdx)
mov    %esi,%eax
sub    %ecx,%eax
lea    (%rcx,%rsi,1),%edx
xorps  %xmm1,%xmm1
xorps  %xmm0,%xmm0
xorps  %xmm3,%xmm3
xorps  %xmm4,%xmm4
xorps  %xmm2,%xmm2
cmp    %edx,%eax
jg     11a0 <parfit+0xa0>
cltq
mov    %ecx,%edx
neg    %edx
lea    (%rdi,%rax,4),%r9
add    %ecx,%ecx
or     $0x1,%rcx
xorps  %xmm6,%xmm6
xorps  %xmm2,%xmm2
xor    %eax,%eax
xorps  %xmm4,%xmm4
xorps  %xmm3,%xmm3
movaps %xmm6,%xmm5
movss  (%r9,%rax,4),%xmm1
lea    (%rdx,%rax,1),%edi
movd   %edi,%xmm6
imul   %edi,%edi
movd   %edi,%xmm7
punpckldq %xmm7,%xmm6
cvtdq2ps %xmm6,%xmm6
addss  %xmm1,%xmm4
shufps $0x0,%xmm1,%xmm1
mulps  %xmm6,%xmm1
shufps $0x55,%xmm6,%xmm6
addss  %xmm6,%xmm2
imul   %edi,%edi
xorps  %xmm6,%xmm6
cvtsi2ss %edi,%xmm6
addss  %xmm6,%xmm3
movaps %xmm1,%xmm6
addps  %xmm5,%xmm6
add    $0x1,%rax
cmp    %eax,%ecx
jne    1150 <parfit+0x50>
addps  %xmm5,%xmm1
xorps  %xmm5,%xmm5
cvtsi2ss %r8d,%xmm5
movaps %xmm2,%xmm6
mulss  %xmm2,%xmm6
mulss  %xmm5,%xmm3
subss  %xmm6,%xmm3
mulss  %xmm1,%xmm3
mulss  %xmm2,%xmm4
shufps $0x55,%xmm1,%xmm1
mulss  %xmm5,%xmm1
subss  %xmm4,%xmm1
mulss  %xmm2,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
cvtss2sd %xmm3,%xmm3
movsd  0xe1e(%rip),%xmm4
mulsd  %xmm4,%xmm3
xorps  %xmm5,%xmm5
cvtss2sd %xmm1,%xmm5
divsd  %xmm5,%xmm3
subsd  %xmm3,%xmm2
addsd  %xmm4,%xmm2
cvttsd2si %xmm2,%ecx
mov    %ecx,%eax
neg    %eax
ucomiss %xmm0,%xmm1
cmovbe %ecx,%eax
ret

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

LLVM IR (-O2):
define i32 @parfit(float* nocapture %0, i32 %1, i32* nocapture %2) {
%4 = load i32, i32* %2
%5 = sdiv i32 %4, 2
%6 = shl nsw i32 %5, 1
%7 = or i32 %6, 1
store i32 %7, i32* %2
%8 = sub i32 %1, %5
%9 = add i32 %5, %1
%10 = icmp sgt i32 %8, %9
br i1 %10, label %42, label %11
11: ; preds = %3
%12 = sext i32 %8 to i64
%13 = add i32 %9, 1
br label %14
14: ; preds = %11, %14
%15 = phi i64 [ %12, %11 ], [ %37, %14 ]
%16 = phi float [ 0.000000e+00, %11 ], [ %29, %14 ]
%17 = phi float [ 0.000000e+00, %11 ], [ %33, %14 ]
%18 = phi float [ 0.000000e+00, %11 ], [ %32, %14 ]
%19 = phi <2 x float> [ zeroinitializer, %11 ], [ %36, %14 ]
%20 = getelementptr inbounds float, float* %0, i64 %15
%21 = load float, float* %20
%22 = trunc i64 %15 to i32
%23 = sub i32 %22, %1
%24 = mul nsw i32 %23, %23
%25 = insertelement <2 x i32> poison, i32 %23, i64 0
%26 = insertelement <2 x i32> %25, i32 %24, i64 1
%27 = sitofp <2 x i32> %26 to <2 x float>
%28 = extractelement <2 x float> %27, i64 1
%29 = fadd float %16, %28
%30 = mul nsw i32 %24, %24
%31 = sitofp i32 %30 to float
%32 = fadd float %18, %31
%33 = fadd float %17, %21
%34 = insertelement <2 x float> poison, float %21, i64 0
%35 = shufflevector <2 x float> %34, <2 x float> poison, <2 x i32> zeroinitializer
%36 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %27, <2 x float> %35, <2 x float> %19)
%37 = add nsw i64 %15, 1
%38 = trunc i64 %37 to i32
%39 = icmp eq i32 %13, %38
br i1 %39, label %40, label %14
40: ; preds = %14
%41 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %27, <2 x float> %35, <2 x float> %19)
br label %42
42: ; preds = %40, %3
%43 = phi float [ 0.000000e+00, %3 ], [ %32, %40 ]
%44 = phi float [ 0.000000e+00, %3 ], [ %33, %40 ]
%45 = phi float [ 0.000000e+00, %3 ], [ %29, %40 ]
%46 = phi <2 x float> [ zeroinitializer, %3 ], [ %41, %40 ]
%47 = sitofp i32 %7 to float
%48 = fneg float %45
%49 = fmul float %45, %48
%50 = tail call float @llvm.fmuladd.f32(float %47, float %43, float %49)
%51 = extractelement <2 x float> %46, i64 0
%52 = fmul float %51, %50
%53 = fmul float %44, %48
%54 = extractelement <2 x float> %46, i64 1
%55 = tail call float @llvm.fmuladd.f32(float %47, float %54, float %53)
%56 = fmul float %45, %55
%57 = sitofp i32 %1 to double
%58 = fpext float %52 to double
%59 = fmul double %58, 5.000000e-01
%60 = fpext float %56 to double
%61 = fdiv double %59, %60
%62 = fsub double %57, %61
%63 = fadd double %62, 5.000000e-01
%64 = fptosi double %63 to i32
%65 = fcmp ogt float %56, 0.000000e+00
%66 = sub nsw i32 0, %64
%67 = select i1 %65, i32 %66, i32 %64
ret i32 %67
}