user: What is the c source code of the assembly code below:

```asm
<LUDecompose_det>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  movslq %esi,%rbx
  mov    %ebx,%edi
  sub    $0x18,%rsp
  call   <dvector@plt>
  mov    %rax,%r12
  test   %rbx,%rbx
  je     L28
  movq   D22(%rip),%xmm3
  movsd  D29(%rip),%xmm5
  xor    %r15d,%r15d
  pxor   %xmm4,%xmm4
  lea    0x0(,%rbx,8),%r13
  lea    D7(%rip),%r14
  nopl   0x0(%rax)
L41:
  mov    0x0(%rbp,%r15,8),%rax
  movapd %xmm4,%xmm1
  lea    0x0(%r13,%rax,1),%rdx
  cs nopw 0x0(%rax,%rax,1)
L68:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  andpd  %xmm3,%xmm0
  maxsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  cmp    %rax,%rdx
  jne    L68
  ucomisd %xmm4,%xmm0
  jp     L95
  je     L18
L95:
  movapd %xmm5,%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%r12,%r15,8)
  add    $0x1,%r15
  cmp    %rbx,%r15
  jne    L41
  xor    %r10d,%r10d
  lea    0x0(%r13,%rbp,1),%r8
  lea    0x8(%rbp),%r9
  xor    %esi,%esi
  lea    0x1(%r10),%rdi
  movsd  D73(%rip),%xmm6
  mov    0x0(%rbp,%rsi,1),%rcx
  movq   D49(%rip),%xmm7
  cmp    %rbx,%rdi
  jae    L13
  nopl   0x0(%rax)
L8:
  mov    %r10,%rdx
  mov    %rdi,%rax
  movapd %xmm4,%xmm1
  nopw   0x0(%rax,%rax,1)
L183:
  mov    0x0(%rbp,%rax,8),%r11
  movsd  (%r11,%rsi,1),%xmm0
  andpd  %xmm3,%xmm0
  mulsd  (%r12,%rax,8),%xmm0
  comisd %xmm1,%xmm0
  maxsd  %xmm1,%xmm0
  cmova  %rax,%rdx
  add    $0x1,%rax
  movapd %xmm0,%xmm1
  cmp    %rbx,%rax
  jne    L183
  cmp    %r10,%rdx
  jne    L111
L72:
  lea    (%rcx,%rsi,1),%r11
  ucomisd (%r11),%xmm4
  jp     L175
  jne    L175
L179:
  movsd  %xmm6,(%r11)
L205:
  cmp    %rbx,%rdi
  je     L0
L175:
  lea    (%r9,%rsi,1),%r10
  nopw   0x0(%rax,%rax,1)
L188:
  mov    (%r10),%rdx
  lea    (%rdx,%rsi,1),%rax
  movsd  (%rax),%xmm2
  divsd  (%r11),%xmm2
  movsd  %xmm2,(%rax)
  mov    %rdi,%rax
  nop
L40:
  movsd  (%rcx,%rax,8),%xmm1
  movsd  (%rdx,%rax,8),%xmm0
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  add    $0x1,%rax
  cmp    %rbx,%rax
  jne    L40
  add    $0x8,%r10
  cmp    %r10,%r8
  jne    L188
  mov    %rdi,%r10
  add    $0x8,%rsi
  lea    0x1(%r10),%rdi
  mov    0x0(%rbp,%rsi,1),%rcx
  cmp    %rbx,%rdi
  jb     L8
L13:
  lea    (%rcx,%rsi,1),%r11
  ucomisd (%r11),%xmm4
  jp     L205
  je     L179
  cmp    %rbx,%rdi
  jne    L175
L0:
  xor    %eax,%eax
  cs nopw 0x0(%rax,%rax,1)
L106:
  mov    0x0(%rbp,%rax,1),%rdx
  mulsd  (%rdx,%rax,1),%xmm5
  add    $0x8,%rax
  cmp    %rax,%r13
  jne    L106
L196:
  mov    %r12,%rdi
  movsd  %xmm5,0x8(%rsp)
  call   <free@plt>
  movsd  0x8(%rsp),%xmm5
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  movapd %xmm5,%xmm0
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L111:
  lea    0x0(,%rdx,8),%r10
  mov    0x0(%rbp,%rdx,8),%rdx
  xor    %eax,%eax
  nopl   0x0(%rax)
L177:
  movsd  (%rdx,%rax,8),%xmm0
  movsd  (%rcx,%rax,8),%xmm1
  movsd  %xmm1,(%rdx,%rax,8)
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %rbx,%rax
  jne    L177
  movsd  (%r12,%rsi,1),%xmm0
  xorpd  %xmm7,%xmm5
  movsd  %xmm0,(%r12,%r10,1)
  jmp    L72
L18:
  mov    %r14,%rdi
  movsd  %xmm0,0x8(%rsp)
  call   <error@plt>
  mov    D29(%rip),%rax
  movsd  0x8(%rsp),%xmm1
  pxor   %xmm4,%xmm4
  movq   D22(%rip),%xmm3
  movq   %rax,%xmm5
  jmp    L95
L28:
  movsd  D29(%rip),%xmm5
  jmp    L196
```

tool_calls: name: parse_data
arguments: data_label: D22
type: qword
name: parse_data
arguments: data_label: D29
type: f64
name: parse_data
arguments: data_label: D7
type: string
name: parse_data
arguments: data_label: D73
type: byte[8]
name: parse_data
arguments: data_label: D49
type: f32[2]

tool: D22: 9223372036854775807

tool: D29: 1

tool: D7: LU decomposition does not like singular matrices

tool: {"D73": [0xd9,0x7d,0xda,0xf5,0xd0,0xf2,0xbe,0x3a]}

tool: {"D49": [0.0,-0.0]}

assistant: ```c
double LUDecompose_det(double **m, const int dim) {
  double big = 0.0;
  double temp = 0.0;
  double *vv = dvector(dim);
  double d = 1.0;
  for (size_t i = 0; i < dim; i++) {
    big = 0.0;
    for (size_t j = 0; j < dim; j++) {
      if ((temp = fabs(m[i][j])) > big)
        big = temp;
    }
    if (big == 0.0) {
      error("LU decomposition does not like singular matrices");
    }
    vv[i] = 1.0 / big;
  }
  size_t i, j;
  size_t imax;
  for (size_t k = 0; k < dim; k++) {
    big = 0.0;
    imax = k;
    for (i = k + 1; i < dim; i++) {
      temp = vv[i] * fabs(m[i][k]);
      if (temp > big) {
        big = temp;
        imax = i;
      }
    }
    if (k != imax) {
      for (j = 0; j < dim; j++) {
        temp = m[imax][j];
        m[imax][j] = m[k][j];
        m[k][j] = temp;
      }
      d = -d;
      vv[imax] = vv[k];
    }
    if (m[k][k] == 0.0)
      m[k][k] = 1.0e-25;
    for (i = k + 1; i < dim; i++) {
      temp = m[i][k] /= m[k][k];
      for (j = k + 1; j < dim; j++) {
        m[i][j] -= temp * m[k][j];
      }
    }
  }
  for (size_t i = 0; i < dim; i++) {
    d *= m[i][i];
  }
  free(vv);
  return d;
}
```