user: Please show me the C implementation that matches the assembly code listed below:

```asm
<vec4_slerp_to>:
  endbr64
  push   %r14
  movaps %xmm0,%xmm1
  pxor   %xmm6,%xmm6
  push   %r12
  mov    %rdi,%r12
  push   %rbp
  mov    %rdx,%rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x28,%rsp
  movss  0xc(%rsi),%xmm0
  movss  (%rsi),%xmm2
  mulss  0xc(%rdx),%xmm0
  mulss  (%rdx),%xmm2
  addss  %xmm2,%xmm0
  movss  0x4(%rsi),%xmm2
  mulss  0x4(%rdx),%xmm2
  addss  %xmm2,%xmm0
  movss  0x8(%rsi),%xmm2
  mulss  0x8(%rdx),%xmm2
  addss  %xmm2,%xmm0
  comiss %xmm0,%xmm6
  ja     L1
  movss  D11(%rip),%xmm6
  movd   %xmm6,%r14d
  movss  %xmm6,0x18(%rsp)
L6:
  movss  %xmm1,0x14(%rsp)
  call   <acosf@plt>
  movss  %xmm0,0xc(%rsp)
  call   <sinf@plt>
  movss  %xmm0,0x10(%rsp)
  movss  0xc(%rsp),%xmm0
  call   <fabs@plt>
  movss  0x14(%rsp),%xmm1
  comiss %xmm1,%xmm0
  jbe    L84
  divss  %xmm0,%xmm1
  movss  0x10(%rsp),%xmm0
  movss  %xmm1,0x14(%rsp)
  call   <fabs@plt>
  movss  0x14(%rsp),%xmm1
  movd   %r14d,%xmm2
  comiss D22(%rip),%xmm0
  subss  %xmm1,%xmm2
  jbe    L60
  comiss D46(%rip),%xmm1
  ja     L35
L60:
  movss  (%rbx),%xmm5
  movss  0x0(%rbp),%xmm0
  pxor   %xmm7,%xmm7
  mulss  %xmm1,%xmm0
  mulss  %xmm2,%xmm5
  addss  %xmm0,%xmm5
  movss  %xmm5,(%r12)
  movss  0x4(%rbx),%xmm4
  movss  0x4(%rbp),%xmm0
  mulss  %xmm2,%xmm4
  mulss  %xmm1,%xmm0
  addss  %xmm0,%xmm4
  movss  %xmm4,0x4(%r12)
  movss  0x8(%rbx),%xmm3
  movss  0x8(%rbp),%xmm0
  mulss  %xmm2,%xmm3
  mulss  %xmm1,%xmm0
  addss  %xmm0,%xmm3
  movaps %xmm5,%xmm0
  mulss  %xmm5,%xmm0
  movss  %xmm3,0x8(%r12)
  mulss  0xc(%rbx),%xmm2
  mulss  0xc(%rbp),%xmm1
  addss  %xmm2,%xmm1
  movaps %xmm4,%xmm2
  mulss  %xmm4,%xmm2
  addss  %xmm2,%xmm0
  movaps %xmm3,%xmm2
  mulss  %xmm3,%xmm2
  addss  %xmm2,%xmm0
  movaps %xmm1,%xmm2
  mulss  %xmm1,%xmm2
  addss  %xmm2,%xmm0
  ucomiss %xmm0,%xmm7
  ja     L42
  sqrtss %xmm0,%xmm0
L2:
  movd   %r14d,%xmm2
  divss  %xmm0,%xmm2
  mulss  %xmm2,%xmm5
  mulss  %xmm2,%xmm4
  mulss  %xmm2,%xmm3
  mulss  %xmm2,%xmm1
  movss  %xmm5,(%r12)
  movss  %xmm4,0x4(%r12)
  movss  %xmm3,0x8(%r12)
  movss  %xmm1,0xc(%r12)
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r14
  ret
  nopw   0x0(%rax,%rax,1)
L84:
  movss  0x10(%rsp),%xmm0
  call   <fabs@plt>
  pxor   %xmm2,%xmm2
  movd   %r14d,%xmm1
  jmp    L60
  nopl   0x0(%rax)
L1:
  movss  D40(%rip),%xmm7
  xorps  D20(%rip),%xmm0
  movss  %xmm7,0x18(%rsp)
  movss  D11(%rip),%xmm7
  movd   %xmm7,%r14d
  jmp    L6
  nopw   0x0(%rax,%rax,1)
L35:
  movd   %r14d,%xmm7
  comiss %xmm1,%xmm7
  jbe    L60
  mulss  0xc(%rsp),%xmm2
  movss  %xmm1,0x1c(%rsp)
  movaps %xmm2,%xmm0
  call   <sinf@plt>
  movss  0x1c(%rsp),%xmm1
  movaps %xmm0,%xmm2
  movss  0xc(%rsp),%xmm0
  divss  0x10(%rsp),%xmm2
  movss  %xmm2,0x14(%rsp)
  mulss  %xmm1,%xmm0
  mulss  0x18(%rsp),%xmm0
  call   <sinf@plt>
  movss  0x14(%rsp),%xmm2
  movaps %xmm0,%xmm1
  divss  0x10(%rsp),%xmm1
  jmp    L60
L42:
  movss  %xmm1,0x18(%rsp)
  movss  %xmm3,0x14(%rsp)
  movss  %xmm4,0x10(%rsp)
  movss  %xmm5,0xc(%rsp)
  call   <sqrtf@plt>
  movss  0x18(%rsp),%xmm1
  movss  0x14(%rsp),%xmm3
  movss  0x10(%rsp),%xmm4
  movss  0xc(%rsp),%xmm5
  jmp    L2
```

tool_calls: name: parse_data
arguments: data_label: D11
type: byte[4]
name: parse_data
arguments: data_label: D22
type: byte[4]
name: parse_data
arguments: data_label: D46
type: byte[4]
name: parse_data
arguments: data_label: D40
type: byte[4]
name: parse_data
arguments: data_label: D20
type: byte[16]

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

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

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

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

tool: {"D20": [0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

assistant: ```c
void vec4_slerp_to(float ret[4], float q1[4], 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;
}
```