user: 请从给定的汇编指令推导出C源代码:

```asm
<poly1305_gmp>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,-0x38(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %rcx,-0x50(%rbp)
  mov    %r8d,-0x54(%rbp)
  mov    -0x20(%rbp),%eax
  mov    %eax,%edi
  call   <mpz_init@plt>
  mov    -0x1c(%rbp),%eax
  mov    %eax,%edi
  call   <mpz_init@plt>
  mov    -0x18(%rbp),%eax
  mov    %eax,%edi
  call   <mpz_init@plt>
  mov    -0x14(%rbp),%eax
  mov    %eax,%edi
  call   <mpz_init@plt>
  mov    -0x20(%rbp),%eax
  sub    $0x8,%rsp
  push   -0x40(%rbp)
  mov    $0x0,%r9d
  mov    $0x0,%r8d
  mov    $0x1,%ecx
  mov    $0xffffffff,%edx
  mov    $0x10,%esi
  mov    %eax,%edi
  call   <mpz_import@plt>
  add    $0x10,%rsp
  mov    -0x1c(%rbp),%eax
  mov    $0x0,%esi
  mov    %eax,%edi
  call   <mpz_set_ui@plt>
  mov    -0x14(%rbp),%eax
  mov    $0x1,%esi
  mov    %eax,%edi
  call   <mpz_set_ui@plt>
  mov    -0x14(%rbp),%ecx
  mov    -0x14(%rbp),%eax
  mov    $0x82,%edx
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <mpz_mul_2exp@plt>
  mov    -0x14(%rbp),%ecx
  mov    -0x14(%rbp),%eax
  mov    $0x5,%edx
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <mpz_sub_ui@plt>
  jmp    L31
L6:
  cmpl   $0xf,-0x54(%rbp)
  ja     L8
  mov    -0x54(%rbp),%eax
  mov    %eax,-0x24(%rbp)
  jmp    L54
L8:
  movl   $0x10,-0x24(%rbp)
L54:
  mov    -0x24(%rbp),%esi
  mov    -0x18(%rbp),%eax
  sub    $0x8,%rsp
  push   -0x50(%rbp)
  mov    $0x0,%r9d
  mov    $0x0,%r8d
  mov    $0x1,%ecx
  mov    $0xffffffff,%edx
  mov    %eax,%edi
  call   <mpz_import@plt>
  add    $0x10,%rsp
  mov    -0x24(%rbp),%eax
  add    %rax,-0x50(%rbp)
  mov    -0x24(%rbp),%eax
  sub    %eax,-0x54(%rbp)
  mov    -0x18(%rbp),%edx
  mov    -0x1c(%rbp),%ecx
  mov    -0x1c(%rbp),%eax
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <mpz_add@plt>
  mov    -0x18(%rbp),%eax
  mov    $0x1,%esi
  mov    %eax,%edi
  call   <mpz_set_ui@plt>
  mov    -0x24(%rbp),%eax
  shl    $0x3,%eax
  mov    %eax,%edx
  mov    -0x18(%rbp),%ecx
  mov    -0x18(%rbp),%eax
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <mpz_mul_2exp@plt>
  mov    -0x18(%rbp),%edx
  mov    -0x1c(%rbp),%ecx
  mov    -0x1c(%rbp),%eax
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <mpz_add@plt>
  mov    -0x20(%rbp),%edx
  mov    -0x1c(%rbp),%ecx
  mov    -0x1c(%rbp),%eax
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <mpz_mul@plt>
  mov    -0x14(%rbp),%edx
  mov    -0x1c(%rbp),%ecx
  mov    -0x1c(%rbp),%eax
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <mpz_tdiv_r@plt>
L31:
  cmpl   $0x0,-0x54(%rbp)
  jne    L6
  mov    -0x18(%rbp),%eax
  sub    $0x8,%rsp
  push   -0x48(%rbp)
  mov    $0x0,%r9d
  mov    $0x0,%r8d
  mov    $0x1,%ecx
  mov    $0xffffffff,%edx
  mov    $0x10,%esi
  mov    %eax,%edi
  call   <mpz_import@plt>
  add    $0x10,%rsp
  mov    -0x18(%rbp),%edx
  mov    -0x1c(%rbp),%ecx
  mov    -0x1c(%rbp),%eax
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <mpz_add@plt>
  movl   $0x0,-0x24(%rbp)
  jmp    L59
L56:
  mov    -0x24(%rbp),%edx
  mov    -0x38(%rbp),%rax
  lea    (%rdx,%rax,1),%rbx
  mov    -0x1c(%rbp),%ecx
  mov    -0x1c(%rbp),%eax
  mov    $0x100,%edx
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <mpz_tdiv_q_ui@plt>
  mov    %al,(%rbx)
  addl   $0x1,-0x24(%rbp)
L59:
  cmpl   $0xf,-0x24(%rbp)
  jbe    L56
  mov    -0x14(%rbp),%eax
  mov    %eax,%edi
  call   <mpz_clear@plt>
  mov    -0x18(%rbp),%eax
  mov    %eax,%edi
  call   <mpz_clear@plt>
  mov    -0x1c(%rbp),%eax
  mov    %eax,%edi
  call   <mpz_clear@plt>
  mov    -0x20(%rbp),%eax
  mov    %eax,%edi
  call   <mpz_clear@plt>
  nop
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

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);
}
```