user: Please show me the C implementation that matches the assembly code listed below:

```asm
<rbcholu>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,0x30(%rsp)
  mov    %esi,0x2c(%rsp)
  mov    %edx,0x3c(%rsp)
  mov    %rcx,0x58(%rsp)
  mov    %r8,0x8(%rsp)
  cmp    $0x20,%esi
  jle    L28
  mov    0x3c(%rsp),%ebx
  lea    0x1(%rbx),%eax
  shl    $0x5,%eax
  lea    -0x20(%rax),%edx
  cltq
  movslq %edx,%rdx
  shl    $0x3,%rax
  lea    0x100(,%rdx,8),%rsi
  movslq 0x2c(%rsp),%rdx
  mov    %rax,0x70(%rsp)
  mov    0x30(%rsp),%rax
  mov    %rsi,0x60(%rsp)
  lea    -0x100(,%rdx,8),%rdi
  mov    %rdx,%r14
  movslq %ebx,%rdx
  mov    %rdi,0x20(%rsp)
  lea    0x0(,%rdx,8),%r12
  mov    %rax,0x18(%rsp)
L87:
  lea    -0x20(%r14),%r13d
  mov    $0x100,%edi
  call   <malloc@plt>
  mov    %r13d,%r15d
  mov    $0x100,%edi
  shl    $0x8,%r15d
  mov    %rax,%rbp
  call   <malloc@plt>
  mov    %r15d,%edi
  mov    %rax,0x10(%rsp)
  call   <malloc@plt>
  mov    %r15d,%edi
  mov    %rax,0x50(%rsp)
  call   <malloc@plt>
  mov    0x18(%rsp),%r15
  mov    0x58(%rsp),%rcx
  mov    $0x20,%esi
  mov    0x3c(%rsp),%edx
  mov    %rax,0x68(%rsp)
  mov    %r15,%rdi
  call   <rcholu@plt>
  mov    0x8(%rsp),%rax
  mov    %r15,%rcx
  xor    %esi,%esi
  pxor   %xmm2,%xmm2
  mov    %rax,%rbx
  mov    %rax,%rdi
  nopw   0x0(%rax,%rax,1)
L13:
  mov    %rdi,%rdx
  xor    %eax,%eax
  nopl   (%rax)
L118:
  movsd  (%rcx,%rax,8),%xmm0
  add    $0x1,%rax
  add    $0x100,%rdx
  movsd  %xmm0,-0x100(%rdx)
  cmp    %eax,%r14d
  jg     L118
  add    $0x1,%esi
  add    $0x8,%rdi
  add    %r12,%rcx
  cmp    $0x20,%esi
  jne    L13
  mov    0x18(%rsp),%rax
  mov    0x8(%rsp),%r8
  mov    %rbp,%r11
  xor    %ecx,%ecx
  lea    0x8(%rbp),%rdi
  mov    0x50(%rsp),%r9
  mov    %r13d,0x38(%rsp)
  lea    0x100(%rbp),%r14
  lea    0x100(%rax),%rsi
  mov    0x10(%rsp),%rax
  mov    %rdi,0x48(%rsp)
  mov    %rbx,0x78(%rsp)
  sub    %rbp,%r9
  mov    %ecx,%ebx
  mov    %rsi,%r13
  lea    -0x8(%rax),%r10
  sub    %rax,%r8
  lea    0x100(%rax),%r15
  sub    %rax,%r11
  mov    $0x100,%eax
  sub    %rbp,%r10
  sub    %rbp,%rax
  mov    %rax,0x40(%rsp)
L8:
  mov    0x10(%rsp),%rax
  nopl   0x0(%rax)
L85:
  movsd  0x2000(%r8,%rax,1),%xmm0
  movsd  %xmm0,(%rax)
  movq   $0x0,(%r11,%rax,1)
  add    $0x8,%rax
  cmp    %r15,%rax
  jne    L85
  mov    0x8(%rsp),%rdi
  mov    0x40(%rsp),%rcx
  xor    %esi,%esi
  movapd %xmm2,%xmm1
  mov    0x48(%rsp),%rdx
  add    %rdi,%rcx
  nopw   0x0(%rax,%rax,1)
L71:
  movsd  (%r10,%rdx,1),%xmm0
  add    $0x1,%esi
  subsd  %xmm1,%xmm0
  divsd  (%rdi),%xmm0
  movsd  %xmm0,-0x8(%rdx)
  cmp    $0x20,%esi
  je     L44
  mov    %rbp,%rax
  movapd %xmm2,%xmm1
  nopw   0x0(%rax,%rax,1)
L67:
  movsd  (%rcx,%rax,1),%xmm0
  mulsd  (%rax),%xmm0
  add    $0x8,%rax
  addsd  %xmm0,%xmm1
  cmp    %rax,%rdx
  jne    L67
  add    $0x108,%rdi
  add    $0x100,%rcx
  add    $0x8,%rdx
  jmp    L71
L44:
  mov    %r13,%rdx
  mov    %rbp,%rax
  nopl   0x0(%rax,%rax,1)
L1:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  movsd  %xmm0,(%rdx)
  add    %r12,%rdx
  movsd  %xmm0,-0x8(%rax,%r9,1)
  cmp    %rax,%r14
  jne    L1
  add    $0x1,%ebx
  add    $0x8,%r13
  add    $0x100,%r9
  add    $0x100,%r8
  cmp    %ebx,0x38(%rsp)
  jg     L8
  mov    0x38(%rsp),%r13d
  mov    0x68(%rsp),%r14
  mov    $0x20,%edx
  mov    0x50(%rsp),%r15
  mov    0x78(%rsp),%rbx
  mov    %r14,%rcx
  mov    %r13d,%esi
  mov    %r15,%rdi
  call   <mtranspose@plt>
  mov    0x8(%rsp),%rdx
  mov    %r13d,%ecx
  mov    %r13d,%r9d
  mov    $0x20,%r8d
  mov    %r14,%rsi
  mov    %r15,%rdi
  call   <mmult@plt>
  mov    %r15,%rdi
  call   <free@plt>
  mov    %r14,%rdi
  call   <free@plt>
  mov    0x10(%rsp),%rdi
  call   <free@plt>
  mov    %rbp,%rdi
  call   <free@plt>
  mov    0x18(%rsp),%rdx
  xor    %ecx,%ecx
  add    0x70(%rsp),%rdx
  cmp    $0x0,%r13d
  jle    L7
L139:
  mov    %rcx,%rax
  nopl   0x0(%rax)
L34:
  movsd  (%rdx,%rax,8),%xmm0
  subsd  (%rbx,%rax,8),%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%r13d
  jg     L34
L7:
  add    $0x1,%rcx
  add    %r12,%rdx
  add    0x20(%rsp),%rbx
  cmp    %ecx,%r13d
  jg     L139
  subq   $0x100,0x20(%rsp)
  mov    0x60(%rsp),%rbx
  add    %rbx,0x18(%rsp)
  cmp    $0x20,%r13d
  jg     L70
  mov    0x2c(%rsp),%esi
  lea    -0x21(%rsi),%eax
  shr    $0x5,%eax
  lea    0x1(%rax),%edx
  neg    %eax
  shl    $0x5,%eax
  imul   %rbx,%rdx
  add    %rdx,0x30(%rsp)
  lea    -0x20(%rsi,%rax,1),%eax
  mov    %eax,0x2c(%rsp)
L28:
  mov    0x58(%rsp),%rcx
  mov    0x3c(%rsp),%edx
  mov    0x2c(%rsp),%esi
  mov    0x30(%rsp),%rdi
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10d0 <rcholu@plt>
L70:
  mov    %r13d,%r14d
  jmp    L87
```

