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

```asm
<sphere_intersect>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdi,%r15
  mov    %rsi,%r14
  mov    %rdx,%rbp
  mov    %rcx,%rbx
  mov    %r8,%r12
  movss  %xmm0,0x8(%rsp)
  mov    %r9d,0xc(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x28(%rsp)
  xor    %eax,%eax
L37:
  movss  0x0(%rbp,%rax,1),%xmm0
  subss  (%r12,%rax,1),%xmm0
  movss  %xmm0,0x1c(%rsp,%rax,1)
  add    $0x4,%rax
  cmp    $0xc,%rax
  jne    L37
  lea    0x1c(%rsp),%r13
  mov    %rbx,%rsi
  mov    %r13,%rdi
  mov    $0x0,%eax
  call   <dot@plt>
  mov    %eax,(%rsp)
  mov    %rbx,%rsi
  mov    %r13,%rdi
  mov    $0x0,%eax
  call   <dot@plt>
  mov    %eax,0x4(%rsp)
  mov    %r13,%rsi
  mov    %r13,%rdi
  mov    $0x0,%eax
  call   <dot@plt>
  mov    %eax,%edx
  mov    (%rsp),%eax
  imul   0x4(%rsp),%eax
  pxor   %xmm2,%xmm2
  cvtsi2ss %eax,%xmm2
  pxor   %xmm1,%xmm1
  cvtsi2ss %edx,%xmm1
  movss  0x8(%rsp),%xmm0
  mulss  %xmm0,%xmm0
  subss  %xmm0,%xmm1
  subss  %xmm1,%xmm2
  pxor   %xmm0,%xmm0
  movss  %xmm2,(%rsp)
  comiss %xmm2,%xmm0
  ja     L29
  cmpl   $0x0,0xc(%rsp)
  je     L9
  mov    %r13,%rdi
  mov    %rbx,%rsi
  mov    $0x0,%eax
  call   <dot@plt>
  neg    %eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  addss  (%rsp),%xmm0
  movd   %xmm0,%r13d
L31:
  pxor   %xmm0,%xmm0
  movd   %r13d,%xmm3
  comiss %xmm3,%xmm0
  jae    L36
  mov    $0x0,%eax
L27:
  movd   %r13d,%xmm0
  mulss  (%rbx,%rax,1),%xmm0
  addss  0x0(%rbp,%rax,1),%xmm0
  movss  %xmm0,(%r15,%rax,1)
  subss  (%r12,%rax,1),%xmm0
  movss  %xmm0,0x10(%rsp,%rax,1)
  add    $0x4,%rax
  cmp    $0xc,%rax
  jne    L27
  lea    0x10(%rsp),%r12
  mov    %rbx,%rsi
  mov    %r12,%rdi
  mov    $0x0,%eax
  call   <dot@plt>
  mov    %eax,%ebp
  mov    %r12,%rsi
  mov    %r12,%rdi
  mov    $0x0,%eax
  call   <dot@plt>
  pxor   %xmm1,%xmm1
  cvtsi2ss %ebp,%xmm1
  addss  %xmm1,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  divss  %xmm0,%xmm1
  mov    $0x0,%eax
  mov    %r12,%rdx
L59:
  movaps %xmm1,%xmm2
  mulss  (%rax,%rdx,1),%xmm2
  movss  (%rbx,%rax,1),%xmm0
  subss  %xmm2,%xmm0
  movss  %xmm0,(%r14,%rax,1)
  add    $0x4,%rax
  cmp    $0xc,%rax
  jne    L59
L10:
  mov    0x28(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L7
  movd   %r13d,%xmm0
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L9:
  lea    0x1c(%rsp),%rdi
  mov    %rbx,%rsi
  mov    $0x0,%eax
  call   <dot@plt>
  neg    %eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  subss  (%rsp),%xmm0
  movd   %xmm0,%r13d
  jmp    L31
L29:
  mov    D17(%rip),%r13d
  jmp    L10
L36:
  mov    D17(%rip),%r13d
  jmp    L10
L7:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D17
datatype: byte[4]

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

assistant: ```c
float sphere_intersect(float *restrict y, float *restrict r, const float *restrict s, const float *restrict d, const float *restrict c, float R, int invert) {
  int i;
  float D, n[3], t, v[3];
  for (i = 0; i < 3; ++i)
    v[i] = s[i] - c[i];
  D = ((dot(v, d)) * (dot(v, d))) - (dot(v, v) - ((R) * (R)));
  if (D < 0)
    return -1;
  if (invert)
    t = -dot(v, d) + D;
  else
    t = -dot(v, d) - D;
  if (t <= 0)
    return -1;
  for (i = 0; i < 3; ++i) {
    y[i] = s[i] + t * d[i];
    n[i] = y[i] - c[i];
  }
  float two_dot_nd_div_sq_n_mag = 2.0f * dot(n, d) / dot(n, n);
  for (i = 0; i < 3; ++i)
    r[i] = d[i] - two_dot_nd_div_sq_n_mag * n[i];
  return t;
}
```