user: What is the c source code of the assembly code below:

```asm
<divide_polynomial>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %rdi,%r12
  push   %rbp
  push   %rbx
  mov    %r8,%rbx
  sub    $0x48,%rsp
  mov    %rsi,(%rsp)
  mov    %rcx,0x38(%rsp)
  cmp    %r8,%r9
  jbe    L51
  test   %r8,%r8
  je     L123
  mov    %rcx,%rdi
  lea    0x1(%r12),%rcx
  lea    -0x1(%r8),%rax
  mov    %rdi,%rdx
  sub    %rcx,%rdx
  cmp    $0xe,%rdx
  jbe    L205
  cmp    $0x6,%rax
  jbe    L205
  cmp    $0xe,%rax
  jbe    L107
  mov    %r8,%rdx
  xor    %eax,%eax
  mov    %rdi,%rcx
  and    $0xfffffffffffffff0,%rdx
  nop
L150:
  movdqu (%r12,%rax,1),%xmm0
  movups %xmm0,(%rcx,%rax,1)
  add    $0x10,%rax
  cmp    %rdx,%rax
  jne    L150
  mov    %rbx,%rax
  and    $0xfffffffffffffff0,%rax
  mov    %eax,%esi
  cmp    %rax,%rbx
  je     L123
  mov    %rbx,%rcx
  sub    %rax,%rcx
  lea    -0x1(%rcx),%rdx
  cmp    $0x6,%rdx
  jbe    L272
L101:
  mov    (%r12,%rax,1),%rdx
  mov    0x38(%rsp),%rdi
  mov    %rdx,(%rdi,%rax,1)
  mov    %rcx,%rdx
  and    $0xfffffffffffffff8,%rdx
  add    %rdx,%rax
  add    %edx,%esi
  cmp    %rdx,%rcx
  je     L123
L272:
  movzbl (%r12,%rax,1),%edx
  mov    0x38(%rsp),%rcx
  mov    %dl,(%rcx,%rax,1)
  lea    0x1(%rsi),%edx
  movslq %edx,%rdx
  cmp    %rdx,%rbx
  jbe    L123
  movzbl (%r12,%rdx,1),%eax
  mov    %al,(%rcx,%rdx,1)
  lea    0x2(%rsi),%edx
  movslq %edx,%rdx
  cmp    %rdx,%rbx
  jbe    L123
  movzbl (%r12,%rdx,1),%eax
  mov    %al,(%rcx,%rdx,1)
  lea    0x3(%rsi),%eax
  cltq
  cmp    %rax,%rbx
  jbe    L123
  movzbl (%r12,%rax,1),%edx
  mov    %dl,(%rcx,%rax,1)
  lea    0x4(%rsi),%eax
  cltq
  cmp    %rax,%rbx
  jbe    L123
  movzbl (%r12,%rax,1),%edx
  mov    %dl,(%rcx,%rax,1)
  lea    0x5(%rsi),%eax
  cltq
  cmp    %rax,%rbx
  jbe    L123
  movzbl (%r12,%rax,1),%edx
  mov    %dl,(%rcx,%rax,1)
  lea    0x6(%rsi),%eax
  cltq
  cmp    %rax,%rbx
  jbe    L123
  movzbl (%r12,%rax,1),%edx
  mov    %dl,(%rcx,%rax,1)
L123:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L51:
  mov    %r8,%rdi
  mov    %r9,0x8(%rsp)
  mov    %rdx,%r14
  call   <malloc@plt>
  mov    0x8(%rsp),%r9
  test   %rax,%rax
  mov    %rax,%r13
  je     L192
L298:
  mov    %rbx,%rdi
  mov    %r9,0x8(%rsp)
  call   <malloc@plt>
  mov    0x8(%rsp),%r9
  test   %rax,%rax
  mov    %rax,%rbp
  je     L99
L159:
  lea    -0x1(%rbx),%r15
  mov    %r9,0x10(%rsp)
  lea    (%r15,%r9,1),%rdi
  imul   %rbx,%rdi
  call   <malloc@plt>
  mov    0x10(%rsp),%r9
  test   %rax,%rax
  mov    %rax,0x8(%rsp)
  je     L184
L234:
  lea    0x1(%rbx),%rax
  xor    %esi,%esi
  mov    %r14,%rdi
  mov    %r9,0x18(%rsp)
  sub    %r9,%rax
  mov    %rax,%rdx
  mov    %rax,0x10(%rsp)
  call   <memset@plt>
  mov    0x18(%rsp),%r9
  mov    0x8(%rsp),%rsi
  mov    %r15,%rcx
  mov    (%rsp),%rdi
  lea    -0x1(%r9),%rax
  mov    %rax,%rdx
  mov    %rax,0x28(%rsp)
  call   <gen_conv_matrix@plt>
  lea    -0x1(%rbx),%eax
  mov    0x18(%rsp),%r9
  cltq
  mov    %rax,0x20(%rsp)
  mov    $0x1,%eax
  sub    %r9,%rax
  mov    %r9d,0x18(%rsp)
  mov    %rax,0x30(%rsp)
  nopl   0x0(%rax)
L57:
  mov    0x10(%rsp),%rcx
  mov    0x8(%rsp),%rsi
  mov    %rbx,%r8
  mov    %r14,%rdx
  mov    %r13,%rdi
  xor    %r15d,%r15d
  call   <inner_product_set@plt>
  test   %rbx,%rbx
  je     L202
L35:
  movzbl 0x0(%r13,%r15,1),%esi
  movzbl (%r12,%r15,1),%edi
  call   <sum@plt>
  mov    %al,0x0(%rbp,%r15,1)
  add    $0x1,%r15
  cmp    %r15,%rbx
  jne    L35
L202:
  mov    0x20(%rsp),%r15
  jmp    L255
  nopl   (%rax)
L94:
  cmpb   $0x0,0x0(%rbp,%r15,1)
  lea    -0x1(%r15),%rax
  jne    L270
  mov    %rax,%r15
L255:
  lea    0x1(%r15),%ecx
  mov    %r15d,%esi
  cmp    $0xffffffff,%r15d
  jne    L94
  mov    $0xffffffffffffffff,%r15
  xor    %ecx,%ecx
L270:
  cmp    0x18(%rsp),%ecx
  jl     L285
  mov    0x28(%rsp),%rdx
  mov    (%rsp),%rax
  movzbl 0x0(%rbp,%r15,1),%edi
  movzbl (%rax,%rdx,1),%esi
  call   <divide@plt>
  mov    0x30(%rsp),%rcx
  lea    (%r14,%r15,1),%rdx
  mov    %al,(%rdx,%rcx,1)
  jmp    L57
  nopl   0x0(%rax)
L285:
  cmp    $0xffffffff,%esi
  je     L134
  lea    0x1(%rsi),%edx
  mov    0x38(%rsp),%rdi
  mov    %rbp,%rsi
  movslq %edx,%rdx
  call   <memcpy@plt>
L134:
  mov    %r13,%rdi
  call   <free@plt>
  mov    %rbp,%rdi
  call   <free@plt>
  mov    0x8(%rsp),%rdi
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10d0 <free@plt>
L205:
  mov    0x38(%rsp),%rcx
  xor    %eax,%eax
  nopw   0x0(%rax,%rax,1)
L62:
  movzbl (%r12,%rax,1),%edx
  mov    %dl,(%rcx,%rax,1)
  add    $0x1,%rax
  cmp    %rax,%rbx
  jne    L62
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L107:
  mov    %r8,%rcx
  xor    %esi,%esi
  xor    %eax,%eax
  jmp    L101
L99:
  xor    %eax,%eax
  call   <abort@plt>
  mov    0x8(%rsp),%r9
  jmp    L159
L184:
  call   <abort@plt>
  mov    0x10(%rsp),%r9
  jmp    L234
L192:
  xor    %eax,%eax
  call   <abort@plt>
  mov    0x8(%rsp),%r9
  jmp    L298
```

