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

```asm
<quaddobl_operate>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x88,%rsp
  mov    %edi,-0x94(%rbp)
  mov    %esi,-0x98(%rbp)
  mov    %rdx,-0xa0(%rbp)
  mov    %rcx,-0xa8(%rbp)
  mov    %r8,-0xb0(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  mov    -0xb0(%rbp),%rdi
  mov    -0xa8(%rbp),%rcx
  mov    -0xa0(%rbp),%rdx
  mov    -0x98(%rbp),%esi
  mov    -0x94(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <numbtrop_quaddobl_initialize@plt>
  mov    %eax,-0x74(%rbp)
  lea    -0x88(%rbp),%rax
  mov    %rax,%rdi
  call   <numbtrop_quaddobl_size@plt>
  mov    %eax,-0x74(%rbp)
  mov    -0x88(%rbp),%eax
  mov    %eax,%esi
  lea    D112(%rip),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  lea    -0x7c(%rbp),%rax
  mov    %rax,%rdi
  call   <numbtrop_quaddobl_dimension@plt>
  mov    %eax,-0x74(%rbp)
  mov    -0x7c(%rbp),%eax
  mov    %eax,%esi
  lea    D4(%rip),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
L62:
  mov    %rsp,%rax
  mov    %rax,%rbx
  mov    -0x98(%rbp),%eax
  shl    $0x2,%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x70(%rbp)
  movslq %eax,%rdx
  mov    %rdx,%r12
  mov    $0x0,%r13d
  movslq %eax,%rdx
  mov    %rdx,%r14
  mov    $0x0,%r15d
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%ecx
  mov    $0x0,%edx
  div    %rcx
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L26:
  cmp    %rdx,%rsp
  je     L6
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L26
L6:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L100
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L100:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x68(%rbp)
  lea    D101(%rip),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  lea    -0x84(%rbp),%rax
  mov    %rax,%rsi
  lea    D93(%rip),%rax
  mov    %rax,%rdi
  call   <scanf@plt>
  mov    -0x84(%rbp),%eax
  test   %eax,%eax
  jg     L23
  mov    %rbx,%rsp
  jmp    L25
L23:
  mov    -0x84(%rbp),%esi
  lea    -0x60(%rbp),%rdi
  mov    -0x68(%rbp),%rcx
  lea    -0x80(%rbp),%rdx
  mov    -0x98(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <numbtrop_quaddobl_retrieve_tropism@plt>
  mov    %eax,-0x74(%rbp)
  mov    -0x80(%rbp),%edx
  mov    -0x84(%rbp),%eax
  mov    %eax,%esi
  lea    D7(%rip),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  lea    D97(%rip),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  movl   $0x0,-0x78(%rbp)
  jmp    L68
L91:
  mov    -0x68(%rbp),%rax
  mov    -0x78(%rbp),%edx
  movslq %edx,%rdx
  mov    (%rax,%rdx,8),%rax
  movq   %rax,%xmm0
  lea    D5(%rip),%rax
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  addl   $0x1,-0x78(%rbp)
L68:
  mov    -0x98(%rbp),%eax
  shl    $0x2,%eax
  cmp    %eax,-0x78(%rbp)
  jl     L91
  movsd  -0x48(%rbp),%xmm2
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mov    -0x60(%rbp),%rax
  movapd %xmm2,%xmm3
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  lea    D8(%rip),%rax
  mov    %rax,%rdi
  mov    $0x4,%eax
  call   <printf@plt>
  lea    -0x80(%rbp),%rax
  mov    %rax,%rsi
  mov    $0x1,%edi
  call   <generate_winding_numbers@plt>
  mov    -0x98(%rbp),%eax
  lea    0x0(,%rax,4),%edx
  mov    -0x68(%rbp),%rax
  mov    %rax,%rsi
  mov    %edx,%edi
  call   <random_doubles@plt>
  lea    -0x60(%rbp),%rax
  mov    %rax,%rsi
  mov    $0x4,%edi
  call   <random_doubles@plt>
  mov    -0x80(%rbp),%edx
  mov    -0x84(%rbp),%eax
  mov    %eax,%esi
  lea    D83(%rip),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  lea    D97(%rip),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  movl   $0x0,-0x78(%rbp)
  jmp    L101
L28:
  mov    -0x68(%rbp),%rax
  mov    -0x78(%rbp),%edx
  movslq %edx,%rdx
  mov    (%rax,%rdx,8),%rax
  movq   %rax,%xmm0
  lea    D5(%rip),%rax
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  addl   $0x1,-0x78(%rbp)
L101:
  mov    -0x98(%rbp),%eax
  shl    $0x2,%eax
  cmp    %eax,-0x78(%rbp)
  jl     L28
  movsd  -0x48(%rbp),%xmm2
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mov    -0x60(%rbp),%rax
  movapd %xmm2,%xmm3
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  lea    D8(%rip),%rax
  mov    %rax,%rdi
  mov    $0x4,%eax
  call   <printf@plt>
  mov    -0x80(%rbp),%edx
  mov    -0x84(%rbp),%esi
  lea    -0x60(%rbp),%rdi
  mov    -0x68(%rbp),%rcx
  mov    -0x98(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <numbtrop_store_quaddobl_tropism@plt>
  mov    %eax,-0x74(%rbp)
  mov    %rbx,%rsp
  mov    -0x84(%rbp),%eax
  test   %eax,%eax
  jg     L62
L25:
  mov    $0x0,%eax
  call   <numbtrop_quaddobl_clear@plt>
  mov    %eax,-0x74(%rbp)
  nop
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L16
  call   <__stack_chk_fail@plt>
L16:
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
```

