user: Translate the below assembly code into C programming language:

```asm
<ludcmp2>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %rdi,0x20(%rsp)
  mov    %esi,%ebp
  mov    %rdx,0x30(%rsp)
  mov    %rcx,%rbx
  mov    %rcx,0x58(%rsp)
  mov    %esi,%edi
  call   <dvector@plt>
  mov    %rax,%r15
  mov    D34(%rip),%rax
  mov    %rax,(%rbx)
  test   %ebp,%ebp
  jle    L50
  mov    %r15,%r13
  movslq %ebp,%rbx
  shl    $0x3,%rbx
  lea    -0x1(%rbp),%eax
  mov    %rax,0x50(%rsp)
  shl    $0x3,%rax
  mov    0x20(%rsp),%rdi
  lea    0x8(%rdi),%rsi
  mov    %rsi,0x28(%rsp)
  lea    (%rsi,%rax,1),%r12
  lea    0x8(%r15,%rax,1),%r14
  mov    %ebp,%eax
  neg    %rax
  shl    $0x3,%rax
  mov    %ebp,0x14(%rsp)
  mov    %rax,%rbp
  jmp    L264
L44:
  movsd  D34(%rip),%xmm0
  divsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x0(%r13)
  add    $0x8,%r13
  add    %rbx,%r12
  cmp    %r14,%r13
  je     L23
L264:
  lea    0x0(%rbp,%r12,1),%rax
  pxor   %xmm4,%xmm4
  movsd  %xmm4,0x8(%rsp)
L253:
  movsd  (%rax),%xmm0
  andpd  D78(%rip),%xmm0
  maxsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x8(%rsp)
  add    $0x8,%rax
  cmp    %r12,%rax
  jne    L253
  pxor   %xmm4,%xmm4
  ucomisd %xmm4,%xmm0
  jp     L44
  jne    L44
  lea    D60(%rip),%rdi
  call   <error@plt>
  jmp    L44
L23:
  mov    0x14(%rsp),%ebp
  lea    0x8(%rbx),%rax
  mov    %rax,0x8(%rsp)
  mov    0x20(%rsp),%r14
  lea    (%r14,%rbx,1),%rdi
  mov    %rdi,0x40(%rsp)
  mov    0x28(%rsp),%rdi
  add    %rbx,%rdi
  mov    %rdi,0x48(%rsp)
  mov    %r14,%r9
  movl   $0xffffffff,0x38(%rsp)
  movl   $0x0,0x1c(%rsp)
  movl   $0x1,0x18(%rsp)
  mov    $0x0,%r8d
  movl   $0x0,0x14(%rsp)
  pxor   %xmm6,%xmm6
  movq   D78(%rip),%xmm3
  movq   D37(%rip),%xmm9
  movapd %xmm6,%xmm5
  movapd %xmm6,%xmm4
  movsd  D47(%rip),%xmm8
  lea    -0x1(%rbp),%eax
  mov    %eax,0x3c(%rsp)
  movsd  D34(%rip),%xmm7
  jmp    L209
L231:
  movsd  (%rdi,%r8,8),%xmm1
  test   %esi,%esi
  jle    L205
  mov    %r9,%rdx
  mov    %rdi,%rax
L123:
  movsd  (%rax),%xmm0
  mulsd  (%rdx),%xmm0
  subsd  %xmm0,%xmm1
  add    $0x8,%rax
  add    %rbx,%rdx
  cmp    %rcx,%rax
  jne    L123
L205:
  add    %rbx,%r10
  add    %rbx,%rdi
  add    0x8(%rsp),%rcx
L138:
  movsd  %xmm1,(%r10)
  add    $0x1,%esi
  cmp    %r8d,%esi
  jne    L231
L79:
  mov    %r8d,%r12d
  cmp    %r8d,%ebp
  jle    L19
  movslq 0x1c(%rsp),%rdx
  lea    (%rdx,%r8,1),%rax
  mov    0x20(%rsp),%r11
  lea    (%r11,%rax,8),%rdi
  mov    0x38(%rsp),%eax
  add    %rdx,%rax
  mov    0x28(%rsp),%rsi
  lea    (%rsi,%rax,8),%rcx
  mov    %r8,%rsi
  movapd %xmm6,%xmm2
  sub    %r9,%r11
  mov    %r11,%r13
  jmp    L143
L237:
  add    $0x1,%rsi
  add    %rbx,%rdi
  add    %rbx,%rcx
  cmp    %esi,%ebp
  jle    L19
L143:
  mov    %esi,%r11d
  mov    %rdi,%r10
  movsd  (%rdi),%xmm1
  test   %r12d,%r12d
  jle    L36
  lea    0x0(%r13,%rdi,1),%rax
  mov    %r9,%rdx
L196:
  movsd  (%rax),%xmm0
  mulsd  (%rdx),%xmm0
  subsd  %xmm0,%xmm1
  add    $0x8,%rax
  add    %rbx,%rdx
  cmp    %rcx,%rax
  jne    L196
L36:
  movsd  %xmm1,(%r10)
  andpd  %xmm3,%xmm1
  mulsd  (%r15,%rsi,8),%xmm1
  comisd %xmm2,%xmm1
  jb     L237
  movapd %xmm1,%xmm2
  mov    %r11d,0x14(%rsp)
  jmp    L237
L19:
  mov    0x14(%rsp),%eax
  cmp    %r12d,%eax
  jne    L230
L16:
  mov    0x30(%rsp),%rax
  mov    0x14(%rsp),%edi
  mov    %edi,(%rax,%r8,4)
  mov    %r14,%rax
  ucomisd (%r14),%xmm5
  jp     L149
  ucomisd (%r14),%xmm4
  jne    L149
  movsd  %xmm8,(%r14)
L149:
  cmp    %r12d,0x3c(%rsp)
  je     L67
  movapd %xmm7,%xmm1
  divsd  (%rax),%xmm1
  mov    0x18(%rsp),%eax
  mov    %eax,%edx
  cmp    %eax,%ebp
  jle    L50
  mov    0x1c(%rsp),%eax
  add    %ebp,%eax
  cltq
  add    %r8,%rax
  mov    0x20(%rsp),%rdi
  lea    (%rdi,%rax,8),%rax
L11:
  movapd %xmm1,%xmm0
  mulsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  add    $0x1,%edx
  add    %rbx,%rax
  cmp    %edx,%ebp
  jne    L11
L88:
  add    $0x1,%r8
  addl   $0x1,0x18(%rsp)
  add    0x8(%rsp),%r14
  add    %ebp,0x1c(%rsp)
  add    $0x8,%r9
  addl   $0x1,0x38(%rsp)
L209:
  mov    0x18(%rsp),%eax
  sub    $0x1,%eax
  test   %eax,%eax
  jle    L79
  movsd  (%r9),%xmm1
  mov    0x48(%rsp),%rcx
  mov    0x40(%rsp),%rdi
  mov    %r9,%r10
  mov    $0x0,%esi
  jmp    L138
L230:
  mov    %eax,%esi
  imul   %ebp,%esi
  movslq %esi,%rsi
  mov    0x20(%rsp),%rax
  lea    (%rax,%rsi,8),%rax
  mov    0x50(%rsp),%rdi
  lea    (%rsi,%rdi,1),%rdx
  mov    0x28(%rsp),%rdi
  lea    (%rdi,%rdx,8),%rcx
  movslq 0x1c(%rsp),%rdx
  sub    %rsi,%rdx
L30:
  movsd  (%rax),%xmm0
  movsd  (%rax,%rdx,8),%xmm1
  movsd  %xmm1,(%rax)
  movsd  %xmm0,(%rax,%rdx,8)
  add    $0x8,%rax
  cmp    %rcx,%rax
  jne    L30
  mov    0x58(%rsp),%rax
  movsd  (%rax),%xmm0
  xorpd  %xmm9,%xmm0
  movsd  %xmm0,(%rax)
  movsd  (%r15,%r8,8),%xmm0
  movslq 0x14(%rsp),%rax
  movsd  %xmm0,(%r15,%rax,8)
  jmp    L16
L67:
  cmp    0x18(%rsp),%ebp
  jg     L88
L50:
  mov    %r15,%rdi
  call   <free@plt>
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: data_label: D34
type: f64
name: read_data
arguments: data_label: D78
type: qword[2]
name: read_data
arguments: data_label: D60
type: string
name: read_data
arguments: data_label: D37
type: f32[2]
name: read_data
arguments: data_label: D47
type: byte[8]

tool: D34: 1

tool: {"D78": [0x7fffffffffffffff,0x0000000000000000]}

tool: D60: Singular matrix in routine ludcmp

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

tool: {"D47": [0x23,0x42,0x92,0x0c,0xa1,0x9c,0xc7,0x3b]}

assistant: ```c
void ludcmp2(double *a, int n, int *indx, double *d) {
  int i, imax, j, k;
  double big, dum, sum, temp;
  double *vv;
  imax = 0;
  vv = dvector(n);
  *d = 1.0;
  for (i = 0; i < n; i++) {
    big = 0.0;
    for (j = 0; j < n; j++)
      if ((temp = fabs(a[i * n + j])) > big)
        big = temp;
    if (big == 0.0)
      error("Singular matrix in routine ludcmp");
    vv[i] = 1.0 / big;
  }
  for (j = 0; j < n; j++) {
    for (i = 0; i < j; i++) {
      sum = a[i * n + j];
      for (k = 0; k < i; k++)
        sum -= a[i * n + k] * a[k * n + j];
      a[i * n + j] = sum;
    }
    big = 0.0;
    for (i = j; i < n; i++) {
      sum = a[i * n + j];
      for (k = 0; k < j; k++)
        sum -= a[i * n + k] * a[k * n + j];
      a[i * n + j] = sum;
      if ((dum = vv[i] * fabs(sum)) >= big) {
        big = dum;
        imax = i;
      }
    }
    if (j != imax) {
      for (k = 0; k < n; k++) {
        dum = a[imax * n + k];
        a[imax * n + k] = a[j * n + k];
        a[j * n + k] = dum;
      }
      *d = -(*d);
      vv[imax] = vv[j];
    }
    indx[j] = imax;
    if (a[j * n + j] == 0.0)
      a[j * n + j] = 1.0e-20;
    if (j != n - 1) {
      dum = 1.0 / (a[j * n + j]);
      for (i = j + 1; i < n; i++)
        a[i * n + j] *= dum;
    }
  }
  free(vv);
}
```