user: What would the C source look like for the assembly code provided here?

```asm
<SolveQuadraticEq>:
  endbr64
  push   %rbx
  sub    $0x20,%rsp
  movsd  %xmm2,0x10(%rsp)
  movsd  %xmm1,(%rsp)
  movsd  %xmm0,0x8(%rsp)
  call   <fabs@plt>
  movsd  (%rsp),%xmm1
  mov    %eax,%ebx
  movapd %xmm1,%xmm0
  call   <fabs@plt>
  pxor   %xmm0,%xmm0
  pxor   %xmm3,%xmm3
  movsd  (%rsp),%xmm1
  cvtsi2sd %ebx,%xmm0
  mulsd  D21(%rip),%xmm0
  movsd  0x8(%rsp),%xmm5
  cvtsi2sd %eax,%xmm3
  movsd  0x10(%rsp),%xmm2
  comisd %xmm0,%xmm3
  pxor   %xmm3,%xmm3
  jae    L266
  ucomisd %xmm3,%xmm5
  jp     L62
  je     L266
L62:
  ucomisd %xmm3,%xmm2
  jp     L228
  je     L22
L228:
  ucomisd %xmm3,%xmm1
  jp     L43
  jne    L43
  movq   D65(%rip),%xmm1
  xorpd  %xmm1,%xmm2
  divsd  %xmm5,%xmm2
  comisd %xmm2,%xmm3
  jbe    L3
L245:
  movsd  D25(%rip),%xmm1
  movapd %xmm1,%xmm4
  jmp    L50
  xchg   %ax,%ax
L266:
  ucomisd %xmm3,%xmm1
  jp     L304
  jne    L304
  ucomisd %xmm3,%xmm2
  jp     L245
  jne    L245
L254:
  movapd %xmm3,%xmm1
  movapd %xmm3,%xmm4
L50:
  movapd %xmm4,%xmm0
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm4,(%rsp)
  call   <fabs@plt>
  movsd  (%rsp),%xmm4
  movsd  0x8(%rsp),%xmm1
  test   %eax,%eax
  movsd  0x10(%rsp),%xmm3
  jg     L172
  movapd %xmm3,%xmm4
L172:
  movapd %xmm1,%xmm0
  movsd  %xmm3,0x8(%rsp)
  movsd  %xmm4,(%rsp)
  movsd  %xmm1,0x10(%rsp)
  call   <fabs@plt>
  movsd  (%rsp),%xmm4
  movsd  0x8(%rsp),%xmm3
  test   %eax,%eax
  jle    L168
  comisd %xmm3,%xmm4
  movsd  0x10(%rsp),%xmm1
  jbe    L105
  comisd %xmm3,%xmm1
  jbe    L105
  add    $0x20,%rsp
  movapd %xmm4,%xmm0
  pop    %rbx
  jmp    10a0 <Min@plt>
  nopl   0x0(%rax)
L304:
  xorpd  D65(%rip),%xmm2
  movapd %xmm2,%xmm4
  divsd  %xmm1,%xmm4
  movapd %xmm4,%xmm1
  jmp    L50
  nopl   0x0(%rax)
L168:
  movapd %xmm3,%xmm1
L105:
  add    $0x20,%rsp
  movapd %xmm4,%xmm0
  pop    %rbx
  jmp    10b0 <Max@plt>
  nopw   0x0(%rax,%rax,1)
L22:
  movapd %xmm1,%xmm4
  xorpd  D65(%rip),%xmm4
  movapd %xmm3,%xmm1
  divsd  %xmm5,%xmm4
  jmp    L50
  nopl   0x0(%rax)
L43:
  movsd  D12(%rip),%xmm4
  movapd %xmm1,%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  %xmm5,%xmm4
  mulsd  %xmm2,%xmm4
  subsd  %xmm4,%xmm0
  comisd %xmm0,%xmm3
  ja     L245
  ucomisd %xmm0,%xmm3
  ja     L42
  sqrtsd %xmm0,%xmm0
L273:
  ucomisd %xmm3,%xmm0
  jp     L208
  je     L191
L208:
  comisd %xmm3,%xmm1
  ja     L9
  comisd %xmm1,%xmm3
  jbe    L254
  subsd  %xmm1,%xmm0
  addsd  %xmm2,%xmm2
  addsd  %xmm5,%xmm5
  divsd  %xmm0,%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm5,%xmm1
  movapd %xmm2,%xmm4
  jmp    L50
  nopl   (%rax)
L3:
  ucomisd %xmm2,%xmm3
  ja     L85
  sqrtsd %xmm2,%xmm2
  movapd %xmm2,%xmm4
L206:
  xorpd  %xmm4,%xmm1
  jmp    L50
  nopl   0x0(%rax,%rax,1)
L191:
  addsd  %xmm5,%xmm5
  movapd %xmm1,%xmm4
  xorpd  D65(%rip),%xmm4
  divsd  %xmm5,%xmm4
  movapd %xmm4,%xmm1
  jmp    L50
  nopl   (%rax)
L9:
  xorpd  D65(%rip),%xmm1
  addsd  %xmm2,%xmm2
  addsd  %xmm5,%xmm5
  movapd %xmm1,%xmm6
  subsd  %xmm0,%xmm6
  divsd  %xmm6,%xmm2
  movapd %xmm6,%xmm4
  divsd  %xmm5,%xmm4
  movapd %xmm2,%xmm1
  jmp    L50
L42:
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm2,0x10(%rsp)
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm5,(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm2
  movsd  0x8(%rsp),%xmm1
  movsd  (%rsp),%xmm5
  jmp    L273
L85:
  movapd %xmm2,%xmm0
  movsd  %xmm3,(%rsp)
  call   <sqrt@plt>
  movq   D65(%rip),%xmm1
  movsd  (%rsp),%xmm3
  movapd %xmm0,%xmm4
  jmp    L206
```

