user: 请从给定的汇编指令推导出C源代码:

```asm
<chol_eps>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  movsd  %xmm0,0x20(%rsp)
  test   %edi,%edi
  jle    L48
  movslq %edi,%rax
  mov    %edi,%r8d
  lea    0x1(%rdi),%r11d
  mov    %rsi,%r9
  mov    %rdx,%rbx
  xor    %r13d,%r13d
  mov    $0x1,%edi
  mov    %r8d,%r15d
  lea    0x8(,%rax,8),%rax
  movl   $0x0,0x38(%rsp)
  pxor   %xmm3,%xmm3
  mov    %rax,0x28(%rsp)
  lea    0x8(%rdx),%rax
  mov    %rax,0x10(%rsp)
  nopl   (%rax)
L0:
  movslq 0x38(%rsp),%rax
  movslq %r13d,%rsi
  xor    %r12d,%r12d
  lea    -0x1(%rdi),%r10d
  shl    $0x3,%rsi
  mov    %r15d,%ebp
  movapd %xmm3,%xmm0
  xor    %r8d,%r8d
  lea    (%r9,%rax,8),%rax
  lea    (%rbx,%rsi,1),%r14
  mov    %rax,0x8(%rsp)
  cmp    %r12d,%r10d
  jle    L79
  nopl   0x0(%rax,%rax,1)
L66:
  movsd  (%r9,%rsi,1),%xmm1
  movslq %r8d,%rax
  subsd  %xmm0,%xmm1
  divsd  (%rbx,%rax,8),%xmm1
L50:
  add    $0x1,%r12d
  movsd  %xmm1,(%rbx,%rsi,1)
  cmp    %r12d,%edi
  je     L107
L11:
  cmp    $0x1,%r12d
  je     L16
  mov    %r12d,%edx
  movslq %ebp,%rax
  movapd %xmm3,%xmm0
  shr    %edx
  lea    (%rbx,%rax,8),%rcx
  xor    %eax,%eax
  shl    $0x4,%rdx
  nopl   0x0(%rax)
L90:
  movupd (%rcx,%rax,1),%xmm2
  movupd (%r14,%rax,1),%xmm4
  add    $0x10,%rax
  mulpd  %xmm4,%xmm2
  addsd  %xmm2,%xmm0
  unpckhpd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  cmp    %rax,%rdx
  jne    L90
  mov    %r12d,%eax
  and    $0xfffffffe,%eax
  test   $0x1,%r12b
  je     L123
L81:
  lea    0x0(%rbp,%rax,1),%edx
  add    %r13d,%eax
  movslq %edx,%rdx
  cltq
  movsd  (%rbx,%rdx,8),%xmm1
  mulsd  (%rbx,%rax,8),%xmm1
  addsd  %xmm1,%xmm0
L123:
  add    %r11d,%r8d
  add    %r15d,%ebp
  add    $0x8,%rsi
  cmp    %r12d,%r10d
  jg     L66
L79:
  mov    0x8(%rsp),%rax
  movsd  0x20(%rsp),%xmm1
  addsd  (%rax),%xmm1
  subsd  %xmm0,%xmm1
  ucomisd %xmm1,%xmm3
  ja     L9
  sqrtsd %xmm1,%xmm1
  add    $0x1,%r12d
  movsd  %xmm1,(%rbx,%rsi,1)
  cmp    %r12d,%edi
  jne    L11
L107:
  add    %r11d,0x38(%rsp)
  lea    0x1(%rdi),%ebp
  add    %r15d,%r13d
  mov    %r9,0x18(%rsp)
  mov    %r11d,0x8(%rsp)
  cmp    %r15d,%edi
  je     L48
  mov    0x10(%rsp),%r14
  mov    %r15d,%eax
  xor    %esi,%esi
  sub    %ebp,%eax
  mov    %r14,%rdi
  lea    0x8(,%rax,8),%rdx
  call   <memset@plt>
  mov    0x28(%rsp),%rdi
  mov    0x18(%rsp),%r9
  pxor   %xmm3,%xmm3
  mov    0x8(%rsp),%r11d
  add    %rdi,%rax
  mov    %ebp,%edi
  mov    %rax,0x10(%rsp)
  jmp    L0
L16:
  xor    %eax,%eax
  movapd %xmm3,%xmm0
  jmp    L81
L48:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L9:
  movapd %xmm1,%xmm0
  mov    %r10d,0x4c(%rsp)
  mov    %r9,0x40(%rsp)
  mov    %edi,0x48(%rsp)
  mov    %r11d,0x3c(%rsp)
  mov    %rsi,0x30(%rsp)
  mov    %r8d,0x18(%rsp)
  call   <sqrt@plt>
  mov    0x4c(%rsp),%r10d
  mov    0x40(%rsp),%r9
  pxor   %xmm3,%xmm3
  mov    0x48(%rsp),%edi
  mov    0x3c(%rsp),%r11d
  movapd %xmm0,%xmm1
  mov    0x30(%rsp),%rsi
  mov    0x18(%rsp),%r8d
  jmp    L50
```

assistant: ```c
void chol_eps(int n, const double A[], double eps, double L[]) {
  for (int i = 0; i < n; ++i) {
    for (int j = 0; j < i + 1; ++j) {
      double s = 0;
      for (int k = 0; k < j; ++k)
        s += L[i * n + k] * L[j * n + k];
      if (j < i)
        L[i * n + j] = (A[i * n + j] - s) / L[j * n + j];
      else
        L[i * n + j] = sqrt(A[i * n + i] + eps - s);
    }
    for (int j = i + 1; j < n; ++j)
      L[i * n + j] = 0;
  }
}
```