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

```asm
<golden>:
  endbr64
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %rdi,%rbx
  mov    %rsi,%rbp
  pxor   %xmm7,%xmm7
  cvtsd2ss %xmm0,%xmm7
  movss  %xmm7,0x10(%rsp)
  pxor   %xmm4,%xmm4
  cvtsd2ss %xmm1,%xmm4
  movss  %xmm4,(%rsp)
  pxor   %xmm5,%xmm5
  cvtsd2ss %xmm2,%xmm5
  movss  %xmm5,0x14(%rsp)
  pxor   %xmm6,%xmm6
  cvtsd2ss %xmm3,%xmm6
  movss  %xmm6,0x18(%rsp)
  subss  %xmm4,%xmm5
  movaps %xmm4,%xmm2
  subss  %xmm7,%xmm2
  movss  D104(%rip),%xmm1
  movaps %xmm5,%xmm3
  andps  %xmm1,%xmm3
  andps  %xmm2,%xmm1
  comiss %xmm1,%xmm3
  jbe    L7
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm5,%xmm0
  mulsd  D21(%rip),%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd %xmm4,%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm7,%xmm7
  cvtsd2ss %xmm0,%xmm7
  movss  %xmm7,0x4(%rsp)
L56:
  pxor   %xmm0,%xmm0
  cvtss2sd (%rsp),%xmm0
  mov    $0x1,%eax
  call   *%rbx
  movss  %xmm0,0x8(%rsp)
  pxor   %xmm0,%xmm0
  cvtss2sd 0x4(%rsp),%xmm0
  mov    $0x1,%eax
  call   *%rbx
  movss  %xmm0,0xc(%rsp)
  pxor   %xmm3,%xmm3
  cvtss2sd 0x18(%rsp),%xmm3
  movsd  %xmm3,0x18(%rsp)
  movss  D104(%rip),%xmm2
  movss  0x14(%rsp),%xmm7
  subss  0x10(%rsp),%xmm7
  movaps %xmm7,%xmm0
  andps  %xmm2,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd %xmm0,%xmm1
  movss  (%rsp),%xmm0
  andps  %xmm2,%xmm0
  cvtss2sd %xmm0,%xmm0
  movss  0x4(%rsp),%xmm5
  andps  %xmm5,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  comisd %xmm0,%xmm1
  ja     L4
L41:
  movss  0xc(%rsp),%xmm3
  comiss 0x8(%rsp),%xmm3
  ja     L20
  movss  0x4(%rsp),%xmm4
  movss  %xmm4,(%rsp)
  movss  %xmm3,0x8(%rsp)
L20:
  movss  (%rsp),%xmm7
  movss  %xmm7,0x0(%rbp)
  movss  0x8(%rsp),%xmm0
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  ret
L7:
  movss  (%rsp),%xmm5
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm5,%xmm0
  cvtss2sd %xmm2,%xmm2
  mulsd  D21(%rip),%xmm2
  subsd  %xmm2,%xmm0
  movss  %xmm5,0x4(%rsp)
  pxor   %xmm6,%xmm6
  cvtsd2ss %xmm0,%xmm6
  movss  %xmm6,(%rsp)
  jmp    L56
L60:
  pxor   %xmm0,%xmm0
  cvtss2sd (%rsp),%xmm0
  mulsd  D93(%rip),%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x10(%rsp),%xmm1
  mulsd  D21(%rip),%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm6,%xmm6
  cvtsd2ss %xmm0,%xmm6
  movss  %xmm6,0xc(%rsp)
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm6,%xmm0
  mov    $0x1,%eax
  call   *%rbx
  movss  0x4(%rsp),%xmm4
  movss  %xmm4,0x14(%rsp)
  movss  (%rsp),%xmm7
  movss  %xmm7,0x4(%rsp)
  movss  0xc(%rsp),%xmm6
  movss  %xmm6,(%rsp)
  movss  0x8(%rsp),%xmm5
  movss  %xmm5,0xc(%rsp)
  movss  %xmm0,0x8(%rsp)
L16:
  movss  0x14(%rsp),%xmm0
  subss  0x10(%rsp),%xmm0
  andps  D104(%rip),%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd %xmm0,%xmm1
  movss  (%rsp),%xmm0
  andps  D104(%rip),%xmm0
  cvtss2sd %xmm0,%xmm0
  movss  0x4(%rsp),%xmm2
  andps  D104(%rip),%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  mulsd  0x18(%rsp),%xmm0
  comisd %xmm0,%xmm1
  jbe    L41
L4:
  movss  0x8(%rsp),%xmm5
  comiss 0xc(%rsp),%xmm5
  jbe    L60
  pxor   %xmm0,%xmm0
  cvtss2sd 0x4(%rsp),%xmm0
  mulsd  D93(%rip),%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x14(%rsp),%xmm1
  mulsd  D21(%rip),%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm3,%xmm3
  cvtsd2ss %xmm0,%xmm3
  movss  %xmm3,0x8(%rsp)
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm3,%xmm0
  mov    $0x1,%eax
  call   *%rbx
  movss  (%rsp),%xmm7
  movss  %xmm7,0x10(%rsp)
  movss  0x4(%rsp),%xmm5
  movss  %xmm5,(%rsp)
  movss  0x8(%rsp),%xmm3
  movss  %xmm3,0x4(%rsp)
  movss  0xc(%rsp),%xmm4
  movss  %xmm4,0x8(%rsp)
  movss  %xmm0,0xc(%rsp)
  jmp    L16
```

tool_calls: name: extract_symbolic_data
arguments: label: D104
datatype: dword[4]
name: extract_symbolic_data
arguments: label: D21
datatype: byte[8]
name: extract_symbolic_data
arguments: label: D93
datatype: f64

tool: {"D104": [0x7fffffff,0x00000000,0x00000000,0x00000000]}

tool: {"D21": [0x0c,0x8d,0x48,0x90,0x21,0x72,0xd8,0x3f]}

tool: D93: 0.61803399

assistant: ```c
float golden(ax, bx, cx, f, tol, xmin)
float ax, bx, cx, tol, *xmin;
float (*f)();
{
  float f0, f1, f2, f3, x0, x1, x2, x3;
  x0 = ax;
  x3 = cx;
  if (fabs(cx - bx) > fabs(bx - ax)) {
    x1 = bx;
    x2 = bx + (1.0 - 0.61803399) * (cx - bx);
  } else {
    x2 = bx;
    x1 = bx - (1.0 - 0.61803399) * (bx - ax);
  }
  f1 = (*f)(x1);
  f2 = (*f)(x2);
  while (fabs(x3 - x0) > tol * (fabs(x1) + fabs(x2))) {
    if (f2 < f1) {
      (x0) = (x1);
      (x1) = (x2);
      (x2) = (0.61803399 * x1 + (1.0 - 0.61803399) * x3);
      (f0) = (f1);
      (f1) = (f2);
      (f2) = ((*f)(x2));
    } else {
      (x3) = (x2);
      (x2) = (x1);
      (x1) = (0.61803399 * x2 + (1.0 - 0.61803399) * x0);
      (f3) = (f2);
      (f2) = (f1);
      (f1) = ((*f)(x1));
    }
  }
  if (f1 < f2) {
    *xmin = x1;
    return f1;
  } else {
    *xmin = x2;
    return f2;
  }
}
```