user: 请将以下的汇编代码转换成C语言代码：

```asm
<vec4_slerp>:
  endbr64
  push   %r12
  movaps %xmm0,%xmm5
  pxor   %xmm7,%xmm7
  mov    %rdi,%r12
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  mov    %rdx,%rbx
  sub    $0x20,%rsp
  movss  0xc(%rsi),%xmm0
  movss  (%rsi),%xmm1
  mulss  0xc(%rdx),%xmm0
  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
  comiss %xmm0,%xmm7
  ja     L65
  movss  D27(%rip),%xmm7
  movss  %xmm7,0xc(%rsp)
L51:
  movss  %xmm7,0x8(%rsp)
  movss  %xmm5,0x4(%rsp)
  call   <acosf@plt>
  movss  %xmm0,0x14(%rsp)
  call   <sinf@plt>
  movss  %xmm0,0x10(%rsp)
  call   <fabs@plt>
  comiss D59(%rip),%xmm0
  movss  0x4(%rsp),%xmm5
  movss  0x8(%rsp),%xmm7
  jbe    L33
  comiss D17(%rip),%xmm5
  ja     L53
L33:
  movaps %xmm7,%xmm2
  subss  %xmm5,%xmm2
L18:
  movss  0x0(%rbp),%xmm4
  movss  (%rbx),%xmm0
  mulss  %xmm2,%xmm4
  mulss  %xmm5,%xmm0
  addss  %xmm0,%xmm4
  movss  %xmm4,(%r12)
  movss  0x4(%rbp),%xmm3
  movss  0x4(%rbx),%xmm0
  mulss  %xmm2,%xmm3
  mulss  %xmm5,%xmm0
  addss  %xmm0,%xmm3
  movss  %xmm3,0x4(%r12)
  movss  0x8(%rbp),%xmm0
  movss  0x8(%rbx),%xmm1
  mulss  %xmm2,%xmm0
  mulss  %xmm5,%xmm1
  movaps %xmm0,%xmm6
  movaps %xmm4,%xmm0
  addss  %xmm1,%xmm6
  mulss  %xmm4,%xmm0
  movaps %xmm3,%xmm1
  mulss  %xmm3,%xmm1
  movss  %xmm6,0x8(%r12)
  mulss  0xc(%rbx),%xmm5
  mulss  0xc(%rbp),%xmm2
  addss  %xmm1,%xmm0
  movaps %xmm6,%xmm1
  mulss  %xmm6,%xmm1
  addss  %xmm5,%xmm2
  pxor   %xmm5,%xmm5
  addss  %xmm1,%xmm0
  movaps %xmm2,%xmm1
  mulss  %xmm2,%xmm1
  addss  %xmm1,%xmm0
  ucomiss %xmm0,%xmm5
  ja     L57
  sqrtss %xmm0,%xmm0
L1:
  divss  %xmm0,%xmm7
  mulss  %xmm7,%xmm4
  mulss  %xmm7,%xmm3
  mulss  %xmm7,%xmm6
  mulss  %xmm7,%xmm2
  movss  %xmm4,(%r12)
  movss  %xmm3,0x4(%r12)
  movss  %xmm6,0x8(%r12)
  movss  %xmm2,0xc(%r12)
  add    $0x20,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
  nopl   0x0(%rax)
L53:
  comiss %xmm5,%xmm7
  jbe    L33
  movaps %xmm7,%xmm0
  movss  0x14(%rsp),%xmm1
  movss  %xmm7,0x1c(%rsp)
  subss  %xmm5,%xmm0
  movss  %xmm5,0x18(%rsp)
  mulss  %xmm1,%xmm0
  call   <sinf@plt>
  movss  0x10(%rsp),%xmm3
  movss  0x18(%rsp),%xmm5
  movaps %xmm0,%xmm2
  movss  0x14(%rsp),%xmm1
  movss  0xc(%rsp),%xmm0
  divss  %xmm3,%xmm2
  movss  %xmm3,0x8(%rsp)
  mulss  %xmm5,%xmm1
  mulss  %xmm1,%xmm0
  movss  %xmm2,0x4(%rsp)
  call   <sinf@plt>
  movss  0x8(%rsp),%xmm3
  movss  0x4(%rsp),%xmm2
  movaps %xmm0,%xmm5
  movss  0x1c(%rsp),%xmm7
  divss  %xmm3,%xmm5
  jmp    L18
  nopl   0x0(%rax)
L65:
  movss  D42(%rip),%xmm6
  xorps  D43(%rip),%xmm0
  movss  D27(%rip),%xmm7
  movss  %xmm6,0xc(%rsp)
  jmp    L51
L57:
  movss  %xmm7,0x14(%rsp)
  movss  %xmm2,0x10(%rsp)
  movss  %xmm6,0xc(%rsp)
  movss  %xmm3,0x8(%rsp)
  movss  %xmm4,0x4(%rsp)
  call   <sqrtf@plt>
  movss  0x14(%rsp),%xmm7
  movss  0x10(%rsp),%xmm2
  movss  0xc(%rsp),%xmm6
  movss  0x8(%rsp),%xmm3
  movss  0x4(%rsp),%xmm4
  jmp    L1
```

tool_calls: name: read_data
arguments: label: D27
datatype: byte[4]
name: read_data
arguments: label: D59
datatype: byte[4]
name: read_data
arguments: label: D17
datatype: byte[4]
name: read_data
arguments: label: D42
datatype: byte[4]
name: read_data
arguments: label: D43
datatype: byte[16]

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

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

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

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

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

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;
}
```