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

```asm
<mp_sqrt_newton>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x40,%rsp
  mov    %edi,-0x14(%rbp)
  mov    %esi,-0x18(%rbp)
  mov    %rdx,-0x20(%rbp)
  mov    %rcx,-0x28(%rbp)
  mov    %r8,-0x30(%rbp)
  mov    %r9,-0x38(%rbp)
  mov    0x10(%rbp),%eax
  sar    %eax
  mov    %eax,-0xc(%rbp)
  mov    -0xc(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x4(%rbp)
  cmpl   $0x1,-0xc(%rbp)
  jg     L9
  movl   $0x2,-0xc(%rbp)
L9:
  mov    -0x14(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  add    $0x1,%eax
  mov    %eax,-0x10(%rbp)
  mov    -0x14(%rbp),%eax
  sub    -0x4(%rbp),%eax
  cmp    %eax,-0x10(%rbp)
  jge    L5
  mov    -0x14(%rbp),%eax
  sub    -0x4(%rbp),%eax
  mov    %eax,-0x10(%rbp)
L5:
  mov    -0xc(%rbp),%eax
  sar    %eax
  lea    0x1(%rax),%edi
  mov    -0x30(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x10(%rbp),%eax
  mov    %rdx,%rcx
  mov    %edi,%edx
  mov    %eax,%edi
  call   <mp_round@plt>
  mov    0x38(%rbp),%rax
  mov    (%rax),%eax
  cmp    %eax,-0xc(%rbp)
  je     L58
  mov    0x18(%rbp),%r8
  mov    -0xc(%rbp),%edi
  mov    -0x38(%rbp),%rcx
  mov    -0x30(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x10(%rbp),%eax
  push   0x30(%rbp)
  push   0x28(%rbp)
  mov    %r8,%r9
  mov    %edi,%r8d
  mov    %eax,%edi
  call   <mp_squh@plt>
  add    $0x10,%rsp
  jmp    L61
L58:
  mov    0x28(%rbp),%r8
  mov    -0xc(%rbp),%edi
  mov    -0x38(%rbp),%rcx
  mov    0x18(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x10(%rbp),%eax
  sub    $0x8,%rsp
  push   0x30(%rbp)
  mov    %r8,%r9
  mov    %edi,%r8d
  mov    %eax,%edi
  call   <mp_squh_use_in1fft@plt>
  add    $0x10,%rsp
L61:
  mov    -0x28(%rbp),%rcx
  mov    -0x4(%rbp),%edx
  mov    -0x18(%rbp),%esi
  mov    -0x14(%rbp),%eax
  mov    %eax,%edi
  call   <mp_round@plt>
  mov    -0x38(%rbp),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x28(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x10(%rbp),%eax
  push   0x30(%rbp)
  push   0x28(%rbp)
  push   0x20(%rbp)
  push   0x18(%rbp)
  mov    0x10(%rbp),%r9d
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <mp_mulh@plt>
  add    $0x20,%rsp
  mov    -0x38(%rbp),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x30(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x10(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <mp_sub@plt>
  mov    -0x30(%rbp),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x30(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x10(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <mp_add@plt>
  mov    0x28(%rbp),%rdi
  mov    -0x38(%rbp),%rcx
  mov    0x18(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x14(%rbp),%eax
  sub    $0x8,%rsp
  push   0x30(%rbp)
  mov    %rdi,%r9
  mov    0x10(%rbp),%r8d
  mov    %eax,%edi
  call   <mp_squh_use_in1fft@plt>
  add    $0x10,%rsp
  mov    -0x38(%rbp),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x20(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x14(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <mp_sub@plt>
  mov    -0x20(%rbp),%rax
  add    $0x4,%rax
  mov    (%rax),%eax
  mov    -0x38(%rbp),%rdx
  add    $0x4,%rdx
  mov    (%rdx),%edx
  sub    %edx,%eax
  mov    %eax,-0x8(%rbp)
  mov    -0x20(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%edx
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%eax
  cmp    %eax,%edx
  jle    L6
  addl   $0x1,-0x8(%rbp)
L6:
  mov    -0x38(%rbp),%rax
  mov    (%rax),%eax
  test   %eax,%eax
  jne    L45
  mov    0x10(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x8(%rbp)
L45:
  mov    -0x30(%rbp),%rcx
  mov    -0x4(%rbp),%edx
  mov    -0x18(%rbp),%esi
  mov    -0x10(%rbp),%eax
  mov    %eax,%edi
  call   <mp_round@plt>
  mov    -0x38(%rbp),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x30(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x10(%rbp),%eax
  push   0x30(%rbp)
  push   0x28(%rbp)
  push   0x20(%rbp)
  push   0x18(%rbp)
  mov    0x10(%rbp),%r9d
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <mp_mulh@plt>
  add    $0x20,%rsp
  mov    0x38(%rbp),%rax
  mov    0x10(%rbp),%edx
  mov    %edx,(%rax)
  mov    -0x38(%rbp),%rcx
  mov    -0x38(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x10(%rbp),%eax
  mov    %eax,%edi
  call   <mp_idiv_2@plt>
  mov    -0x38(%rbp),%rax
  mov    (%rax),%eax
  test   %eax,%eax
  je     L4
  mov    -0x28(%rbp),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x28(%rbp),%rdx
  mov    -0x18(%rbp),%esi
  mov    -0x14(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <mp_add@plt>
L4:
  mov    -0x8(%rbp),%eax
  leave
  ret
```

