user: 请从给定的汇编指令推导出C源代码:

```asm
<gen_brune_stf>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x50,%rsp
  mov    %rdi,-0x28(%rbp)
  mov    %rsi,-0x30(%rbp)
  mov    %rdx,-0x38(%rbp)
  mov    %ecx,-0x3c(%rbp)
  mov    %r8,-0x48(%rbp)
  mov    -0x3c(%rbp),%edx
  mov    -0x38(%rbp),%rax
  mov    %edx,%esi
  mov    %rax,%rdi
  call   <zapit@plt>
  mov    $0x1,%edi
  call   <exp@plt>
  pxor   %xmm1,%xmm1
  cvtss2sd %xmm0,%xmm1
  movsd  D14(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x30(%rbp),%rax
  movss  (%rax),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,-0x10(%rbp)
  movss  -0x10(%rbp),%xmm1
  movss  D4(%rip),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0xc(%rbp)
  mov    -0x48(%rbp),%rax
  movss  (%rax),%xmm1
  movss  -0xc(%rbp),%xmm0
  divss  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd %xmm0,%xmm1
  movsd  D13(%rip),%xmm0
  addsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x18(%rbp)
  mov    -0x18(%rbp),%eax
  cmp    -0x3c(%rbp),%eax
  jle    L17
  mov    -0x3c(%rbp),%eax
  mov    %eax,-0x18(%rbp)
L17:
  cmpl   $0x0,-0x18(%rbp)
  jne    L14
  mov    $0x0,%eax
  jmp    L13
L14:
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm0
  mov    -0x30(%rbp),%rax
  movss  (%rax),%xmm1
  divss  %xmm1,%xmm0
  movss  %xmm0,-0x8(%rbp)
  mov    -0x48(%rbp),%rax
  movss  (%rax),%xmm0
  mov    -0x30(%rbp),%rax
  movss  (%rax),%xmm1
  divss  %xmm1,%xmm0
  movss  %xmm0,-0x4(%rbp)
  movl   $0x0,-0x1c(%rbp)
  jmp    L6
L74:
  pxor   %xmm0,%xmm0
  cvtsi2ssl -0x1c(%rbp),%xmm0
  mulss  -0x4(%rbp),%xmm0
  mulss  -0x8(%rbp),%xmm0
  movss  %xmm0,-0x4c(%rbp)
  mov    -0x1c(%rbp),%eax
  neg    %eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  mulss  -0x4(%rbp),%xmm0
  cvttss2si %xmm0,%eax
  mov    %eax,%edi
  call   <exp@plt>
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mulss  -0x4c(%rbp),%xmm0
  movss  %xmm0,(%rax)
  addl   $0x1,-0x1c(%rbp)
L6:
  mov    -0x1c(%rbp),%eax
  cmp    -0x18(%rbp),%eax
  jl     L74
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x14(%rbp)
  movl   $0x0,-0x1c(%rbp)
  jmp    L25
L59:
  mov    -0x48(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  -0x14(%rbp),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x14(%rbp)
  addl   $0x1,-0x1c(%rbp)
L25:
  mov    -0x1c(%rbp),%eax
  cmp    -0x18(%rbp),%eax
  jl     L59
  pxor   %xmm0,%xmm0
  comiss -0x14(%rbp),%xmm0
  jb     L86
  mov    $0x0,%eax
  jmp    L13
L86:
  mov    -0x18(%rbp),%eax
L13:
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_tag: D14
datatype: f64
name: parse_data
arguments: data_tag: D4
datatype: f32
name: parse_data
arguments: data_tag: D13
datatype: f64

tool: D14: 1.745

tool: D4: 3

tool: D13: 0.5

assistant: ```c
int gen_brune_stf(float *slip, float *t0, float *stf, int nt, float *dt) {
  int it, nstf;
  float t95, tend, sfac, tfac;
  float sum;
  zapit(stf, nt);
  t95 = 1.745 * exp(1.0) * (*t0);
  tend = 3.0 * t95;
  nstf = (int)((tend) / (*dt) + 0.5);
  if (nstf > nt)
    nstf = nt;
  if (nstf == 0)
    return (0);
  sfac = (*slip) / (*t0);
  tfac = (*dt) / (*t0);
  for (it = 0; it < nstf; it++)
    stf[it] = sfac * (it * tfac) * exp(-it * tfac);
  sum = 0.0;
  for (it = 0; it < nstf; it++)
    sum = sum + (*dt) * stf[it];
  if (sum <= 0.0)
    return (0);
  return (nstf);
}
```