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

```asm
<get_delaunay_shortest_vectors_2D>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %rdi,%r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x118,%rsp
  movdqu (%rdi),%xmm2
  movdqu 0x10(%rdi),%xmm3
  movsd  %xmm0,0x20(%rsp)
  movsd  (%rdi),%xmm0
  addsd  0x18(%rdi),%xmm0
  lea    0x30(%rsp),%r15
  movdqu 0x20(%rdi),%xmm4
  movdqu 0x30(%rdi),%xmm5
  mov    %rsi,0x28(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x108(%rsp)
  xor    %eax,%eax
  mov    0x40(%rdi),%rax
  movsd  %xmm0,0xe8(%rsp)
  movsd  0x8(%rdi),%xmm0
  addsd  0x20(%rdi),%xmm0
  movl   $0x3,0x14(%rsp)
  mov    %rax,0xe0(%rsp)
  lea    0xa0(%rsp),%rax
  mov    %rax,0x18(%rsp)
  lea    0xe8(%rsp),%rax
  movsd  %xmm0,0xf0(%rsp)
  movsd  0x10(%rdi),%xmm0
  addsd  0x28(%rdi),%xmm0
  movaps %xmm2,0xa0(%rsp)
  movaps %xmm3,0xb0(%rsp)
  movaps %xmm4,0xc0(%rsp)
  movaps %xmm5,0xd0(%rsp)
  movsd  %xmm0,0xf8(%rsp)
  mov    %rax,0x8(%rsp)
L66:
  mov    0x18(%rsp),%r14
  mov    %r14,%rbp
L19:
  mov    %r14,%rdi
  mov    %r14,%r12
  add    $0x18,%r14
  call   <mat_norm_squared_d3@plt>
  mov    %r14,%rdi
  mov    %rax,%rbx
  call   <mat_norm_squared_d3@plt>
  cmp    %rax,%rbx
  jg     L79
L28:
  cmp    %r14,0x8(%rsp)
  jne    L19
  subl   $0x1,0x14(%rsp)
  jne    L66
  mov    0x28(%rsp),%rax
  mov    $0x1,%ebx
  lea    0x50(%rsp),%r12
  movsd  0xa0(%rsp),%xmm0
  movsd  %xmm0,0x50(%rsp)
  movsd  (%rax),%xmm0
  movsd  %xmm0,0x58(%rsp)
  movsd  0xa8(%rsp),%xmm0
  movsd  %xmm0,0x68(%rsp)
  movsd  0x8(%rax),%xmm0
  movsd  %xmm0,0x70(%rsp)
  movsd  0xb0(%rsp),%xmm0
  movsd  %xmm0,0x80(%rsp)
  movsd  0x10(%rax),%xmm0
  movsd  %xmm0,0x88(%rsp)
L34:
  movsd  0x18(%rbp),%xmm0
  mov    %r12,%rdi
  movsd  %xmm0,0x60(%rsp)
  movsd  0x20(%rbp),%xmm0
  movsd  %xmm0,0x78(%rsp)
  movsd  0x28(%rbp),%xmm0
  movsd  %xmm0,0x90(%rsp)
  call   <mat_get_determinant_d3@plt>
  mov    %eax,%edi
  call   <mat_Dabs@plt>
  comisd 0x20(%rsp),%xmm0
  ja     L72
  add    $0x1,%ebx
  add    $0x18,%rbp
  cmp    $0x4,%ebx
  jne    L34
L4:
  mov    0x108(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L42
  add    $0x118,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L79:
  mov    %r12,%rsi
  mov    %r15,%rdi
  call   <mat_copy_vector_d3@plt>
  mov    %r14,%rsi
  mov    %r12,%rdi
  call   <mat_copy_vector_d3@plt>
  mov    %r15,%rsi
  mov    %r14,%rdi
  call   <mat_copy_vector_d3@plt>
  jmp    L28
  cs nopw 0x0(%rax,%rax,1)
L72:
  movsd  0xa0(%rsp),%xmm0
  movslq %ebx,%rbx
  lea    (%rbx,%rbx,2),%rax
  shl    $0x3,%rax
  movsd  %xmm0,0x0(%r13)
  movsd  0xa0(%rsp,%rax,1),%xmm0
  movsd  %xmm0,0x18(%r13)
  movsd  0xa8(%rsp),%xmm0
  movsd  %xmm0,0x8(%r13)
  movsd  0xa8(%rsp,%rax,1),%xmm0
  movsd  %xmm0,0x20(%r13)
  movsd  0xb0(%rsp),%xmm0
  movsd  %xmm0,0x10(%r13)
  movsd  0xb0(%rsp,%rax,1),%xmm0
  movsd  %xmm0,0x28(%r13)
  jmp    L4
L42:
  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;
    }
  }
}
```