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

```asm
<__emul>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %rdx,%r13
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  mov    %rdi,%rbp
  sub    $0x40,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x38(%rsp)
  xor    %eax,%eax
  call   <__eisnan@plt>
  test   %rax,%rax
  jne    L10
  mov    %r12,%rdi
  call   <__eisnan@plt>
  test   %rax,%rax
  jne    L125
  mov    %rbp,%rdi
  call   <__eisinf@plt>
  test   %rax,%rax
  je     L128
  mov    %r12,%rdi
  call   <__eiiszero@plt>
  test   %rax,%rax
  jne    L181
L128:
  mov    %r12,%rdi
  call   <__eisinf@plt>
  test   %rax,%rax
  je     L93
  mov    %rbp,%rdi
  call   <__eiiszero@plt>
  test   %rax,%rax
  je     L93
L181:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L53
  add    $0x40,%rsp
  mov    %r13,%rdi
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1190 <__enan_NBITS@plt>
  nopw   0x0(%rax,%rax,1)
L125:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L53
  mov    %r13,%rsi
  mov    %r12,%rdi
L158:
  add    $0x40,%rsp
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1170 <__emov@plt>
  nopl   0x0(%rax)
L93:
  mov    %rbp,%rdi
  call   <__eisinf@plt>
  test   %rax,%rax
  jne    L172
  mov    %r12,%rdi
  call   <__eisinf@plt>
  test   %rax,%rax
  jne    L172
  mov    %rsp,%r15
  mov    %rbp,%rdi
  lea    0x20(%rsp),%r14
  mov    %r15,%rsi
  call   <__emovi@plt>
  mov    %r12,%rdi
  mov    %r14,%rsi
  call   <__emovi@plt>
  movzwl 0x2(%rsp),%eax
  movzwl 0x22(%rsp),%r12d
  test   %ax,%ax
  je     L147
  movzwl %ax,%ebp
  mov    %r12d,%ecx
L40:
  lea    0x24(%rsp),%rax
  lea    0x30(%rsp),%rdx
  test   %cx,%cx
  jne    L135
L169:
  cmp    %rax,%rdx
  je     L68
  add    $0x2,%rax
  cmpw   $0x0,-0x2(%rax)
  je     L169
  mov    %r14,%rdi
  call   <__enormlz@plt>
  sub    %rax,%r12
L135:
  mov    %r14,%rsi
  mov    %r15,%rdi
  call   <__emulm@plt>
  xor    %edx,%edx
  mov    %r14,%rdi
  mov    $0x50,%r9d
  mov    %eax,%esi
  lea    -0x3ffe(%rbp,%r12,1),%rcx
  mov    $0x40,%r8d
  call   <__emdnorm@plt>
  movzwl 0x20(%rsp),%eax
  mov    %r13,%rsi
  mov    %r14,%rdi
  cmp    %ax,(%rsp)
  setne  %al
  movzbl %al,%eax
  neg    %eax
  mov    %ax,0x20(%rsp)
  call   <__emovo@plt>
L48:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L53
  add    $0x40,%rsp
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  xchg   %ax,%ax
L172:
  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    L53
  add    $0x40,%rsp
  mov    %r13,%rdi
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    11e0 <__einfin@plt>
  nopl   0x0(%rax,%rax,1)
L10:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L53
  mov    %r13,%rsi
  mov    %rbp,%rdi
  jmp    L158
  nopl   0x0(%rax,%rax,1)
L147:
  lea    0x4(%rsp),%rax
  lea    0x10(%rsp),%rdx
L67:
  cmp    %rdx,%rax
  je     L68
  add    $0x2,%rax
  cmpw   $0x0,-0x2(%rax)
  je     L67
  mov    %r15,%rdi
  call   <__enormlz@plt>
  movzwl 0x22(%rsp),%ecx
  neg    %rax
  mov    %rax,%rbp
  jmp    L40
  nopw   0x0(%rax,%rax,1)
L68:
  mov    %r13,%rdi
  call   <__eclear@plt>
  jmp    L48
L53:
  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);
}
```