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

```asm
<simplex>:
  endbr64
  push   %rbp
  mov    %rcx,%r11
  mov    %rsp,%rbp
  push   %r15
  mov    %r9,%r15
  push   %r14
  lea    (%rsi,%rdi,1),%r14d
  push   %r13
  push   %r12
  mov    %edi,%r12d
  lea    0x1(%r14),%edi
  push   %rbx
  mov    %esi,%ebx
  sub    $0x38,%rsp
  mov    %rdx,-0x50(%rbp)
  lea    0x1(%r12),%edx
  imul   %rdi,%rdx
  mov    %r8,-0x58(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  lea    0xf(,%rdx,4),%rdx
  mov    %esi,%eax
  mov    %rsp,%rsi
  mov    %rdx,%rcx
  and    $0xfffffffffffff000,%rdx
  mov    %rax,-0x48(%rbp)
  mov    %rdi,%rax
  sub    %rdx,%rsi
  and    $0xfffffffffffffff0,%rcx
  cmp    %rsi,%rsp
  je     L148
L45:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rsi,%rsp
  jne    L45
L148:
  and    $0xfff,%ecx
  sub    %rcx,%rsp
  test   %rcx,%rcx
  jne    L29
L212:
  mov    -0x48(%rbp),%rsi
  mov    %rsp,%r13
  xor    %r9d,%r9d
  xor    %r8d,%r8d
  lea    0x0(,%rsi,4),%r10
  test   %eax,%eax
  je     L50
  movss  D14(%rip),%xmm1
  movss  D5(%rip),%xmm2
  cmp    %r8d,%ebx
  jbe    L236
  nopl   0x0(%rax,%rax,1)
L125:
  movss  (%r11,%r9,1),%xmm0
  lea    (%r15,%r9,1),%rcx
  xor    %edx,%edx
  mov    %r8d,%esi
  xorps  %xmm2,%xmm0
  movss  %xmm0,0x0(%r13,%r9,1)
  test   %r12d,%r12d
  je     L3
  xchg   %ax,%ax
L78:
  lea    0x1(%rdx),%eax
  movss  (%rcx),%xmm0
  add    %r10,%rcx
  mov    %rax,%rdx
  imul   %rdi,%rax
  add    %rsi,%rax
  movss  %xmm0,0x0(%r13,%rax,4)
  cmp    %edx,%r12d
  jne    L78
L3:
  lea    0x1(%r8),%eax
  add    $0x4,%r9
  cmp    %r8d,%r14d
  je     L50
  mov    %eax,%r8d
  cmp    %r8d,%ebx
  ja     L125
L236:
  movl   $0x0,0x0(%r13,%r9,1)
  test   %r12d,%r12d
  je     L3
  cmp    %r8d,%r14d
  je     L51
  mov    %r8d,%ecx
  xor    %eax,%eax
  mov    %r8d,%esi
  sub    %ebx,%ecx
  lea    0x1(%rax),%edx
  cmp    %ecx,%eax
  mov    %edx,%eax
  je     L238
  nop
L211:
  imul   %rdi,%rax
  add    %rsi,%rax
  movl   $0x0,0x0(%r13,%rax,4)
  cmp    %r12d,%edx
  je     L3
  mov    %edx,%eax
L69:
  lea    0x1(%rax),%edx
  cmp    %ecx,%eax
  mov    %edx,%eax
  jne    L211
L238:
  imul   %rdi,%rax
  add    %rsi,%rax
  movss  %xmm1,0x0(%r13,%rax,4)
  cmp    %r12d,%edx
  je     L3
  mov    %edx,%eax
  jmp    L69
  nopw   0x0(%rax,%rax,1)
L50:
  mov    %r13,%rdx
  mov    %r14d,%esi
  mov    %r12d,%edi
  call   <simplex2@plt>
  mov    %r14d,%eax
  mov    %rsp,%rdx
  lea    0xf(,%rax,4),%rax
  shr    $0x4,%rax
  shl    $0x4,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  sub    %rcx,%rdx
  cmp    %rdx,%rsp
  je     L226
L64:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rdx,%rsp
  jne    L64
L226:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  jne    L109
L117:
  mov    %rsp,%r15
  mov    %r13,%rcx
  mov    %r14d,%esi
  mov    %r12d,%edi
  mov    %r15,%rdx
  call   <solution@plt>
  test   %ebx,%ebx
  je     L53
  mov    -0x48(%rbp),%rdx
  mov    -0x50(%rbp),%rdi
  mov    %r15,%rsi
  shl    $0x2,%rdx
  call   <memcpy@plt>
L53:
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L153
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
  cs nopw 0x0(%rax,%rax,1)
L51:
  mov    -0x58(%rbp),%rcx
  xor    %edx,%edx
  mov    %r14d,%esi
  nopl   0x0(%rax)
L179:
  lea    0x1(%rdx),%eax
  movss  (%rcx),%xmm0
  add    $0x4,%rcx
  mov    %rax,%rdx
  imul   %rdi,%rax
  add    %rsi,%rax
  movss  %xmm0,0x0(%r13,%rax,4)
  cmp    %r12d,%edx
  jne    L179
  jmp    L3
L109:
  orq    $0x0,-0x8(%rsp,%rax,1)
  jmp    L117
L29:
  orq    $0x0,-0x8(%rsp,%rcx,1)
  jmp    L212
L153:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: data_label: D14
type_infered: f32
name: parse_data
arguments: data_label: D5
type_infered: byte[4]

tool: D14: 1

tool: {"D5": [0x00,0x00,0x00,0x80]}

assistant: ```c
void(simplex)(unsigned int D, unsigned int N, float x[N], const float c[N], const float b[D], const float A[D][N]) {
  float A2[D + 1][N + D + 1];
  for (unsigned int i = 0; i < N + D + 1; i++) {
    A2[0][i] = (i < N) ? -c[i] : 0.;
    for (unsigned int j = 0; j < D; j++) {
      if (i < N)
        A2[1 + j][i] = A[j][i];
      else if (i == N + D)
        A2[1 + j][i] = b[j];
      else
        A2[1 + j][i] = (i - N == j) ? 1. : 0.;
    }
  }
  simplex2(D, N + D, A2);
  float x2[D + N];
  solution(D, D + N, x2, A2);
  for (unsigned int i = 0; i < N; i++)
    x[i] = x2[i];
}
```