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

```asm
<newviewFlexGamma_Ancestral>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x78,%rsp
  mov    %edi,%ebx
  mov    %edi,0x38(%rsp)
  mov    %rsi,0x48(%rsp)
  mov    %rdx,0x40(%rsp)
  mov    %rcx,(%rsp)
  mov    %r8,0x20(%rsp)
  mov    %r9,0x18(%rsp)
  mov    0xd0(%rsp),%r12d
  movslq %r12d,%r15
  mov    %r12d,%edi
  shl    $0x5,%edi
  call   <rax_malloc@plt>
  mov    %rax,0x8(%rsp)
  mov    %r12d,%eax
  imul   %r12d,%eax
  mov    %eax,0x28(%rsp)
  lea    0x0(,%r12,4),%eax
  mov    %eax,0x2c(%rsp)
  cmp    $0x1,%ebx
  je     L105
  cmp    $0x2,%ebx
  je     L485
  test   %ebx,%ebx
  je     L481
  mov    $0x0,%edi
  call   <assert@plt>
L273:
  mov    0x8(%rsp),%rdi
  call   <rax_free@plt>
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L485:
  cmpl   $0x0,0xb8(%rsp)
  jle    L273
  mov    0x2c(%rsp),%eax
  movslq %eax,%rdi
  mov    %rdi,0x38(%rsp)
  shl    $0x3,%r15
  lea    -0x1(%r12),%ebx
  mov    0x8(%rsp),%rdi
  lea    0x8(%rdi,%rbx,8),%rsi
  mov    %rsi,0x30(%rsp)
  lea    -0x1(%rax),%eax
  lea    0x8(%rdi,%rax,8),%rax
  mov    %rax,0x50(%rsp)
  movq   $0x0,0x20(%rsp)
  movl   $0x0,0x18(%rsp)
  jmp    L164
L105:
  cmpl   $0x0,0xb8(%rsp)
  jle    L273
  mov    0x2c(%rsp),%eax
  movslq %eax,%rdi
  mov    %rdi,0x58(%rsp)
  mov    0xb8(%rsp),%edi
  mov    %rdi,0x60(%rsp)
  lea    0x0(,%r15,8),%r13
  lea    -0x1(%r12),%ebx
  mov    0x8(%rsp),%rdi
  lea    0x8(%rdi,%rbx,8),%rsi
  mov    %rsi,0x50(%rsp)
  lea    -0x1(%rax),%eax
  lea    0x8(%rdi,%rax,8),%rax
  mov    %rax,0x68(%rsp)
  movq   $0x0,0x38(%rsp)
  movq   $0x0,0x10(%rsp)
  jmp    L78
L481:
  cmpl   $0x0,0xb8(%rsp)
  jle    L273
  mov    0xb8(%rsp),%eax
  mov    %rax,0x48(%rsp)
  lea    0x0(,%r15,8),%r14
  lea    -0x1(%r12),%ebx
  mov    0x8(%rsp),%rax
  lea    0x8(%rax,%rbx,8),%rax
  mov    %rax,0x50(%rsp)
  mov    $0x0,%r13d
  jmp    L270
L223:
  mov    0x2c(%rsp),%r9d
  mov    %r12d,%r8d
  mov    0x40(%rsp),%ecx
  mov    0xd8(%rsp),%rdx
  mov    0xe0(%rsp),%rsi
  mov    0x8(%rsp),%rdi
  call   <ancestralGamma@plt>
  add    $0x1,%r13
  cmp    %r13,0x48(%rsp)
  je     L273
L270:
  mov    %r13d,0x40(%rsp)
  mov    0x50(%rsp),%r11
  mov    0x8(%rsp),%rbp
  movl   $0x0,0x10(%rsp)
  movl   $0x4,0x30(%rsp)
L168:
  mov    0x18(%rsp),%rax
  movzbl (%rax,%r13,1),%r8d
  imul   %r12d,%r8d
  movslq %r8d,%r8
  mov    0xb0(%rsp),%rax
  movzbl (%rax,%r13,1),%eax
  imul   %r12d,%eax
  cltq
  mov    %rbp,%rdx
  test   %r12d,%r12d
  jle    L331
  mov    %rbp,%rcx
L50:
  movq   $0x0,(%rcx)
  add    $0x8,%rcx
  cmp    %r11,%rcx
  jne    L50
  mov    (%rsp),%rdi
  movslq 0x10(%rsp),%rcx
  shl    $0x3,%rcx
  mov    %rcx,%rsi
  add    0xc0(%rsp),%rsi
  add    0xc8(%rsp),%rcx
  mov    0x20(%rsp),%r10
  lea    (%r10,%r8,8),%r9
  lea    (%r10,%rax,8),%r8
  mov    0x38(%rsp),%r10d
L52:
  mov    $0x0,%eax
  pxor   %xmm2,%xmm2
  movapd %xmm2,%xmm1
L384:
  movsd  (%r9,%rax,8),%xmm0
  mulsd  (%rsi,%rax,8),%xmm0
  addsd  %xmm0,%xmm1
  movsd  (%r8,%rax,8),%xmm0
  mulsd  (%rcx,%rax,8),%xmm0
  addsd  %xmm0,%xmm2
  mov    %rax,%r15
  add    $0x1,%rax
  cmp    %rbx,%r15
  jne    L384
  mulsd  %xmm2,%xmm1
  mov    $0x0,%eax
L401:
  movapd %xmm1,%xmm0
  mulsd  (%rdi,%rax,8),%xmm0
  addsd  (%rdx,%rax,8),%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  mov    %rax,%r15
  add    $0x1,%rax
  cmp    %rbx,%r15
  jne    L401
  add    $0x1,%r10d
  add    %r14,%rdi
  add    %r14,%rsi
  add    %r14,%rcx
  cmp    %r10d,%r12d
  jne    L52
L331:
  add    %r14,%rbp
  mov    0x28(%rsp),%edi
  add    %edi,0x10(%rsp)
  add    %r14,%r11
  subl   $0x1,0x30(%rsp)
  jne    L168
  jmp    L223
L153:
  comisd D58(%rip),%xmm0
  seta   %cl
  movzbl %cl,%ecx
L410:
  add    $0x1,%edx
  add    $0x8,%rax
  test   %ecx,%ecx
  je     L288
  cmp    %edx,%esi
  jle    L288
L333:
  movsd  (%rax),%xmm0
  pxor   %xmm3,%xmm3
  comisd %xmm0,%xmm3
  ja     L153
  movsd  D17(%rip),%xmm5
  comisd %xmm0,%xmm5
  seta   %cl
  movzbl %cl,%ecx
  jmp    L410
L47:
  mov    0x8(%rsp),%rax
  mov    0x68(%rsp),%rdx
L422:
  movsd  (%rax),%xmm0
  mulsd  D14(%rip),%xmm0
  movsd  %xmm0,(%rax)
  add    $0x8,%rax
  cmp    %rax,%rdx
  jne    L422
  jmp    L266
L328:
  comisd D58(%rip),%xmm0
  seta   %cl
  movzbl %cl,%ecx
L344:
  add    $0x1,%edx
  add    $0x8,%rax
  test   %ecx,%ecx
  je     L459
  cmp    %edx,%esi
  jle    L459
L182:
  movsd  (%rax),%xmm0
  pxor   %xmm4,%xmm4
  comisd %xmm0,%xmm4
  ja     L328
  movsd  D17(%rip),%xmm6
  comisd %xmm0,%xmm6
  seta   %cl
  movzbl %cl,%ecx
  jmp    L344
L172:
  mov    0x8(%rsp),%rax
  mov    0x50(%rsp),%rdx
L318:
  movsd  (%rax),%xmm0
  mulsd  D14(%rip),%xmm0
  movsd  %xmm0,(%rax)
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L318
  jmp    L225
L288:
  test   %ecx,%ecx
  jne    L47
L266:
  mov    0x2c(%rsp),%r9d
  mov    %r12d,%r8d
  mov    0x48(%rsp),%ecx
  mov    0xd8(%rsp),%rdx
  mov    0xe0(%rsp),%rsi
  mov    0x8(%rsp),%rdi
  call   <ancestralGamma@plt>
  addq   $0x1,0x10(%rsp)
  mov    0x10(%rsp),%rax
  mov    0x58(%rsp),%rsi
  add    %rsi,0x38(%rsp)
  cmp    0x60(%rsp),%rax
  je     L273
L78:
  mov    0x10(%rsp),%eax
  mov    %eax,0x48(%rsp)
  mov    0x40(%rsp),%rax
  mov    0x38(%rsp),%rdi
  lea    (%rax,%rdi,8),%r8
  mov    0x50(%rsp),%r11
  mov    0x8(%rsp),%rbp
  mov    $0x0,%r14d
  movl   $0x4,0x30(%rsp)
L195:
  mov    0x18(%rsp),%rax
  mov    0x10(%rsp),%rdi
  movzbl (%rax,%rdi,1),%eax
  imul   %r12d,%eax
  cltq
  mov    %rbp,%rdx
  test   %r12d,%r12d
  jle    L365
  mov    %rbp,%rcx
L488:
  movq   $0x0,(%rcx)
  add    $0x8,%rcx
  cmp    %r11,%rcx
  jne    L488
  mov    (%rsp),%rdi
  movslq %r14d,%rcx
  shl    $0x3,%rcx
  mov    %rcx,%rsi
  add    0xc0(%rsp),%rsi
  add    0xc8(%rsp),%rcx
  mov    0x20(%rsp),%r10
  lea    (%r10,%rax,8),%r9
  mov    $0x0,%r10d
L283:
  mov    $0x0,%eax
  pxor   %xmm2,%xmm2
  movapd %xmm2,%xmm1
L152:
  movsd  (%r9,%rax,8),%xmm0
  mulsd  (%rsi,%rax,8),%xmm0
  addsd  %xmm0,%xmm1
  movsd  (%r8,%rax,8),%xmm0
  mulsd  (%rcx,%rax,8),%xmm0
  addsd  %xmm0,%xmm2
  mov    %rax,%r15
  add    $0x1,%rax
  cmp    %rbx,%r15
  jne    L152
  mulsd  %xmm2,%xmm1
  mov    $0x0,%eax
L366:
  movapd %xmm1,%xmm0
  mulsd  (%rdi,%rax,8),%xmm0
  addsd  (%rdx,%rax,8),%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  mov    %rax,%r15
  add    $0x1,%rax
  cmp    %rbx,%r15
  jne    L366
  add    $0x1,%r10d
  add    %r13,%rdi
  add    %r13,%rsi
  add    %r13,%rcx
  cmp    %r10d,%r12d
  jne    L283
L365:
  add    %r13,%rbp
  add    %r13,%r8
  add    0x28(%rsp),%r14d
  add    %r13,%r11
  subl   $0x1,0x30(%rsp)
  jne    L195
  mov    0x2c(%rsp),%esi
  test   %esi,%esi
  jle    L266
  mov    0x8(%rsp),%rax
  mov    $0x0,%edx
  jmp    L333
L459:
  test   %ecx,%ecx
  jne    L172
L225:
  mov    0x2c(%rsp),%r9d
  mov    %r12d,%r8d
  mov    0x18(%rsp),%r14d
  mov    %r14d,%ecx
  mov    0xd8(%rsp),%rdx
  mov    0xe0(%rsp),%rsi
  mov    0x8(%rsp),%rdi
  call   <ancestralGamma@plt>
  mov    %r14d,%eax
  add    $0x1,%eax
  mov    %eax,0x18(%rsp)
  mov    0x38(%rsp),%rsi
  add    %rsi,0x20(%rsp)
  cmp    %eax,0xb8(%rsp)
  je     L273
L164:
  mov    0x20(%rsp),%rax
  lea    0x0(,%rax,8),%r8
  mov    0x48(%rsp),%rax
  lea    (%rax,%r8,1),%r9
  add    0x40(%rsp),%r8
  mov    0x30(%rsp),%r11
  mov    0x8(%rsp),%rbp
  mov    $0x0,%r13d
  movl   $0x4,0x10(%rsp)
L258:
  mov    %rbp,%rdx
  test   %r12d,%r12d
  jle    L340
  mov    %rbp,%rax
L96:
  movq   $0x0,(%rax)
  add    $0x8,%rax
  cmp    %r11,%rax
  jne    L96
  mov    (%rsp),%rdi
  movslq %r13d,%rcx
  shl    $0x3,%rcx
  mov    %rcx,%rsi
  add    0xc0(%rsp),%rsi
  add    0xc8(%rsp),%rcx
  mov    $0x0,%r10d
L218:
  mov    $0x0,%eax
  pxor   %xmm2,%xmm2
  movapd %xmm2,%xmm1
L101:
  movsd  (%r9,%rax,8),%xmm0
  mulsd  (%rsi,%rax,8),%xmm0
  addsd  %xmm0,%xmm1
  movsd  (%r8,%rax,8),%xmm0
  mulsd  (%rcx,%rax,8),%xmm0
  addsd  %xmm0,%xmm2
  mov    %rax,%r14
  add    $0x1,%rax
  cmp    %rbx,%r14
  jne    L101
  mulsd  %xmm2,%xmm1
  mov    $0x0,%eax
L58:
  movapd %xmm1,%xmm0
  mulsd  (%rdi,%rax,8),%xmm0
  addsd  (%rdx,%rax,8),%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  mov    %rax,%r14
  add    $0x1,%rax
  cmp    %rbx,%r14
  jne    L58
  add    $0x1,%r10d
  add    %r15,%rdi
  add    %r15,%rsi
  add    %r15,%rcx
  cmp    %r10d,%r12d
  jne    L218
L340:
  add    %r15,%rbp
  add    %r15,%r9
  add    %r15,%r8
  add    0x28(%rsp),%r13d
  add    %r15,%r11
  subl   $0x1,0x10(%rsp)
  jne    L258
  mov    0x2c(%rsp),%esi
  test   %esi,%esi
  jle    L225
  mov    0x8(%rsp),%rax
  mov    $0x0,%edx
  jmp    L182
```