assistant: ```c
int rbcholu(double *A, int N, int stride, double *UB, double *UT) {
  int bs, bb, i, j, Nb, t, k, u, v, w, sc;
  double *b, *x, *U12, *U12T;
  double sum;
  bs = (int)32;
  bb = bs * bs;
  if (N <= 32) {
    sc = rcholu(A, N, stride, UB);
    if (sc == -1) {
      return -1;
    }
  } else {
    Nb = N - bs;
    x = (double *)malloc(sizeof(double) * bs);
    b = (double *)malloc(sizeof(double) * bs);
    U12T = (double *)malloc(sizeof(double) * Nb * bs);
    U12 = (double *)malloc(sizeof(double) * Nb * bs);
    rcholu(A, bs, stride, UB);
    for (i = 0; i < bs; ++i) {
      t = i * stride;
      u = 0;
      for (j = 0; j < N; ++j) {
        UT[u + i] = A[j + t];
        u += bs;
      }
    }
    for (k = 0; k < Nb; ++k) {
      u = k * bs;
      for (i = 0; i < bs; ++i) {
        b[i] = UT[bb + u + i];
        x[i] = 0.;
      }
      for (i = 0; i < bs; ++i) {
        t = i * bs;
        sum = 0;
        for (j = 0; j < i; ++j) {
          sum += UT[t + j] * x[j];
        }
        x[i] = (b[i] - sum) / UT[t + i];
      }
      v = bs + k;
      for (i = 0; i < bs; ++i) {
        A[v] = x[i];
        U12T[u + i] = x[i];
        v += stride;
      }
    }
    mtranspose(U12T, Nb, bs, U12);
    mmult(U12T, U12, UT, Nb, bs, Nb);
    free(U12T);
    free(U12);
    free(b);
    free(x);
    for (i = 0; i < Nb; ++i) {
      u = bs * stride + bs + i * stride;
      w = i * Nb;
      for (j = i; j < Nb; j++) {
        A[j + u] -= UT[j + w];
      }
    }
    sc = rbcholu(A + bs * stride + bs, Nb, stride, UB, UT);
    if (sc == -1) {
      return -1;
    }
  }
  return sc;
}
```