user: What is the c source code of the assembly code below:

```asm
<simplx>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %rdi,0x10(%rsp)
  mov    %rsi,%rbp
  mov    %edx,%r13d
  mov    %rcx,%r12
  mov    %r8d,0x3c(%rsp)
  movsd  %xmm0,0x28(%rsp)
  movsd  %xmm1,0x30(%rsp)
  mov    %r9,0x20(%rsp)
  mov    0xc0(%rsp),%r15
  movslq %r8d,%rax
  lea    0x0(,%rax,8),%r14
  mov    %r14,%rdi
  call   <malloc@plt>
  test   %rax,%rax
  je     L627
  mov    %rax,%rbx
  mov    %r14,%rdi
  call   <malloc@plt>
  mov    %rax,0x18(%rsp)
  test   %rax,%rax
  je     L627
  lea    0x1(%r13),%eax
  mov    %eax,0x38(%rsp)
  mov    0xb0(%rsp),%rax
  movl   $0x0,(%rax)
  lea    -0x1(%r13),%r14d
  lea    0x8(%rbx,%r14,8),%rax
  mov    %rax,0x40(%rsp)
  jmp    L621
L627:
  lea    D67(%rip),%rdi
  call   <puts@plt>
  mov    $0x0,%eax
  jmp    L380
L377:
  movslq %r8d,%rcx
  comisd 0x0(%rbp,%rcx,8),%xmm0
  jbe    L478
  cmp    %eax,%edx
  cmovne %esi,%r8d
L478:
  lea    0x1(%rax),%rdx
  cmp    %rdi,%rax
  je     L405
  mov    %rdx,%rax
L162:
  mov    %eax,%esi
  movslq (%r9),%rdx
  movsd  0x0(%rbp,%rdx,8),%xmm0
  comisd 0x0(%rbp,%rax,8),%xmm0
  jbe    L18
  mov    %eax,(%r9)
L18:
  movsd  0x0(%rbp,%rax,8),%xmm0
  mov    (%r15),%edx
  movslq %edx,%rcx
  comisd 0x0(%rbp,%rcx,8),%xmm0
  jbe    L377
  mov    %esi,(%r15)
  mov    %edx,%r8d
  jmp    L478
L405:
  mov    %r8d,0x4(%rsp)
L97:
  movslq (%r15),%rax
  movsd  0x0(%rbp,%rax,8),%xmm0
  mov    0xb8(%rsp),%rax
  movslq (%rax),%rax
  movsd  0x0(%rbp,%rax,8),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  andpd  D48(%rip),%xmm1
  addsd  %xmm1,%xmm1
  andpd  D48(%rip),%xmm0
  movapd %xmm2,%xmm3
  andpd  D48(%rip),%xmm3
  addsd  %xmm3,%xmm0
  mulsd  0x28(%rsp),%xmm0
  comisd %xmm1,%xmm0
  jae    L114
  movsd  0x30(%rsp),%xmm5
  comisd %xmm2,%xmm5
  ja     L178
  mov    0xb0(%rsp),%rax
  mov    (%rax),%eax
  cmp    0xa8(%rsp),%eax
  je     L161
  add    $0x1,%eax
  mov    0xb0(%rsp),%rdi
  mov    %eax,(%rdi)
  test   %r13d,%r13d
  jle    L156
  mov    %rbx,%rax
  mov    0x40(%rsp),%rdx
L496:
  movq   $0x0,(%rax)
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L496
  cmpl   $0x0,0x38(%rsp)
  jle    L524
L319:
  mov    (%r15),%r8d
  mov    $0x0,%edi
  mov    $0x0,%ecx
  mov    0x10(%rsp),%r9
  jmp    L299
L114:
  mov    %rbx,%rdi
  call   <free@plt>
  mov    0x18(%rsp),%rdi
  call   <free@plt>
  mov    $0x1,%eax
L380:
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L178:
  mov    %rbx,%rdi
  call   <free@plt>
  mov    0x18(%rsp),%rdi
  call   <free@plt>
  mov    $0x2,%eax
  jmp    L380
L161:
  mov    %rbx,%rdi
  call   <free@plt>
  mov    0x18(%rsp),%rdi
  call   <free@plt>
  mov    $0x3,%eax
  jmp    L380
L539:
  lea    0x1(%rcx),%eax
  add    %r13d,%edi
  cmp    %ecx,%r13d
  jle    L484
  mov    %eax,%ecx
L299:
  cmp    %ecx,%r8d
  je     L539
  test   %r13d,%r13d
  jle    L539
  movslq %edi,%rax
  lea    (%r9,%rax,8),%rsi
  mov    $0x0,%eax
L19:
  movsd  (%rbx,%rax,8),%xmm0
  addsd  (%rsi,%rax,8),%xmm0
  movsd  %xmm0,(%rbx,%rax,8)
  mov    %rax,%rdx
  add    $0x1,%rax
  cmp    %r14,%rdx
  jne    L19
  jmp    L539
L484:
  test   %r13d,%r13d
  jle    L50
L524:
  mov    $0x0,%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %r13d,%xmm1
  mov    0x10(%rsp),%rsi
  mov    0x20(%rsp),%rdi
L378:
  movsd  (%rbx,%rax,8),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rbx,%rax,8)
  movslq (%rdi,%rax,4),%rcx
  addsd  %xmm0,%xmm0
  mov    %r13d,%edx
  imul   (%r15),%edx
  add    %eax,%edx
  movslq %edx,%rdx
  subsd  (%rsi,%rdx,8),%xmm0
  movsd  %xmm0,(%r12,%rcx,8)
  add    $0x1,%rax
  cmp    %eax,%r13d
  jg     L378
L50:
  mov    0x3c(%rsp),%eax
  test   %eax,%eax
  jle    L90
  mov    %eax,%edx
  mov    $0x0,%eax
  mov    0x18(%rsp),%rcx
L43:
  movsd  (%r12,%rax,8),%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L43
L90:
  mov    %r12,%rdi
  mov    $0x0,%eax
  call   *0xa0(%rsp)
  movsd  %xmm0,0x8(%rsp)
  mov    0xb8(%rsp),%rax
  movslq (%rax),%rax
  movsd  0x0(%rbp,%rax,8),%xmm0
  comisd 0x8(%rsp),%xmm0
  ja     L164
  movslq 0x4(%rsp),%rax
  movsd  0x8(%rsp),%xmm7
  comisd 0x0(%rbp,%rax,8),%xmm7
  jae    L460
  mov    $0x0,%edx
  test   %r13d,%r13d
  jle    L495
  mov    0x10(%rsp),%rcx
  mov    0x20(%rsp),%rsi
L515:
  movslq (%rsi,%rdx,4),%rax
  movsd  (%r12,%rax,8),%xmm0
  mov    %r13d,%eax
  imul   (%r15),%eax
  add    %edx,%eax
  cltq
  movsd  %xmm0,(%rcx,%rax,8)
  mov    %rdx,%rax
  add    $0x1,%rdx
  cmp    %r14,%rax
  jne    L515
L495:
  movslq (%r15),%rax
  movsd  0x8(%rsp),%xmm6
  movsd  %xmm6,0x0(%rbp,%rax,8)
L621:
  mov    0xb8(%rsp),%rax
  movl   $0x0,(%rax)
  movsd  0x0(%rbp),%xmm0
  comisd 0x8(%rbp),%xmm0
  setbe  %al
  movzbl %al,%eax
  seta   %dil
  movzbl %dil,%edi
  mov    %edi,0x4(%rsp)
  mov    %eax,(%r15)
  mov    %r13d,%edi
  mov    $0x0,%eax
  cmpl   $0x0,0x38(%rsp)
  jle    L97
  mov    0x4(%rsp),%r8d
  mov    0xb8(%rsp),%r9
  jmp    L162
L164:
  test   %r13d,%r13d
  jle    L113
  mov    $0x0,%eax
  mov    0x18(%rsp),%rcx
  mov    0x20(%rsp),%rsi
L576:
  movslq (%rsi,%rax,4),%rdx
  movsd  (%r12,%rdx,8),%xmm0
  addsd  %xmm0,%xmm0
  subsd  (%rbx,%rax,8),%xmm0
  movsd  %xmm0,(%rcx,%rdx,8)
  mov    %rax,%rdx
  add    $0x1,%rax
  cmp    %r14,%rdx
  jne    L576
  mov    0x18(%rsp),%rdi
  mov    $0x0,%eax
  call   *0xa0(%rsp)
  movapd %xmm0,%xmm1
  mov    0xb8(%rsp),%rax
  movslq (%rax),%rax
  movsd  0x0(%rbp,%rax,8),%xmm0
  comisd %xmm1,%xmm0
  ja     L448
  mov    $0x0,%edx
  mov    0x10(%rsp),%rcx
  mov    0x20(%rsp),%rsi
L445:
  movslq (%rsi,%rdx,4),%rax
  movsd  (%r12,%rax,8),%xmm0
  mov    %r13d,%eax
  imul   (%r15),%eax
  add    %edx,%eax
  cltq
  movsd  %xmm0,(%rcx,%rax,8)
  mov    %rdx,%rax
  add    $0x1,%rdx
  cmp    %r14,%rax
  jne    L445
L499:
  movslq (%r15),%rax
  movsd  0x8(%rsp),%xmm7
  movsd  %xmm7,0x0(%rbp,%rax,8)
  jmp    L621
L448:
  mov    $0x0,%edx
  mov    0x18(%rsp),%rsi
  mov    0x10(%rsp),%rdi
  mov    0x20(%rsp),%r8
L476:
  mov    %r13d,%eax
  imul   (%r15),%eax
  add    %edx,%eax
  cltq
  movslq (%r8,%rdx,4),%rcx
  movsd  (%rsi,%rcx,8),%xmm0
  movsd  %xmm0,(%rdi,%rax,8)
  mov    %rdx,%rax
  add    $0x1,%rdx
  cmp    %r14,%rax
  jne    L476
L76:
  movslq (%r15),%rax
  movsd  %xmm1,0x0(%rbp,%rax,8)
  jmp    L621
L460:
  movslq (%r15),%rax
  movsd  0x0(%rbp,%rax,8),%xmm0
  comisd 0x8(%rsp),%xmm0
  ja     L277
L47:
  test   %r13d,%r13d
  jle    L477
  mov    %r13d,%eax
  imul   (%r15),%eax
  cltq
  mov    0x10(%rsp),%rdi
  lea    (%rdi,%rax,8),%rcx
  mov    $0x0,%eax
  mov    0x18(%rsp),%rsi
  mov    0x20(%rsp),%rdi
L349:
  movslq (%rdi,%rax,4),%rdx
  movsd  (%rcx,%rax,8),%xmm0
  mulsd  D14(%rip),%xmm0
  movsd  (%rbx,%rax,8),%xmm1
  mulsd  D14(%rip),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rsi,%rdx,8)
  mov    %rax,%rdx
  add    $0x1,%rax
  cmp    %r14,%rdx
  jne    L349
  mov    0x18(%rsp),%rdi
  mov    $0x0,%eax
  call   *0xa0(%rsp)
  movapd %xmm0,%xmm1
  movslq (%r15),%rax
  movsd  0x0(%rbp,%rax,8),%xmm0
  comisd %xmm1,%xmm0
  ja     L149
L280:
  mov    %r13d,%eax
  mov    %rax,0x8(%rsp)
  movl   $0x0,0x4(%rsp)
  mov    $0x0,%edi
  cmpl   $0x0,0x38(%rsp)
  jle    L621
  mov    %rbx,0x58(%rsp)
  mov    %r14,0x50(%rsp)
  mov    %rdi,%r14
  mov    %rbp,0x48(%rsp)
  mov    0x20(%rsp),%rbx
  mov    0xb8(%rsp),%rbp
  mov    %r15,0xc0(%rsp)
  mov    0x10(%rsp),%r15
  jmp    L119
L277:
  test   %r13d,%r13d
  jle    L159
  mov    $0x0,%edx
  mov    0x10(%rsp),%rcx
  mov    0x20(%rsp),%rsi
L248:
  movslq (%rsi,%rdx,4),%rax
  movsd  (%r12,%rax,8),%xmm0
  mov    %r13d,%eax
  imul   (%r15),%eax
  add    %edx,%eax
  cltq
  movsd  %xmm0,(%rcx,%rax,8)
  mov    %rdx,%rax
  add    $0x1,%rdx
  cmp    %r14,%rax
  jne    L248
L159:
  movslq (%r15),%rax
  movsd  0x8(%rsp),%xmm5
  movsd  %xmm5,0x0(%rbp,%rax,8)
  jmp    L47
L149:
  mov    $0x0,%edx
  mov    0x18(%rsp),%rsi
  mov    0x10(%rsp),%rdi
  mov    0x20(%rsp),%r8
L185:
  mov    %r13d,%eax
  imul   (%r15),%eax
  add    %edx,%eax
  cltq
  movslq (%r8,%rdx,4),%rcx
  movsd  (%rsi,%rcx,8),%xmm0
  movsd  %xmm0,(%rdi,%rax,8)
  mov    %rdx,%rax
  add    $0x1,%rdx
  cmp    %r14,%rax
  jne    L185
L238:
  movslq (%r15),%rax
  movsd  %xmm1,0x0(%rbp,%rax,8)
  jmp    L621
L300:
  mov    %r12,%rdi
  mov    $0x0,%eax
  call   *0xa0(%rsp)
  mov    0x48(%rsp),%rax
  movsd  %xmm0,(%rax,%r14,8)
L52:
  lea    0x1(%r14),%rax
  add    %r13d,0x4(%rsp)
  cmp    0x8(%rsp),%r14
  je     L382
  mov    %rax,%r14
L119:
  cmp    %r14d,0x0(%rbp)
  je     L52
  test   %r13d,%r13d
  jle    L300
  movslq 0x4(%rsp),%rax
  lea    (%r15,%rax,8),%rcx
  mov    $0x0,%eax
  mov    0x50(%rsp),%rdi
L288:
  movslq (%rbx,%rax,4),%rsi
  mov    %r13d,%edx
  imul   0x0(%rbp),%edx
  add    %eax,%edx
  movslq %edx,%rdx
  movsd  (%r15,%rdx,8),%xmm0
  addsd  (%rcx,%rax,8),%xmm0
  mulsd  D14(%rip),%xmm0
  movsd  %xmm0,(%r12,%rsi,8)
  movslq (%rbx,%rax,4),%rdx
  movsd  (%r12,%rdx,8),%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  mov    %rax,%rdx
  add    $0x1,%rax
  cmp    %rdi,%rdx
  jne    L288
  jmp    L300
L382:
  mov    0x58(%rsp),%rbx
  mov    0x50(%rsp),%r14
  mov    0x48(%rsp),%rbp
  mov    0xc0(%rsp),%r15
  jmp    L621
L113:
  mov    0x18(%rsp),%rdi
  mov    $0x0,%eax
  call   *0xa0(%rsp)
  movapd %xmm0,%xmm1
  mov    0xb8(%rsp),%rax
  movslq (%rax),%rax
  movsd  0x0(%rbp,%rax,8),%xmm0
  comisd %xmm1,%xmm0
  ja     L76
  jmp    L499
L477:
  mov    0x18(%rsp),%rdi
  mov    $0x0,%eax
  call   *0xa0(%rsp)
  movapd %xmm0,%xmm1
  movslq (%r15),%rax
  movsd  0x0(%rbp,%rax,8),%xmm0
  comisd %xmm1,%xmm0
  ja     L238
  jmp    L280
L156:
  cmpl   $0x0,0x38(%rsp)
  jg     L319
  jmp    L50
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D67
type: string
name: extract_symbolic_data
arguments: data_label: D48
type: qword[2]
name: extract_symbolic_data
arguments: data_label: D14
type: f64

tool: D67: SIMPLX: Error - not enough memory

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

tool: D14: 0.5

assistant: ```c
int simplx(double *p, double *y, int ndim, double *pr, int na, double ftol, double atol, int *lista, double (*funk)(), int itmax, int *iter, int *ilo, int *ihi) {
  static double alpha = 1.0, beta = 0.5, gamma = 2.0;
  double *pbar, *prr;
  double yprr, ypr;
  int inhi, mpts;
  int i, j;
  if (!(pbar = (double *)malloc(na * sizeof(double))) || !(prr = (double *)malloc(na * sizeof(double)))) {
    printf("SIMPLX: Error - not enough memory\n");
    return (0);
  }
  mpts = ndim + 1;
  *iter = 0;
  while (1) {
    *ilo = 0;
    if (y[0] > y[1]) {
      *ihi = 0;
      inhi = 1;
    } else {
      *ihi = 1;
      inhi = 0;
    }
    for (i = 0; i < mpts; i++) {
      if (y[i] < y[*ilo]) {
        *ilo = i;
      }
      if (y[i] > y[*ihi]) {
        inhi = *ihi;
        *ihi = i;
      } else if (y[i] > y[inhi]) {
        if (i != *ihi) {
          inhi = i;
        }
      }
    }
    if (fabs(y[*ihi] - y[*ilo]) * (double)2.0 <= ftol * (fabs(y[*ihi]) + fabs(y[*ilo]))) {
      free(pbar);
      free(prr);
      return (1);
    }
    if (y[*ilo] < atol) {
      free(pbar);
      free(prr);
      return (2);
    }
    if (*iter == itmax) {
      free(pbar);
      free(prr);
      return (3);
    }
    ++(*iter);
    for (j = 0; j < ndim; j++) {
      *(pbar + j) = (double)0.0;
    }
    for (i = 0; i < mpts; i++) {
      if (i != *ihi) {
        for (j = 0; j < ndim; j++) {
          *(pbar + j) += p[i * ndim + j];
        }
      }
    }
    for (j = 0; j < ndim; j++) {
      pbar[j] /= (double)ndim;
      pr[lista[j]] = *(pbar + j) * ((double)1.0 + alpha) - alpha * p[*ihi * ndim + j];
    }
    for (j = 0; j < na; j++) {
      *(prr + j) = pr[j];
    }
    ypr = (*funk)(pr);
    if (ypr < y[*ilo]) {
      for (j = 0; j < ndim; j++) {
        *(prr + lista[j]) = pr[lista[j]] * gamma + *(pbar + j) * ((double)1.0 - gamma);
      }
      yprr = (*funk)(prr);
      if (yprr < y[*ilo]) {
        for (j = 0; j < ndim; j++) {
          p[*ihi * ndim + j] = *(prr + lista[j]);
        }
        y[*ihi] = yprr;
      } else {
        for (j = 0; j < ndim; j++) {
          p[*ihi * ndim + j] = pr[lista[j]];
        }
        y[*ihi] = ypr;
      }
    } else if (ypr >= y[inhi]) {
      if (ypr < y[*ihi]) {
        for (j = 0; j < ndim; j++) {
          p[*ihi * ndim + j] = pr[lista[j]];
        }
        y[*ihi] = ypr;
      }
      for (j = 0; j < ndim; j++) {
        *(prr + lista[j]) = p[*ihi * ndim + j] * beta + *(pbar + j) * ((double)1.0 - beta);
      }
      yprr = (*funk)(prr);
      if (yprr < y[*ihi]) {
        for (j = 0; j < ndim; j++) {
          p[*ihi * ndim + j] = *(prr + lista[j]);
        }
        y[*ihi] = yprr;
      } else {
        for (i = 0; i < mpts; i++) {
          if (i != *ilo) {
            for (j = 0; j < ndim; j++) {
              pr[lista[j]] = (p[i * ndim + j] + p[*ilo * ndim + j]) * (double)0.5;
              p[i * ndim + j] = pr[lista[j]];
            }
            y[i] = (*funk)(pr);
          }
        }
      }
    } else {
      for (j = 0; j < ndim; j++) {
        p[*ihi * ndim + j] = pr[lista[j]];
      }
      y[*ihi] = ypr;
    }
  }
}
```