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

```asm
<apprgrdn>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  movsd  %xmm0,0x10(%rsp)
  mov    %r8,0x18(%rsp)
  test   %di,%di
  je     L152
  mov    %rsi,%r14
  mov    %rdx,%rbp
  mov    %rcx,%r13
  mov    %r9d,%r15d
  movzwl %di,%eax
  shl    $0x3,%rax
  mov    %rax,0x20(%rsp)
  mov    $0x0,%ebx
  jmp    L103
L37:
  movsd  0x8(%rsp),%xmm0
  andpd  D3(%rip),%xmm0
  movsd  D20(%rip),%xmm2
  comisd %xmm0,%xmm2
  jbe    L145
  movsd  D64(%rip),%xmm3
  movsd  %xmm3,0x8(%rsp)
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  ja     L145
  movsd  %xmm2,0x8(%rsp)
L145:
  movq   %r12,%xmm0
  addsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x0(%rbp,%rbx,1)
  mov    %rbp,%rdi
  mov    $0x0,%eax
  call   *%r13
  subsd  0x10(%rsp),%xmm0
  divsd  0x8(%rsp),%xmm0
  movsd  %xmm0,(%r14,%rbx,1)
L127:
  mov    %r12,0x0(%rbp,%rbx,1)
  add    $0x8,%rbx
  cmp    0x20(%rsp),%rbx
  je     L152
L103:
  mov    0x0(%rbp,%rbx,1),%r12
  movq   %r12,%xmm0
  andpd  D3(%rip),%xmm0
  movsd  D20(%rip),%xmm3
  comisd %xmm0,%xmm3
  jbe    L153
  movapd %xmm3,%xmm0
L153:
  mov    0x18(%rsp),%rax
  movsd  (%rax,%rbx,1),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,0x8(%rsp)
  test   %r15w,%r15w
  je     L37
  andpd  D3(%rip),%xmm0
  movsd  D10(%rip),%xmm2
  comisd %xmm0,%xmm2
  jbe    L146
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  ja     L75
  movw   $0x1,0x2c(%rsp)
  movsd  %xmm2,0x8(%rsp)
  jmp    L22
L146:
  movw   $0x0,0x2c(%rsp)
L22:
  movq   %r12,%xmm0
  addsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x0(%rbp,%rbx,1)
  mov    %rbp,%rdi
  mov    $0x0,%eax
  call   *%r13
  movsd  0x10(%rsp),%xmm7
  ucomisd %xmm7,%xmm0
  jp     L94
  je     L143
L94:
  subsd  0x10(%rsp),%xmm0
  movsd  0x8(%rsp),%xmm6
  divsd  %xmm6,%xmm0
  movsd  %xmm0,(%r14,%rbx,1)
  cmpw   $0x0,0x2c(%rsp)
  je     L127
  movq   %r12,%xmm4
  subsd  %xmm6,%xmm4
  movsd  %xmm4,0x0(%rbp,%rbx,1)
  mov    %rbp,%rdi
  mov    $0x0,%eax
  call   *%r13
  movapd %xmm0,%xmm1
  movsd  0x10(%rsp),%xmm0
  subsd  %xmm1,%xmm0
  divsd  0x8(%rsp),%xmm0
  addsd  (%r14,%rbx,1),%xmm0
  mulsd  D19(%rip),%xmm0
  movsd  %xmm0,(%r14,%rbx,1)
  jmp    L127
L75:
  movw   $0x1,0x2c(%rsp)
  movsd  D1(%rip),%xmm7
  movsd  %xmm7,0x8(%rsp)
  jmp    L22
L143:
  movw   $0x3,0x2e(%rsp)
L123:
  movsd  D96(%rip),%xmm5
  mulsd  0x8(%rsp),%xmm5
  movsd  %xmm5,0x8(%rsp)
  movq   %r12,%xmm0
  addsd  %xmm5,%xmm0
  movsd  %xmm0,0x0(%rbp,%rbx,1)
  mov    %rbp,%rdi
  mov    $0x0,%eax
  call   *%r13
  movsd  0x10(%rsp),%xmm2
  ucomisd %xmm0,%xmm2
  jp     L94
  jne    L94
  subw   $0x1,0x2e(%rsp)
  jne    L123
  jmp    L94
L152:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D3
type_infered: qword[2]
name: extract_symbolic_data
arguments: data_tag: D20
type_infered: byte[8]
name: extract_symbolic_data
arguments: data_tag: D64
type_infered: byte[8]
name: extract_symbolic_data
arguments: data_tag: D10
type_infered: byte[8]
name: extract_symbolic_data
arguments: data_tag: D19
type_infered: f64
name: extract_symbolic_data
arguments: data_tag: D1
type_infered: byte[8]
name: extract_symbolic_data
arguments: data_tag: D96
type_infered: f64

tool: {"D3": [0x7fffffffffffffff,0x0000000000000000]}

tool: {"D20": [0x9b,0x2b,0xa1,0x86,0x9b,0x84,0xf6,0x3c]}

tool: {"D64": [0x9b,0x2b,0xa1,0x86,0x9b,0x84,0xf6,0xbc]}

tool: {"D10": [0xbb,0xbd,0xd7,0xd9,0xdf,0x7c,0xeb,0x3d]}

tool: D19: 0.5

tool: {"D1": [0xbb,0xbd,0xd7,0xd9,0xdf,0x7c,0xeb,0xbd]}

tool: D96: 10

assistant: ```c
void apprgrdn(unsigned short n, double g[], double x[], double f, double fun(), double deltax[], unsigned short obj) {
  double const lowbndobj = 2.0e-10, lowbndcnt = 5.0e-15, ten = 10.0, half = 0.5;
  double d, y, fi;
  unsigned short i, j, center = 0;
  for (i = 0; i < n; i++) {
    y = x[i];
    d = ((lowbndcnt) > (fabs(y)) ? (lowbndcnt) : (fabs(y)));
    d *= deltax[i];
    if (obj) {
      if (fabs(d) < lowbndobj) {
        if (deltax[i] < 0.0)
          d = -lowbndobj;
        else
          d = lowbndobj;
        center = 1;
      } else
        center = 0;
    } else if (fabs(d) < lowbndcnt) {
      if (deltax[i] < 0.0)
        d = -lowbndcnt;
      else
        d = lowbndcnt;
    }
    x[i] = y + d;
    fi = fun(x);
    if (obj) {
      if (fi == f) {
        for (j = 1; j <= 3; j++) {
          d *= ten;
          x[i] = y + d;
          fi = fun(x);
          if (fi != f)
            break;
        }
      }
    }
    g[i] = (fi - f) / d;
    if (obj) {
      if (center) {
        x[i] = y - d;
        fi = fun(x);
        g[i] = half * (g[i] + (f - fi) / d);
      }
    }
    x[i] = y;
  }
}
```