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

```asm
<xdurbin>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x1e8,%rsp
  movss  (%rdi),%xmm1
  movss  0x4(%rdi),%xmm0
  xorps  D10(%rip),%xmm0
  mov    %ecx,0x2c(%rsp)
  mov    %r8,0x38(%rsp)
  movss  D22(%rip),%xmm2
  divss  %xmm1,%xmm0
  mov    %fs:0x28,%rax
  mov    %rax,0x1d8(%rsp)
  xor    %eax,%eax
  movss  %xmm0,(%rsi)
  movss  %xmm0,(%rdx)
  mulss  %xmm0,%xmm0
  subss  %xmm0,%xmm2
  mulss  %xmm1,%xmm2
  cmp    $0x1,%ecx
  jle    L5
  mov    %rdi,0x30(%rsp)
  mov    %rdx,%rbp
  lea    0x4(%rsi),%r14
  lea    -0x8(%rdi),%rcx
  lea    0x40(%rsp),%r10
  lea    0x34(%rsp),%rax
  mov    $0x4,%r12d
  mov    $0x1,%ebx
  mov    %r10,%r9
  mov    %rax,%r10
  nopl   (%rax)
L87:
  lea    -0x1(%rbx),%r11d
  mov    %ebx,%r15d
  cmp    $0x2,%r11d
  jbe    L17
  mov    %ebx,%edx
  lea    0x10(%rbp),%r8
  pxor   %xmm1,%xmm1
  mov    %rbp,%rax
  shr    $0x2,%edx
  lea    -0x1(%rdx),%esi
  mov    %rcx,%rdx
  shl    $0x4,%rsi
  add    %r8,%rsi
  cs nopw 0x0(%rax,%rax,1)
L68:
  movups (%rdx),%xmm0
  movups (%rax),%xmm5
  add    $0x10,%rax
  sub    $0x10,%rdx
  shufps $0x1b,%xmm0,%xmm0
  mulps  %xmm5,%xmm0
  subss  %xmm0,%xmm1
  movaps %xmm0,%xmm3
  shufps $0x55,%xmm0,%xmm3
  subss  %xmm3,%xmm1
  movaps %xmm0,%xmm3
  unpckhps %xmm0,%xmm3
  shufps $0xff,%xmm0,%xmm0
  subss  %xmm3,%xmm1
  subss  %xmm0,%xmm1
  cmp    %rax,%rsi
  jne    L68
  mov    %ebx,%eax
  and    $0xfffffffc,%eax
  cmp    %ebx,%eax
  je     L107
L206:
  mov    %ebx,%edx
  mov    0x30(%rsp),%r8
  movslq %eax,%rsi
  sub    %eax,%edx
  lea    0x0(,%rsi,4),%rdi
  movslq %edx,%rdx
  movss  (%r8,%rdx,4),%xmm0
  mulss  0x0(%rbp,%rsi,4),%xmm0
  lea    0x1(%rax),%esi
  subss  %xmm0,%xmm1
  cmp    %esi,%ebx
  jle    L59
  mov    %ebx,%edx
  add    $0x2,%eax
  sub    %esi,%edx
  movslq %edx,%rdx
  movss  (%r8,%rdx,4),%xmm0
  mulss  0x4(%rbp,%rdi,1),%xmm0
  subss  %xmm0,%xmm1
  cmp    %eax,%ebx
  jle    L59
  mov    %ebx,%edx
  sub    %eax,%edx
  movslq %edx,%rax
  movss  (%r8,%rax,4),%xmm0
  mulss  0x8(%rbp,%rdi,1),%xmm0
  subss  %xmm0,%xmm1
L59:
  subss  0x10(%rcx),%xmm1
  lea    0x4(%r12),%rdx
  mov    %r9,%rdi
  mov    %rbp,%rsi
  mov    %r10,0x10(%rsp)
  mov    %r14,%r13
  add    $0x4,%r14
  mov    %rcx,0x8(%rsp)
  divss  %xmm2,%xmm1
  mov    %rdx,(%rsp)
  movss  %xmm2,0x28(%rsp)
  mov    %r11d,0x18(%rsp)
  movss  %xmm1,-0x4(%r14)
  movss  %xmm1,0x0(%rbp,%r12,1)
  movss  %xmm1,0x20(%rsp)
  call   <memcpy@plt>
  lea    (%r12,%rbp,1),%rsi
  mov    (%rsp),%rdx
  movss  0x28(%rsp),%xmm2
  cmp    %rsi,%r13
  mov    %rax,%r9
  mov    0x8(%rsp),%rcx
  mov    0x10(%rsp),%r10
  setae  %r8b
  cmp    %r14,%rbp
  mov    %rbp,%rax
  setae  %dil
  or     %r8b,%dil
  je     L155
  cmp    $0x1,%ebx
  je     L155
  mov    0x18(%rsp),%r11d
  cmp    $0x2,%r11d
  jbe    L32
  movss  0x20(%rsp),%xmm1
  lea    0x10(%rbp),%r8
L44:
  mov    %ebx,%esi
  shufps $0x0,%xmm1,%xmm1
  shr    $0x2,%esi
  lea    -0x1(%rsi),%edi
  mov    %r10,%rsi
  shl    $0x4,%rdi
  add    %r8,%rdi
  nopl   0x0(%rax)
L21:
  movups (%rsi),%xmm0
  movups (%rax),%xmm7
  add    $0x10,%rax
  sub    $0x10,%rsi
  shufps $0x1b,%xmm0,%xmm0
  mulps  %xmm1,%xmm0
  addps  %xmm7,%xmm0
  movups %xmm0,-0x10(%rax)
  cmp    %rdi,%rax
  jne    L21
  mov    %ebx,%eax
  and    $0xfffffffc,%eax
  mov    %eax,%edi
  cmp    %ebx,%eax
  je     L3
  mov    %ebx,%r15d
  sub    %eax,%r15d
  cmp    $0x1,%r15d
  je     L46
L83:
  lea    0x0(,%rax,4),%rsi
  sub    $0x8,%r12
  movss  0x0(%r13),%xmm1
  sub    %rsi,%r12
  lea    0x0(%rbp,%rsi,1),%rax
  movq   0x40(%rsp,%r12,1),%xmm0
  shufps $0xe0,%xmm1,%xmm1
  shufps $0xe1,%xmm0,%xmm0
  mulps  %xmm1,%xmm0
  movq   (%rax),%xmm1
  addps  %xmm0,%xmm1
  movlps %xmm1,(%rax)
  mov    %r15d,%eax
  and    $0xfffffffe,%eax
  add    %eax,%edi
  cmp    %eax,%r15d
  je     L3
L46:
  movslq %edi,%rax
  lea    0x0(%rbp,%rax,4),%rsi
  mov    %ebx,%eax
  sub    %edi,%eax
  cltq
  movss  0x3c(%rsp,%rax,4),%xmm0
  mulss  0x0(%r13),%xmm0
  addss  (%rsi),%xmm0
  movss  %xmm0,(%rsi)
L3:
  mov    %rdx,%r12
L41:
  movss  0x0(%r13),%xmm0
  add    $0x1,%ebx
  add    $0x4,%rcx
  add    $0x4,%r10
  movss  D22(%rip),%xmm1
  mulss  %xmm0,%xmm0
  subss  %xmm0,%xmm1
  mulss  %xmm1,%xmm2
  cmp    %ebx,0x2c(%rsp)
  jne    L87
L5:
  mov    0x38(%rsp),%rax
  movss  %xmm2,(%rax)
  mov    0x1d8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L77
  add    $0x1e8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L107:
  subss  0x10(%rcx),%xmm1
  lea    0x4(%r12),%rdx
  mov    %rbp,%rsi
  mov    %r9,%rdi
  mov    %r14,%r13
  add    $0x4,%r14
  mov    %r10,0x20(%rsp)
  mov    %r8,0x18(%rsp)
  divss  %xmm2,%xmm1
  mov    %rcx,0x10(%rsp)
  mov    %rdx,(%rsp)
  movss  %xmm2,0x8(%rsp)
  movss  %xmm1,-0x4(%r14)
  movss  %xmm1,0x0(%rbp,%r12,1)
  movss  %xmm1,0x28(%rsp)
  call   <memcpy@plt>
  cmp    %r14,%rbp
  mov    (%rsp),%rdx
  lea    (%r12,%rbp,1),%rsi
  mov    %rax,%r9
  movss  0x28(%rsp),%xmm1
  movss  0x8(%rsp),%xmm2
  mov    %rbp,%rax
  mov    0x10(%rsp),%rcx
  mov    0x18(%rsp),%r8
  mov    0x20(%rsp),%r10
  jae    L44
  cmp    %rsi,%r13
  jae    L44
L155:
  mov    %rdx,%r12
  lea    0xc(%r10),%rdx
  xchg   %ax,%ax
L102:
  movss  (%rdx),%xmm0
  mulss  0x0(%r13),%xmm0
  add    $0x4,%rax
  sub    $0x4,%rdx
  addss  -0x4(%rax),%xmm0
  movss  %xmm0,-0x4(%rax)
  cmp    %rax,%rsi
  jne    L102
  jmp    L41
  cs nopw 0x0(%rax,%rax,1)
L17:
  pxor   %xmm1,%xmm1
  xor    %eax,%eax
  jmp    L206
L32:
  xor    %edi,%edi
  xor    %eax,%eax
  jmp    L83
L77:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: data_label: D10
type_guessed: byte[16]
name: parse_data
arguments: data_label: D22
type_guessed: f32

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

tool: D22: 1

assistant: ```c
void xdurbin(r, k, a, p, ex) register int p;
register float *r, *k, *a, *ex;
{
  float bb[100];
  register int i, j;
  register float e, s, *b = bb;
  e = *r;
  *k = -r[1] / e;
  *a = *k;
  e *= (float)(1. - (*k) * (*k));
  for (i = 1; i < p; i++) {
    s = 0;
    for (j = 0; j < i; j++) {
      s -= a[j] * r[i - j];
    }
    k[i] = (s - r[i + 1]) / e;
    a[i] = k[i];
    for (j = 0; j <= i; j++) {
      b[j] = a[j];
    }
    for (j = 0; j < i; j++) {
      a[j] += k[i] * b[i - j - 1];
    }
    e *= (float)(1. - (k[i] * k[i]));
  }
  *ex = e;
}
```