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

```asm
<__emul>:
  endbr64
  push   %r14
  push   %r13
  mov    %rdx,%r13
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  sub    $0x40,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x38(%rsp)
  xor    %eax,%eax
  call   <__eisnan@plt>
  test   %rax,%rax
  jne    L109
  mov    %r12,%rdi
  call   <__eisnan@plt>
  test   %rax,%rax
  jne    L133
  mov    %rbp,%rdi
  call   <__eisinf@plt>
  test   %rax,%rax
  je     L131
  mov    %r12,%rdi
  call   <__eiiszero@plt>
  test   %rax,%rax
  jne    L172
L131:
  mov    %r12,%rdi
  call   <__eisinf@plt>
  test   %rax,%rax
  je     L110
  mov    %rbp,%rdi
  call   <__eiiszero@plt>
  test   %rax,%rax
  je     L110
L172:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L199
  add    $0x40,%rsp
  mov    %r13,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  jmp    1190 <__enan_NBITS@plt>
L133:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L199
  mov    %r13,%rsi
  mov    %r12,%rdi
L39:
  add    $0x40,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  jmp    1170 <__emov@plt>
  nopl   0x0(%rax,%rax,1)
L110:
  mov    %rbp,%rdi
  call   <__eisinf@plt>
  test   %rax,%rax
  jne    L188
  mov    %r12,%rdi
  call   <__eisinf@plt>
  test   %rax,%rax
  jne    L188
  mov    %rsp,%r14
  mov    %rbp,%rdi
  lea    0x20(%rsp),%rbp
  mov    %r14,%rsi
  call   <__emovi@plt>
  mov    %rbp,%rsi
  mov    %r12,%rdi
  call   <__emovi@plt>
  movzwl 0x2(%rsp),%eax
  movzwl 0x22(%rsp),%ebx
  test   %ax,%ax
  je     L97
  movzwl %ax,%r12d
  mov    %ebx,%eax
L5:
  test   %ax,%ax
  jne    L17
  movzwl 0x26(%rsp),%eax
  or     0x24(%rsp),%ax
  or     0x28(%rsp),%ax
  or     0x2a(%rsp),%ax
  or     0x2c(%rsp),%ax
  or     0x2e(%rsp),%ax
  jne    L67
L167:
  mov    %r13,%rdi
  call   <__eclear@plt>
L78:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L199
  add    $0x40,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nop
L188:
  mov    %rbp,%rdi
  call   <__eisneg@plt>
  mov    %r12,%rdi
  mov    %eax,%ebp
  call   <__eisneg@plt>
  cmp    %eax,%ebp
  setne  %al
  movzbl %al,%eax
  shl    $0xf,%eax
  mov    %ax,0xa(%r13)
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L199
  add    $0x40,%rsp
  mov    %r13,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  jmp    11e0 <__einfin@plt>
  nopw   0x0(%rax,%rax,1)
L109:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L199
  mov    %r13,%rsi
  mov    %rbp,%rdi
  jmp    L39
  nop
L97:
  movzwl 0x4(%rsp),%eax
  or     0x6(%rsp),%ax
  or     0x8(%rsp),%ax
  or     0xa(%rsp),%ax
  or     0xc(%rsp),%ax
  or     0xe(%rsp),%ax
  je     L167
  mov    %r14,%rdi
  call   <__enormlz@plt>
  neg    %rax
  mov    %rax,%r12
  movzwl 0x22(%rsp),%eax
  jmp    L5
L67:
  mov    %rbp,%rdi
  call   <__enormlz@plt>
  sub    %rax,%rbx
L17:
  mov    %rbp,%rsi
  mov    %r14,%rdi
  call   <__emulm@plt>
  xor    %edx,%edx
  mov    %rbp,%rdi
  mov    $0x50,%r9d
  mov    %eax,%esi
  lea    -0x3ffe(%r12,%rbx,1),%rcx
  mov    $0x40,%r8d
  call   <__emdnorm@plt>
  movzwl 0x20(%rsp),%eax
  mov    %r13,%rsi
  mov    %rbp,%rdi
  cmp    %ax,(%rsp)
  setne  %al
  movzbl %al,%eax
  neg    %eax
  mov    %ax,0x20(%rsp)
  call   <__emovo@plt>
  jmp    L78
L199:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void __emul(const short unsigned int *a, const short unsigned int *b, short unsigned int *c) {
  unsigned short ai[(6 + 3)], bi[(6 + 3)];
  int i, j;
  long lt, lta, ltb;
  if (__eisnan(a)) {
    __emov(a, c);
    return;
  }
  if (__eisnan(b)) {
    __emov(b, c);
    return;
  }
  if ((__eisinf(a) && __eiiszero(b)) || (__eisinf(b) && __eiiszero(a))) {
    ;
    __enan_NBITS(c);
    return;
  }
  if (__eisinf(a) || __eisinf(b)) {
    if (__eisneg(a) ^ __eisneg(b))
      *(c + (6 - 1)) = 0x8000;
    else
      *(c + (6 - 1)) = 0;
    __einfin(c);
    return;
  }
  __emovi(a, ai);
  __emovi(b, bi);
  lta = ai[1];
  ltb = bi[1];
  if (ai[1] == 0) {
    for (i = 1; i < (6 + 3) - 1; i++) {
      if (ai[i] != 0) {
        lta -= __enormlz(ai);
        goto mnzer1;
      }
    }
    __eclear(c);
    return;
  }
mnzer1:
  if (bi[1] == 0) {
    for (i = 1; i < (6 + 3) - 1; i++) {
      if (bi[i] != 0) {
        ltb -= __enormlz(bi);
        goto mnzer2;
      }
    }
    __eclear(c);
    return;
  }
mnzer2:
  j = __emulm(ai, bi);
  lt = lta + ltb - ((0x3fff) - 1);
  __emdnorm(bi, j, 0, lt, 64, (((6 + 3) - 4) * 16));
  if (ai[0] == bi[0])
    bi[0] = 0;
  else
    bi[0] = 0xffff;
  __emovo(bi, c);
}
```