void chsone(float bins[], float ebins[], int nbins, int knstrn, float *df,
            float *chsq, float *prob) {
  float gammq(float a, float x);
  void nrerror(char error_text[]);
  int j;
  float temp;
  *df = nbins - knstrn;
  *chsq = 0.0;
  for (j = 1; j <= nbins; j++) {
    if (ebins[j] <= 0.0)
      nrerror("Bad expected number in chsone");
    temp = bins[j] - ebins[j];
    *chsq += temp * temp / ebins[j];
  }
  *prob = gammq(0.5 * (*df), 0.5 * (*chsq));
}


x86 assembly (clang, -O0):
<chsone>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x18(%rbp),%eax
cvtsi2ss %eax,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     11ec <chsone+0xcc>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jb     119c <chsone+0x7c>
lea    0xe71(%rip),%rdi
call   1040 <nrerror@plt>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
mulss  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
mov    -0x28(%rbp),%rax
addss  (%rax),%xmm0
movss  %xmm0,(%rax)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1166 <chsone+0x46>
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe00(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x28(%rbp),%rax
movss  (%rax),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  0xde4(%rip),%xmm1
mulsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm1
call   1030 <gammq@plt>
mov    0x10(%rbp),%rax
movss  %xmm0,(%rax)
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<chsone>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %edx,%eax
sub    %ecx,%eax
cvtsi2ss %eax,%xmm0
mov    0x40(%rsp),%rbp
movss  %xmm0,(%r8)
movl   $0x0,(%r9)
test   %edx,%edx
jle    11e0 <chsone+0xc0>
mov    %r9,%r12
mov    %edx,%r15d
mov    %rsi,%rbx
mov    %rdi,%r13
mov    %r8,(%rsp)
add    $0x1,%r15d
mov    $0x1,%r14d
xorps  %xmm3,%xmm3
lea    0xe94(%rip),%rbp
xorps  %xmm1,%xmm1
jmp    11a9 <chsone+0x89>
cs nopw 0x0(%rax,%rax,1)
nop
movaps %xmm1,%xmm2
movss  0x0(%r13,%r14,4),%xmm1
subss  %xmm0,%xmm1
mulss  %xmm1,%xmm1
divss  %xmm0,%xmm1
addss  %xmm2,%xmm1
movss  %xmm1,(%r12)
add    $0x1,%r14
cmp    %r14,%r15
je     11cd <chsone+0xad>
movss  (%rbx,%r14,4),%xmm0
ucomiss %xmm0,%xmm3
jb     1180 <chsone+0x60>
mov    %rbp,%rdi
call   1040 <nrerror@plt>
xorps  %xmm3,%xmm3
movss  (%rbx,%r14,4),%xmm0
movss  (%r12),%xmm1
jmp    1180 <chsone+0x60>
mulss  0xe2b(%rip),%xmm1
mov    0x40(%rsp),%rbp
mov    (%rsp),%r8
jmp    11e3 <chsone+0xc3>
xorps  %xmm1,%xmm1
movss  (%r8),%xmm0
mulss  0xe10(%rip),%xmm0
call   1030 <gammq@plt>
movss  %xmm0,0x0(%rbp)
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [30 x i8] c"Bad expected number in chsone\00"
define void @chsone(float* %0, float* %1, i32 %2, i32 %3, float* %4, float* %5, float* %6) {
%8 = alloca float*
%9 = alloca float*
%10 = alloca i32
%11 = alloca i32
%12 = alloca float*
%13 = alloca float*
%14 = alloca float*
%15 = alloca i32
%16 = alloca float
store float* %0, float** %8
store float* %1, float** %9
store i32 %2, i32* %10
store i32 %3, i32* %11
store float* %4, float** %12
store float* %5, float** %13
store float* %6, float** %14
%17 = load i32, i32* %10
%18 = load i32, i32* %11
%19 = sub nsw i32 %17, %18
%20 = sitofp i32 %19 to float
%21 = load float*, float** %12
store float %20, float* %21
%22 = load float*, float** %13
store float 0.000000e+00, float* %22
store i32 1, i32* %15
br label %23
23: ; preds = %60, %7
%24 = load i32, i32* %15
%25 = load i32, i32* %10
%26 = icmp sle i32 %24, %25
br i1 %26, label %27, label %63
27: ; preds = %23
%28 = load float*, float** %9
%29 = load i32, i32* %15
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float, float* %28, i64 %30
%32 = load float, float* %31
%33 = fpext float %32 to double
%34 = fcmp ole double %33, 0.000000e+00
br i1 %34, label %35, label %36
35: ; preds = %27
br label %36
36: ; preds = %35, %27
%37 = load float*, float** %8
%38 = load i32, i32* %15
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds float, float* %37, i64 %39
%41 = load float, float* %40
%42 = load float*, float** %9
%43 = load i32, i32* %15
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds float, float* %42, i64 %44
%46 = load float, float* %45
%47 = fsub float %41, %46
store float %47, float* %16
%48 = load float, float* %16
%49 = load float, float* %16
%50 = fmul float %48, %49
%51 = load float*, float** %9
%52 = load i32, i32* %15
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds float, float* %51, i64 %53
%55 = load float, float* %54
%56 = fdiv float %50, %55
%57 = load float*, float** %13
%58 = load float, float* %57
%59 = fadd float %58, %56
store float %59, float* %57
br label %60
60: ; preds = %36
%61 = load i32, i32* %15
%62 = add nsw i32 %61, 1
store i32 %62, i32* %15
br label %23
63: ; preds = %23
%64 = load float*, float** %12
%65 = load float, float* %64
%66 = fpext float %65 to double
%67 = fmul double 5.000000e-01, %66
%68 = fptrunc double %67 to float
%69 = load float*, float** %13
%70 = load float, float* %69
%71 = fpext float %70 to double
%72 = fmul double 5.000000e-01, %71
%73 = fptrunc double %72 to float
%74 = call float @gammq(float %68, float %73)
%75 = load float*, float** %14
store float %74, float* %75
ret void
}

LLVM IR (-O2):
@.str = [30 x i8] c"Bad expected number in chsone\00"
define void @chsone(float* nocapture %0, float* nocapture %1, i32 %2, i32 %3, float* nocapture %4, float* nocapture %5, float* nocapture writeonly %6) {
%8 = sub nsw i32 %2, %3
%9 = sitofp i32 %8 to float
store float %9, float* %4
store float 0.000000e+00, float* %5
%10 = icmp slt i32 %2, 1
br i1 %10, label %36, label %11
11: ; preds = %7
%12 = add nuw i32 %2, 1
%13 = zext i32 %12 to i64
br label %14
14: ; preds = %11, %23
%15 = phi float [ 0.000000e+00, %11 ], [ %31, %23 ]
%16 = phi i64 [ 1, %11 ], [ %32, %23 ]
%17 = getelementptr inbounds float, float* %1, i64 %16
%18 = load float, float* %17
%19 = fcmp ugt float %18, 0.000000e+00
br i1 %19, label %23, label %20
20: ; preds = %14
tail call void @nrerror(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str, i64 0, i64 0))
%21 = load float, float* %17
%22 = load float, float* %5
br label %23
23: ; preds = %20, %14
%24 = phi float [ %22, %20 ], [ %15, %14 ]
%25 = phi float [ %21, %20 ], [ %18, %14 ]
%26 = getelementptr inbounds float, float* %0, i64 %16
%27 = load float, float* %26
%28 = fsub float %27, %25
%29 = fmul float %28, %28
%30 = fdiv float %29, %25
%31 = fadd float %24, %30
store float %31, float* %5
%32 = add nuw nsw i64 %16, 1
%33 = icmp eq i64 %32, %13
br i1 %33, label %34, label %14
34: ; preds = %23
%35 = fmul float %31, 5.000000e-01
br label %36
36: ; preds = %34, %7
%37 = phi float [ %35, %34 ], [ 0.000000e+00, %7 ]
%38 = load float, float* %4
%39 = fmul float %38, 5.000000e-01
%40 = tail call float @gammq(float %39, float %37)
store float %40, float* %6
ret void
}