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

```asm
<solve_quartic_2>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  movss  %xmm0,-0x74(%rbp)
  movss  %xmm1,-0x78(%rbp)
  movss  %xmm2,-0x7c(%rbp)
  movss  %xmm3,-0x80(%rbp)
  mov    %rdi,-0x88(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movss  -0x74(%rbp),%xmm0
  movaps %xmm0,%xmm1
  mulss  %xmm0,%xmm1
  movss  D182(%rip),%xmm0
  mulss  %xmm1,%xmm0
  movss  D184(%rip),%xmm2
  movaps %xmm0,%xmm1
  divss  %xmm2,%xmm1
  movss  -0x78(%rbp),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,-0x6c(%rbp)
  movss  -0x74(%rbp),%xmm0
  mulss  %xmm0,%xmm0
  mulss  -0x74(%rbp),%xmm0
  movss  D184(%rip),%xmm1
  divss  %xmm1,%xmm0
  movss  -0x74(%rbp),%xmm1
  mulss  -0x78(%rbp),%xmm1
  movss  D169(%rip),%xmm2
  divss  %xmm2,%xmm1
  subss  %xmm1,%xmm0
  movss  -0x7c(%rbp),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x68(%rbp)
  movss  -0x74(%rbp),%xmm0
  mulss  %xmm0,%xmm0
  mulss  -0x74(%rbp),%xmm0
  mulss  -0x74(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd %xmm0,%xmm1
  movsd  D256(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movss  -0x74(%rbp),%xmm0
  mulss  %xmm0,%xmm0
  mulss  -0x78(%rbp),%xmm0
  movss  D147(%rip),%xmm2
  divss  %xmm2,%xmm0
  cvtss2sd %xmm0,%xmm0
  addsd  %xmm1,%xmm0
  movss  -0x74(%rbp),%xmm1
  mulss  -0x7c(%rbp),%xmm1
  movss  D43(%rip),%xmm2
  divss  %xmm2,%xmm1
  pxor   %xmm2,%xmm2
  cvtss2sd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  pxor   %xmm0,%xmm0
  cvtss2sd -0x80(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,-0x64(%rbp)
  pxor   %xmm0,%xmm0
  ucomiss -0x68(%rbp),%xmm0
  jp     L17
  pxor   %xmm0,%xmm0
  ucomiss -0x68(%rbp),%xmm0
  jne    L17
  lea    -0x20(%rbp),%rdx
  movss  -0x64(%rbp),%xmm0
  mov    -0x6c(%rbp),%eax
  mov    %rdx,%rsi
  movaps %xmm0,%xmm1
  movd   %eax,%xmm0
  mov    $0x1,%edi
  call   <solve_quadratic_2@plt>
  mov    %eax,-0x44(%rbp)
  cmpl   $0x0,-0x44(%rbp)
  jne    L25
  mov    $0x0,%eax
  jmp    L41
L25:
  cmpl   $0x2,-0x44(%rbp)
  sete   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  movl   $0x0,-0x70(%rbp)
  movss  -0x20(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comiss %xmm1,%xmm0
  jb     L55
  mov    -0x20(%rbp),%eax
  movd   %eax,%xmm0
  call   <sqrt@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x40(%rbp)
  movss  -0x40(%rbp),%xmm0
  movss  D181(%rip),%xmm1
  xorps  %xmm1,%xmm0
  movss  %xmm0,-0x3c(%rbp)
  movss  -0x74(%rbp),%xmm0
  movss  D43(%rip),%xmm2
  movaps %xmm0,%xmm1
  divss  %xmm2,%xmm1
  movss  -0x40(%rbp),%xmm0
  subss  %xmm1,%xmm0
  mov    -0x88(%rbp),%rax
  movss  %xmm0,(%rax)
  movss  -0x74(%rbp),%xmm0
  movss  D43(%rip),%xmm2
  movaps %xmm0,%xmm1
  divss  %xmm2,%xmm1
  mov    -0x88(%rbp),%rax
  add    $0x4,%rax
  movss  -0x3c(%rbp),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  addl   $0x2,-0x70(%rbp)
L55:
  movss  -0x1c(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comiss %xmm1,%xmm0
  jb     L10
  mov    -0x1c(%rbp),%eax
  movd   %eax,%xmm0
  call   <sqrt@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x38(%rbp)
  movss  -0x38(%rbp),%xmm0
  movss  D181(%rip),%xmm1
  xorps  %xmm1,%xmm0
  movss  %xmm0,-0x34(%rbp)
  movss  -0x74(%rbp),%xmm0
  movss  D43(%rip),%xmm2
  movaps %xmm0,%xmm1
  divss  %xmm2,%xmm1
  mov    -0x70(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  -0x38(%rbp),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x74(%rbp),%xmm0
  movss  D43(%rip),%xmm2
  movaps %xmm0,%xmm1
  divss  %xmm2,%xmm1
  mov    -0x70(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  -0x34(%rbp),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  addl   $0x2,-0x70(%rbp)
L10:
  mov    -0x70(%rbp),%eax
  jmp    L41
L17:
  movss  -0x68(%rbp),%xmm0
  movss  D181(%rip),%xmm1
  xorps  %xmm1,%xmm0
  movaps %xmm0,%xmm3
  mulss  -0x68(%rbp),%xmm3
  movss  -0x6c(%rbp),%xmm0
  mulss  %xmm0,%xmm0
  movss  -0x64(%rbp),%xmm2
  movss  D43(%rip),%xmm1
  mulss  %xmm2,%xmm1
  subss  %xmm1,%xmm0
  movd   %xmm0,%edx
  movss  -0x6c(%rbp),%xmm0
  addss  %xmm0,%xmm0
  cvttss2si %xmm0,%eax
  lea    -0x20(%rbp),%rcx
  mov    %rcx,%rsi
  movaps %xmm3,%xmm1
  movd   %edx,%xmm0
  mov    %eax,%edi
  call   <solve_cubic@plt>
  movss  -0x20(%rbp),%xmm0
  movss  %xmm0,-0x60(%rbp)
  mov    -0x60(%rbp),%eax
  movd   %eax,%xmm0
  call   <sqrt@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x5c(%rbp)
  movss  -0x68(%rbp),%xmm0
  movss  D181(%rip),%xmm1
  xorps  %xmm1,%xmm0
  divss  -0x5c(%rbp),%xmm0
  movss  %xmm0,-0x58(%rbp)
  movss  -0x6c(%rbp),%xmm0
  addss  -0x60(%rbp),%xmm0
  addss  -0x58(%rbp),%xmm0
  movss  D169(%rip),%xmm1
  divss  %xmm1,%xmm0
  movss  %xmm0,-0x54(%rbp)
  movss  -0x6c(%rbp),%xmm0
  addss  -0x60(%rbp),%xmm0
  subss  -0x58(%rbp),%xmm0
  movss  D169(%rip),%xmm1
  divss  %xmm1,%xmm0
  movss  %xmm0,-0x50(%rbp)
  movss  D225(%rip),%xmm0
  movss  %xmm0,-0x30(%rbp)
  movss  D225(%rip),%xmm0
  movss  %xmm0,-0x2c(%rbp)
  movss  D225(%rip),%xmm0
  movss  %xmm0,-0x28(%rbp)
  movss  D225(%rip),%xmm0
  movss  %xmm0,-0x24(%rbp)
  lea    -0x30(%rbp),%rdx
  movss  -0x54(%rbp),%xmm0
  mov    -0x5c(%rbp),%eax
  mov    %rdx,%rsi
  movaps %xmm0,%xmm1
  movd   %eax,%xmm0
  mov    $0x1,%edi
  call   <solve_quadratic_2@plt>
  mov    %eax,-0x4c(%rbp)
  movss  -0x5c(%rbp),%xmm0
  movss  D181(%rip),%xmm1
  xorps  %xmm1,%xmm0
  movd   %xmm0,%eax
  lea    -0x28(%rbp),%rdx
  movss  -0x50(%rbp),%xmm0
  mov    %rdx,%rsi
  movaps %xmm0,%xmm1
  movd   %eax,%xmm0
  mov    $0x1,%edi
  call   <solve_quadratic_2@plt>
  mov    %eax,-0x48(%rbp)
  movss  -0x30(%rbp),%xmm0
  movss  -0x74(%rbp),%xmm1
  movss  D43(%rip),%xmm2
  divss  %xmm2,%xmm1
  subss  %xmm1,%xmm0
  mov    -0x88(%rbp),%rax
  movss  %xmm0,(%rax)
  movss  -0x2c(%rbp),%xmm0
  movss  -0x74(%rbp),%xmm1
  movss  D43(%rip),%xmm2
  divss  %xmm2,%xmm1
  mov    -0x88(%rbp),%rax
  add    $0x4,%rax
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x28(%rbp),%xmm0
  movss  -0x74(%rbp),%xmm1
  movss  D43(%rip),%xmm2
  divss  %xmm2,%xmm1
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x24(%rbp),%xmm0
  movss  -0x74(%rbp),%xmm1
  movss  D43(%rip),%xmm2
  divss  %xmm2,%xmm1
  mov    -0x4c(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x4c(%rbp),%edx
  mov    -0x48(%rbp),%eax
  add    %edx,%eax
L41:
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L13
  call   <__stack_chk_fail@plt>
L13:
  leave
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D182
datatype: f32
name: unpack_data
arguments: data_label: D184
datatype: f32
name: unpack_data
arguments: data_label: D169
datatype: f32
name: unpack_data
arguments: data_label: D256
datatype: f32[2]
name: unpack_data
arguments: data_label: D147
datatype: f32
name: unpack_data
arguments: data_label: D43
datatype: f32
name: unpack_data
arguments: data_label: D181
datatype: byte[4]
name: unpack_data
arguments: data_label: D225
datatype: byte[4]

tool: D182: 3

tool: D184: 8

tool: D169: 2

tool: {"D256": [0.0,-1.0625]}

tool: D147: 16

tool: D43: 4

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

tool: {"D225": [0x00,0x00,0x80,0xbf]}

assistant: ```c
int solve_quartic_2(float p, float q, float r, float s, float t[4]) {
  float a = q - p * p * 3 / 8;
  float b = p * p * p / 8 - p * q / 2 + r;
  float c = p * p * p * p * (-3.0 / 256.0) + p * p * q / 16 - p * r / 4 + s;
  if (b == 0.0) {
    float z[2];
    int count = solve_quadratic_2(1, a, c, z);
    if (count == 0)
      return 0;
    assert(count == 2);
    int solutions = 0;
    if (z[0] >= 0) {
      float x1 = sqrt(z[0]);
      float x2 = -x1;
      t[0] = x1 - p / 4;
      t[1] = x2 - p / 4;
      solutions += 2;
    }
    if (z[1] >= 0) {
      float x1 = sqrt(z[1]);
      float x2 = -x1;
      t[solutions + 0] = x1 - p / 4;
      t[solutions + 1] = x2 - p / 4;
      solutions += 2;
    }
    return solutions;
  }
  float t3[4];
  solve_cubic(2 * a, a * a - 4 * c, -b * b, t3);
  float z = t3[0];
  float alpha = sqrt(z);
  float rho = -b / alpha;
  float beta = (a + z + rho) / 2;
  float delta = (a + z - rho) / 2;
  float t1[2] = {-1, -1}, t2[2] = {-1, -1};
  int n1 = solve_quadratic_2(1, alpha, beta, t1);
  int n2 = solve_quadratic_2(1, -alpha, delta, t2);
  t[0] = t1[0] - p / 4;
  t[1] = t1[1] - p / 4;
  t[n1 + 0] = t2[0] - p / 4;
  t[n1 + 1] = t2[1] - p / 4;
  return n1 + n2;
}
```