user: 请从给定的汇编指令推导出C源代码:

```asm
<student_t_p_2tail>:
  endbr64
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  movss  %xmm0,0x4(%rsp)
  lea    0x14(%rsp),%rsi
  movaps %xmm1,%xmm0
  lea    0x10(%rsp),%rdi
  mov    %fs:0x28,%rax
  mov    %rax,0x18(%rsp)
  xor    %eax,%eax
  call   <df_to_index@plt>
  movss  0x4(%rsp),%xmm4
  pxor   %xmm7,%xmm7
  comiss %xmm4,%xmm7
  jbe    L107
  xorps  D9(%rip),%xmm4
  movss  %xmm4,0x4(%rsp)
L107:
  movss  0x14(%rsp),%xmm0
  mov    0x10(%rsp),%esi
  xor    %edi,%edi
  call   <get_t@plt>
  movss  0x4(%rsp),%xmm6
  movaps %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  comiss %xmm1,%xmm6
  jbe    L1
L47:
  mov    0x18(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L54
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
  nopl   0x0(%rax)
L1:
  movss  0x14(%rsp),%xmm0
  mov    0x10(%rsp),%esi
  mov    $0xbd,%edi
  call   <get_t@plt>
  movaps %xmm0,%xmm1
  comiss 0x4(%rsp),%xmm1
  movss  D1(%rip),%xmm0
  ja     L47
  mov    $0xbd,%ebp
  xor    %r12d,%r12d
  jmp    L66
  nopl   0x0(%rax,%rax,1)
L28:
  ucomiss %xmm0,%xmm3
  cmovae %ebx,%ebp
L111:
  mov    %ebp,%eax
  sub    %r12d,%eax
  cmp    $0x1,%eax
  jle    L26
L66:
  lea    (%r12,%rbp,1),%eax
  movss  0x14(%rsp),%xmm0
  mov    0x10(%rsp),%esi
  mov    %eax,%ebx
  shr    $0x1f,%ebx
  add    %eax,%ebx
  sar    %ebx
  mov    %ebx,%edi
  mov    %ebx,%r13d
  call   <get_t@plt>
  movss  0x4(%rsp),%xmm3
  comiss %xmm3,%xmm0
  jb     L28
  movss  0x4(%rsp),%xmm5
  ucomiss %xmm0,%xmm5
  jp     L9
  je     L97
L9:
  mov    %ebx,%r12d
  jmp    L111
  nopw   0x0(%rax,%rax,1)
L97:
  mov    %r13d,%ebp
L79:
  mov    0x10(%rsp),%esi
  movss  0x14(%rsp),%xmm0
  mov    %r13d,%edi
  call   <get_t@plt>
  mov    0x10(%rsp),%esi
  mov    %ebp,%edi
  movss  %xmm0,0x8(%rsp)
  movss  0x14(%rsp),%xmm0
  call   <get_t@plt>
  movss  0x4(%rsp),%xmm7
  movss  %xmm0,0xc(%rsp)
  comiss %xmm0,%xmm7
  setae  %dil
  movzbl %dil,%edi
  call   <assert@plt>
  movss  0x8(%rsp),%xmm1
  xor    %edi,%edi
  comiss 0x4(%rsp),%xmm1
  setae  %dil
  call   <assert@plt>
  mov    %ebp,%edi
  call   <index_to_p@plt>
  mov    %r13d,%edi
  mov    %eax,%ebp
  call   <index_to_p@plt>
  movss  0xc(%rsp),%xmm2
  movss  0x8(%rsp),%xmm1
  mov    %ebp,%esi
  movss  0x4(%rsp),%xmm0
  mov    %eax,%edi
  call   <linear@plt>
  jmp    L47
L54:
  call   <__stack_chk_fail@plt>
L26:
  mov    %r12d,%r13d
  jmp    L79
```

tool_calls: name: read_data
arguments: data_label: D9
type_infered: byte[16]
name: read_data
arguments: data_label: D1
type_infered: f32

tool: {"D9": [0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: D1: 1

assistant: ```c
float student_t_p_2tail(float t, float df) {
  int df_index1, p_index_min, p_index_mid, p_index_max;
  float df_weight1, t_max, t_mid, t_min;
  p_index_min = 0;
  p_index_max = 189;
  df_to_index(df, &df_index1, &df_weight1);
  if (t < 0)
    t = -t;
  if (t > get_t(p_index_min, df_index1, df_weight1))
    return 0.000;
  if (t < get_t(p_index_max, df_index1, df_weight1))
    return 1.000;
  while (p_index_max - p_index_min > 1) {
    p_index_mid = (p_index_min + p_index_max) / 2;
    t_mid = get_t(p_index_mid, df_index1, df_weight1);
    if (t_mid >= t)
      p_index_min = p_index_mid;
    if (t_mid <= t)
      p_index_max = p_index_mid;
  }
  t_max = get_t(p_index_min, df_index1, df_weight1);
  t_min = get_t(p_index_max, df_index1, df_weight1);
  assert(t >= t_min);
  assert(t <= t_max);
  return linear(t, t_max, index_to_p(p_index_min), t_min, index_to_p(p_index_max));
}
```