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

```asm
<main>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x6f0,%rsp
  mov    %edi,-0x6e4(%rbp)
  mov    %rsi,-0x6f0(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  cmpl   $0x1,-0x6e4(%rbp)
  jle    L18
  mov    -0x6f0(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%rax
  mov    %rax,%rdi
  call   <atoi@plt>
  mov    %eax,%edi
  call   <srandom@plt>
L18:
  movl   $0x0,-0x6dc(%rbp)
  jmp    L113
L105:
  mov    $0x0,%eax
  call   <random@plt>
  movq   %xmm0,%rax
  movsd  D1(%rip),%xmm1
  movq   %rax,%xmm0
  divsd  %xmm1,%xmm0
  mov    -0x6dc(%rbp),%eax
  cltq
  movsd  %xmm0,-0x6d0(%rbp,%rax,8)
  mov    -0x6dc(%rbp),%eax
  cltq
  movsd  -0x6d0(%rbp,%rax,8),%xmm0
  mov    -0x6dc(%rbp),%eax
  cltq
  movsd  %xmm0,-0x490(%rbp,%rax,8)
  addl   $0x1,-0x6dc(%rbp)
L113:
  cmpl   $0x47,-0x6dc(%rbp)
  jle    L105
  lea    -0x6d0(%rbp),%rdx
  lea    -0x250(%rbp),%rax
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <cl_invert@plt>
  lea    -0x490(%rbp),%rdx
  lea    -0x250(%rbp),%rcx
  lea    -0x6d0(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <mul@plt>
  movl   $0x0,-0x6e0(%rbp)
  jmp    L14
L53:
  movl   $0x0,-0x6dc(%rbp)
  jmp    L31
L50:
  movl   $0x0,-0x6d8(%rbp)
  jmp    L20
L51:
  movl   $0x0,-0x6d4(%rbp)
  jmp    L7
L66:
  mov    -0x6e0(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x6dc(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x6d8(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x6d4(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  lea    -0x6d0(%rbp),%rdx
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm2
  mov    -0x6e0(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x6dc(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x6d8(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x6d4(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  movsd  -0x6d0(%rbp,%rax,8),%xmm1
  mov    -0x6e0(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x6dc(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x6d8(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x6d4(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  lea    -0x250(%rbp),%rdx
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mov    -0x6e0(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x6dc(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x6d8(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x6d4(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  mov    -0x250(%rbp,%rax,8),%rsi
  mov    -0x6d4(%rbp),%edi
  mov    -0x6d8(%rbp),%ecx
  mov    -0x6dc(%rbp),%edx
  mov    -0x6e0(%rbp),%eax
  movapd %xmm2,%xmm3
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rsi,%xmm0
  mov    %edi,%r8d
  mov    %eax,%esi
  lea    D18(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  addl   $0x1,-0x6d4(%rbp)
L7:
  cmpl   $0x1,-0x6d4(%rbp)
  jle    L66
  addl   $0x1,-0x6d8(%rbp)
L20:
  cmpl   $0x2,-0x6d8(%rbp)
  jle    L51
  addl   $0x1,-0x6dc(%rbp)
L31:
  cmpl   $0x1,-0x6dc(%rbp)
  jle    L50
  addl   $0x1,-0x6e0(%rbp)
L14:
  cmpl   $0x2,-0x6e0(%rbp)
  jle    L53
  mov    $0x0,%eax
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L17
  call   <__stack_chk_fail@plt>
L17:
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: tag: D1
type_guessed: byte[8]
name: retrieve_data
arguments: tag: D18
type_guessed: string

tool: {"D1": [0x00,0x00,0xc0,0xff,0xff,0xff,0xdf,0x41]}

tool: D18:  [%d,%d][%d,%d]  %15f %15f   %15f %15f


assistant: ```c
int main(int argc, char *argv[]) {
  double cl[(2 * 2 * 3 * 3 * 2)];
  double v[(2 * 2 * 3 * 3 * 2)];
  double ci[(2 * 2 * 3 * 3 * 2)];
  int a, i, b, j;
  if (argc > 1)
    srandom(atoi(argv[1]));
  for (i = 0; i < (2 * 2 * 3 * 3 * 2); i++)
    v[i] = cl[i] = random() / (double)0x7fffffff;
  cl_invert(ci, cl);
  mul(cl, ci, v);
  for (a = 0; a < 3; a++)
    for (i = 0; i < 2; i++) {
      for (b = 0; b < 3; b++)
        for (j = 0; j < 2; j++) {
          printf(" [%d,%d][%d,%d]  %15f %15f   %15f %15f\n", a, i, b, j, (&((ci)[2 * ((j) + 2 * ((b) + 3 * ((i) + 2 * (a))))]))[0], (&((ci)[2 * ((j) + 2 * ((b) + 3 * ((i) + 2 * (a))))]))[1], (&((cl)[2 * ((j) + 2 * ((b) + 3 * ((i) + 2 * (a))))]))[0], (&((cl)[2 * ((j) + 2 * ((b) + 3 * ((i) + 2 * (a))))]))[1]);
        }
    }
  return 0;
}
```