user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<bn_mul_normal>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %rdi,-0x38(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    %edx,-0x44(%rbp)
  mov    %rcx,-0x50(%rbp)
  mov    %r8d,-0x48(%rbp)
  mov    -0x44(%rbp),%eax
  cmp    -0x48(%rbp),%eax
  jge    L20
  mov    -0x44(%rbp),%eax
  mov    %eax,-0x24(%rbp)
  mov    -0x48(%rbp),%eax
  mov    %eax,-0x44(%rbp)
  mov    -0x24(%rbp),%eax
  mov    %eax,-0x48(%rbp)
  mov    -0x40(%rbp),%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x50(%rbp),%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x18(%rbp),%rax
  mov    %rax,-0x50(%rbp)
L20:
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x20(%rbp)
  mov    -0x50(%rbp),%rax
  mov    (%rax),%rcx
  mov    -0x44(%rbp),%edx
  mov    -0x40(%rbp),%rsi
  mov    -0x38(%rbp),%rax
  mov    %rax,%rdi
  call   <bn_mul_words@plt>
  mov    -0x20(%rbp),%rdx
  mov    %rax,(%rdx)
L89:
  subl   $0x1,-0x48(%rbp)
  cmpl   $0x0,-0x48(%rbp)
  jle    L78
  mov    -0x50(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%rcx
  mov    -0x38(%rbp),%rax
  lea    0x8(%rax),%rdi
  mov    -0x20(%rbp),%rax
  lea    0x8(%rax),%rbx
  mov    -0x44(%rbp),%edx
  mov    -0x40(%rbp),%rax
  mov    %rax,%rsi
  call   <bn_mul_add_words@plt>
  mov    %rax,(%rbx)
  subl   $0x1,-0x48(%rbp)
  cmpl   $0x0,-0x48(%rbp)
  jle    L68
  mov    -0x50(%rbp),%rax
  add    $0x10,%rax
  mov    (%rax),%rcx
  mov    -0x38(%rbp),%rax
  lea    0x10(%rax),%rdi
  mov    -0x20(%rbp),%rax
  lea    0x10(%rax),%rbx
  mov    -0x44(%rbp),%edx
  mov    -0x40(%rbp),%rax
  mov    %rax,%rsi
  call   <bn_mul_add_words@plt>
  mov    %rax,(%rbx)
  subl   $0x1,-0x48(%rbp)
  cmpl   $0x0,-0x48(%rbp)
  jle    L6
  mov    -0x50(%rbp),%rax
  add    $0x18,%rax
  mov    (%rax),%rcx
  mov    -0x38(%rbp),%rax
  lea    0x18(%rax),%rdi
  mov    -0x20(%rbp),%rax
  lea    0x18(%rax),%rbx
  mov    -0x44(%rbp),%edx
  mov    -0x40(%rbp),%rax
  mov    %rax,%rsi
  call   <bn_mul_add_words@plt>
  mov    %rax,(%rbx)
  subl   $0x1,-0x48(%rbp)
  cmpl   $0x0,-0x48(%rbp)
  jle    L5
  mov    -0x50(%rbp),%rax
  add    $0x20,%rax
  mov    (%rax),%rcx
  mov    -0x38(%rbp),%rax
  lea    0x20(%rax),%rdi
  mov    -0x20(%rbp),%rax
  lea    0x20(%rax),%rbx
  mov    -0x44(%rbp),%edx
  mov    -0x40(%rbp),%rax
  mov    %rax,%rsi
  call   <bn_mul_add_words@plt>
  mov    %rax,(%rbx)
  addq   $0x20,-0x20(%rbp)
  addq   $0x20,-0x38(%rbp)
  addq   $0x20,-0x50(%rbp)
  jmp    L89
L78:
  nop
  jmp    L14
L68:
  nop
  jmp    L14
L6:
  nop
  jmp    L14
L5:
  nop
L14:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

assistant: ```c
void bn_mul_normal(unsigned long long *r, unsigned long long *a, int na, unsigned long long *b, int nb) {
  unsigned long long *rr;
  if (na < nb) {
    int itmp;
    unsigned long long *ltmp;
    itmp = na;
    na = nb;
    nb = itmp;
    ltmp = a;
    a = b;
    b = ltmp;
  }
  rr = &(r[na]);
  rr[0] = bn_mul_words(r, a, na, b[0]);
  for (;;) {
    if (--nb <= 0)
      return;
    rr[1] = bn_mul_add_words(&(r[1]), a, na, b[1]);
    if (--nb <= 0)
      return;
    rr[2] = bn_mul_add_words(&(r[2]), a, na, b[2]);
    if (--nb <= 0)
      return;
    rr[3] = bn_mul_add_words(&(r[3]), a, na, b[3]);
    if (--nb <= 0)
      return;
    rr[4] = bn_mul_add_words(&(r[4]), a, na, b[4]);
    rr += 4;
    r += 4;
    b += 4;
  }
}
```