user: Translate the below assembly code into C programming language:

```asm
<draw_ellipse>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %edi,%ebp
  mov    %edi,0x24(%rsp)
  mov    %esi,%r12d
  mov    %esi,0x10(%rsp)
  mov    %edx,0xc(%rsp)
  mov    %ecx,%r15d
  mov    %ecx,0x8(%rsp)
  mov    %r8d,%r13d
  mov    %edx,%ebx
  imul   %edx,%ebx
  mov    %ebx,0x1c(%rsp)
  mov    %ecx,%r14d
  imul   %ecx,%r14d
  mov    %r14d,0x20(%rsp)
  mov    %ecx,%eax
  add    %esi,%eax
  mov    %eax,%esi
  mov    %r8d,%edx
  call   <draw_pixel@plt>
  mov    %r12d,%esi
  sub    %r15d,%esi
  mov    %r13d,%edx
  mov    %ebp,%edi
  call   <draw_pixel@plt>
  mov    %ebx,%eax
  imul   %r15d,%eax
  lea    (%rax,%rax,1),%r15d
  mov    %ebx,%ecx
  lea    0x3(%rbx),%ebx
  test   %ecx,%ecx
  cmovns %ecx,%ebx
  sar    $0x2,%ebx
  sub    %eax,%ebx
  lea    (%rcx,%rcx,1),%eax
  mov    %eax,0x2c(%rsp)
  lea    (%r14,%r14,1),%eax
  mov    %eax,0x18(%rsp)
  mov    %ebp,%eax
  lea    0x1(%rbp),%ebp
  add    %eax,%eax
  mov    %eax,0x28(%rsp)
  mov    $0x0,%r12d
  mov    %r15d,%r14d
  jmp    L14
L19:
  add    0x18(%rsp),%r12d
  cmp    %r12d,%r14d
  jle    L64
  mov    0x10(%rsp),%r15d
  sub    0x8(%rsp),%r15d
  mov    %r13d,%edx
  mov    %r15d,%esi
  mov    %ebp,%edi
  call   <draw_pixel@plt>
  mov    0x28(%rsp),%eax
  sub    %ebp,%eax
  mov    %r13d,%edx
  mov    %r15d,%esi
  mov    %eax,0x14(%rsp)
  mov    %eax,%edi
  call   <draw_pixel@plt>
  mov    0x8(%rsp),%ecx
  add    0x10(%rsp),%ecx
  mov    %ecx,%r15d
  mov    %r13d,%edx
  mov    %ecx,%esi
  mov    %ebp,%edi
  call   <draw_pixel@plt>
  mov    %r13d,%edx
  mov    %r15d,%esi
  mov    0x14(%rsp),%edi
  call   <draw_pixel@plt>
  add    $0x1,%ebp
L14:
  mov    0x20(%rsp),%eax
  add    %r12d,%eax
  add    %eax,%ebx
  js     L19
  sub    0x2c(%rsp),%r14d
  sub    %r14d,%ebx
  subl   $0x1,0x8(%rsp)
  jmp    L19
L64:
  mov    0xc(%rsp),%ebx
  mov    %ebx,%edi
  mov    0x24(%rsp),%r15d
  add    %r15d,%edi
  mov    %r13d,%edx
  mov    0x10(%rsp),%r14d
  mov    %r14d,%esi
  call   <draw_pixel@plt>
  mov    %r15d,%edi
  sub    %ebx,%edi
  mov    %r13d,%edx
  mov    %r14d,%esi
  call   <draw_pixel@plt>
  mov    0x20(%rsp),%ecx
  imul   %ecx,%ebx
  mov    %ebx,%eax
  lea    (%rbx,%rbx,1),%ebx
  mov    %ebx,0x8(%rsp)
  lea    0x3(%rcx),%ebx
  test   %ecx,%ecx
  cmovns %ecx,%ebx
  sar    $0x2,%ebx
  sub    %eax,%ebx
  mov    0x1c(%rsp),%eax
  add    %eax,%eax
  mov    %eax,0x14(%rsp)
  lea    -0x1(%r14),%ebp
  lea    (%r14,%r14,1),%eax
  mov    %eax,0x20(%rsp)
  mov    $0x0,%r12d
  jmp    L31
L70:
  mov    0x18(%rsp),%ecx
  sub    %ecx,0x8(%rsp)
  mov    0x8(%rsp),%eax
  sub    %eax,%ebx
  subl   $0x1,0xc(%rsp)
L59:
  add    0x14(%rsp),%r12d
  cmp    %r12d,0x8(%rsp)
  jl     L24
  mov    0xc(%rsp),%r15d
  mov    0x24(%rsp),%r14d
  add    %r14d,%r15d
  mov    %r13d,%edx
  mov    %ebp,%esi
  mov    %r15d,%edi
  call   <draw_pixel@plt>
  sub    0xc(%rsp),%r14d
  mov    %r13d,%edx
  mov    %ebp,%esi
  mov    %r14d,0x10(%rsp)
  mov    %r14d,%edi
  call   <draw_pixel@plt>
  mov    0x20(%rsp),%eax
  sub    %ebp,%eax
  mov    %eax,%r14d
  mov    %r13d,%edx
  mov    %eax,%esi
  mov    %r15d,%edi
  call   <draw_pixel@plt>
  mov    %r13d,%edx
  mov    %r14d,%esi
  mov    0x10(%rsp),%edi
  call   <draw_pixel@plt>
  sub    $0x1,%ebp
L31:
  mov    0x1c(%rsp),%eax
  add    %r12d,%eax
  add    %eax,%ebx
  jns    L70
  jmp    L59
L24:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void draw_ellipse(int x, int y, int a, int b, int color) {
  int wx, wy;
  int thresh;
  int asq = a * a;
  int bsq = b * b;
  int xa, ya;
  draw_pixel(x, y + b, color);
  draw_pixel(x, y - b, color);
  wx = 0;
  wy = b;
  xa = 0;
  ya = asq * 2 * b;
  thresh = asq / 4 - asq * b;
  for (;;) {
    thresh += xa + bsq;
    if (thresh >= 0) {
      ya -= asq * 2;
      thresh -= ya;
      wy--;
    }
    xa += bsq * 2;
    wx++;
    if (xa >= ya)
      break;
    draw_pixel(x + wx, y - wy, color);
    draw_pixel(x - wx, y - wy, color);
    draw_pixel(x + wx, y + wy, color);
    draw_pixel(x - wx, y + wy, color);
  }
  draw_pixel(x + a, y, color);
  draw_pixel(x - a, y, color);
  wx = a;
  wy = 0;
  xa = bsq * 2 * a;
  ya = 0;
  thresh = bsq / 4 - bsq * a;
  for (;;) {
    thresh += ya + asq;
    if (thresh >= 0) {
      xa -= bsq * 2;
      thresh = thresh - xa;
      wx--;
    }
    ya += asq * 2;
    wy++;
    if (ya > xa)
      break;
    draw_pixel(x + wx, y - wy, color);
    draw_pixel(x - wx, y - wy, color);
    draw_pixel(x + wx, y + wy, color);
    draw_pixel(x - wx, y + wy, color);
  }
}
```