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

```asm
<Mutation_DoubleSkeletalChange>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,%r14
  mov    %rsi,0x38(%rsp)
  mov    %rdx,%r10
  mov    %rdx,0x28(%rsp)
  mov    %rcx,%r15
  mov    %r8d,%ebx
  mov    %r8d,0xc(%rsp)
  mov    %r9,0x18(%rsp)
  mov    0xc0(%rsp),%rdi
  mov    %rdi,0x20(%rsp)
  mov    0xc8(%rsp),%rdi
  mov    %rdi,0x40(%rsp)
  mov    0xd0(%rsp),%rax
  mov    %rax,%rsi
  mov    %rax,0x48(%rsp)
  movsd  %xmm0,0x50(%rsp)
  mov    0xe0(%rsp),%rdi
  mov    %rdi,0x58(%rsp)
  mov    0xe8(%rsp),%r13
  mov    0xf0(%rsp),%rbp
  mov    %rbp,0x10(%rsp)
  mov    0xf8(%rsp),%r12
  mov    %r12,0x30(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x78(%rsp)
  xor    %eax,%eax
  lea    0x70(%rsp),%rdi
  mov    0xd8(%rsp),%edx
  movsd  (%r10),%xmm0
  call   <self_adj_index_search@plt>
  test   %ebx,%ebx
  jle    L133
  mov    0x38(%rsp),%rax
  lea    0x8(%rax),%rdi
  lea    0x8(%rbp),%rsi
  lea    0x1(%rbx),%r10d
  lea    0x0(,%r10,4),%r8
  mov    $0x1,%r9d
  mov    %r12,%r11
L179:
  mov    (%r14,%r9,4),%eax
  mov    %eax,0x0(%r13,%r9,4)
  movsd  (%r15,%r9,8),%xmm0
  movsd  %xmm0,(%r11,%r9,8)
  mov    $0x4,%eax
L159:
  mov    (%rdi),%rdx
  mov    (%rdx,%rax,1),%ecx
  mov    (%rsi),%rdx
  mov    %ecx,(%rdx,%rax,1)
  add    $0x4,%rax
  cmp    %rax,%r8
  jne    L159
  add    $0x1,%r9
  add    $0x8,%rdi
  add    $0x8,%rsi
  cmp    %r10,%r9
  jne    L179
L133:
  mov    0xc(%rsp),%eax
  lea    -0x1(%rax),%ebx
  mov    %r14,0x60(%rsp)
  mov    %ebx,%r14d
  mov    %r15,0x68(%rsp)
  mov    %eax,%r15d
  jmp    L131
L62:
  cmp    %ebx,%eax
  jne    L176
  cmp    %ebp,%r12d
  jne    L176
L131:
  mov    %r15d,%esi
  mov    $0x1,%edi
  call   <integerrng@plt>
  mov    %eax,%ebp
  mov    %r14d,%esi
  mov    $0x1,%edi
  call   <integerrng@plt>
  mov    %eax,%ebx
  cmp    %eax,%ebp
  jg     L157
  add    $0x1,%ebx
  cmp    %ebx,%ebp
  jg     L157
  mov    %ebx,%eax
  mov    %ebp,%ebx
  mov    %eax,%ebp
L157:
  mov    %r15d,%esi
  mov    $0x1,%edi
  call   <integerrng@plt>
  mov    %eax,%r12d
  mov    %r14d,%esi
  mov    $0x1,%edi
  call   <integerrng@plt>
  cmp    %eax,%r12d
  jg     L62
  add    $0x1,%eax
  cmp    %eax,%r12d
  jg     L62
  mov    %eax,%edx
  mov    %r12d,%eax
  mov    %edx,%r12d
  jmp    L62
L176:
  mov    0x60(%rsp),%r14
  mov    0x68(%rsp),%r15
  movslq %ebx,%rbx
  movslq %ebp,%rdx
  mov    0x10(%rsp),%rdi
  mov    (%rdi,%rbx,8),%rcx
  lea    (%rcx,%rdx,4),%rcx
  mov    $0x1,%edx
  mov    %edx,%esi
  sub    (%rcx),%esi
  mov    %esi,(%rcx)
  cltq
  movslq %r12d,%rcx
  mov    (%rdi,%rax,8),%rsi
  lea    (%rsi,%rcx,4),%rax
  sub    (%rax),%edx
  mov    %edx,(%rax)
  cmp    %ebp,%r12d
  je     L127
  mov    0x18(%rsp),%rax
  movl   $0x2,(%rax)
  mov    0x20(%rsp),%rax
  mov    %ebp,0x4(%rax)
  mov    %r12d,0x8(%rax)
L123:
  mov    0x18(%rsp),%rax
  mov    (%rax),%r9d
  mov    0x20(%rsp),%r8
  mov    0x30(%rsp),%rcx
  mov    0xc(%rsp),%edx
  mov    0x10(%rsp),%rsi
  mov    %r13,%rdi
  call   <cost@plt>
  movq   %xmm0,%rbx
  lea    0x74(%rsp),%rdi
  mov    0xd8(%rsp),%edx
  mov    0x48(%rsp),%rsi
  call   <self_adj_index_search@plt>
  movslq 0x74(%rsp),%rax
  mov    0x40(%rsp),%rsi
  movsd  (%rsi,%rax,8),%xmm0
  xorpd  D19(%rip),%xmm0
  movq   %rbx,%xmm1
  movsd  0x50(%rsp),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movslq 0x70(%rsp),%rax
  addsd  (%rsi,%rax,8),%xmm0
  mov    0x28(%rsp),%rax
  movsd  (%rax),%xmm1
  divsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  comisd D18(%rip),%xmm0
  jbe    L26
  pxor   %xmm0,%xmm0
L26:
  call   <exp@plt>
  mov    0x58(%rsp),%rbp
  movsd  %xmm0,0x0(%rbp)
  mov    $0x0,%eax
  call   <uniformrng@plt>
  movsd  0x0(%rbp),%xmm1
  comisd %xmm0,%xmm1
  jbe    L57
  mov    0x28(%rsp),%rax
  mov    %rbx,(%rax)
  mov    0xc(%rsp),%eax
  test   %eax,%eax
  jle    L57
  mov    0x10(%rsp),%rdi
  add    $0x8,%rdi
  mov    0x38(%rsp),%rsi
  add    $0x8,%rsi
  lea    0x1(%rax),%r10d
  lea    0x0(,%r10,4),%r8
  mov    $0x1,%r9d
  mov    0x30(%rsp),%r11
L85:
  mov    0x0(%r13,%r9,4),%eax
  mov    %eax,(%r14,%r9,4)
  movsd  (%r11,%r9,8),%xmm0
  movsd  %xmm0,(%r15,%r9,8)
  mov    $0x4,%eax
L39:
  mov    (%rdi),%rdx
  mov    (%rdx,%rax,1),%ecx
  mov    (%rsi),%rdx
  mov    %ecx,(%rdx,%rax,1)
  add    $0x4,%rax
  cmp    %r8,%rax
  jne    L39
  add    $0x1,%r9
  add    $0x8,%rdi
  add    $0x8,%rsi
  cmp    %r10,%r9
  jne    L85
L57:
  mov    0x78(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L135
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L127:
  mov    0x18(%rsp),%rax
  movl   $0x1,(%rax)
  mov    0x20(%rsp),%rax
  mov    %ebp,0x4(%rax)
  jmp    L123
L135:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_tag: D19
type_infered: f32[4]
name: retrieve_data
arguments: data_tag: D18
type_infered: f32

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

tool: D18: 0

assistant: ```c
void Mutation_DoubleSkeletalChange(int *z, int **zmat, double *fz, double *fz_cond, int N_node, int *changelength, int *changelist, double *theta, double *grid_points, int grid_size, double temp, double *accpr_pop, int *z_new, int **zmat_new, double *fz_cond_new) {
  int k_old;
  int k_new;
  int j;
  int i;
  int nodeA1;
  int nodeA2;
  int nodeB1;
  int nodeB2;
  double fz_new;
  double rat;
  double un;
  self_adj_index_search(&k_old, *fz, grid_points, grid_size);
  for (i = 1; i <= N_node; i++) {
    z_new[i] = z[i];
    fz_cond_new[i] = fz_cond[i];
    for (j = 1; j <= N_node; j++)
      zmat_new[i][j] = zmat[i][j];
  }
  nodeA1 = 0;
  nodeA2 = 0;
  nodeB1 = 0;
  nodeB2 = 0;
  while (nodeA1 == nodeB1 && nodeA2 == nodeB2) {
    nodeA1 = integerrng(1, N_node);
    nodeA2 = integerrng(1, N_node - 1);
    if (nodeA2 >= nodeA1)
      nodeA2++;
    if (nodeA1 > nodeA2) {
      i = nodeA1;
      nodeA1 = nodeA2;
      nodeA2 = i;
    }
    nodeB1 = integerrng(1, N_node);
    nodeB2 = integerrng(1, N_node - 1);
    if (nodeB2 >= nodeB1)
      nodeB2++;
    if (nodeB1 > nodeB2) {
      i = nodeB1;
      nodeB1 = nodeB2;
      nodeB2 = i;
    }
  }
  zmat_new[nodeA1][nodeA2] = 1 - zmat_new[nodeA1][nodeA2];
  zmat_new[nodeB1][nodeB2] = 1 - zmat_new[nodeB1][nodeB2];
  if (nodeA2 == nodeB2) {
    *changelength = 1;
    changelist[1] = nodeA2;
  } else {
    *changelength = 2;
    changelist[1] = nodeA2;
    changelist[2] = nodeB2;
  }
  fz_new = cost(z_new, zmat_new, N_node, fz_cond_new, changelist, *changelength);
  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 = exp((((0.0) < (rat)) ? (0.0) : (rat)));
  un = uniformrng();
  if (*accpr_pop > un) {
    *fz = fz_new;
    for (i = 1; i <= N_node; i++) {
      z[i] = z_new[i];
      fz_cond[i] = fz_cond_new[i];
      for (j = 1; j <= N_node; j++)
        zmat[i][j] = zmat_new[i][j];
    }
  }
}
```