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

```asm
<cholmod>:
  endbr64
  push   %r15
  push   %r14
  mov    %rdi,%r14
  push   %r13
  push   %r12
  push   %rbp
  movslq %esi,%rbp
  push   %rbx
  mov    %ebp,%r13d
  mov    %rdx,%rbx
  mov    %rbp,%r12
  imul   %ebp,%r13d
  sub    $0x98,%rsp
  mov    %rdi,0x30(%rsp)
  lea    0x0(,%r13,8),%edi
  movsd  %xmm0,0x8(%rsp)
  call   <malloc@plt>
  movsd  0x8(%rsp),%xmm5
  movsd  D11(%rip),%xmm8
  mov    %rax,0x60(%rsp)
  mov    %rax,%rdi
  mulsd  %xmm5,%xmm8
  test   %ebp,%ebp
  jle    L35
  mov    %ebp,%ecx
  lea    0x8(,%rbp,8),%rsi
  mov    %r14,%rdx
  lea    (%rax,%rcx,8),%rcx
  nopl   (%rax)
L25:
  movsd  (%rdx),%xmm0
  add    $0x8,%rax
  add    %rsi,%rdx
  movsd  %xmm0,-0x8(%rax)
  cmp    %rcx,%rax
  jne    L25
  lea    -0x1(%r13),%eax
  test   %r13d,%r13d
  mov    %rbx,%rdi
  movsd  %xmm5,0x10(%rsp)
  lea    0x8(,%rax,8),%rdx
  mov    $0x8,%eax
  movsd  %xmm8,0x8(%rsp)
  cmovle %rax,%rdx
  xor    %esi,%esi
  call   <memset@plt>
  mov    0x60(%rsp),%rdi
  mov    %r12d,%esi
  call   <array_max_abs@plt>
  movsd  0x10(%rsp),%xmm5
  pxor   %xmm3,%xmm3
  movsd  0x8(%rsp),%xmm8
  ucomisd %xmm3,%xmm5
  jp     L362
  je     L89
L362:
  movsd  D8(%rip),%xmm7
  mulsd  %xmm5,%xmm7
L197:
  xor    %r14d,%r14d
  mov    0x30(%rsp),%r13
  xor    %r15d,%r15d
  movapd %xmm3,%xmm6
  lea    0x0(,%rbp,8),%r11
  movq   D49(%rip),%xmm4
  shl    $0x4,%rbp
  movapd %xmm3,%xmm0
  lea    0x8(%r11),%rax
  mov    %rbp,0x58(%rsp)
  xor    %ebp,%ebp
  mov    %rax,0x48(%rsp)
  lea    (%rbx,%r11,1),%rax
  mov    %rax,0x40(%rsp)
  mov    %rax,0x18(%rsp)
  lea    -0x1(%r12),%eax
  mov    %rax,0x38(%rsp)
  lea    -0x2(%r12),%eax
  movq   $0x1,0x20(%rsp)
  mov    %eax,0x54(%rsp)
  mov    %r11,%rax
  mov    %r14,%r11
  mov    %rax,%r14
  nopl   0x0(%rax)
L255:
  mov    0x30(%rsp),%rax
  mov    0x20(%rsp),%r9
  mov    %r15d,0x8(%rsp)
  mov    %r15d,0x10(%rsp)
  movsd  (%rax,%r11,1),%xmm1
  mov    %r9d,0x50(%rsp)
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rbx,%r11,1)
  cmp    %r9d,%r12d
  jle    L317
  lea    0x0(%rbp,%r12,1),%esi
  mov    0x18(%rsp),%r8
  mov    %r15d,%r10d
  mov    %r11,0x28(%rsp)
  movslq %esi,%rax
  and    $0xfffffffe,%r10d
  movapd %xmm3,%xmm9
  mov    %r15d,%r11d
  lea    (%rbx,%rax,8),%rcx
  mov    %r15d,%eax
  sub    %r14,%r8
  shr    %eax
  lea    -0x1(%rax),%edi
  add    $0x1,%rdi
  shl    $0x4,%rdi
  nopw   0x0(%rax,%rax,1)
L158:
  test   %r15,%r15
  je     L211
  cmp    $0x1,%r11d
  je     L5
  xor    %eax,%eax
  movapd %xmm3,%xmm0
  nopl   0x0(%rax)
L144:
  movupd (%r8,%rax,1),%xmm2
  movupd (%rcx,%rax,1),%xmm1
  add    $0x10,%rax
  mulpd  %xmm1,%xmm2
  addsd  %xmm2,%xmm0
  unpckhpd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  cmp    %rdi,%rax
  jne    L144
  mov    %r10d,%eax
  cmp    0x8(%rsp),%r10d
  je     L203
L57:
  lea    0x0(%rbp,%rax,1),%edx
  add    %esi,%eax
  movslq %edx,%rdx
  cltq
  movsd  (%rbx,%rdx,8),%xmm1
  mulsd  (%rbx,%rax,8),%xmm1
  addsd  %xmm1,%xmm0
L203:
  movsd  0x0(%r13,%r9,8),%xmm1
  add    $0x1,%r9
  add    %r12d,%esi
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rcx,%r15,8)
  andpd  %xmm4,%xmm1
  add    %r14,%rcx
  maxsd  %xmm9,%xmm1
  movapd %xmm1,%xmm9
  cmp    %r9d,%r12d
  jg     L158
  mov    0x28(%rsp),%r11
  movsd  (%rbx,%r11,1),%xmm1
L27:
  divsd  %xmm5,%xmm9
  maxsd  %xmm8,%xmm9
  movapd %xmm9,%xmm0
  mulsd  %xmm9,%xmm0
  comisd %xmm0,%xmm1
  ja     L277
  comisd %xmm9,%xmm7
  jbe    L271
  movapd %xmm7,%xmm0
  movapd %xmm7,%xmm9
  mulsd  %xmm7,%xmm0
L271:
  subsd  %xmm1,%xmm0
  maxsd  %xmm6,%xmm0
  movapd %xmm0,%xmm6
L150:
  movsd  %xmm9,(%rbx,%r11,1)
  cmp    0x38(%rsp),%r15
  je     L66
L347:
  mov    0x40(%rsp),%rax
  mov    0x8(%rsp),%ecx
  lea    0x0(%rbp,%r12,1),%esi
  mov    %esi,%ebp
  movsd  (%rax,%r11,1),%xmm0
  lea    0x2(%rcx),%edx
  divsd  %xmm9,%xmm0
  movsd  %xmm0,(%rax,%r11,1)
  cmp    %edx,%r12d
  jle    L344
  mov    0x54(%rsp),%r10d
  lea    (%r12,%rsi,1),%eax
  xor    %r9d,%r9d
  cltq
  mov    %r9,%r8
  mov    %r10d,%edi
  add    %r15,%rax
  sub    %ecx,%edi
  shl    $0x3,%rax
  mov    %edi,%ecx
  sub    $0x1,%rcx
  imul   %r14,%rcx
  test   %r12d,%r12d
  cmovns %rcx,%r8
  add    %rax,%r8
  cmp    %r11,%r8
  setl   %r8b
  test   %r12d,%r12d
  cmovns %r9,%rcx
  add    %rax,%rcx
  cmp    %r11,%rcx
  setg   %cl
  or     %r8b,%cl
  je     L48
  mov    %r10d,%ecx
  sub    0x50(%rsp),%ecx
  cmp    $0x2,%ecx
  jbe    L48
  movsd  (%rbx,%r11,1),%xmm1
  mov    %edi,%ecx
  mov    0x58(%rsp),%r9
  add    %rbx,%rax
  shr    %ecx
  xor    %r8d,%r8d
  unpcklpd %xmm1,%xmm1
  nopw   0x0(%rax,%rax,1)
L117:
  movsd  (%rax),%xmm0
  add    $0x1,%r8d
  movhpd (%rax,%r14,1),%xmm0
  divpd  %xmm1,%xmm0
  movlpd %xmm0,(%rax)
  movhpd %xmm0,(%rax,%r14,1)
  add    %r9,%rax
  cmp    %ecx,%r8d
  jne    L117
  mov    %edi,%eax
  and    $0xfffffffe,%eax
  add    %eax,%edx
  cmp    %edi,%eax
  je     L344
  imul   %r12d,%edx
  add    0x10(%rsp),%edx
  movslq %edx,%rdx
  lea    (%rbx,%rdx,8),%rax
  movsd  (%rax),%xmm0
  divsd  (%rbx,%r11,1),%xmm0
  movsd  %xmm0,(%rax)
L344:
  mov    %r15d,%eax
  test   %r15d,%r15d
  je     L97
L68:
  mov    0x50(%rsp),%eax
  mov    0x18(%rsp),%rdi
  movapd %xmm3,%xmm0
  shr    %eax
  sub    $0x1,%eax
  shl    $0x4,%rax
  lea    0x10(%rdi,%rax,1),%rdx
  mov    %rdi,%rax
  cs nopw 0x0(%rax,%rax,1)
L280:
  movupd (%rax),%xmm2
  add    $0x10,%rax
  mulpd  %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  unpckhpd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  cmp    %rdx,%rax
  jne    L280
  mov    0x50(%rsp),%edi
  mov    %edi,%eax
  and    $0xfffffffe,%eax
  and    $0x1,%edi
  je     L219
L315:
  add    %esi,%eax
  cltq
  movsd  (%rbx,%rax,8),%xmm1
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm0
L219:
  addq   $0x1,0x20(%rsp)
  add    $0x1,%r15
  add    %r14,%r13
  add    0x48(%rsp),%r11
  add    %r14,0x18(%rsp)
  jmp    L255
  nopl   0x0(%rax,%rax,1)
L211:
  movapd %xmm3,%xmm0
  jmp    L203
L5:
  movapd %xmm3,%xmm0
  xor    %eax,%eax
  jmp    L57
L277:
  ucomisd %xmm1,%xmm3
  ja     L275
  sqrtsd %xmm1,%xmm1
  movapd %xmm1,%xmm9
  movsd  %xmm9,(%rbx,%r11,1)
  cmp    0x38(%rsp),%r15
  jne    L347
L66:
  mov    0x60(%rsp),%rdi
  movsd  %xmm6,0x8(%rsp)
  call   <free@plt>
  movsd  0x8(%rsp),%xmm6
  add    $0x98,%rsp
  pop    %rbx
  pop    %rbp
  movapd %xmm6,%xmm0
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L48:
  add    %rbx,%rax
  nopw   0x0(%rax,%rax,1)
L316:
  add    $0x1,%edx
  movsd  (%rax),%xmm0
  divsd  (%rbx,%r11,1),%xmm0
  movsd  %xmm0,(%rax)
  add    %r14,%rax
  cmp    %edx,%r12d
  jne    L316
  mov    %r15d,%eax
  test   %r15d,%r15d
  jne    L68
L97:
  movapd %xmm3,%xmm0
  jmp    L315
L317:
  movapd %xmm3,%xmm9
  jmp    L27
L35:
  test   %r13d,%r13d
  jne    L206
  movsd  %xmm5,0x10(%rsp)
  mov    %ebp,%esi
  movsd  %xmm8,0x8(%rsp)
L355:
  call   <array_max_abs@plt>
  movsd  0x10(%rsp),%xmm5
  pxor   %xmm3,%xmm3
  ucomisd %xmm3,%xmm5
  jp     L41
  movsd  0x8(%rsp),%xmm8
  jne    L41
L89:
  ucomisd %xmm0,%xmm3
  ja     L11
  movapd %xmm0,%xmm5
  sqrtsd %xmm5,%xmm5
L243:
  movsd  D8(%rip),%xmm7
  mulsd  %xmm5,%xmm7
  test   %r12d,%r12d
  jg     L197
L41:
  movapd %xmm3,%xmm6
  jmp    L66
L206:
  lea    -0x1(%r13),%eax
  mov    %rbx,%rdi
  movsd  %xmm5,0x10(%rsp)
  lea    0x8(,%rax,8),%rdx
  mov    $0x8,%eax
  movsd  %xmm8,0x8(%rsp)
  cmovle %rax,%rdx
  xor    %esi,%esi
  call   <memset@plt>
  mov    0x60(%rsp),%rdi
  mov    %ebp,%esi
  jmp    L355
  nopw   0x0(%rax,%rax,1)
L275:
  movapd %xmm1,%xmm0
  mov    %r11,0x88(%rsp)
  movsd  %xmm3,0x80(%rsp)
  movsd  %xmm5,0x78(%rsp)
  movsd  %xmm6,0x70(%rsp)
  movsd  %xmm7,0x68(%rsp)
  movsd  %xmm8,0x28(%rsp)
  call   <sqrt@plt>
  movsd  0x28(%rsp),%xmm8
  movsd  0x68(%rsp),%xmm7
  movsd  0x70(%rsp),%xmm6
  movsd  0x78(%rsp),%xmm5
  movapd %xmm0,%xmm9
  movq   D49(%rip),%xmm4
  mov    0x88(%rsp),%r11
  movsd  0x80(%rsp),%xmm3
  jmp    L150
L11:
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm8,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x10(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm8
  movapd %xmm0,%xmm5
  jmp    L243
```