assistant: ```c
int mp_sqrt_newton(int n, int radix, int in[], int inout[], int inout_rev[], int tmp[], int nfft, double tmp1fft[], double tmp2fft[], int ip[], double w[], int *n_tmp1fft) {
  void mp_round(int n, int radix, int m, int inout[]);
  void mp_add(int n, int radix, int in1[], int in2[], int out[]);
  void mp_sub(int n, int radix, int in1[], int in2[], int out[]);
  void mp_idiv_2(int n, int radix, int in[], int out[]);
  void mp_mulh(int n, int radix, int in1[], int in2[], int out[], int nfft, double in1fft[], double outfft[], int ip[], double w[]);
  void mp_squh(int n, int radix, int in[], int out[], int nfft, double inoutfft[], int ip[], double w[]);
  void mp_squh_use_in1fft(int n, int radix, double inoutfft[], int out[], int nfft, int ip[], double w[]);
  int n_h, nfft_h, shift, prc;
  nfft_h = nfft >> 1;
  shift = nfft_h + 1;
  if (nfft_h < 2) {
    nfft_h = 2;
  }
  n_h = n / 2 + 1;
  if (n_h < n - shift) {
    n_h = n - shift;
  }
  mp_round(n_h, radix, (nfft_h >> 1) + 1, inout_rev);
  if (*n_tmp1fft != nfft_h) {
    mp_squh(n_h, radix, inout_rev, tmp, nfft_h, tmp1fft, ip, w);
  } else {
    mp_squh_use_in1fft(n_h, radix, tmp1fft, tmp, nfft_h, ip, w);
  }
  mp_round(n, radix, shift, inout);
  mp_mulh(n_h, radix, inout, tmp, tmp, nfft, tmp1fft, tmp2fft, ip, w);
  mp_sub(n_h, radix, inout_rev, tmp, tmp);
  mp_add(n_h, radix, inout_rev, tmp, inout_rev);
  mp_squh_use_in1fft(n, radix, tmp1fft, tmp, nfft, ip, w);
  mp_sub(n, radix, in, tmp, tmp);
  prc = in[1] - tmp[1];
  if (in[2] > tmp[2]) {
    prc++;
  }
  if (tmp[0] == 0) {
    prc = nfft + 1;
  }
  mp_round(n_h, radix, shift, inout_rev);
  mp_mulh(n_h, radix, inout_rev, tmp, tmp, nfft, tmp1fft, tmp2fft, ip, w);
  *n_tmp1fft = nfft;
  mp_idiv_2(n_h, radix, tmp, tmp);
  if (tmp[0] != 0) {
    mp_add(n, radix, inout, tmp, inout);
  }
  return prc;
}
```