tool_calls: name: extract_symbolic_data
arguments: label: D21
datatype: f64
name: extract_symbolic_data
arguments: label: D65
datatype: f32[4]
name: extract_symbolic_data
arguments: label: D25
datatype: f32[2]
name: extract_symbolic_data
arguments: label: D12
datatype: f64

tool: D21: 1000000000000

tool: {"D65": [0.0,-0.0,0.0,0.0]}

tool: {"D25": [0.0,-6.1103515625]}

tool: D12: 4

assistant: ```c
double SolveQuadraticEq(double AA, double BB, double CC) {
  double X1 = 0.0, X2 = 0.0;
  double DD, TEMP;
  double Time = 0.0;
  if ((fabs(AA) * 1.0e12) <= fabs(BB) || (AA == 0.0)) {
    if (BB != 0.0) {
      X1 = -CC / BB;
      X2 = X1;
      goto ok15;
    } else {
      if (CC == 0.0) {
        X1 = 0.0;
        X2 = 0.0;
        goto ok15;
      }
      X1 = -1.0e4;
      X2 = X1;
      goto ok15;
    }
  }
  if (CC == 0) {
    X1 = -BB / AA;
    X2 = 0.0;
    goto ok15;
  }
  if (BB == 0.0) {
    TEMP = -CC / AA;
    if (TEMP < 0.0) {
      X1 = -1.0e4;
      X2 = -1.0e4;
      goto ok15;
    }
    X1 = sqrt(TEMP);
    X2 = -X1;
    goto ok15;
  }
  DD = BB * BB - 4.0 * AA * CC;
  if (DD < 0.0) {
    X1 = -1.0e4;
    X2 = -1.0e4;
    goto ok15;
  }
  DD = sqrt(DD);
  if (DD == 0.0) {
    X1 = -BB / (2.0 * AA);
    X2 = X1;
    goto ok15;
  }
  if (BB > 0.0) {
    X1 = (-BB - DD) / (2.0 * AA);
    X2 = (2.0 * CC) / (-BB - DD);
    goto ok15;
  }
  if (BB < 0.0) {
    X2 = (-BB + DD) / (2.0 * AA);
    X1 = (2.0 * CC) / (-BB + DD);
    goto ok15;
  }
ok15:
  if (fabs(X1) < 4e-10)
    X1 = 0.0;
  if (fabs(X2) < 4e-10)
    X2 = 0.0;
  if ((X1 > 0.0) && (X2 > 0.0)) {
    Time = Min(X1, X2);
  } else {
    Time = Max(X1, X2);
  }
  return Time;
}
```