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

```asm
<poly1305_gmp>:
  endbr64
  push   %r15
  xor    %r11d,%r11d
  xor    %r15d,%r15d
  push   %r14
  mov    %rsi,%r14
  push   %r13
  mov    %rcx,%r13
  push   %r12
  mov    %r8d,%r12d
  push   %rbp
  xor    %ebp,%ebp
  push   %rbx
  mov    %rdi,%rbx
  mov    %r11d,%edi
  sub    $0x18,%rsp
  mov    %rdx,(%rsp)
  mov    %r11d,0xc(%rsp)
  call   <mpz_init@plt>
  mov    %ebp,%edi
  call   <mpz_init@plt>
  mov    %r15d,%edi
  call   <mpz_init@plt>
  mov    %r15d,%edi
  call   <mpz_init@plt>
  sub    $0x8,%rsp
  xor    %r9d,%r9d
  xor    %r8d,%r8d
  push   %r14
  mov    0x1c(%rsp),%edi
  mov    $0x1,%ecx
  mov    $0xffffffff,%edx
  mov    $0x10,%esi
  call   <mpz_import@plt>
  xor    %esi,%esi
  mov    %ebp,%edi
  call   <mpz_set_ui@plt>
  mov    $0x1,%esi
  mov    %r15d,%edi
  call   <mpz_set_ui@plt>
  mov    $0x82,%edx
  mov    %r15d,%esi
  mov    %r15d,%edi
  call   <mpz_mul_2exp@plt>
  mov    $0x5,%edx
  mov    %r15d,%esi
  mov    %r15d,%edi
  call   <mpz_sub_ui@plt>
  pop    %r9
  pop    %r10
  test   %r12d,%r12d
  je     L26
  mov    %r12d,%r14d
  and    $0xf,%r14d
  nopl   0x0(%rax)
L16:
  cmp    %r14d,%r12d
  jne    L42
  sub    $0x8,%rsp
  xor    %r15d,%r15d
  xor    %ebp,%ebp
  xor    %r9d,%r9d
  push   %r13
  xor    %r8d,%r8d
  mov    $0x1,%ecx
  mov    $0xffffffff,%edx
  mov    %r12d,%esi
  mov    %r15d,%edi
  call   <mpz_import@plt>
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebp,%edi
  call   <mpz_add@plt>
  mov    $0x1,%esi
  mov    %r15d,%edi
  call   <mpz_set_ui@plt>
  lea    0x0(,%r12,8),%edx
  mov    %r15d,%esi
  mov    %r15d,%edi
  call   <mpz_mul_2exp@plt>
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebp,%edi
  call   <mpz_add@plt>
  xor    %edx,%edx
  mov    %ebp,%esi
  mov    %ebp,%edi
  call   <mpz_mul@plt>
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebp,%edi
  call   <mpz_tdiv_r@plt>
  pop    %rax
  pop    %rdx
L26:
  sub    $0x8,%rsp
  xor    %r15d,%r15d
  mov    $0x1,%ecx
  xor    %r9d,%r9d
  push   0x8(%rsp)
  xor    %r8d,%r8d
  mov    $0xffffffff,%edx
  mov    %r15d,%edi
  mov    $0x10,%esi
  xor    %ebp,%ebp
  lea    0x10(%rbx),%r12
  call   <mpz_import@plt>
  mov    %ebp,%esi
  mov    %r15d,%edx
  mov    %ebp,%edi
  call   <mpz_add@plt>
  pop    %rcx
  pop    %rsi
  nopl   0x0(%rax)
L47:
  xor    %ebp,%ebp
  mov    $0x100,%edx
  add    $0x1,%rbx
  mov    %ebp,%esi
  mov    %ebp,%edi
  call   <mpz_tdiv_q_ui@plt>
  mov    %al,-0x1(%rbx)
  cmp    %r12,%rbx
  jne    L47
  xor    %r15d,%r15d
  mov    %r15d,%edi
  call   <mpz_clear@plt>
  mov    %r15d,%edi
  call   <mpz_clear@plt>
  mov    %ebp,%edi
  call   <mpz_clear@plt>
  add    $0x18,%rsp
  xor    %edi,%edi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1110 <mpz_clear@plt>
  nopl   0x0(%rax,%rax,1)
L42:
  sub    $0x8,%rsp
  xor    %r15d,%r15d
  xor    %r8d,%r8d
  xor    %ebp,%ebp
  push   %r13
  xor    %r9d,%r9d
  mov    $0x1,%ecx
  mov    $0xffffffff,%edx
  mov    $0x10,%esi
  mov    %r15d,%edi
  sub    $0x10,%r12d
  add    $0x10,%r13
  call   <mpz_import@plt>
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebp,%edi
  call   <mpz_add@plt>
  mov    $0x1,%esi
  mov    %r15d,%edi
  call   <mpz_set_ui@plt>
  mov    $0x80,%edx
  mov    %r15d,%esi
  mov    %r15d,%edi
  call   <mpz_mul_2exp@plt>
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebp,%edi
  call   <mpz_add@plt>
  xor    %edx,%edx
  mov    %ebp,%esi
  mov    %ebp,%edi
  call   <mpz_mul@plt>
  mov    %ebp,%edi
  mov    %r15d,%edx
  mov    %ebp,%esi
  call   <mpz_tdiv_r@plt>
  pop    %rdi
  pop    %r8
  test   %r12d,%r12d
  jne    L16
  jmp    L26
```

assistant: ```c
void poly1305_gmp(unsigned char out[16], const unsigned char r[16], const unsigned char s[16], const unsigned char m[], unsigned int l) {
  unsigned int j;
  mpz_t rbar;
  mpz_t h;
  mpz_t c;
  mpz_t p;
  mpz_init(rbar);
  mpz_init(h);
  mpz_init(c);
  mpz_init(p);
  mpz_import(rbar, 16, -1, 1, 0, 0, r);
  mpz_set_ui(h, 0);
  mpz_set_ui(p, 1);
  mpz_mul_2exp(p, p, 130);
  mpz_sub_ui(p, p, 5);
  while (l > 0) {
    if (l < 16)
      j = l;
    else
      j = 16;
    mpz_import(c, j, -1, 1, 0, 0, m);
    m += j;
    l -= j;
    mpz_add(h, h, c);
    mpz_set_ui(c, 1);
    mpz_mul_2exp(c, c, 8 * j);
    mpz_add(h, h, c);
    mpz_mul(h, h, rbar);
    mpz_tdiv_r(h, h, p);
  }
  mpz_import(c, 16, -1, 1, 0, 0, s);
  mpz_add(h, h, c);
  for (j = 0; j < 16; ++j)
    out[j] = mpz_tdiv_q_ui(h, h, 256);
  mpz_clear(p);
  mpz_clear(c);
  mpz_clear(h);
  mpz_clear(rbar);
}
```