user: 请从给定的汇编指令推导出C源代码:

```asm
<Mutation_Kpoint>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %rsi,-0x50(%rbp)
  mov    %edx,-0x54(%rbp)
  mov    %rcx,-0x60(%rbp)
  mov    %r8,-0x68(%rbp)
  mov    %r9d,-0x58(%rbp)
  movsd  %xmm0,-0x70(%rbp)
  movsd  %xmm1,-0x78(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0x80(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0x88(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    -0x50(%rbp),%rax
  mov    (%rax),%rcx
  mov    -0x58(%rbp),%edx
  mov    -0x68(%rbp),%rsi
  lea    -0x34(%rbp),%rax
  movq   %rcx,%xmm0
  mov    %rax,%rdi
  call   <self_adj_index_search@plt>
  movl   $0x1,-0x2c(%rbp)
  jmp    L28
L93:
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rax,%rdx
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x88(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x2c(%rbp)
L28:
  mov    -0x2c(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  jle    L93
  mov    -0x54(%rbp),%eax
  mov    %eax,%esi
  mov    $0x1,%edi
  call   <integerrng@plt>
  mov    %eax,-0x24(%rbp)
  mov    $0x0,%eax
  call   <normalrng@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm3
  mulsd  -0x78(%rbp),%xmm3
  movsd  %xmm3,-0x90(%rbp)
  mov    -0x34(%rbp),%eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  pxor   %xmm2,%xmm2
  cvtsi2sdl -0x58(%rbp),%xmm2
  movsd  D46(%rip),%xmm1
  addsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  mov    D0(%rip),%rdx
  movq   %rdx,%xmm0
  mov    %eax,%edi
  call   <pow@plt>
  mulsd  -0x90(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    $0x0,%eax
  call   <uniformrng@plt>
  movsd  D0(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L18
  cmpl   $0x1,-0x54(%rbp)
  jle    L18
  mov    -0x54(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,%esi
  mov    $0x1,%edi
  call   <integerrng@plt>
  mov    %eax,-0x28(%rbp)
  mov    -0x28(%rbp),%eax
  cmp    -0x24(%rbp),%eax
  jl     L85
  addl   $0x1,-0x28(%rbp)
L85:
  mov    $0x0,%eax
  call   <normalrng@plt>
  movsd  -0x78(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,(%rax)
L18:
  mov    -0x54(%rbp),%edx
  mov    -0x88(%rbp),%rax
  mov    %edx,%esi
  mov    %rax,%rdi
  call   <cost@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x58(%rbp),%edx
  mov    -0x68(%rbp),%rsi
  mov    -0x18(%rbp),%rcx
  lea    -0x30(%rbp),%rax
  movq   %rcx,%xmm0
  mov    %rax,%rdi
  call   <self_adj_index_search@plt>
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D8(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  movapd %xmm1,%xmm2
  divsd  -0x70(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm0,%xmm1
  mov    -0x50(%rbp),%rax
  movsd  (%rax),%xmm0
  divsd  -0x70(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L50
  mov    D46(%rip),%rax
  jmp    L90
L50:
  mov    -0x10(%rbp),%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  movq   %xmm0,%rax
L90:
  mov    -0x80(%rbp),%rdx
  mov    %rax,(%rdx)
  mov    $0x0,%eax
  call   <uniformrng@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x20(%rbp)
  mov    -0x80(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd -0x20(%rbp),%xmm0
  ja     L58
  jmp    L115
L58:
  mov    -0x50(%rbp),%rax
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movl   $0x1,-0x2c(%rbp)
  jmp    L89
L22:
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x88(%rbp),%rax
  add    %rax,%rdx
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x2c(%rbp)
L89:
  mov    -0x2c(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  jle    L22
L115:
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L74
  call   <__stack_chk_fail@plt>
L74:
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D46
type_infered: f64
name: extract_symbolic_data
arguments: data_label: D0
type_infered: f64
name: extract_symbolic_data
arguments: data_label: D8
type_infered: f32[2]

tool: D46: 1

tool: D0: 0.5

tool: {"D8": [0.0,-0.0]}

assistant: ```c
void Mutation_Kpoint(double *z, double *fz, int N_dimension, double *theta, double *grid_points, int grid_size, double temp, double scl, double *accpr_pop, double *z_new) {
  int k_old;
  int k_new;
  int i;
  int lab_1;
  int lab_2;
  double fz_new;
  double rat;
  double un;
  self_adj_index_search(&k_old, *fz, grid_points, grid_size);
  for (i = 1; i <= N_dimension; i++)
    z_new[i] = z[i];
  lab_1 = integerrng(1, N_dimension);
  un = normalrng() * scl * pow(k_old / (grid_size + 1.0), 0.5);
  z_new[lab_1] += un;
  if ((uniformrng() < 0.5) && (N_dimension >= 2)) {
    lab_2 = integerrng(1, N_dimension - 1);
    if (lab_2 >= lab_1)
      lab_2++;
    un = normalrng() * scl;
    z_new[lab_2] += un;
  }
  fz_new = cost(z_new, N_dimension);
  self_adj_index_search(&k_new, fz_new, grid_points, grid_size);
  rat = -theta[k_new] - fz_new / temp + theta[k_old] + *fz / temp;
  *accpr_pop = ((rat > 0.0) ? 1.0 : exp(rat));
  un = uniformrng();
  if (*accpr_pop > un) {
    *fz = fz_new;
    for (i = 1; i <= N_dimension; i++)
      z[i] = z_new[i];
  }
}
```