user: Please show me the C implementation that matches the assembly code listed below:

```asm
<simplex2>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x48,%rsp
  mov    %edi,%r12d
  mov    %esi,%ebx
  mov    %rdx,%r13
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  lea    0x1(%rsi),%eax
  mov    %eax,-0x54(%rbp)
  mov    %eax,%eax
  mov    %rax,-0x60(%rbp)
  lea    0x1(%rdi),%r14d
  mov    %esi,%edi
  mov    %rdi,-0x68(%rbp)
  shl    $0x2,%rdi
  mov    %rdi,-0x50(%rbp)
  lea    0x0(,%rax,4),%r15
  jmp    L146
L50:
  cmp    %r8d,%ebx
  je     L66
  cmp    $0x1,%r14d
  jbe    L66
  mov    %r8d,%ecx
  mov    -0x60(%rbp),%rax
  add    %rcx,%rax
  lea    0x0(%r13,%rax,4),%rax
  mov    $0x1,%edx
  pxor   %xmm2,%xmm2
  mov    $0xffffffff,%esi
  mov    -0x68(%rbp),%rdi
  sub    %rcx,%rdi
  jmp    L127
L49:
  mov    %edx,%esi
  movaps %xmm1,%xmm2
L4:
  lea    0x1(%rdx),%ecx
  add    %r15,%rax
  cmp    %edx,%r12d
  je     L33
  mov    %ecx,%edx
L127:
  movss  (%rax),%xmm0
  comiss D7(%rip),%xmm0
  jbe    L4
  movss  (%rax,%rdi,4),%xmm1
  divss  %xmm0,%xmm1
  cmp    $0xffffffff,%esi
  je     L49
  comiss %xmm1,%xmm2
  jbe    L4
  jmp    L49
L33:
  cmp    $0xffffffff,%esi
  je     L66
  mov    %esi,%ecx
  mov    %r13,%rdx
  mov    -0x54(%rbp),%esi
  mov    %r14d,%edi
  call   <trafo@plt>
  mov    -0x50(%rbp),%rax
  add    $0xf,%rax
  shr    $0x4,%rax
  shl    $0x4,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L137:
  cmp    %rdx,%rsp
  je     L89
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L137
L89:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  je     L19
  orq    $0x0,-0x8(%rsp,%rdx,1)
L19:
  mov    %rsp,%rax
  mov    %r13,%rcx
  mov    %rax,-0x70(%rbp)
  mov    %rax,%rdx
  mov    %ebx,%esi
  mov    %r12d,%edi
  call   <solution@plt>
  mov    %r13,%rcx
  mov    -0x70(%rbp),%rdx
  mov    %ebx,%esi
  mov    %r12d,%edi
  call   <feasible_p@plt>
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0x48(%rbp),%rsp
L146:
  mov    %rsp,-0x48(%rbp)
  test   %ebx,%ebx
  je     L66
  mov    %r13,%rax
  mov    $0x0,%r8d
L36:
  pxor   %xmm3,%xmm3
  comiss (%rax),%xmm3
  ja     L50
  add    $0x1,%r8d
  add    $0x4,%rax
  cmp    %r8d,%ebx
  jne    L36
L66:
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L116
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
L116:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: tag: D7
datatype: f32

tool: D7: 0

assistant: ```c
void simplex2(unsigned int D, unsigned int N, float A[D + 1][N + 1]) {
  while (true) {
    unsigned int i = 0;
    for (i = 0; i < N; i++)
      if (A[0][i] < 0.)
        break;
    if (i == N)
      break;
    int pivot_index = -1;
    float pivot_value = 0.;
    for (unsigned int j = 1; j < D + 1; j++) {
      if (0. < A[j][i]) {
        float nval = A[j][N] / A[j][i];
        if ((-1 == pivot_index) || (nval < pivot_value)) {
          pivot_value = nval;
          pivot_index = j;
        }
      }
    }
    if (-1 == pivot_index)
      break;
    trafo(D + 1, N + 1, A, pivot_index, i);
    float x[N];
    solution(D, N, x, A);
    assert(feasible_p(D, N, x, A));
  }
}
```