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

```asm
<gen_ucsbT_stf>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x38(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %ecx,-0x4c(%rbp)
  mov    %r8,-0x58(%rbp)
  mov    %r9,-0x60(%rbp)
  movss  D72(%rip),%xmm0
  movss  %xmm0,-0x20(%rbp)
  mov    -0x4c(%rbp),%edx
  mov    -0x48(%rbp),%rax
  mov    %edx,%esi
  mov    %rax,%rdi
  call   <zapit@plt>
  mov    -0x60(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x1c(%rbp)
  pxor   %xmm1,%xmm1
  cvtss2sd -0x1c(%rbp),%xmm1
  movsd  D83(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0x60(%rbp),%rax
  movss  (%rax),%xmm1
  cvtss2sd %xmm1,%xmm1
  divsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,-0x18(%rbp)
  movss  -0x1c(%rbp),%xmm0
  subss  -0x18(%rbp),%xmm0
  movss  %xmm0,-0x14(%rbp)
  movss  -0x18(%rbp),%xmm0
  addss  %xmm0,%xmm0
  movss  %xmm0,-0x10(%rbp)
  mov    -0x58(%rbp),%rax
  movss  (%rax),%xmm1
  movss  -0x1c(%rbp),%xmm0
  divss  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd %xmm0,%xmm1
  movsd  D41(%rip),%xmm0
  addsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x2c(%rbp),%eax
  cmp    -0x4c(%rbp),%eax
  jle    L0
  mov    -0x4c(%rbp),%eax
  mov    %eax,-0x2c(%rbp)
L0:
  cmpl   $0x0,-0x2c(%rbp)
  jne    L68
  mov    $0x0,%eax
  jmp    L59
L68:
  movl   $0x0,-0x30(%rbp)
  jmp    L4
L105:
  pxor   %xmm1,%xmm1
  cvtsi2ssl -0x30(%rbp),%xmm1
  mov    -0x58(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0xc(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x24(%rbp)
  movss  -0x18(%rbp),%xmm0
  comiss -0xc(%rbp),%xmm0
  jbe    L140
  movss  -0x20(%rbp),%xmm0
  mulss  -0xc(%rbp),%xmm0
  divss  -0x18(%rbp),%xmm0
  movss  %xmm0,-0x8(%rbp)
  movss  -0x8(%rbp),%xmm1
  movss  D85(%rip),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x4(%rbp)
  mov    -0x8(%rbp),%eax
  movd   %eax,%xmm0
  call   <cos@plt>
  movsd  D99(%rip),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  D99(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  mov    -0x4(%rbp),%eax
  movd   %eax,%xmm0
  call   <sin@plt>
  movsd  D106(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  addsd  -0x68(%rbp),%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,-0x24(%rbp)
  jmp    L126
L140:
  movss  -0x10(%rbp),%xmm0
  comiss -0xc(%rbp),%xmm0
  jbe    L87
  movss  -0x20(%rbp),%xmm0
  mulss  -0xc(%rbp),%xmm0
  divss  -0x18(%rbp),%xmm0
  movss  %xmm0,-0x8(%rbp)
  movss  -0xc(%rbp),%xmm0
  subss  -0x18(%rbp),%xmm0
  mulss  -0x20(%rbp),%xmm0
  divss  -0x14(%rbp),%xmm0
  movss  %xmm0,-0x4(%rbp)
  mov    -0x8(%rbp),%eax
  movd   %eax,%xmm0
  call   <cos@plt>
  movsd  D99(%rip),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  D71(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  mov    -0x4(%rbp),%eax
  movd   %eax,%xmm0
  call   <cos@plt>
  movsd  D114(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  addsd  -0x68(%rbp),%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,-0x24(%rbp)
  jmp    L126
L87:
  movss  -0x1c(%rbp),%xmm0
  comiss -0xc(%rbp),%xmm0
  jbe    L126
  movss  -0xc(%rbp),%xmm0
  subss  -0x18(%rbp),%xmm0
  mulss  -0x20(%rbp),%xmm0
  divss  -0x14(%rbp),%xmm0
  movss  %xmm0,-0x8(%rbp)
  mov    -0x8(%rbp),%eax
  movd   %eax,%xmm0
  call   <cos@plt>
  movsd  D114(%rip),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  D114(%rip),%xmm0
  addsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,-0x24(%rbp)
L126:
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  -0x24(%rbp),%xmm0
  movss  %xmm0,(%rax)
  addl   $0x1,-0x30(%rbp)
L4:
  mov    -0x30(%rbp),%eax
  cmp    -0x2c(%rbp),%eax
  jl     L105
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x28(%rbp)
  movl   $0x0,-0x30(%rbp)
  jmp    L23
L53:
  mov    -0x58(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  -0x28(%rbp),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x28(%rbp)
  addl   $0x1,-0x30(%rbp)
L23:
  mov    -0x30(%rbp),%eax
  cmp    -0x2c(%rbp),%eax
  jl     L53
  pxor   %xmm0,%xmm0
  comiss -0x28(%rbp),%xmm0
  jb     L15
  mov    $0x0,%eax
  jmp    L59
L15:
  mov    -0x38(%rbp),%rax
  movss  (%rax),%xmm0
  divss  -0x28(%rbp),%xmm0
  movss  %xmm0,-0x28(%rbp)
  movl   $0x0,-0x30(%rbp)
  jmp    L14
L71:
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mulss  -0x28(%rbp),%xmm0
  movss  %xmm0,(%rax)
  addl   $0x1,-0x30(%rbp)
L14:
  mov    -0x30(%rbp),%eax
  cmp    -0x2c(%rbp),%eax
  jl     L71
  mov    -0x2c(%rbp),%eax
L59:
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: tag: D72
data_type: f32
name: retrieve_data
arguments: tag: D83
data_type: f64
name: retrieve_data
arguments: tag: D41
data_type: f64
name: retrieve_data
arguments: tag: D85
data_type: byte[4]
name: retrieve_data
arguments: tag: D99
data_type: f64
name: retrieve_data
arguments: tag: D106
data_type: f64
name: retrieve_data
arguments: tag: D71
data_type: f64
name: retrieve_data
arguments: tag: D114
data_type: f64

tool: D72: 3.141592654

tool: D83: 0.13

tool: D41: 0.5

tool: {"D85": [0x00,0x00,0x00,0x3f]}

tool: D99: 0.7

tool: D106: 0.6

tool: D71: 1

tool: D114: 0.3

assistant: ```c
int gen_ucsbT_stf(float *slip, float *t0, float *stf, int nt, float *dt, float *beta) {
  int it, nstf;
  float tau, tau1, tau2, tau1x2, arg1, arg2;
  float sum, t, alpha;
  float pi = 3.141592654;
  zapit(stf, nt);
  tau = (*beta) * (*t0);
  tau1 = 0.13 * tau / (*beta);
  tau2 = tau - tau1;
  tau1x2 = 2.0 * tau1;
  nstf = (int)((tau) / (*dt) + 0.5);
  if (nstf > nt)
    nstf = nt;
  if (nstf == 0)
    return (0);
  for (it = 0; it < nstf; it++) {
    t = it * (*dt);
    alpha = 0.0;
    if (t < tau1) {
      arg1 = pi * t / tau1;
      arg2 = 0.5 * arg1;
      alpha = 0.7 - 0.7 * cos(arg1) + 0.6 * sin(arg2);
    } else if (t < tau1x2) {
      arg1 = pi * t / tau1;
      arg2 = pi * (t - tau1) / tau2;
      alpha = 1.0 - 0.7 * cos(arg1) + 0.3 * cos(arg2);
    } else if (t < tau) {
      arg1 = pi * (t - tau1) / tau2;
      alpha = 0.3 + 0.3 * cos(arg1);
    }
    stf[it] = alpha;
  }
  sum = 0.0;
  for (it = 0; it < nstf; it++)
    sum = sum + (*dt) * stf[it];
  if (sum <= 0.0)
    return (0);
  sum = (*slip) / sum;
  for (it = 0; it < nstf; it++)
    stf[it] = stf[it] * sum;
  return (nstf);
}
```