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

```asm
<vec4_slerp>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x18,%rsp
  mov    %rdi,%r12
  mov    %rsi,%rbp
  mov    %rdx,%rbx
  movss  %xmm0,0x4(%rsp)
  movss  0xc(%rsi),%xmm0
  mulss  0xc(%rdx),%xmm0
  movss  (%rsi),%xmm1
  mulss  (%rdx),%xmm1
  addss  %xmm1,%xmm0
  movss  0x4(%rsi),%xmm1
  mulss  0x4(%rdx),%xmm1
  addss  %xmm1,%xmm0
  movss  0x8(%rsi),%xmm1
  mulss  0x8(%rdx),%xmm1
  addss  %xmm1,%xmm0
  mov    D8(%rip),%r13d
  pxor   %xmm1,%xmm1
  comiss %xmm0,%xmm1
  ja     L0
  mov    D31(%rip),%r13d
L0:
  movd   %r13d,%xmm7
  mulss  %xmm7,%xmm0
  call   <acosf@plt>
  movss  %xmm0,0x8(%rsp)
  call   <sinf@plt>
  movss  %xmm0,0xc(%rsp)
  call   <fabs@plt>
  comiss D45(%rip),%xmm0
  jbe    L30
  movss  0x4(%rsp),%xmm7
  comiss D57(%rip),%xmm7
  jbe    L30
  movss  D31(%rip),%xmm0
  comiss %xmm7,%xmm0
  ja     L59
L30:
  movss  D31(%rip),%xmm0
  subss  0x4(%rsp),%xmm0
  movd   %xmm0,%r14d
L12:
  movd   %r14d,%xmm2
  movd   %r14d,%xmm1
  mulss  0x0(%rbp),%xmm1
  movss  0x4(%rsp),%xmm3
  movaps %xmm3,%xmm0
  mulss  (%rbx),%xmm0
  addss  %xmm0,%xmm1
  movd   %xmm1,%r14d
  movss  %xmm1,(%r12)
  movaps %xmm2,%xmm1
  mulss  0x4(%rbp),%xmm1
  movaps %xmm3,%xmm0
  mulss  0x4(%rbx),%xmm0
  addss  %xmm0,%xmm1
  movd   %xmm1,%r13d
  movss  %xmm1,0x4(%r12)
  movaps %xmm2,%xmm1
  mulss  0x8(%rbp),%xmm1
  movaps %xmm3,%xmm0
  mulss  0x8(%rbx),%xmm0
  addss  %xmm0,%xmm1
  movd   %xmm1,%r15d
  movss  %xmm1,0x8(%r12)
  mulss  0xc(%rbp),%xmm2
  mulss  0xc(%rbx),%xmm3
  movaps %xmm2,%xmm4
  addss  %xmm3,%xmm4
  movss  %xmm4,0x4(%rsp)
  movd   %r14d,%xmm0
  mulss  %xmm0,%xmm0
  movd   %r13d,%xmm1
  mulss  %xmm1,%xmm1
  addss  %xmm1,%xmm0
  movd   %r15d,%xmm1
  mulss  %xmm1,%xmm1
  addss  %xmm1,%xmm0
  movaps %xmm4,%xmm1
  mulss  %xmm4,%xmm1
  addss  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  ucomiss %xmm0,%xmm1
  ja     L14
  sqrtss %xmm0,%xmm0
L52:
  movss  D31(%rip),%xmm1
  divss  %xmm0,%xmm1
  movd   %r14d,%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,(%r12)
  movd   %r13d,%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,0x4(%r12)
  movd   %r15d,%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,0x8(%r12)
  mulss  0x4(%rsp),%xmm1
  movss  %xmm1,0xc(%r12)
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L59:
  subss  %xmm7,%xmm0
  mulss  0x8(%rsp),%xmm0
  call   <sinf@plt>
  movaps %xmm0,%xmm4
  divss  0xc(%rsp),%xmm4
  movd   %xmm4,%r14d
  movss  0x8(%rsp),%xmm0
  mulss  0x4(%rsp),%xmm0
  movd   %r13d,%xmm5
  mulss  %xmm5,%xmm0
  call   <sinf@plt>
  divss  0xc(%rsp),%xmm0
  movss  %xmm0,0x4(%rsp)
  jmp    L12
L14:
  call   <sqrtf@plt>
  jmp    L52
```

tool_calls: name: disassemble_data
arguments: tag: D8
type_infered: byte[4]
name: disassemble_data
arguments: tag: D31
type_infered: byte[4]
name: disassemble_data
arguments: tag: D45
type_infered: byte[4]
name: disassemble_data
arguments: tag: D57
type_infered: byte[4]

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

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

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

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

assistant: ```c
void vec4_slerp(float ret[4], float q1[4], float q2[4], float t) {
  float cos_fi, sin_fi, fi, k1, k2, sign;
  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);
  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;
}
```