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

```asm
<bandet>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,-0x48(%rsp)
  mov    %edx,-0x10(%rsp)
  test   %esi,%esi
  jle    L101
  test   %edx,%edx
  jle    L101
  lea    (%rsi,%rsi,1),%eax
  movslq %esi,%rdi
  mov    %esi,%r10d
  mov    -0x48(%rsp),%rsi
  lea    0x1(%rax),%ebx
  mov    %rdi,-0x40(%rsp)
  cltq
  xor    %r14d,%r14d
  mov    %ebx,-0x20(%rsp)
  movslq %ebx,%rbx
  lea    (%rsi,%rdi,8),%r8
  neg    %rax
  add    %rbx,%rdi
  mov    %rbx,-0x30(%rsp)
  movsd  (%r8),%xmm3
  lea    0x0(,%rbx,8),%rcx
  lea    (%rsi,%rdi,8),%r15
  mov    -0x10(%rsp),%edi
  lea    -0x1(%rbx),%r13
  mov    %rcx,-0x8(%rsp)
  mov    %rbx,-0x38(%rsp)
  mov    $0x1,%r11d
  add    $0x1,%edi
  mov    %edi,-0xc(%rsp)
  lea    0x0(,%rax,8),%rdi
L64:
  mov    -0x10(%rsp),%eax
  lea    -0x1(%r11),%ebx
  sub    %r11d,%eax
  cmp    %r10d,%eax
  cmovg  %r10d,%eax
  test   %eax,%eax
  jne    L108
L79:
  lea    0x1(%r11),%esi
  cmp    -0xc(%rsp),%esi
  je     L101
  movslq %r14d,%rdx
  mov    -0x48(%rsp),%rbx
  lea    -0x10(%r15),%rcx
  lea    -0x8(%r15),%rax
  add    -0x40(%rsp),%rdx
  cmp    %r11d,%r10d
  movsd  (%r15),%xmm3
  cmovle %r10d,%r11d
  lea    (%rbx,%rdx,8),%rdx
  lea    -0x1(%r11),%r9d
  shl    $0x3,%r9
  sub    %r9,%rcx
L62:
  movsd  (%rax),%xmm0
  mulsd  0x8(%rdx),%xmm0
  sub    $0x8,%rax
  add    %rdi,%rdx
  subsd  %xmm0,%xmm3
  cmp    %rax,%rcx
  jne    L62
  mov    -0x8(%rsp),%rax
  mov    -0x30(%rsp),%rbx
  movsd  %xmm3,(%r15)
  mov    %esi,%r11d
  add    %rbx,-0x38(%rsp)
  add    -0x20(%rsp),%r14d
  add    %rax,%r15
  add    %rax,%r8
  jmp    L64
  nopl   0x0(%rax,%rax,1)
L108:
  mov    -0x38(%rsp),%rsi
  mov    -0x30(%rsp),%rcx
  lea    0x1(%rax),%ebp
  mov    %r14d,%eax
  sub    -0x20(%rsp),%eax
  sub    %rsi,%rcx
  movslq %esi,%rdx
  cltq
  add    -0x40(%rsp),%rax
  lea    -0x1(%rcx,%rdx,1),%rsi
  mov    -0x48(%rsp),%rcx
  lea    (%rcx,%rax,8),%rax
  mov    $0x1,%ecx
  mov    %rax,-0x28(%rsp)
  lea    (%r10,%r14,1),%eax
  mov    %eax,-0x1c(%rsp)
  lea    -0x10(%r8),%rax
  mov    %rax,-0x18(%rsp)
  jmp    L100
  nopl   0x0(%rax)
L102:
  divsd  %xmm3,%xmm1
  add    $0x1,%rcx
  movsd  %xmm1,(%r8,%rsi,8)
  add    %r13,%rsi
  cmp    %rbp,%rcx
  je     L79
L100:
  mov    %r10d,%edx
  movsd  (%r8,%rsi,8),%xmm1
  sub    %ecx,%edx
  cmp    %ebx,%edx
  cmovg  %ebx,%edx
  test   %edx,%edx
  jle    L102
  mov    -0x1c(%rsp),%eax
  mov    -0x48(%rsp),%r9
  sub    $0x1,%edx
  shl    $0x3,%rdx
  add    %ecx,%eax
  cltq
  lea    (%r9,%rax,8),%r12
  mov    -0x18(%rsp),%r9
  lea    -0x8(%r8),%rax
  movsd  (%r12),%xmm2
  sub    %rdx,%r9
  mov    -0x28(%rsp),%rdx
  nopl   (%rax)
L85:
  movsd  (%rax),%xmm0
  mulsd  0x8(%rdx,%rcx,8),%xmm0
  subsd  %xmm0,%xmm2
  movsd  (%rax,%rsi,8),%xmm0
  mulsd  0x8(%rdx),%xmm0
  add    %rdi,%rdx
  sub    $0x8,%rax
  subsd  %xmm0,%xmm1
  cmp    %rax,%r9
  jne    L85
  movsd  %xmm2,(%r12)
  jmp    L102
  nopl   0x0(%rax,%rax,1)
L101:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void bandet(double *e, int m, int n) {
  int i, k, l, mi, lm, km, m2p1;
  double di, du, dl;
  m2p1 = 2 * m + 1;
  if (m <= 0) {
    return;
  }
  for (i = 1; i <= n; i++) {
    di = e[(i - 1) * m2p1 + m];
    mi = (m < (i - 1)) ? m : i - 1;
    if (mi >= 1) {
      for (k = 1; k <= mi; k++) {
        di -= e[(i - 1) * m2p1 - k + m] * e[(i - k - 1) * m2p1 + k + m];
      }
      e[(i - 1) * m2p1 + m] = di;
    }
    lm = (m < (n - i)) ? m : n - i;
    if (lm >= 1) {
      for (l = 1; l <= lm; l++) {
        dl = e[(i + l - 1) * m2p1 - l + m];
        km = ((m - l) < (i - 1)) ? m - l : i - 1;
        if (km >= 1) {
          du = e[(i - 1) * m2p1 + l + m];
          for (k = 1; k <= km; k++) {
            du -= e[(i - 1) * m2p1 - k + m] * e[(i - k - 1) * m2p1 + l + k + m];
            dl -= e[(l + i - 1) * m2p1 - l - k + m] * e[(i - k - 1) * m2p1 + k + m];
          }
          e[(i - 1) * m2p1 + l + m] = du;
        }
        e[(i + l - 1) * m2p1 - l + m] = dl / di;
      }
    }
  }
}
```