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

```asm
<project_v>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %edi,-0x54(%rbp)
  mov    %esi,-0x58(%rbp)
  mov    %edx,-0x5c(%rbp)
  mov    %ecx,-0x60(%rbp)
  mov    %r8,-0x68(%rbp)
  mov    %r9,-0x70(%rbp)
  mov    -0x54(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x54(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x54(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x54(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x5c(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x60(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x54(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x5c(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x60(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x54(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movl   $0x0,-0x48(%rbp)
  jmp    L19
L23:
  movl   $0x0,-0x44(%rbp)
  jmp    L21
L38:
  mov    -0x54(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x48(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x44(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x58(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x48(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x44(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x5c(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x60(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x48(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x44(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x5c(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x60(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x48(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x44(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x40(%rbp),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movsd  -0x38(%rbp),%xmm1
  mulsd  -0x8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x30(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x20(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  mov    -0x5c(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x60(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x48(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x44(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x8(%rbp),%xmm1
  movsd  -0x38(%rbp),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x30(%rbp),%xmm1
  mulsd  -0x18(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x28(%rbp),%xmm1
  mulsd  -0x20(%rbp),%xmm1
  mov    -0x5c(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x60(%rbp),%eax
  add    %eax,%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x48(%rbp),%eax
  add    %edx,%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x44(%rbp),%eax
  add    %edx,%eax
  add    %eax,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x44(%rbp)
L21:
  cmpl   $0x1,-0x44(%rbp)
  jle    L38
  addl   $0x1,-0x48(%rbp)
L19:
  cmpl   $0x2,-0x48(%rbp)
  jle    L23
  nop
  nop
  pop    %rbp
  ret
```

assistant: ```c
void project_v(int a, int i, int b, int j, double r[], double v[]) {
  int c, k;
  double xi, yi, xj, yj;
  xi = (&((r)[2 * ((i) + 2 * ((a) + 3 * ((i) + 2 * (a))))]))[0];
  yi = (&((r)[2 * ((i) + 2 * ((a) + 3 * ((i) + 2 * (a))))]))[1];
  xj = (&((r)[2 * ((i) + 2 * ((a) + 3 * ((j) + 2 * (b))))]))[0];
  yj = (&((r)[2 * ((i) + 2 * ((a) + 3 * ((j) + 2 * (b))))]))[1];
  for (c = 0; c < 3; c++)
    for (k = 0; k < 2; k++) {
      double zi = (&((v)[2 * ((k) + 2 * ((c) + 3 * ((i) + 2 * (a))))]))[0];
      double ti = (&((v)[2 * ((k) + 2 * ((c) + 3 * ((i) + 2 * (a))))]))[1];
      double zj = (&((v)[2 * ((k) + 2 * ((c) + 3 * ((j) + 2 * (b))))]))[0];
      double tj = (&((v)[2 * ((k) + 2 * ((c) + 3 * ((j) + 2 * (b))))]))[1];
      (&((v)[2 * ((k) + 2 * ((c) + 3 * ((j) + 2 * (b))))]))[0] = xi * zj - yi * tj - xj * zi + yj * ti;
      (&((v)[2 * ((k) + 2 * ((c) + 3 * ((j) + 2 * (b))))]))[1] = xi * tj + yi * zj - xj * ti - yj * zi;
    }
}
```