user: Please show me the C implementation that matches the assembly code listed below:

```asm
<bcuint>:
  endbr64
  push   %r13
  movaps %xmm1,%xmm7
  movaps %xmm3,%xmm6
  mov    %rsi,%r13
  push   %r12
  subss  %xmm0,%xmm7
  subss  %xmm2,%xmm6
  mov    %rdi,%r12
  push   %rbp
  mov    %rcx,%rbp
  push   %rbx
  mov    %rdx,%rbx
  mov    $0x0,%edx
  sub    $0x28,%rsp
  ucomiss %xmm0,%xmm1
  setnp  %al
  cmovne %edx,%eax
  test   %al,%al
  jne    L19
  ucomiss %xmm2,%xmm3
  setnp  %dl
  cmove  %edx,%eax
  test   %al,%al
  jne    L19
L18:
  subss  %xmm2,%xmm5
  subss  %xmm0,%xmm4
  mov    0x20(%r12),%rcx
  movl   $0x0,(%rbx)
  movl   $0x0,0x0(%rbp)
  pxor   %xmm2,%xmm2
  mov    0x18(%r12),%rdx
  pxor   %xmm3,%xmm3
  mov    0x10(%r12),%rax
  movl   $0x0,0x0(%r13)
  pxor   %xmm8,%xmm8
  pxor   %xmm1,%xmm1
  divss  %xmm6,%xmm5
  movss  0x10(%rcx),%xmm0
  divss  %xmm7,%xmm4
  mulss  %xmm5,%xmm0
  cvtss2sd %xmm5,%xmm3
  addss  0xc(%rcx),%xmm0
  mulss  %xmm5,%xmm0
  addss  0x8(%rcx),%xmm0
  mulss  %xmm5,%xmm0
  mulss  %xmm4,%xmm2
  cvtss2sd %xmm4,%xmm1
  addss  %xmm2,%xmm0
  addss  0x4(%rcx),%xmm0
  movsd  D4(%rip),%xmm2
  movss  %xmm0,0x0(%r13)
  pxor   %xmm0,%xmm0
  cvtss2sd 0x10(%rcx),%xmm0
  mulsd  %xmm2,%xmm0
  cvtss2sd 0xc(%rcx),%xmm8
  addsd  %xmm8,%xmm8
  mulsd  %xmm3,%xmm0
  addsd  %xmm8,%xmm0
  movss  0x0(%rbp),%xmm8
  mulss  %xmm4,%xmm8
  mulsd  %xmm3,%xmm0
  cvtss2sd %xmm8,%xmm8
  addsd  %xmm8,%xmm0
  pxor   %xmm8,%xmm8
  cvtss2sd 0x8(%rcx),%xmm8
  addsd  %xmm8,%xmm0
  pxor   %xmm8,%xmm8
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,0x0(%rbp)
  pxor   %xmm0,%xmm0
  cvtss2sd 0x10(%rcx),%xmm0
  mulsd  %xmm2,%xmm0
  cvtss2sd 0x10(%rdx),%xmm8
  addsd  %xmm8,%xmm8
  mulsd  %xmm1,%xmm0
  addsd  %xmm8,%xmm0
  movss  (%rbx),%xmm8
  mulss  %xmm5,%xmm8
  mulsd  %xmm1,%xmm0
  cvtss2sd %xmm8,%xmm8
  addsd  %xmm8,%xmm0
  pxor   %xmm8,%xmm8
  cvtss2sd 0x10(%rax),%xmm8
  addsd  %xmm8,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rbx)
  movss  0x10(%rdx),%xmm0
  movss  0x0(%r13),%xmm8
  mulss  %xmm5,%xmm0
  addss  0xc(%rdx),%xmm0
  mulss  %xmm4,%xmm8
  mulss  %xmm5,%xmm0
  addss  0x8(%rdx),%xmm0
  mulss  %xmm5,%xmm0
  addss  %xmm8,%xmm0
  addss  0x4(%rdx),%xmm0
  pxor   %xmm8,%xmm8
  movss  %xmm0,0x0(%r13)
  pxor   %xmm0,%xmm0
  cvtss2sd 0x10(%rdx),%xmm0
  mulsd  %xmm2,%xmm0
  cvtss2sd 0xc(%rdx),%xmm8
  addsd  %xmm8,%xmm8
  mulsd  %xmm3,%xmm0
  addsd  %xmm8,%xmm0
  movss  0x0(%rbp),%xmm8
  mulss  %xmm4,%xmm8
  mulsd  %xmm3,%xmm0
  cvtss2sd %xmm8,%xmm8
  addsd  %xmm8,%xmm0
  pxor   %xmm8,%xmm8
  cvtss2sd 0x8(%rdx),%xmm8
  addsd  %xmm8,%xmm0
  pxor   %xmm8,%xmm8
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,0x0(%rbp)
  pxor   %xmm0,%xmm0
  cvtss2sd 0xc(%rdx),%xmm8
  cvtss2sd 0xc(%rcx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm8,%xmm8
  mov    0x8(%r12),%rsi
  mulsd  %xmm1,%xmm0
  addsd  %xmm8,%xmm0
  movss  (%rbx),%xmm8
  mulss  %xmm5,%xmm8
  mulsd  %xmm1,%xmm0
  cvtss2sd %xmm8,%xmm8
  addsd  %xmm8,%xmm0
  pxor   %xmm8,%xmm8
  cvtss2sd 0xc(%rax),%xmm8
  addsd  %xmm8,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rbx)
  movss  0x10(%rax),%xmm0
  movss  0x0(%r13),%xmm8
  mulss  %xmm5,%xmm0
  addss  0xc(%rax),%xmm0
  mulss  %xmm4,%xmm8
  mulss  %xmm5,%xmm0
  addss  0x8(%rax),%xmm0
  mulss  %xmm5,%xmm0
  addss  %xmm8,%xmm0
  addss  0x4(%rax),%xmm0
  pxor   %xmm8,%xmm8
  movss  %xmm0,0x0(%r13)
  pxor   %xmm0,%xmm0
  cvtss2sd 0x10(%rax),%xmm0
  mulsd  %xmm2,%xmm0
  cvtss2sd 0xc(%rax),%xmm8
  addsd  %xmm8,%xmm8
  mulsd  %xmm3,%xmm0
  addsd  %xmm8,%xmm0
  movss  0x0(%rbp),%xmm8
  mulss  %xmm4,%xmm8
  mulsd  %xmm3,%xmm0
  cvtss2sd %xmm8,%xmm8
  addsd  %xmm8,%xmm0
  pxor   %xmm8,%xmm8
  cvtss2sd 0x8(%rax),%xmm8
  addsd  %xmm8,%xmm0
  pxor   %xmm8,%xmm8
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,0x0(%rbp)
  pxor   %xmm0,%xmm0
  cvtss2sd 0x8(%rcx),%xmm0
  mulsd  %xmm2,%xmm0
  cvtss2sd 0x8(%rdx),%xmm8
  addsd  %xmm8,%xmm8
  mulsd  %xmm1,%xmm0
  addsd  %xmm8,%xmm0
  movss  (%rbx),%xmm8
  mulss  %xmm5,%xmm8
  mulsd  %xmm1,%xmm0
  cvtss2sd %xmm8,%xmm8
  addsd  %xmm8,%xmm0
  pxor   %xmm8,%xmm8
  cvtss2sd 0x8(%rax),%xmm8
  addsd  %xmm8,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rbx)
  movss  0x10(%rsi),%xmm0
  movss  0x0(%r13),%xmm8
  mulss  %xmm5,%xmm0
  addss  0xc(%rsi),%xmm0
  mulss  %xmm4,%xmm8
  mulss  %xmm5,%xmm0
  addss  0x8(%rsi),%xmm0
  mulss  %xmm5,%xmm0
  addss  %xmm8,%xmm0
  addss  0x4(%rsi),%xmm0
  pxor   %xmm8,%xmm8
  movss  %xmm0,0x0(%r13)
  mulss  0x0(%rbp),%xmm4
  pxor   %xmm0,%xmm0
  cvtss2sd 0x10(%rsi),%xmm0
  mulsd  %xmm2,%xmm0
  cvtss2sd 0xc(%rsi),%xmm8
  addsd  %xmm8,%xmm8
  cvtss2sd %xmm4,%xmm4
  mulsd  %xmm3,%xmm0
  addsd  %xmm8,%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm0,%xmm4
  pxor   %xmm0,%xmm0
  cvtss2sd 0x8(%rsi),%xmm0
  addsd  %xmm0,%xmm4
  pxor   %xmm0,%xmm0
  cvtsd2ss %xmm4,%xmm4
  movss  %xmm4,0x0(%rbp)
  cvtss2sd 0x4(%rcx),%xmm0
  mulsd  %xmm2,%xmm0
  mulss  (%rbx),%xmm5
  pxor   %xmm2,%xmm2
  cvtss2sd 0x4(%rdx),%xmm2
  addsd  %xmm2,%xmm2
  mulsd  %xmm1,%xmm0
  cvtss2sd %xmm5,%xmm5
  addsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x4(%rax),%xmm1
  addsd  %xmm5,%xmm0
  addsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  divss  %xmm7,%xmm0
  movss  %xmm0,(%rbx)
  movss  0x0(%rbp),%xmm0
  divss  %xmm6,%xmm0
  movss  %xmm0,0x0(%rbp)
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
  nop
L19:
  lea    D11(%rip),%rdi
  movss  %xmm5,0x1c(%rsp)
  movss  %xmm4,0x18(%rsp)
  movss  %xmm2,0x14(%rsp)
  movss  %xmm0,0x10(%rsp)
  movss  %xmm6,0xc(%rsp)
  movss  %xmm7,0x8(%rsp)
  call   <nrerror@plt>
  movss  0x1c(%rsp),%xmm5
  movss  0x18(%rsp),%xmm4
  movss  0x14(%rsp),%xmm2
  movss  0x10(%rsp),%xmm0
  movss  0xc(%rsp),%xmm6
  movss  0x8(%rsp),%xmm7
  jmp    L18
```