assistant: ```c
void divide_polynomial(unsigned char *dividend, unsigned char *divisor, unsigned char *quotient, unsigned char *remainder, size_t dividend_len, size_t divisor_len) {
  if (divisor_len > dividend_len) {
    for (int i = 0; i < dividend_len; i++) {
      remainder[i] = dividend[i];
    }
    return;
  }
  unsigned char *quotient_divisor_product = malloc(dividend_len);
  if (!quotient_divisor_product)
    abort();
  unsigned char *product_dividend_difference = malloc(dividend_len);
  if (!product_dividend_difference)
    abort();
  unsigned char *divisor_conv = malloc(dividend_len * (dividend_len + divisor_len - 1));
  if (!divisor_conv)
    abort();
  size_t quotient_len = dividend_len - (divisor_len - 1);
  memset(quotient, 0, quotient_len);
  gen_conv_matrix(divisor, divisor_conv, divisor_len - 1, dividend_len - 1);
  for (;;) {
    inner_product_set(quotient_divisor_product, divisor_conv, quotient, quotient_len, dividend_len);
    for (int i = 0; i < dividend_len; i++) {
      product_dividend_difference[i] = sum(dividend[i], quotient_divisor_product[i]);
    }
    int diff_greatest_nonzero = dividend_len;
    for (;;) {
      diff_greatest_nonzero--;
      if (diff_greatest_nonzero == -1)
        break;
      if (product_dividend_difference[diff_greatest_nonzero])
        break;
    }
    if (diff_greatest_nonzero < (int)divisor_len - 1) {
      for (int i = 0; i <= diff_greatest_nonzero; i++) {
        remainder[i] = product_dividend_difference[i];
      }
      free(quotient_divisor_product);
      free(product_dividend_difference);
      free(divisor_conv);
      return;
    }
    unsigned char factor = divide(product_dividend_difference[diff_greatest_nonzero], divisor[divisor_len - 1]);
    quotient[diff_greatest_nonzero - (divisor_len - 1)] = factor;
  }
}
```