tool_calls: name: retrieve_data
arguments: data_tag: D11
type_infered: byte[8]
name: retrieve_data
arguments: data_tag: D8
type_infered: byte[8]
name: retrieve_data
arguments: data_tag: D49
type_infered: qword

tool: {"D11": [0x23,0x26,0xc1,0x61,0xea,0x4e,0x27,0x3f]}

tool: {"D8": [0x40,0xeb,0xd6,0x89,0x33,0xfa,0x60,0x3e]}

tool: D49: 9223372036854775807

assistant: ```c
double cholmod(double *A, int N, double *L, double maxinp) {
  int i, j, k, step, step2;
  double beta, d1, ls;
  double *U22;
  double maxadd, maxdiag, minl, minl2, minljj;
  U22 = (double *)malloc(sizeof(double) * N * N);
  minl = sqrt(sqrt((double)1e-15)) * maxinp;
  maxadd = 0.0;
  for (i = 0; i < N; ++i) {
    U22[i] = A[i + i * N];
  }
  for (i = 0; i < N * N; ++i) {
    L[i] = 0.0;
  }
  maxdiag = array_max_abs(U22, N);
  if (maxinp == 0.0) {
    maxinp = sqrt(maxdiag);
  }
  minl2 = sqrt((double)1e-15) * maxinp;
  for (j = 0; j < N; ++j) {
    step = j * N;
    ls = 0.0;
    for (i = 0; i < j; ++i) {
      ls += L[step + i] * L[step + i];
    }
    L[step + j] = A[step + j] - ls;
    minljj = 0.0;
    for (i = j + 1; i < N; ++i) {
      ls = 0.0;
      step2 = i * N;
      for (k = 0; k < j; ++k) {
        ls += L[step + k] * L[step2 + k];
      }
      L[step2 + j] = A[step + i] - ls;
      if (fabs(L[step2 + j]) > minljj) {
        minljj = fabs(L[step2 + j]);
      }
    }
    if (minljj / maxinp > minl) {
      minljj = minljj / maxinp;
    } else {
      minljj = minl;
    }
    if (L[step + j] > minljj * minljj) {
      L[step + j] = sqrt(L[step + j]);
    } else {
      if (minljj < minl2) {
        minljj = minl2;
      }
      if (maxadd < (minljj * minljj - L[step + j])) {
        maxadd = minljj * minljj - L[step + j];
      }
      L[step + j] = minljj;
    }
    for (i = j + 1; i < N; ++i) {
      L[i * N + j] /= L[step + j];
    }
  }
  free(U22);
  return maxadd;
}
```