user: 请将以下的汇编代码转换成C语言代码：

```asm
<get_delaunay_shortest_vectors_2D>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x108,%rsp
  mov    %rdi,%r15
  mov    %rsi,0x18(%rsp)
  movsd  %xmm0,0x10(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0xf8(%rsp)
  xor    %eax,%eax
  mov    %rdi,%rdx
  lea    0x90(%rsp),%r12
  lea    0xd8(%rsp),%r14
  mov    %r12,%rax
L12:
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  movsd  0x8(%rdx),%xmm0
  movsd  %xmm0,0x8(%rax)
  movsd  0x10(%rdx),%xmm0
  movsd  %xmm0,0x10(%rax)
  add    $0x18,%rdx
  add    $0x18,%rax
  cmp    %r14,%rax
  jne    L12
  movsd  (%r15),%xmm0
  addsd  0x18(%r15),%xmm0
  movsd  %xmm0,0xd8(%rsp)
  movsd  0x8(%r15),%xmm0
  addsd  0x20(%r15),%xmm0
  movsd  %xmm0,0xe0(%rsp)
  movsd  0x10(%r15),%xmm0
  addsd  0x28(%r15),%xmm0
  movsd  %xmm0,0xe8(%rsp)
  movl   $0x3,0x4(%rsp)
  lea    0x20(%rsp),%rax
  mov    %rax,0x8(%rsp)
  jmp    L83
L86:
  cmp    %r14,%rbx
  je     L92
L33:
  mov    %rbx,%r13
  mov    %rbx,%rdi
  call   <mat_norm_squared_d3@plt>
  mov    %rax,%rbp
  add    $0x18,%rbx
  mov    %rbx,%rdi
  call   <mat_norm_squared_d3@plt>
  cmp    %rax,%rbp
  jle    L86
  mov    %r13,%rsi
  mov    0x8(%rsp),%rbp
  mov    %rbp,%rdi
  call   <mat_copy_vector_d3@plt>
  mov    %rbx,%rsi
  mov    %r13,%rdi
  call   <mat_copy_vector_d3@plt>
  mov    %rbp,%rsi
  mov    %rbx,%rdi
  call   <mat_copy_vector_d3@plt>
  jmp    L86
L92:
  subl   $0x1,0x4(%rsp)
  je     L98
L83:
  mov    %r12,%rbx
  jmp    L33
L98:
  movsd  0x90(%rsp),%xmm0
  movsd  %xmm0,0x40(%rsp)
  mov    0x18(%rsp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,0x48(%rsp)
  movsd  0x98(%rsp),%xmm0
  movsd  %xmm0,0x58(%rsp)
  movsd  0x8(%rax),%xmm0
  movsd  %xmm0,0x60(%rsp)
  movsd  0xa0(%rsp),%xmm0
  movsd  %xmm0,0x70(%rsp)
  movsd  0x10(%rax),%xmm0
  movsd  %xmm0,0x78(%rsp)
  mov    $0x1,%ebx
  lea    0x40(%rsp),%rbp
L16:
  movsd  0x18(%r12),%xmm0
  movsd  %xmm0,0x50(%rsp)
  movsd  0x20(%r12),%xmm0
  movsd  %xmm0,0x68(%rsp)
  movsd  0x28(%r12),%xmm0
  movsd  %xmm0,0x80(%rsp)
  mov    %rbp,%rdi
  call   <mat_get_determinant_d3@plt>
  mov    %eax,%edi
  call   <mat_Dabs@plt>
  comisd 0x10(%rsp),%xmm0
  ja     L101
  add    $0x1,%ebx
  add    $0x18,%r12
  cmp    $0x4,%ebx
  jne    L16
L40:
  mov    0xf8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L21
  add    $0x108,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L101:
  movsd  0x90(%rsp),%xmm0
  movsd  %xmm0,(%r15)
  movslq %ebx,%rbx
  lea    (%rbx,%rbx,1),%rax
  lea    (%rax,%rbx,1),%rdx
  movsd  0x90(%rsp,%rdx,8),%xmm0
  movsd  %xmm0,0x18(%r15)
  movsd  0x98(%rsp),%xmm0
  movsd  %xmm0,0x8(%r15)
  movsd  0x98(%rsp,%rdx,8),%xmm0
  movsd  %xmm0,0x20(%r15)
  movsd  0xa0(%rsp),%xmm0
  movsd  %xmm0,0x10(%r15)
  movsd  0xa0(%rsp,%rdx,8),%xmm0
  movsd  %xmm0,0x28(%r15)
  jmp    L40
L21:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void get_delaunay_shortest_vectors_2D(double basis[3][3], const double unique_vec[3], const double symprec) {
  int i, j;
  double b[4][3], tmpmat[3][3];
  double tmpvec[3];
  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      b[i][j] = basis[i][j];
    }
  }
  for (i = 0; i < 3; i++) {
    b[3][i] = basis[0][i] + basis[1][i];
  }
  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      if (mat_norm_squared_d3(b[j]) > mat_norm_squared_d3(b[j + 1])) {
        mat_copy_vector_d3(tmpvec, b[j]);
        mat_copy_vector_d3(b[j], b[j + 1]);
        mat_copy_vector_d3(b[j + 1], tmpvec);
      }
    }
  }
  for (i = 0; i < 3; i++) {
    tmpmat[i][0] = b[0][i];
    tmpmat[i][1] = unique_vec[i];
  }
  for (i = 1; i < 4; i++) {
    for (j = 0; j < 3; j++) {
      tmpmat[j][2] = b[i][j];
    }
    if (mat_Dabs(mat_get_determinant_d3(tmpmat)) > symprec) {
      for (j = 0; j < 3; j++) {
        basis[0][j] = b[0][j];
        basis[1][j] = b[i][j];
      }
      break;
    }
  }
}
```