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

```asm
<k_cholsl>:
  endbr64
  push   %r15
  mov    %esi,%eax
  sub    $0x1,%esi
  mov    %rdi,%r11
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %rdx,%rbp
  push   %rbx
  mov    %esi,-0xc(%rsp)
  test   %eax,%eax
  jle    L114
  movss  (%rcx),%xmm0
  mov    %rcx,%r14
  divss  (%rdx),%xmm0
  movss  %xmm0,(%r8)
  cmp    $0x1,%eax
  je     L56
  movslq %eax,%rbx
  lea    -0xc(%r8),%rdx
  mov    %eax,%r10d
  mov    $0x1,%esi
  mov    %rbx,-0x18(%rsp)
  lea    0x4(,%rbx,4),%r15
  xor    %r9d,%r9d
  lea    -0x10(%r11,%r15,1),%rcx
  nopw   0x0(%rax,%rax,1)
L77:
  movss  (%r14,%rsi,4),%xmm1
  mov    %r9d,%edi
  mov    %esi,%ebx
  cmp    $0x2,%r9d
  jbe    L91
  mov    %esi,%r12d
  xor    %edi,%edi
  shr    $0x2,%r12d
  neg    %r12
  shl    $0x4,%r12
  nopl   0x0(%rax)
L3:
  movups (%rdx,%rdi,1),%xmm0
  movups (%rcx,%rdi,1),%xmm2
  sub    $0x10,%rdi
  shufps $0x1b,%xmm2,%xmm2
  shufps $0x1b,%xmm0,%xmm0
  mulps  %xmm2,%xmm0
  subss  %xmm0,%xmm1
  movaps %xmm0,%xmm2
  shufps $0x55,%xmm0,%xmm2
  subss  %xmm2,%xmm1
  movaps %xmm0,%xmm2
  unpckhps %xmm0,%xmm2
  shufps $0xff,%xmm0,%xmm0
  subss  %xmm2,%xmm1
  subss  %xmm0,%xmm1
  cmp    %r12,%rdi
  jne    L3
  mov    %ebx,%r12d
  mov    %r9d,%edi
  and    $0xfffffffc,%r12d
  sub    %r12d,%edi
  cmp    %ebx,%r12d
  je     L28
L91:
  lea    (%r10,%rdi,1),%ebx
  movslq %edi,%r12
  movslq %ebx,%rbx
  lea    0x0(,%r12,4),%r13
  movss  (%r11,%rbx,4),%xmm0
  mulss  (%r8,%r12,4),%xmm0
  lea    -0x1(%rdi),%ebx
  subss  %xmm0,%xmm1
  test   %edi,%edi
  je     L28
  add    %r10d,%ebx
  movslq %ebx,%rbx
  movss  (%r11,%rbx,4),%xmm0
  mulss  -0x4(%r8,%r13,1),%xmm0
  lea    -0x2(%rdi),%ebx
  subss  %xmm0,%xmm1
  cmp    $0x1,%edi
  je     L28
  add    %r10d,%ebx
  movslq %ebx,%rbx
  movss  (%r11,%rbx,4),%xmm0
  mulss  -0x8(%r8,%r13,1),%xmm0
  subss  %xmm0,%xmm1
L28:
  divss  0x0(%rbp,%rsi,4),%xmm1
  add    $0x1,%r9d
  movss  %xmm1,0x10(%rdx)
  add    %eax,%r10d
  add    $0x1,%rsi
  add    %r15,%rcx
  add    $0x4,%rdx
  cmp    -0x18(%rsp),%rsi
  jne    L77
L104:
  mov    %eax,%r10d
  mov    -0x18(%rsp),%rsi
  movslq -0xc(%rsp),%rdi
  mov    %eax,%r9d
  imul   %eax,%r10d
  shl    $0x2,%rsi
  xchg   %ax,%ax
L121:
  movss  (%r8,%rdi,4),%xmm1
  cmp    %r9d,%eax
  jle    L96
  movslq %r10d,%rcx
  movslq %r9d,%rdx
  add    %rdi,%rcx
  lea    (%r11,%rcx,4),%rcx
  nopl   0x0(%rax,%rax,1)
L63:
  movss  (%rcx),%xmm0
  mulss  (%r8,%rdx,4),%xmm0
  add    $0x1,%rdx
  add    %rsi,%rcx
  subss  %xmm0,%xmm1
  cmp    %edx,%eax
  jg     L63
L96:
  divss  0x0(%rbp,%rdi,4),%xmm1
  movss  %xmm1,(%r8,%rdi,4)
  sub    $0x1,%rdi
  sub    $0x1,%r9d
  sub    %eax,%r10d
  test   %edi,%edi
  jns    L121
L135:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L114:
  mov    -0xc(%rsp),%edx
  test   %edx,%edx
  js     L135
  movslq %eax,%rsi
  mov    %rsi,-0x18(%rsp)
  jmp    L104
  nopl   0x0(%rax)
L56:
  movq   $0x1,-0x18(%rsp)
  jmp    L104
```

assistant: ```c
void k_cholsl(float *a, int n, float p[], float b[], float x[]) {
  int i, k;
  float sum;
  for (i = 0; i < n; i++) {
    for (sum = b[i], k = i - 1; k >= 0; k--)
      sum -= a[i * n + k] * x[k];
    x[i] = sum / p[i];
  }
  for (i = n - 1; i >= 0; i--) {
    for (sum = x[i], k = i + 1; k < n; k++)
      sum -= a[k * n + i] * x[k];
    x[i] = sum / p[i];
  }
}
```