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

```asm
<draw_ellipse>:
  endbr64
  push   %r15
  mov    %ecx,%eax
  lea    (%rcx,%rsi,1),%r15d
  push   %r14
  mov    %edx,%r14d
  imul   %ecx,%eax
  push   %r13
  imul   %edx,%r14d
  mov    %ecx,%r13d
  push   %r12
  mov    %edx,%r12d
  mov    %r8d,%edx
  push   %rbp
  mov    %edi,%ebp
  push   %rbx
  mov    %r8d,%ebx
  sub    $0x38,%rsp
  mov    %esi,0x10(%rsp)
  mov    %r15d,%esi
  mov    %edi,0x18(%rsp)
  mov    %r14d,0x20(%rsp)
  mov    %eax,0x1c(%rsp)
  call   <draw_pixel@plt>
  mov    0x10(%rsp),%esi
  mov    %ebx,%edx
  mov    %ebp,%edi
  sub    %r13d,%esi
  call   <draw_pixel@plt>
  mov    %r14d,%eax
  mov    %r14d,%ecx
  lea    (%r14,%r14,1),%edx
  imul   %r13d,%eax
  sar    $0x2,%ecx
  xor    %r14d,%r14d
  mov    %r12d,0x2c(%rsp)
  mov    %r13d,%r11d
  mov    %r14d,%r12d
  mov    %edx,0x28(%rsp)
  lea    (%rax,%rax,1),%r10d
  sub    %eax,%ecx
  mov    0x1c(%rsp),%eax
  mov    %ecx,%r14d
  mov    %r10d,%r13d
  add    %eax,%eax
  mov    %eax,0x14(%rsp)
  mov    %ebp,%eax
  lea    0x1(%rbp),%ebp
  add    %eax,%eax
  mov    %eax,0x24(%rsp)
  jmp    L43
  nop
L79:
  sub    0x28(%rsp),%r13d
  add    0x14(%rsp),%r12d
  sub    $0x1,%r11d
  sub    %r13d,%r14d
  cmp    %r12d,%r13d
  jle    L64
  mov    0x10(%rsp),%eax
  lea    (%r11,%rax,1),%r15d
L71:
  mov    0x10(%rsp),%esi
  mov    %ebp,%edi
  mov    %ebx,%edx
  mov    %r11d,0xc(%rsp)
  sub    %r11d,%esi
  mov    %esi,0x8(%rsp)
  call   <draw_pixel@plt>
  mov    0x24(%rsp),%r9d
  mov    0x8(%rsp),%esi
  mov    %ebx,%edx
  sub    %ebp,%r9d
  mov    %r9d,%edi
  mov    %r9d,0x8(%rsp)
  call   <draw_pixel@plt>
  mov    %ebp,%edi
  mov    %ebx,%edx
  mov    %r15d,%esi
  call   <draw_pixel@plt>
  mov    0x8(%rsp),%edi
  mov    %ebx,%edx
  mov    %r15d,%esi
  add    $0x1,%ebp
  call   <draw_pixel@plt>
  mov    0xc(%rsp),%r11d
L43:
  mov    0x1c(%rsp),%eax
  add    %r12d,%eax
  add    %eax,%r14d
  jns    L79
  add    0x14(%rsp),%r12d
  cmp    %r12d,%r13d
  jg     L71
L64:
  mov    0x2c(%rsp),%r12d
  mov    0x18(%rsp),%ebp
  mov    %ebx,%edx
  xor    %r13d,%r13d
  mov    0x10(%rsp),%r15d
  lea    (%r12,%rbp,1),%r14d
  mov    %r15d,%esi
  mov    %r14d,%edi
  call   <draw_pixel@plt>
  mov    %ebp,%edi
  mov    %r15d,%esi
  mov    %ebx,%edx
  sub    %r12d,%edi
  call   <draw_pixel@plt>
  mov    0x1c(%rsp),%ecx
  mov    %r12d,0x8(%rsp)
  mov    %ecx,%eax
  sar    $0x2,%ecx
  imul   %r12d,%eax
  mov    %ecx,%ebp
  mov    %r14d,%r12d
  sub    %eax,%ebp
  lea    (%rax,%rax,1),%r9d
  mov    0x20(%rsp),%eax
  mov    %ebp,%r14d
  mov    %r13d,%ebp
  mov    %r9d,%r13d
  add    %eax,%eax
  mov    %eax,0x1c(%rsp)
  mov    %r15d,%eax
  lea    -0x1(%r15),%r15d
  add    %eax,%eax
  mov    %eax,0x24(%rsp)
  jmp    L7
  nopw   0x0(%rax,%rax,1)
L37:
  sub    0x14(%rsp),%r13d
  add    0x1c(%rsp),%ebp
  subl   $0x1,0x8(%rsp)
  sub    %r13d,%r14d
  mov    0x8(%rsp),%eax
  cmp    %ebp,%r13d
  jl     L49
  mov    0x18(%rsp),%ecx
  lea    (%rax,%rcx,1),%r12d
L61:
  mov    %r15d,%esi
  mov    %ebx,%edx
  mov    %r12d,%edi
  call   <draw_pixel@plt>
  mov    0x18(%rsp),%r10d
  mov    %r15d,%esi
  mov    %ebx,%edx
  sub    0x8(%rsp),%r10d
  mov    %r10d,%edi
  mov    %r10d,0x10(%rsp)
  call   <draw_pixel@plt>
  mov    0x24(%rsp),%esi
  mov    %ebx,%edx
  mov    %r12d,%edi
  sub    %r15d,%esi
  sub    $0x1,%r15d
  mov    %esi,0xc(%rsp)
  call   <draw_pixel@plt>
  mov    0xc(%rsp),%esi
  mov    0x10(%rsp),%edi
  mov    %ebx,%edx
  call   <draw_pixel@plt>
L7:
  mov    0x20(%rsp),%eax
  add    %ebp,%eax
  add    %eax,%r14d
  jns    L37
  add    0x1c(%rsp),%ebp
  cmp    %ebp,%r13d
  jge    L61
L49:
  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);
  }
}
```