tool_calls: name: disassemble_data
arguments: label: D4
type_guessed: f64
name: disassemble_data
arguments: label: D11
type_guessed: string

tool: D4: 3

tool: D11: Bad input in routine bcuint

assistant: ```c
void bcuint(float **c, float x1l, float x1u, float x2l, float x2u, float x1, float x2, float *ansy, float *ansy1, float *ansy2) {
  int i;
  float t, u, d1, d2;
  d1 = x1u - x1l;
  d2 = x2u - x2l;
  if (x1u == x1l || x2u == x2l)
    nrerror("Bad input in routine bcuint");
  t = (x1 - x1l) / d1;
  u = (x2 - x2l) / d2;
  *ansy = (*ansy2) = (*ansy1) = 0.0;
  for (i = 4; i >= 1; i--) {
    *ansy = t * (*ansy) + ((c[i][4] * u + c[i][3]) * u + c[i][2]) * u + c[i][1];
    *ansy2 = t * (*ansy2) + (3.0 * c[i][4] * u + 2.0 * c[i][3]) * u + c[i][2];
    *ansy1 = u * (*ansy1) + (3.0 * c[4][i] * t + 2.0 * c[3][i]) * t + c[2][i];
  }
  *ansy1 /= d1;
  *ansy2 /= d2;
}
```