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

```asm
<simplex_contract>:
  endbr64
  test   %edx,%edx
  js     L83
  imul   %edx,%ecx
  mov    %edx,%r11d
  mov    %ecx,%r8d
  test   %edx,%edx
  je     L83
  push   %r15
  movslq %edx,%rax
  mov    %rdi,%rsi
  movss  D0(%rip),%xmm1
  push   %r14
  movaps D0(%rip),%xmm2
  xor    %r14d,%r14d
  movq   D0(%rip),%xmm3
  push   %r13
  xor    %r13d,%r13d
  push   %r12
  push   %rbp
  movslq %r8d,%rbp
  push   %rbx
  mov    %ecx,%ebx
  mov    %rdi,%rcx
  lea    -0x1(%rdx),%edi
  not    %ebx
  lea    (%rsi,%rbp,4),%r10
  mov    %rax,-0x38(%rsp)
  shl    $0x2,%rax
  mov    %rax,-0x30(%rsp)
  lea    0x1(%rdx),%eax
  mov    %eax,-0x24(%rsp)
  mov    %edx,%eax
  shr    $0x2,%eax
  mov    %edi,-0x1c(%rsp)
  lea    -0x1(%rax),%r9d
  mov    %edx,%eax
  mov    %rbp,-0x40(%rsp)
  and    $0xfffffffc,%eax
  add    $0x1,%r9
  mov    %eax,-0x18(%rsp)
  mov    %edi,%eax
  lea    0x1(%r8),%edi
  shl    $0x4,%r9
  mov    %edi,-0x20(%rsp)
  mov    %edx,%edi
  lea    0x4(,%rax,4),%rax
  and    $0x3,%edi
  mov    %rax,-0x10(%rsp)
  mov    %edi,-0x14(%rsp)
  nopl   0x0(%rax,%rax,1)
L127:
  mov    -0x20(%rsp),%eax
  lea    (%rax,%rbx,1),%edx
  cmp    $0x2,%ebx
  jbe    L84
  cmp    $0x1,%r11d
  je     L84
  cmpl   $0x2,-0x1c(%rsp)
  jbe    L1
  xor    %eax,%eax
  nop
L98:
  movups (%rcx,%rax,1),%xmm0
  movups (%r10,%rax,1),%xmm6
  addps  %xmm0,%xmm0
  addps  %xmm6,%xmm0
  mulps  %xmm2,%xmm0
  movups %xmm0,(%rcx,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%r9
  jne    L98
  mov    -0x18(%rsp),%eax
  cmp    %eax,%r11d
  je     L119
  mov    -0x14(%rsp),%edi
  mov    %edi,%r15d
  cmp    $0x1,%edi
  je     L125
  mov    %eax,%edi
L108:
  lea    0x0(%r13,%rdi,1),%r12
  add    -0x40(%rsp),%rdi
  lea    (%rsi,%r12,4),%r12
  movq   (%rsi,%rdi,4),%xmm0
  mov    %r15d,%edi
  movq   (%r12),%xmm4
  and    $0xfffffffe,%edi
  add    %edi,%eax
  addps  %xmm4,%xmm4
  addps  %xmm4,%xmm0
  mulps  %xmm3,%xmm0
  movlps %xmm0,(%r12)
  cmp    %r15d,%edi
  je     L119
L125:
  add    %eax,%edx
  add    %r8d,%eax
  movslq %edx,%rdx
  cltq
  movss  (%rsi,%rdx,4),%xmm0
  addss  %xmm0,%xmm0
  addss  (%rsi,%rax,4),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,(%rsi,%rdx,4)
L119:
  mov    -0x38(%rsp),%rax
  add    $0x1,%r14d
  add    %r11d,%ebx
  add    -0x30(%rsp),%rcx
  sub    %rax,%rbp
  add    %rax,%r13
  cmp    -0x24(%rsp),%r14d
  jne    L127
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L84:
  mov    -0x10(%rsp),%rax
  lea    (%rax,%rcx,1),%rdx
  mov    %rcx,%rax
  nopl   0x0(%rax)
L99:
  movss  (%rax),%xmm0
  addss  %xmm0,%xmm0
  addss  (%rax,%rbp,4),%xmm0
  add    $0x4,%rax
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x4(%rax)
  cmp    %rdx,%rax
  jne    L99
  jmp    L119
L1:
  mov    %r11d,%r15d
  xor    %edi,%edi
  xor    %eax,%eax
  jmp    L108
L83:
  ret
```

tool_calls: name: extract_symbolic_data
arguments: label: D0
datatype: byte[16]

tool: {"D0": [0xab,0xaa,0xaa,0x3e,0xab,0xaa,0xaa,0x3e,0xab,0xaa,0xaa,0x3e,0xab,0xaa,0xaa,0x3e]}

assistant: ```c
simplex_contract(pt, rf, n, best, worst, howmuch) int n, best, worst;
float howmuch;
float (*pt)[], (*rf)[];
{
  int i, j;
  float r1, r2, r3, randf();
  for (i = 0; i < n + 1; i++) {
    r1 = 2.;
    r2 = 1.;
    r3 = 1. / (r1 + r2);
    for (j = 0; j < n; j++)
      (*pt)[i * n + j] = r3 * (r1 * (*pt)[i * n + j] + r2 * (*pt)[n * best + j]);
  }
}
```