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

```asm
<main>:
  endbr64
  push   %r15
  mov    $0x10,%esi
  mov    $0x1,%edi
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x28(%rsp)
  xor    %eax,%eax
  call   <vector@plt>
  mov    %rax,%r13
  lea    0x24(%rsp),%rax
  mov    %rax,0x18(%rsp)
L1:
  mov    $0x10,%esi
  lea    D54(%rip),%rdi
  xor    %eax,%eax
  call   <printf@plt>
  mov    0x18(%rsp),%rsi
  lea    D88(%rip),%rdi
  call   <scanf@plt>
  movss  0x24(%rsp),%xmm1
  pxor   %xmm6,%xmm6
  comiss %xmm1,%xmm6
  jae    L28
  xor    %ebx,%ebx
  jmp    L24
  nopl   0x0(%rax)
L118:
  movss  0x24(%rsp),%xmm1
L24:
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm1,%xmm1
  cvtsi2sd %ebx,%xmm0
  mulsd  D25(%rip),%xmm0
  divsd  %xmm1,%xmm0
  call   <sin@plt>
  movss  %xmm0,0x4(%r13,%rbx,4)
  add    $0x1,%rbx
  cmp    $0x10,%rbx
  jne    L118
  mov    $0x10,%esi
  mov    %r13,%rdi
  lea    0x4(%r13),%r12
  call   <sinft@plt>
  movl   $0x501502f9,0x10(%rsp)
  lea    0x44(%r13),%rax
  movss  D119(%rip),%xmm1
  mov    %rax,0x8(%rsp)
  mov    %r12,%rax
  nopl   0x0(%rax)
L47:
  movss  (%rax),%xmm0
  add    $0x4,%rax
  movaps %xmm0,%xmm2
  maxss  %xmm1,%xmm0
  minss  0x10(%rsp),%xmm2
  movaps %xmm0,%xmm1
  movss  %xmm2,0x10(%rsp)
  cmp    %rax,0x8(%rsp)
  jne    L47
  subss  %xmm2,%xmm1
  movss  D84(%rip),%xmm6
  mov    $0x1,%ebp
  lea    D16(%rip),%r14
  divss  %xmm1,%xmm6
  movss  %xmm6,0x14(%rsp)
  nop
L129:
  movss  0x0(%r13,%rbp,4),%xmm0
  mov    %ebp,%esi
  mov    %r14,%rdi
  xor    %eax,%eax
  subss  0x10(%rsp),%xmm0
  mulss  0x14(%rsp),%xmm0
  xor    %r15d,%r15d
  cvtss2sd %xmm0,%xmm0
  addsd  D5(%rip),%xmm0
  addsd  D106(%rip),%xmm0
  cvttsd2si %xmm0,%ebx
  call   <printf@plt>
  test   %ebx,%ebx
  js     L87
  xchg   %ax,%ax
L7:
  mov    $0x2a,%edi
  add    $0x1,%r15d
  call   <putchar@plt>
  cmp    %r15d,%ebx
  jge    L7
L87:
  mov    $0xa,%edi
  add    $0x1,%rbp
  call   <putchar@plt>
  cmp    $0x11,%rbp
  jne    L129
  lea    D82(%rip),%rdi
  call   <puts@plt>
  xor    %eax,%eax
  call   <getchar@plt>
  mov    $0x10,%esi
  mov    %r13,%rdi
  call   <sinft@plt>
  movl   $0x501502f9,0x10(%rsp)
  movss  D119(%rip),%xmm1
  nopw   0x0(%rax,%rax,1)
L62:
  movss  (%r12),%xmm0
  add    $0x4,%r12
  movaps %xmm0,%xmm3
  maxss  %xmm1,%xmm0
  minss  0x10(%rsp),%xmm3
  movaps %xmm0,%xmm1
  movss  %xmm3,0x10(%rsp)
  cmp    %r12,0x8(%rsp)
  jne    L62
  subss  %xmm3,%xmm1
  movss  D84(%rip),%xmm7
  mov    $0x1,%ebp
  divss  %xmm1,%xmm7
  movss  %xmm7,0x8(%rsp)
  nopw   0x0(%rax,%rax,1)
L31:
  movss  0x0(%r13,%rbp,4),%xmm0
  mov    %ebp,%esi
  mov    %r14,%rdi
  xor    %eax,%eax
  subss  0x10(%rsp),%xmm0
  mulss  0x8(%rsp),%xmm0
  xor    %r12d,%r12d
  cvtss2sd %xmm0,%xmm0
  addsd  D5(%rip),%xmm0
  addsd  D106(%rip),%xmm0
  cvttsd2si %xmm0,%ebx
  call   <printf@plt>
  test   %ebx,%ebx
  js     L43
  cs nopw 0x0(%rax,%rax,1)
L55:
  mov    $0x2a,%edi
  add    $0x1,%r12d
  call   <putchar@plt>
  cmp    %r12d,%ebx
  jge    L55
L43:
  mov    $0xa,%edi
  add    $0x1,%rbp
  call   <putchar@plt>
  cmp    $0x11,%rbp
  jne    L31
  jmp    L1
L28:
  mov    $0x10,%edx
  mov    $0x1,%esi
  mov    %r13,%rdi
  call   <free_vector@plt>
  mov    0x28(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L13
  add    $0x38,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L13:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: tag: D54
type_infered: string
name: extract_symbolic_data
arguments: tag: D88
type_infered: string
name: extract_symbolic_data
arguments: tag: D25
type_infered: byte[8]
name: extract_symbolic_data
arguments: tag: D119
type_infered: byte[4]
name: extract_symbolic_data
arguments: tag: D84
type_infered: f32
name: extract_symbolic_data
arguments: tag: D16
type_infered: string
name: extract_symbolic_data
arguments: tag: D5
type_infered: byte[8]
name: extract_symbolic_data
arguments: tag: D106
type_infered: f64
name: extract_symbolic_data
arguments: tag: D82
type_infered: string

tool: D54: 
Period of sinusoid in channels (3-%2d)


tool: D88: %f

tool: {"D25": [0x4a,0xd8,0x12,0x4d,0xfb,0x21,0x19,0x40]}

tool: {"D119": [0xf9,0x02,0x15,0xd0]}

tool: D84: 30

tool: D16: %4d  

tool: {"D5": [0xfc,0xa9,0xf1,0xd2,0x4d,0x62,0x50,0x3f]}

tool: D106: 0.5

tool: D82: press RETURN to continue ...

assistant: ```c
int main(void) {
  float big, per, scal, small, *data;
  int i, j, nlim;
  data = vector(1, 16);
  for (;;) {
    printf("\nPeriod of sinusoid in channels (3-%2d)\n", 16);
    scanf("%f", &per);
    if (per <= 0.0)
      break;
    for (i = 1; i <= 16; i++)
      data[i] = sin(2.0 * 3.1415926 * (i - 1) / per);
    sinft(data, 16);
    big = -1.0e10;
    small = 1.0e10;
    for (i = 1; i <= 16; i++) {
      if (data[i] < small)
        small = data[i];
      if (data[i] > big)
        big = data[i];
    }
    scal = 30.0 / (big - small);
    for (i = 1; i <= 16; i++) {
      nlim = (int)(scal * (data[i] - small) + 1.0e-3 + 0.5);
      printf("%4d  ", i);
      for (j = 1; j <= nlim + 1; j++)
        printf("*");
      printf("\n");
    }
    printf("press RETURN to continue ...\n");
    (void)getchar();
    sinft(data, 16);
    big = -1.0e10;
    small = 1.0e10;
    for (i = 1; i <= 16; i++) {
      if (data[i] < small)
        small = data[i];
      if (data[i] > big)
        big = data[i];
    }
    scal = 30.0 / (big - small);
    for (i = 1; i <= 16; i++) {
      nlim = (int)(scal * (data[i] - small) + 1.0e-3 + 0.5);
      printf("%4d  ", i);
      for (j = 1; j <= nlim + 1; j++)
        printf("*");
      printf("\n");
    }
  }
  free_vector(data, 1, 16);
  return 0;
}
```