user: What is the c source code of the assembly code below:

```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  D32(%rip),%xmm0
  movsd  0x8(%rsp),%xmm5
  cvtsi2sd %eax,%xmm3
  movsd  0x10(%rsp),%xmm2
  comisd %xmm0,%xmm3
  pxor   %xmm3,%xmm3
  jae    L291
  ucomisd %xmm3,%xmm5
  jp     L236
  je     L291
L236:
  ucomisd %xmm3,%xmm2
  jp     L204
  je     L100
L204:
  ucomisd %xmm3,%xmm1
  jp     L133
  jne    L133
  movq   D74(%rip),%xmm1
  xorpd  %xmm1,%xmm2
  divsd  %xmm5,%xmm2
  comisd %xmm2,%xmm3
  jbe    L155
L131:
  movsd  D63(%rip),%xmm1
  movapd %xmm1,%xmm4
  jmp    L7
  xchg   %ax,%ax
L291:
  ucomisd %xmm3,%xmm1
  jp     L51
  jne    L51
  ucomisd %xmm3,%xmm2
  jp     L131
  jne    L131
L267:
  movapd %xmm3,%xmm1
  movapd %xmm3,%xmm4
L7:
  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     L77
  movapd %xmm3,%xmm4
L77:
  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    L255
  comisd %xmm3,%xmm4
  movsd  0x10(%rsp),%xmm1
  jbe    L270
  comisd %xmm3,%xmm1
  jbe    L270
  add    $0x20,%rsp
  movapd %xmm4,%xmm0
  pop    %rbx
  jmp    10a0 <Min@plt>
  nopl   0x0(%rax)
L51:
  xorpd  D74(%rip),%xmm2
  movapd %xmm2,%xmm4
  divsd  %xmm1,%xmm4
  movapd %xmm4,%xmm1
  jmp    L7
  nopl   0x0(%rax)
L255:
  movapd %xmm3,%xmm1
L270:
  add    $0x20,%rsp
  movapd %xmm4,%xmm0
  pop    %rbx
  jmp    10b0 <Max@plt>
  nopw   0x0(%rax,%rax,1)
L100:
  movapd %xmm1,%xmm4
  xorpd  D74(%rip),%xmm4
  movapd %xmm3,%xmm1
  divsd  %xmm5,%xmm4
  jmp    L7
  nopl   0x0(%rax)
L133:
  movsd  D51(%rip),%xmm4
  movapd %xmm1,%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  %xmm5,%xmm4
  mulsd  %xmm2,%xmm4
  subsd  %xmm4,%xmm0
  comisd %xmm0,%xmm3
  ja     L131
  ucomisd %xmm0,%xmm3
  ja     L251
  sqrtsd %xmm0,%xmm0
L31:
  ucomisd %xmm3,%xmm0
  jp     L83
  je     L257
L83:
  comisd %xmm3,%xmm1
  ja     L152
  comisd %xmm1,%xmm3
  jbe    L267
  subsd  %xmm1,%xmm0
  addsd  %xmm2,%xmm2
  addsd  %xmm5,%xmm5
  divsd  %xmm0,%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm5,%xmm1
  movapd %xmm2,%xmm4
  jmp    L7
  nopl   (%rax)
L155:
  ucomisd %xmm2,%xmm3
  ja     L180
  sqrtsd %xmm2,%xmm2
  movapd %xmm2,%xmm4
L130:
  xorpd  %xmm4,%xmm1
  jmp    L7
  nopl   0x0(%rax,%rax,1)
L257:
  addsd  %xmm5,%xmm5
  movapd %xmm1,%xmm4
  xorpd  D74(%rip),%xmm4
  divsd  %xmm5,%xmm4
  movapd %xmm4,%xmm1
  jmp    L7
  nopl   (%rax)
L152:
  xorpd  D74(%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    L7
L251:
  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    L31
L180:
  movapd %xmm2,%xmm0
  movsd  %xmm3,(%rsp)
  call   <sqrt@plt>
  movq   D74(%rip),%xmm1
  movsd  (%rsp),%xmm3
  movapd %xmm0,%xmm4
  jmp    L130
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D32
datatype: f64
name: extract_symbolic_data
arguments: data_label: D74
datatype: f32[4]
name: extract_symbolic_data
arguments: data_label: D63
datatype: f32[2]
name: extract_symbolic_data
arguments: data_label: D51
datatype: f64

tool: D32: 1000000000000

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

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

tool: D51: 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;
}
```