tool_calls: name: disassemble_data
arguments: data_label: D58
type_infered: f32[4]
name: disassemble_data
arguments: data_label: D17
type_infered: f32[2]
name: disassemble_data
arguments: data_label: D14
type_infered: f64

tool: {"D58": [0.0,-4.3655745685100555e-10,0.0,4.3655745685100555e-10]}

tool: {"D17": [0.0,4.3655745685100555e-10]}

tool: D14: 1.157920892373162e+77

assistant: ```c
void newviewFlexGamma_Ancestral(int tipCase, double *x1, double *x2, double *extEV, double *tipVector, unsigned char *tipX1, unsigned char *tipX2, int n, double *left, double *right, const int numStates, double *diagptable, double *sumBuffer) {
  double *v, *x3 = (double *)rax_malloc(sizeof(double) * 4 * numStates);
  double x1px2;
  int i, j, l, k, scale;
  double *vl, *vr, al, ar;
  const int statesSquare = numStates * numStates, gammaStates = 4 * numStates;
  switch (tipCase) {
  case 0: {
    for (i = 0; i < n; i++) {
      for (k = 0; k < 4; k++) {
        vl = &(tipVector[numStates * tipX1[i]]);
        vr = &(tipVector[numStates * tipX2[i]]);
        v = &(x3[numStates * k]);
        for (l = 0; l < numStates; l++)
          v[l] = 0;
        for (l = 0; l < numStates; l++) {
          al = 0.0;
          ar = 0.0;
          for (j = 0; j < numStates; j++) {
            al += vl[j] * left[k * statesSquare + l * numStates + j];
            ar += vr[j] * right[k * statesSquare + l * numStates + j];
          }
          x1px2 = al * ar;
          for (j = 0; j < numStates; j++)
            v[j] += x1px2 * extEV[numStates * l + j];
        }
      }
      ancestralGamma(x3, sumBuffer, diagptable, i, numStates, gammaStates);
    }
  } break;
  case 1: {
    for (i = 0; i < n; i++) {
      for (k = 0; k < 4; k++) {
        vl = &(tipVector[numStates * tipX1[i]]);
        vr = &(x2[gammaStates * i + numStates * k]);
        v = &(x3[numStates * k]);
        for (l = 0; l < numStates; l++)
          v[l] = 0;
        for (l = 0; l < numStates; l++) {
          al = 0.0;
          ar = 0.0;
          for (j = 0; j < numStates; j++) {
            al += vl[j] * left[k * statesSquare + l * numStates + j];
            ar += vr[j] * right[k * statesSquare + l * numStates + j];
          }
          x1px2 = al * ar;
          for (j = 0; j < numStates; j++)
            v[j] += x1px2 * extEV[numStates * l + j];
        }
      }
      v = x3;
      scale = 1;
      for (l = 0; scale && (l < gammaStates); l++)
        scale = ((((v[l]) < 0) ? (-(v[l])) : (v[l])) < (1.0 / 115792089237316195423570985008687907853269984665640564039457584007913129639936.0));
      if (scale) {
        for (l = 0; l < gammaStates; l++)
          v[l] *= 115792089237316195423570985008687907853269984665640564039457584007913129639936.0;
      }
      ancestralGamma(x3, sumBuffer, diagptable, i, numStates, gammaStates);
    }
  } break;
  case 2:
    for (i = 0; i < n; i++) {
      for (k = 0; k < 4; k++) {
        vl = &(x1[gammaStates * i + numStates * k]);
        vr = &(x2[gammaStates * i + numStates * k]);
        v = &(x3[numStates * k]);
        for (l = 0; l < numStates; l++)
          v[l] = 0;
        for (l = 0; l < numStates; l++) {
          al = 0.0;
          ar = 0.0;
          for (j = 0; j < numStates; j++) {
            al += vl[j] * left[k * statesSquare + l * numStates + j];
            ar += vr[j] * right[k * statesSquare + l * numStates + j];
          }
          x1px2 = al * ar;
          for (j = 0; j < numStates; j++)
            v[j] += x1px2 * extEV[numStates * l + j];
        }
      }
      v = x3;
      scale = 1;
      for (l = 0; scale && (l < gammaStates); l++)
        scale = (((((v[l]) < 0) ? (-(v[l])) : (v[l])) < (1.0 / 115792089237316195423570985008687907853269984665640564039457584007913129639936.0)));
      if (scale) {
        for (l = 0; l < gammaStates; l++)
          v[l] *= 115792089237316195423570985008687907853269984665640564039457584007913129639936.0;
      }
      ancestralGamma(x3, sumBuffer, diagptable, i, numStates, gammaStates);
    }
    break;
  default:
    assert(0);
  }
  rax_free(x3);
}
```