user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<vec4_slerp_to>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x40,%rsp
  mov    %rdi,-0x28(%rbp)
  mov    %rsi,-0x30(%rbp)
  mov    %rdx,-0x38(%rbp)
  movss  %xmm0,-0x3c(%rbp)
  movss  D30(%rip),%xmm0
  movss  %xmm0,-0x14(%rbp)
  mov    -0x30(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm1
  mov    -0x38(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x30(%rbp),%rax
  movss  (%rax),%xmm2
  mov    -0x38(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x30(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x30(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x10(%rbp)
  pxor   %xmm0,%xmm0
  comiss -0x10(%rbp),%xmm0
  jbe    L3
  movss  D42(%rip),%xmm0
  jmp    L40
L3:
  movss  D30(%rip),%xmm0
L40:
  movss  %xmm0,-0xc(%rbp)
  movss  -0xc(%rbp),%xmm0
  mulss  -0x10(%rbp),%xmm0
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  call   <acosf@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x8(%rbp)
  mov    -0x8(%rbp),%eax
  movd   %eax,%xmm0
  call   <sinf@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x4(%rbp)
  mov    -0x8(%rbp),%eax
  movd   %eax,%xmm0
  call   <fabs@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x14(%rbp)
  movss  -0x14(%rbp),%xmm0
  comiss -0x3c(%rbp),%xmm0
  jbe    L63
  movss  -0x3c(%rbp),%xmm0
  divss  -0x14(%rbp),%xmm0
  movss  %xmm0,-0x14(%rbp)
  jmp    L0
L63:
  movss  D30(%rip),%xmm0
  movss  %xmm0,-0x14(%rbp)
L0:
  mov    -0x4(%rbp),%eax
  movd   %eax,%xmm0
  call   <fabs@plt>
  movd   %xmm0,%eax
  movd   %eax,%xmm3
  comiss D89(%rip),%xmm3
  jbe    L32
  movss  -0x14(%rbp),%xmm0
  comiss D20(%rip),%xmm0
  jbe    L32
  movss  D30(%rip),%xmm0
  comiss -0x14(%rbp),%xmm0
  jbe    L32
  movss  D30(%rip),%xmm0
  subss  -0x14(%rbp),%xmm0
  mulss  -0x8(%rbp),%xmm0
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  call   <sinf@plt>
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  divss  -0x4(%rbp),%xmm0
  movss  %xmm0,-0x1c(%rbp)
  movss  -0x8(%rbp),%xmm0
  mulss  -0x14(%rbp),%xmm0
  mulss  -0xc(%rbp),%xmm0
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  call   <sinf@plt>
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  divss  -0x4(%rbp),%xmm0
  movss  %xmm0,-0x18(%rbp)
  jmp    L76
L32:
  movss  D30(%rip),%xmm0
  subss  -0x14(%rbp),%xmm0
  movss  %xmm0,-0x1c(%rbp)
  movss  -0x14(%rbp),%xmm0
  movss  %xmm0,-0x18(%rbp)
L76:
  mov    -0x30(%rbp),%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x1c(%rbp),%xmm1
  mov    -0x38(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  -0x18(%rbp),%xmm0
  addss  %xmm1,%xmm0
  mov    -0x28(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x1c(%rbp),%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  -0x18(%rbp),%xmm0
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x1c(%rbp),%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  -0x18(%rbp),%xmm0
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x1c(%rbp),%xmm1
  mov    -0x38(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  -0x18(%rbp),%xmm0
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  movd   %xmm1,%eax
  movd   %eax,%xmm0
  call   <sqrtf@plt>
  movaps %xmm0,%xmm1
  movss  D30(%rip),%xmm0
  divss  %xmm1,%xmm0
  movss  %xmm0,-0x8(%rbp)
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  -0x8(%rbp),%xmm0
  mov    -0x28(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  mulss  -0x8(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  mulss  -0x8(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  mulss  -0x8(%rbp),%xmm0
  movss  %xmm0,(%rax)
  movss  -0x14(%rbp),%xmm0
  leave
  ret
```

tool_calls: name: unpack_data
arguments: tag: D30
type: byte[4]
name: unpack_data
arguments: tag: D42
type: byte[4]
name: unpack_data
arguments: tag: D89
type: byte[4]
name: unpack_data
arguments: tag: D20
type: byte[4]

tool: {"D30": [0x00,0x00,0x80,0x3f]}

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

tool: {"D89": [0xac,0xc5,0x27,0x37]}

tool: {"D20": [0x17,0xb7,0xd1,0x38]}

assistant: ```c
float vec4_slerp_to(float ret[4], const float q1[4], const float q2[4], float max_step_rad) {
  float cos_fi, sin_fi, fi, k1, k2, sign;
  float t = 1.0f;
  cos_fi = q1[3] * q2[3] + q1[0] * q2[0] + q1[1] * q2[1] + q1[2] * q2[2];
  sign = (cos_fi < 0.0f) ? (-1.0f) : (1.0f);
  fi = acosf(sign * cos_fi);
  sin_fi = sinf(fi);
  t = fabs(fi);
  if (t > max_step_rad) {
    t = max_step_rad / t;
  } else {
    t = 1.0f;
  }
  if ((fabs(sin_fi) > 0.00001f) && (t > 0.0001f) && (t < 1.0f)) {
    k1 = sinf(fi * (1.0f - t)) / sin_fi;
    k2 = sinf(fi * t * sign) / sin_fi;
  } else {
    k1 = 1.0f - t;
    k2 = t;
  }
  ret[0] = k1 * q1[0] + k2 * q2[0];
  ret[1] = k1 * q1[1] + k2 * q2[1];
  ret[2] = k1 * q1[2] + k2 * q2[2];
  ret[3] = k1 * q1[3] + k2 * q2[3];
  fi = 1.0f / (sqrtf((ret)[0] * (ret)[0] + (ret)[1] * (ret)[1] + (ret)[2] * (ret)[2] + (ret)[3] * (ret)[3]));
  ret[0] *= fi;
  ret[1] *= fi;
  ret[2] *= fi;
  ret[3] *= fi;
  return t;
}
```