tool_calls: name: disassemble_data
arguments: tag: D112
data_type: string
name: disassemble_data
arguments: tag: D4
data_type: string
name: disassemble_data
arguments: tag: D101
data_type: string
name: disassemble_data
arguments: tag: D93
data_type: string
name: disassemble_data
arguments: tag: D7
data_type: string
name: disassemble_data
arguments: tag: D97
data_type: string
name: disassemble_data
arguments: tag: D5
data_type: string
name: disassemble_data
arguments: tag: D8
data_type: string
name: disassemble_data
arguments: tag: D83
data_type: string

tool: D112: The number of tropisms : %d


tool: D4: The retrieved dimension : %d


tool: D101: Give an index (<= 0 to exit): 

tool: D93: %d

tool: D7: The tropism %d has winding number %d,


tool: D97: with coordinates :

tool: D5:  %.3e

tool: D8: 
and error : %.3e %.3e %.3e %.3e


tool: D83: Changing tropism %d to winding number %d,


assistant: ```c
void quaddobl_operate(int nbt, int dim, int *wind, double *dir, double *err) {
  int fail, size, k, idx, wnd, retdim;
  fail = numbtrop_quaddobl_initialize(nbt, dim, wind, dir, err);
  fail = numbtrop_quaddobl_size(&size);
  printf("The number of tropisms : %d\n", size);
  fail = numbtrop_quaddobl_dimension(&retdim);
  printf("The retrieved dimension : %d\n", retdim);
  do {
    double diridx[4 * dim];
    double erridx[4];
    printf("Give an index (<= 0 to exit): ");
    scanf("%d", &idx);
    if (idx <= 0)
      break;
    fail = numbtrop_quaddobl_retrieve_tropism(dim, idx, &wnd, diridx, erridx);
    printf("The tropism %d has winding number %d,\n", idx, wnd);
    printf("with coordinates :");
    for (k = 0; k < 4 * dim; k++)
      printf(" %.3e", diridx[k]);
    printf("\nand error : %.3e %.3e %.3e %.3e\n", erridx[0], erridx[1], erridx[2], erridx[3]);
    generate_winding_numbers(1, &wnd);
    random_doubles(4 * dim, diridx);
    random_doubles(4, erridx);
    printf("Changing tropism %d to winding number %d,\n", idx, wnd);
    printf("with coordinates :");
    for (k = 0; k < 4 * dim; k++)
      printf(" %.3e", diridx[k]);
    printf("\nand error : %.3e %.3e %.3e %.3e\n", erridx[0], erridx[1], erridx[2], erridx[3]);
    fail = numbtrop_store_quaddobl_tropism(dim, idx, wnd, diridx, erridx);
  } while (idx > 0);
  fail = numbtrop_quaddobl_clear();
}
```