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

```asm
<init_array>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  test   %edi,%edi
  jle    L41
  pxor   %xmm2,%xmm2
  lea    -0x1(%rdi),%r15d
  mov    %rsi,0x48(%rsp)
  mov    %edi,%ebx
  cvtsi2sd %edi,%xmm2
  mov    %r15d,0xc(%rsp)
  mov    %rsi,%r12
  lea    0x8(%rsi),%rcx
  movsd  D19(%rip),%xmm1
  lea    0x2(%r15),%rbp
  mov    %rsi,%r14
  mov    $0x1,%r13d
  xchg   %ax,%ax
L11:
  xor    %esi,%esi
  nopw   0x0(%rax,%rax,1)
L217:
  mov    %esi,%eax
  pxor   %xmm0,%xmm0
  neg    %eax
  cltd
  idiv   %ebx
  cvtsi2sd %edx,%xmm0
  divsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%r14,%rsi,8)
  add    $0x1,%rsi
  cmp    %r13,%rsi
  jne    L217
  add    $0x1,%r13
  movsd  %xmm2,0x10(%rsp)
  add    $0x3e80,%r14
  cmp    %rbp,%r13
  je     L71
  mov    %ebx,%eax
  mov    %rcx,%rdi
  xor    %esi,%esi
  sub    %r13d,%eax
  lea    0x8(,%rax,8),%rdx
  call   <memset@plt>
  movsd  0x10(%rsp),%xmm2
  mov    %rax,%rcx
  mov    D19(%rip),%rax
  add    $0x3e88,%rcx
  mov    %rax,-0x3e90(%rcx)
  movq   %rax,%xmm1
  jmp    L11
L71:
  movslq 0xc(%rsp),%rax
  mov    $0x8,%esi
  mov    $0x3d0900,%edi
  imul   $0x3e80,%r15,%r15
  lea    (%r12,%rax,8),%rax
  movsd  %xmm1,(%rax,%r15,1)
  call   <polybench_alloc_data@plt>
  mov    %rax,%rbp
  mov    %rax,0x40(%rsp)
  lea    0x3e80(%rax,%r15,1),%r14
  mov    %ebx,%eax
  mov    %rax,0x28(%rsp)
  lea    0x0(,%rax,8),%r13
  mov    %rbp,%rcx
  nopw   0x0(%rax,%rax,1)
L101:
  mov    %rcx,%rdi
  mov    %r13,%rdx
  xor    %esi,%esi
  call   <memset@plt>
  mov    %rax,%rcx
  add    $0x3e80,%rcx
  cmp    %r14,%rcx
  jne    L101
  movslq %ebx,%rax
  lea    0x3e80(%r12,%r15,1),%r15
  movq   $0x0,0x20(%rsp)
  imul   $0x3e80,%rax,%rdx
  lea    0x0(%rbp,%rax,8),%rax
  mov    %rax,0x30(%rsp)
  mov    %ebx,%eax
  shr    %eax
  sub    $0x1,%eax
  lea    -0x3e78(%r12,%rdx,1),%r8
  shl    $0x4,%rax
  mov    %rax,%rdi
  mov    %rax,0x38(%rsp)
  mov    %ebx,%eax
  and    $0xfffffffe,%eax
  add    $0x10,%rdi
  mov    %eax,%r11d
  mov    %eax,0x10(%rsp)
  movslq %eax,%r13
  imul   $0x3e80,%r11,%r11
  mov    %rdi,0x18(%rsp)
  add    %r12,%r11
  xchg   %ax,%ax
L47:
  mov    0x30(%rsp),%r14
  movslq 0x20(%rsp),%r10
  mov    %rbp,%rdi
  mov    %r12,%rcx
  xor    %esi,%esi
  nopw   0x0(%rax,%rax,1)
L65:
  cmp    %r8,%rdi
  mov    %r12,%rdx
  setae  %r9b
  cmp    %r12,%r14
  setbe  %al
  or     %r9d,%eax
  cmpl   $0x1,0xc(%rsp)
  seta   %r9b
  test   %r9b,%al
  je     L113
  lea    0x8(%rcx),%rax
  cmp    %rax,%rdi
  setae  %r9b
  cmp    %rcx,%r14
  setbe  %al
  or     %r9b,%al
  je     L113
  mov    0x18(%rsp),%rax
  movsd  (%rcx),%xmm1
  lea    (%rax,%rdi,1),%r9
  unpcklpd %xmm1,%xmm1
  mov    %rdi,%rax
  nopl   0x0(%rax,%rax,1)
L151:
  movsd  (%rdx),%xmm0
  movupd (%rax),%xmm3
  add    $0x10,%rax
  add    $0x7d00,%rdx
  movhpd -0x3e80(%rdx),%xmm0
  mulpd  %xmm1,%xmm0
  addpd  %xmm3,%xmm0
  movups %xmm0,-0x10(%rax)
  cmp    %r9,%rax
  jne    L151
  cmp    %ebx,0x10(%rsp)
  je     L134
  movslq %esi,%rax
  movsd  (%rcx),%xmm0
  imul   $0x7d0,%rax,%rax
  mulsd  (%r11,%r10,8),%xmm0
  add    %r13,%rax
  addsd  0x0(%rbp,%rax,8),%xmm0
  movsd  %xmm0,0x0(%rbp,%rax,8)
L134:
  lea    0x1(%rsi),%eax
  add    $0x3e80,%rcx
  add    $0x3e80,%rdi
  add    $0x3e80,%r14
  cmp    %eax,%ebx
  je     L195
  mov    %eax,%esi
  jmp    L65
  nopl   0x0(%rax,%rax,1)
L113:
  mov    %r12,%rdx
  mov    %rdi,%rax
  cs nopw 0x0(%rax,%rax,1)
L147:
  movsd  (%rcx),%xmm0
  mulsd  (%rdx),%xmm0
  add    $0x3e80,%rdx
  add    $0x8,%rax
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %r15,%rdx
  jne    L147
  jmp    L134
  nopl   0x0(%rax)
L195:
  addq   $0x1,0x20(%rsp)
  add    $0x8,%r12
  add    $0x8,%r15
  add    $0x8,%r8
  mov    0x20(%rsp),%rax
  cmp    %rax,0x28(%rsp)
  jne    L47
  mov    %rbp,%r10
  mov    $0x8,%r11d
  mov    0x38(%rsp),%r12
  xor    %ecx,%ecx
  neg    %r10
  mov    %r11,%rax
  mov    0x48(%rsp),%rdi
  movslq 0x10(%rsp),%rdx
  sub    %r10,%rax
  mov    0x28(%rsp),%r14
  mov    0x40(%rsp),%r13
  add    $0x10,%r12
  mov    0xc(%rsp),%r15d
  cmp    %rax,%rdi
  je     L63
  nopl   (%rax)
L162:
  lea    (%r12,%r13,1),%r9
  mov    %r13,%rax
  lea    (%rdi,%r10,1),%r8
  cmp    $0x2,%r15d
  jbe    L63
  nopl   0x0(%rax)
L146:
  movupd (%rax),%xmm4
  movups %xmm4,(%r8,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%r9
  jne    L146
  cmp    %ebx,0x10(%rsp)
  je     L215
  movslq %ecx,%rax
  imul   $0x7d0,%rax,%rax
  add    %rdx,%rax
  movsd  0x0(%rbp,%rax,8),%xmm0
  movsd  %xmm0,(%rdi,%rdx,8)
L215:
  lea    0x1(%rcx),%eax
  add    $0x3e80,%rdi
  add    $0x3e80,%r13
  sub    $0x3e80,%r10
  cmp    %esi,%ecx
  je     L174
  mov    %eax,%ecx
  mov    %r11,%rax
  sub    %r10,%rax
  cmp    %rax,%rdi
  jne    L162
L63:
  xor    %eax,%eax
  nopl   0x0(%rax,%rax,1)
L130:
  movsd  0x0(%r13,%rax,8),%xmm0
  movsd  %xmm0,(%rdi,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%r14
  jne    L130
  jmp    L215
L41:
  mov    $0x8,%esi
  mov    $0x3d0900,%edi
  call   <polybench_alloc_data@plt>
  mov    %rax,%rbp
L174:
  add    $0x58,%rsp
  mov    %rbp,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1070 <free@plt>
```

tool_calls: name: unpack_data
arguments: data_label: D19
data_type: byte[8]

tool: {"D19": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

assistant: ```c
void init_array(int n, double A[2000 + 0][2000 + 0]) {
  int i, j;
  for (i = 0; i < n; i++) {
    for (j = 0; j <= i; j++)
      A[i][j] = (double)(-j % n) / n + 1;
    for (j = i + 1; j < n; j++) {
      A[i][j] = 0;
    }
    A[i][i] = 1;
  }
  int r, s, t;
  double(*B)[2000 + 0][2000 + 0];
  B = (double(*)[2000 + 0][2000 + 0]) polybench_alloc_data((2000 + 0) * (2000 + 0), sizeof(double));
  ;
  for (r = 0; r < n; ++r)
    for (s = 0; s < n; ++s)
      (*B)[r][s] = 0;
  for (t = 0; t < n; ++t)
    for (r = 0; r < n; ++r)
      for (s = 0; s < n; ++s)
        (*B)[r][s] += A[r][t] * A[s][t];
  for (r = 0; r < n; ++r)
    for (s = 0; s < n; ++s)
      A[r][s] = (*B)[r][s];
  free((void *)B);
  ;
}
```