user: 提供的汇编代码用C语言表示会是什么样子:

```asm
<Mutation_int_Gibbs>:
  endbr64
  push   %r15
  mov    %rdi,%r10
  mov    %r9d,%r15d
  push   %r14
  mov    %r8,%r14
  push   %r13
  mov    %edx,%r13d
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %rcx,0x28(%rsp)
  mov    0xa8(%rsp),%rbp
  movsd  %xmm0,0x30(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
  test   %edx,%edx
  jle    L14
  lea    0x1(%rdx),%ecx
  mov    $0x1,%eax
L13:
  mov    (%r10,%rax,4),%edx
  mov    %edx,0x0(%rbp,%rax,4)
  add    $0x1,%rax
  cmp    %rax,%rcx
  jne    L13
L14:
  mov    0xa0(%rsp),%eax
  test   %eax,%eax
  jle    L21
  lea    0x50(%rsp),%rax
  mov    %r10,0x48(%rsp)
  mov    $0x1,%r8d
  mov    %rax,0x38(%rsp)
  lea    0x54(%rsp),%rax
  mov    %rax,0x40(%rsp)
  jmp    L51
  nopl   0x0(%rax)
L1:
  movsd  0x18(%rsp),%xmm2
  mov    (%rbx),%eax
  add    $0x1,%r8d
  mov    0x48(%rsp),%rsi
  mov    0x20(%rsp),%r9
  movsd  %xmm2,(%r12)
  mov    %eax,(%rsi,%r9,4)
  cmp    %r8d,0xa0(%rsp)
  jl     L21
L51:
  mov    %r13d,%esi
  mov    $0x1,%edi
  mov    %r8d,0x14(%rsp)
  call   <integerrng@plt>
  mov    0x38(%rsp),%rdi
  mov    %r15d,%edx
  mov    %r14,%rsi
  movsd  (%r12),%xmm0
  mov    %eax,%ebx
  call   <self_adj_index_search@plt>
  movslq %ebx,%r9
  mov    $0x1,%eax
  mov    %r13d,%esi
  lea    0x0(%rbp,%r9,4),%rbx
  mov    %rbp,%rdi
  mov    %r9,0x20(%rsp)
  sub    (%rbx),%eax
  mov    %eax,(%rbx)
  call   <cost@plt>
  mov    0x40(%rsp),%rdi
  mov    %r14,%rsi
  mov    %r15d,%edx
  movsd  %xmm0,0x8(%rsp)
  call   <self_adj_index_search@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm4
  movsd  (%r12),%xmm1
  mov    0x28(%rsp),%rdx
  movapd %xmm2,%xmm0
  movslq 0x54(%rsp),%rax
  movsd  %xmm2,0x18(%rsp)
  divsd  %xmm4,%xmm0
  divsd  %xmm4,%xmm1
  addsd  (%rdx,%rax,8),%xmm0
  movslq 0x50(%rsp),%rax
  subsd  (%rdx,%rax,8),%xmm0
  subsd  %xmm1,%xmm0
  call   <exp@plt>
  xor    %eax,%eax
  movsd  %xmm0,0x8(%rsp)
  call   <uniformrng@plt>
  movsd  0x8(%rsp),%xmm1
  mov    0x14(%rsp),%r8d
  mov    $0x1,%ecx
  movapd %xmm0,%xmm3
  movsd  D12(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  D12(%rip),%xmm1
  divsd  %xmm0,%xmm1
  comisd %xmm3,%xmm1
  jae    L1
  sub    (%rbx),%ecx
  add    $0x1,%r8d
  mov    %ecx,(%rbx)
  cmp    %r8d,0xa0(%rsp)
  jge    L51
L21:
  mov    0x58(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L27
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L27:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D12
type_guessed: f64

tool: D12: 1

assistant: ```c
void Mutation_int_Gibbs(int *z, double *fz, int N_dimension, double *theta, double *grid_points, int grid_size, double temp, int Krep, int *z_new) {
  int k_old;
  int k_new;
  int i;
  double Pr;
  int lab;
  int gg;
  double fz_new;
  double rat;
  double un;
  for (i = 1; i <= N_dimension; i++)
    z_new[i] = z[i];
  for (gg = 1; gg <= Krep; gg++) {
    lab = integerrng(1, N_dimension);
    self_adj_index_search(&k_old, *fz, grid_points, grid_size);
    z_new[lab] = 1 - z_new[lab];
    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;
    Pr = 1.0 / (1.0 + exp(rat));
    un = uniformrng();
    if (Pr >= un) {
      *fz = fz_new;
      z[lab] = z_new[lab];
    } else {
      z_new[lab] = 1 - z_new[lab];
    }
  }
}
```