user: 下面的汇编指令对应的C源代码是什么:

```asm
<ComputeHueQuadrature>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xa0,%rsp
  movsd  %xmm0,-0x88(%rbp)
  mov    %rdi,-0x90(%rbp)
  mov    %rsi,-0x98(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  D118(%rip),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  D16(%rip),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  D0(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  D44(%rip),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  D113(%rip),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  D132(%rip),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  D35(%rip),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x88(%rbp),%xmm0
  comisd D52(%rip),%xmm0
  jb     L79
  movsd  D75(%rip),%xmm0
  comisd -0x88(%rbp),%xmm0
  jbe    L79
  movl   $0x0,-0x78(%rbp)
  movl   $0x1,-0x74(%rbp)
  jmp    L88
L79:
  movsd  -0x88(%rbp),%xmm0
  comisd D75(%rip),%xmm0
  jb     L45
  movsd  D39(%rip),%xmm0
  comisd -0x88(%rbp),%xmm0
  jbe    L45
  movl   $0x1,-0x78(%rbp)
  movl   $0x2,-0x74(%rbp)
  jmp    L88
L45:
  movsd  -0x88(%rbp),%xmm0
  comisd D39(%rip),%xmm0
  jb     L83
  movsd  D117(%rip),%xmm0
  comisd -0x88(%rbp),%xmm0
  jbe    L83
  movl   $0x2,-0x78(%rbp)
  movl   $0x3,-0x74(%rbp)
  jmp    L88
L83:
  movl   $0x3,-0x78(%rbp)
  movl   $0x0,-0x74(%rbp)
L88:
  mov    -0x78(%rbp),%eax
  cltq
  movsd  -0x50(%rbp,%rax,8),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0x74(%rbp),%eax
  cltq
  movsd  -0x50(%rbp,%rax,8),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  mov    -0x78(%rbp),%eax
  cltq
  movsd  -0x30(%rbp,%rax,8),%xmm0
  movsd  %xmm0,-0x60(%rbp)
  mov    -0x74(%rbp),%eax
  cltq
  movsd  -0x30(%rbp,%rax,8),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x68(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x70(%rbp),%xmm1
  movsd  -0x88(%rbp),%xmm0
  subsd  -0x60(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x58(%rbp),%xmm1
  subsd  -0x60(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  addsd  -0x70(%rbp),%xmm0
  mov    -0x98(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  -0x88(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x60(%rbp),%xmm1
  movsd  D113(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  divsd  -0x70(%rbp),%xmm0
  movsd  -0x88(%rbp),%xmm1
  subsd  -0x60(%rbp),%xmm1
  movapd %xmm1,%xmm2
  divsd  -0x70(%rbp),%xmm2
  movsd  -0x58(%rbp),%xmm1
  subsd  -0x88(%rbp),%xmm1
  divsd  -0x68(%rbp),%xmm1
  addsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  addsd  -0x60(%rbp),%xmm0
  mov    -0x90(%rbp),%rax
  movsd  %xmm0,(%rax)
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L62
  call   <__stack_chk_fail@plt>
L62:
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: data_label: D118
type: f64
name: retrieve_data
arguments: data_label: D16
type: f64
name: retrieve_data
arguments: data_label: D0
type: f64
name: retrieve_data
arguments: data_label: D44
type: f64
name: retrieve_data
arguments: data_label: D113
type: f64
name: retrieve_data
arguments: data_label: D132
type: byte[8]
name: retrieve_data
arguments: data_label: D35
type: byte[8]
name: retrieve_data
arguments: data_label: D52
type: f64[2]
name: retrieve_data
arguments: data_label: D75
type: f64[2]
name: retrieve_data
arguments: data_label: D39
type: f64[2]
name: retrieve_data
arguments: data_label: D117
type: f64

tool: D118: 0.8

tool: D16: 0.7

tool: D0: 1

tool: D44: 1.2

tool: D113: 100

tool: {"D132": [0x00,0x00,0x00,0x00,0x00,0x00,0x69,0x40]}

tool: {"D35": [0x00,0x00,0x00,0x00,0x00,0xc0,0x72,0x40]}

tool: {"D52": [20.14,90.0]}

tool: {"D75": [90.0,164.25]}

tool: {"D39": [164.25,237.53]}

tool: D117: 237.53

assistant: ```c
void ComputeHueQuadrature(double h, double *H, double *e) {
  double e_tab[] = {0.8, 0.7, 1.0, 1.2};
  double H_tab[] = {0, 100, 200, 300};
  int p1, p2;
  double e1, e2, h1, h2;
  if (h >= 20.14 && h < 90.0) {
    p1 = 0;
    p2 = 1;
  } else if (h >= 90.0 && h < 164.25) {
    p1 = 1;
    p2 = 2;
  } else if (h >= 164.25 && h < 237.53) {
    p1 = 2;
    p2 = 3;
  } else {
    p1 = 3;
    p2 = 0;
  }
  e1 = e_tab[p1];
  e2 = e_tab[p2];
  h1 = H_tab[p1];
  h2 = H_tab[p2];
  *e = e1 + ((e2 - e1) * (h - h1) / (h2 - h1));
  *H = h1 + (100. * (h - h1) / e1) / ((h - h1) / e1 + (h2 - h) / e2);
}
```