user: What would the C source look like for the assembly code provided here?

```asm
<smath>:
  endbr64
  push   %r15
  push   %r14
  mov    %rsi,%r14
  lea    0x1(%rdx),%esi
  push   %r13
  mov    %rdi,%r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x418,%rsp
  mov    %edx,0x14(%rsp)
  movss  %xmm0,0x18(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x408(%rsp)
  mov    %edx,%eax
  mov    %esi,0x1c(%rsp)
  add    %eax,%eax
  mov    %eax,%ebx
  js     L112
  lea    0x2(%rax),%eax
  pxor   %xmm7,%xmm7
  mov    %ebx,%r15d
  mov    $0x1,%ebp
  cvtsi2ss %eax,%xmm7
  add    $0x2,%r15
  lea    0x20(%rsp),%r12
  movss  %xmm7,0x10(%rsp)
  movss  D35(%rip),%xmm7
  movss  %xmm7,0xc(%rsp)
  nopl   (%rax)
L113:
  pxor   %xmm0,%xmm0
  cvtsi2ss %ebp,%xmm0
  mulss  D69(%rip),%xmm0
  divss  0x10(%rsp),%xmm0
  cvttss2si %xmm0,%edi
  call   <cos@plt>
  pxor   %xmm1,%xmm1
  movss  0xc(%rsp),%xmm0
  cvtsi2ss %rax,%xmm1
  subss  %xmm1,%xmm0
  mulss  D42(%rip),%xmm0
  movss  %xmm0,-0x4(%r12,%rbp,4)
  add    $0x1,%rbp
  cmp    %rbp,%r15
  jne    L113
L112:
  movslq 0x14(%rsp),%rax
  lea    0x40(%rsp),%r15
  movl   $0x0,0x20(%rsp,%rax,4)
  movss  0x20(%rsp),%xmm1
  xor    %eax,%eax
  nopl   0x0(%rax)
L81:
  movss  (%r14,%rax,1),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,(%r15,%rax,1)
  add    $0x4,%rax
  cmp    $0x140,%rax
  jne    L81
  cmpl   $0x1,0x14(%rsp)
  jle    L189
  mov    0x14(%rsp),%eax
  lea    0x20(%rsp),%rsi
  lea    0x140(%r14),%rdi
  lea    0x180(%rsp),%rcx
  sub    $0x2,%eax
  lea    0x24(%rsp,%rax,4),%r8
L89:
  movss  0x4(%rsi),%xmm1
  mov    %r15,%rax
  mov    %rdi,%rdx
  xchg   %ax,%ax
L147:
  movss  (%rdx),%xmm0
  add    $0x4,%rax
  add    $0x4,%rdx
  mulss  %xmm1,%xmm0
  addss  -0x4(%rax),%xmm0
  movss  %xmm0,-0x4(%rax)
  cmp    %rax,%rcx
  jne    L147
  add    $0x4,%rsi
  add    $0x140,%rdi
  cmp    %rsi,%r8
  jne    L89
L189:
  cmp    0x1c(%rsp),%ebx
  jl     L157
  movslq 0x1c(%rsp),%rdi
  lea    0x20(%rsp),%r12
  lea    0x180(%rsp),%rcx
  lea    (%rdi,%rdi,4),%r8d
  shl    $0x4,%r8d
  movslq %r8d,%r8
  nopw   0x0(%rax,%rax,1)
L3:
  movss  (%r12,%rdi,4),%xmm1
  mov    %r15,%rax
  lea    (%r14,%r8,4),%rdx
  nopl   (%rax)
L65:
  movss  (%rdx),%xmm0
  add    $0x4,%rax
  add    $0x4,%rdx
  mulss  %xmm1,%xmm0
  addss  -0x4(%rax),%xmm0
  movss  %xmm0,-0x4(%rax)
  cmp    %rax,%rcx
  jne    L65
  add    $0x1,%rdi
  add    $0x50,%r8
  cmp    %edi,%ebx
  jge    L3
L157:
  mov    0x14(%rsp),%eax
  pxor   %xmm2,%xmm2
  movaps %xmm2,%xmm4
  movaps %xmm2,%xmm3
  lea    (%rax,%rax,4),%eax
  shl    $0x4,%eax
  cltq
  lea    (%r14,%rax,4),%rbx
  xor    %eax,%eax
  nop
L1:
  movss  (%rbx,%rax,1),%xmm0
  movaps %xmm0,%xmm1
  mulss  %xmm0,%xmm1
  addss  %xmm1,%xmm3
  movss  (%r15,%rax,1),%xmm1
  add    $0x4,%rax
  mulss  %xmm1,%xmm0
  movaps %xmm1,%xmm5
  mulss  %xmm1,%xmm5
  addss  %xmm0,%xmm4
  addss  %xmm5,%xmm2
  cmp    $0x140,%rax
  jne    L1
  movaps %xmm2,%xmm0
  movss  %xmm4,0x14(%rsp)
  movss  %xmm3,0x10(%rsp)
  movss  %xmm2,0xc(%rsp)
  call   <fabs@plt>
  movsd  D48(%rip),%xmm1
  movss  0x10(%rsp),%xmm3
  movss  0x14(%rsp),%xmm4
  movss  D35(%rip),%xmm2
  comisd %xmm0,%xmm1
  movaps %xmm3,%xmm0
  ja     L52
  movss  0xc(%rsp),%xmm2
  divss  %xmm2,%xmm0
L52:
  movss  %xmm4,0x14(%rsp)
  movss  %xmm2,0x10(%rsp)
  movss  %xmm3,0xc(%rsp)
  call   <sqrt@plt>
  pxor   %xmm6,%xmm6
  pxor   %xmm5,%xmm5
  movss  0x14(%rsp),%xmm4
  movss  0x10(%rsp),%xmm2
  movss  0xc(%rsp),%xmm3
  cvtsi2ss %rax,%xmm6
  xor    %eax,%eax
  cs nopw 0x0(%rax,%rax,1)
L156:
  movss  (%r15,%rax,1),%xmm1
  mulss  %xmm6,%xmm1
  movss  %xmm1,0x0(%r13,%rax,1)
  movss  (%rbx,%rax,1),%xmm0
  add    $0x4,%rax
  subss  %xmm1,%xmm0
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm5
  cmp    $0x140,%rax
  jne    L156
  movss  0x18(%rsp),%xmm0
  mulss  %xmm3,%xmm0
  comiss %xmm0,%xmm5
  jbe    L54
  movss  D35(%rip),%xmm7
  comiss %xmm3,%xmm7
  movss  %xmm7,0xc(%rsp)
  ja     L154
  mulss  %xmm3,%xmm2
  movaps %xmm3,%xmm0
  mulss  %xmm3,%xmm0
L117:
  movaps %xmm4,%xmm1
  mulss  %xmm4,%xmm1
  subss  %xmm1,%xmm2
  divss  %xmm0,%xmm2
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm2,%xmm0
  comisd D33(%rip),%xmm0
  ja     L83
  movss  0xc(%rsp),%xmm2
  pxor   %xmm5,%xmm5
L19:
  xor    %eax,%eax
  nopl   0x0(%rax)
L101:
  movss  (%rbx,%rax,1),%xmm1
  movss  (%r15,%rax,1),%xmm0
  mulss  %xmm2,%xmm1
  mulss  %xmm5,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,0x0(%r13,%rax,1)
  add    $0x4,%rax
  cmp    $0x140,%rax
  jne    L101
L54:
  mov    0x408(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L14
  add    $0x418,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L154:
  movss  0xc(%rsp),%xmm3
  movaps %xmm3,%xmm0
  jmp    L117
L83:
  movss  0x18(%rsp),%xmm7
  movss  %xmm4,0x14(%rsp)
  movss  %xmm3,0x10(%rsp)
  movaps %xmm7,%xmm1
  movaps %xmm7,%xmm0
  mulss  %xmm7,%xmm1
  mulss  D41(%rip),%xmm1
  subss  %xmm1,%xmm0
  divss  %xmm2,%xmm0
  call   <sqrt@plt>
  pxor   %xmm5,%xmm5
  movss  0x14(%rsp),%xmm4
  movss  0x10(%rsp),%xmm3
  cvtsi2ss %rax,%xmm5
  movss  0x18(%rsp),%xmm0
  xorps  D90(%rip),%xmm0
  mulss  D42(%rip),%xmm0
  movss  0xc(%rsp),%xmm2
  mulss  %xmm5,%xmm4
  divss  %xmm3,%xmm4
  subss  %xmm4,%xmm0
  addss  %xmm0,%xmm2
  jmp    L19
L14:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: label: D35
type: f32
name: parse_data
arguments: label: D69
type: byte[4]
name: parse_data
arguments: label: D42
type: f32
name: parse_data
arguments: label: D48
type: f32[2]
name: parse_data
arguments: label: D33
type: f64[2]
name: parse_data
arguments: label: D41
type: byte[4]
name: parse_data
arguments: label: D90
type: byte[16]

tool: D35: 1

tool: {"D69": [0xdb,0x0f,0xc9,0x40]}

tool: D42: 0.5

tool: {"D48": [0.0,1.875]}

tool: {"D33": [0.0001,1.0609978955e-314]}

tool: {"D41": [0x00,0x00,0x80,0x3e]}

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

assistant: ```c
void smath(float *odata, float *sseq, int hl, float alpha0) {
  int i, k;
  float w00, w10, w11, A, B, C, *psseq, err, errs;
  float surround[240];
  float wt[2 * 3 + 1];
  float denom;
  for (i = 1; i <= 2 * hl + 1; i++) {
    wt[i - 1] = (float)0.5 * (1 - (float)cos(2 * (float)3.14159265358979323846 * i / (2 * hl + 2)));
  }
  wt[hl] = 0.0;
  for (i = 0; i < 80; i++) {
    surround[i] = sseq[i] * wt[0];
  }
  for (k = 1; k < hl; k++) {
    psseq = sseq + k * 80;
    for (i = 0; i < 80; i++) {
      surround[i] += psseq[i] * wt[k];
    }
  }
  for (k = hl + 1; k <= 2 * hl; k++) {
    psseq = sseq + k * 80;
    for (i = 0; i < 80; i++) {
      surround[i] += psseq[i] * wt[k];
    }
  }
  w00 = w10 = w11 = 0.0;
  psseq = sseq + hl * 80;
  for (i = 0; i < 80; i++) {
    w00 += psseq[i] * psseq[i];
    w11 += surround[i] * surround[i];
    w10 += surround[i] * psseq[i];
  }
  if (fabs(w11) < 1.0) {
    w11 = 1.0;
  }
  C = (float)sqrt(w00 / w11);
  errs = 0.0;
  psseq = sseq + hl * 80;
  for (i = 0; i < 80; i++) {
    odata[i] = C * surround[i];
    err = psseq[i] - odata[i];
    errs += err * err;
  }
  if (errs > alpha0 * w00) {
    if (w00 < 1) {
      w00 = 1;
    }
    denom = (w11 * w00 - w10 * w10) / (w00 * w00);
    if (denom > 0.0001) {
      A = (float)sqrt((alpha0 - alpha0 * alpha0 / 4) / denom);
      B = -alpha0 / 2 - A * w10 / w00;
      B = B + 1;
    } else {
      A = 0.0;
      B = 1.0;
    }
    psseq = sseq + hl * 80;
    for (i = 0; i < 80; i++) {
      odata[i] = A * surround[i] + B * psseq[i];
    }
  }
}
```