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

```asm
<mp_inv_newton>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %edi,%ebp
  mov    %esi,%r12d
  mov    %rdx,0x8(%rsp)
  mov    %r8,%r13
  mov    %r9,%rbx
  mov    0x60(%rsp),%r14d
  mov    %r14d,%r15d
  sar    %r15d
  add    $0x1,%r15d
  mov    %edi,%eax
  shr    $0x1f,%eax
  add    %edi,%eax
  sar    %eax
  add    $0x1,%eax
  sub    %r15d,%edi
  cmp    %edi,%eax
  mov    %edi,%esi
  cmovge %eax,%esi
  mov    %esi,0x1c(%rsp)
  mov    %rcx,0x10(%rsp)
  mov    %r15d,%edx
  mov    %r12d,%esi
  mov    %ebp,%edi
  call   <mp_round@plt>
  push   0x70(%rsp)
  push   0x70(%rsp)
  mov    %r14d,%r9d
  mov    %r13,%r8
  mov    0x18(%rsp),%rcx
  mov    0x20(%rsp),%rdx
  mov    %r12d,%esi
  mov    %ebp,%edi
  call   <mp_mulh@plt>
  mov    %rbx,%rdx
  mov    %r12d,%esi
  mov    %ebp,%edi
  call   <mp_load_1@plt>
  mov    %rbx,%r8
  mov    %r13,%rcx
  mov    %rbx,%rdx
  mov    %r12d,%esi
  mov    %ebp,%edi
  call   <mp_sub@plt>
  push   0x80(%rsp)
  push   %r14
  mov    %r13,%r9
  mov    0x28(%rsp),%r8
  mov    %r15d,%ecx
  mov    0x88(%rsp),%rdx
  mov    %r12d,%esi
  mov    %ebp,%edi
  call   <mp_mulh_use_in1fft@plt>
  add    $0x20,%rsp
  mov    %rbx,%r8
  mov    %r13,%rcx
  mov    %rbx,%rdx
  mov    %r12d,%esi
  mov    0x1c(%rsp),%r15d
  mov    %r15d,%edi
  call   <mp_sub@plt>
  mov    0x4(%rbx),%eax
  neg    %eax
  lea    0x1(%r14),%r13d
  cmpl   $0x0,(%rbx)
  cmovne %eax,%r13d
  push   0x70(%rsp)
  push   %r14
  mov    %rbx,%r9
  mov    %rbx,%r8
  mov    $0x0,%ecx
  mov    0x78(%rsp),%rdx
  mov    %r12d,%esi
  mov    %r15d,%edi
  call   <mp_mulh_use_in1fft@plt>
  add    $0x10,%rsp
  cmpl   $0x0,(%rbx)
  jne    L15
L13:
  mov    %r13d,%eax
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L15:
  mov    0x10(%rsp),%rdx
  mov    %rdx,%r8
  mov    %rbx,%rcx
  mov    %r12d,%esi
  mov    %ebp,%edi
  call   <mp_add@plt>
  jmp    L13
```

assistant: ```c
int mp_inv_newton(int n, int radix, int in[], int inout[], int tmp1[], int tmp2[], int nfft, double tmp1fft[], double tmp2fft[]) {
  void mp_load_1(int n, int radix, int out[]);
  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_mulh(int n, int radix, int in1[], int in2[], int out[], int nfft, double in1fft[], double outfft[]);
  void mp_mulh_use_in1fft(int n, int radix, double in1fft[], int shift, int in2[], int out[], int nfft, double outfft[]);
  int n_h, shift, prc;
  shift = (nfft >> 1) + 1;
  n_h = n / 2 + 1;
  if (n_h < n - shift) {
    n_h = n - shift;
  }
  mp_round(n, radix, shift, inout);
  mp_mulh(n, radix, inout, in, tmp1, nfft, tmp1fft, tmp2fft);
  mp_load_1(n, radix, tmp2);
  mp_sub(n, radix, tmp2, tmp1, tmp2);
  mp_mulh_use_in1fft(n, radix, tmp1fft, shift, in, tmp1, nfft, tmp2fft);
  mp_sub(n_h, radix, tmp2, tmp1, tmp2);
  prc = -tmp2[1];
  if (tmp2[0] == 0) {
    prc = nfft + 1;
  }
  mp_mulh_use_in1fft(n_h, radix, tmp1fft, 0, tmp2, tmp2, nfft, tmp2fft);
  if (tmp2[0] != 0) {
    mp_add(n, radix, inout, tmp2, inout);
  }
  return prc;
}
```