user: 下面的汇编指令对应的C源代码是什么:

```asm
<get_Delaunay_reduction>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xb0,%rsp
  mov    %rdi,-0x98(%rbp)
  mov    %rsi,-0xa0(%rbp)
  movsd  %xmm0,-0xa8(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    -0xa0(%rbp),%rdx
  lea    -0x70(%rbp),%rax
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <get_exteneded_basis@plt>
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movl   $0x0,-0x88(%rbp)
  jmp    L138
L50:
  movl   $0x0,-0x84(%rbp)
  jmp    L169
L124:
  mov    -0x84(%rbp),%eax
  movslq %eax,%rcx
  mov    -0x88(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  add    %rcx,%rax
  movsd  -0x70(%rbp,%rax,8),%xmm1
  mov    -0x84(%rbp),%eax
  movslq %eax,%rcx
  mov    -0x88(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  add    %rcx,%rax
  movsd  -0x70(%rbp,%rax,8),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x80(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  addl   $0x1,-0x84(%rbp)
L169:
  cmpl   $0x2,-0x84(%rbp)
  jle    L124
  addl   $0x1,-0x88(%rbp)
L138:
  cmpl   $0x3,-0x88(%rbp)
  jle    L50
L151:
  mov    -0xa8(%rbp),%rdx
  lea    -0x70(%rbp),%rax
  movq   %rdx,%xmm0
  mov    %rax,%rdi
  call   <get_Delaunay_reduction_basis@plt>
  test   %rax,%rax
  jne    L103
  jmp    L151
L103:
  nop
  mov    -0xa8(%rbp),%rdx
  lea    -0x70(%rbp),%rax
  movq   %rdx,%xmm0
  mov    %rax,%rdi
  call   <get_Delaunay_shortest_vectors@plt>
  movl   $0x0,-0x88(%rbp)
  jmp    L133
L146:
  movl   $0x0,-0x84(%rbp)
  jmp    L156
L15:
  mov    -0x88(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x98(%rbp),%rax
  lea    (%rdx,%rax,1),%rcx
  mov    -0x88(%rbp),%eax
  movslq %eax,%rsi
  mov    -0x84(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  add    %rsi,%rax
  movsd  -0x70(%rbp,%rax,8),%xmm0
  mov    -0x84(%rbp),%eax
  cltq
  movsd  %xmm0,(%rcx,%rax,8)
  addl   $0x1,-0x84(%rbp)
L156:
  cmpl   $0x2,-0x84(%rbp)
  jle    L15
  addl   $0x1,-0x88(%rbp)
L133:
  cmpl   $0x2,-0x88(%rbp)
  jle    L146
  mov    -0x98(%rbp),%rax
  mov    %rax,%rdi
  call   <mat_get_determinant_d3@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x78(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <mat_Dabs@plt>
  movsd  -0xa8(%rbp),%xmm1
  comisd %xmm0,%xmm1
  ja     L148
  pxor   %xmm0,%xmm0
  comisd -0x78(%rbp),%xmm0
  jbe    L31
  movl   $0x0,-0x88(%rbp)
  jmp    L82
L70:
  movl   $0x0,-0x84(%rbp)
  jmp    L26
L163:
  mov    -0x88(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x98(%rbp),%rax
  add    %rax,%rdx
  mov    -0x84(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm0
  mov    -0x88(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x98(%rbp),%rax
  add    %rax,%rdx
  movq   D9(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mov    -0x84(%rbp),%eax
  cltq
  movsd  %xmm0,(%rdx,%rax,8)
  addl   $0x1,-0x84(%rbp)
L26:
  cmpl   $0x2,-0x84(%rbp)
  jle    L163
  addl   $0x1,-0x88(%rbp)
L82:
  cmpl   $0x2,-0x88(%rbp)
  jle    L70
L31:
  mov    $0x1,%eax
  jmp    L81
L148:
  nop
  mov    $0x0,%eax
L81:
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L40
  call   <__stack_chk_fail@plt>
L40:
  leave
  ret
```

tool_calls: name: parse_data
arguments: tag: D9
type: byte[8]

tool: {"D9": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

assistant: ```c
int get_Delaunay_reduction(double red_lattice[3][3], double lattice[3][3], const double symprec) {
  int i, j;
  double volume, sum;
  double basis[4][3];
  get_exteneded_basis(basis, lattice);
  sum = 0;
  for (i = 0; i < 4; i++) {
    for (j = 0; j < 3; j++) {
      sum += basis[i][j] * basis[i][j];
    }
  }
  while (1) {
    if (get_Delaunay_reduction_basis(basis, symprec)) {
      break;
    }
  }
  get_Delaunay_shortest_vectors(basis, symprec);
  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      red_lattice[i][j] = basis[j][i];
    }
  }
  volume = mat_get_determinant_d3(red_lattice);
  if (mat_Dabs(volume) < symprec) {
    ;
    goto err;
  }
  if (volume < 0) {
    for (i = 0; i < 3; i++) {
      for (j = 0; j < 3; j++) {
        red_lattice[i][j] = -red_lattice[i][j];
      }
    }
  }
  return 1;
err:
  return 0;
}
```