user: Please show me the C implementation that matches the assembly code listed below:

```asm
<drawEllipse>:
  endbr64
  push   %r15
  mov    %edx,%r11d
  mov    $0x1,%r10d
  push   %r14
  push   %r13
  mov    %esi,%r13d
  push   %r12
  push   %rbp
  push   %rbx
  mov    %ecx,%ebx
  sub    $0x28,%rsp
  neg    %r11d
  cmovs  %edx,%r11d
  neg    %ebx
  mov    %edi,0x10(%rsp)
  cmovs  %ecx,%ebx
  xor    %r15d,%r15d
  xor    %r9d,%r9d
  imul   %edx,%edx
  xor    %r14d,%r14d
  mov    %r11d,%r12d
  imul   %ecx,%ecx
  mov    %ebx,0x1c(%rsp)
  mov    %r9d,%ebx
  lea    (%rdx,%rdx,1),%eax
  mov    %edx,0xc(%rsp)
  lea    (%rcx,%rcx,1),%ebp
  mov    %eax,0x4(%rsp)
  lea    (%r11,%r11,1),%eax
  sub    %eax,%r10d
  mov    %ebp,%eax
  mov    %ecx,0x14(%rsp)
  imul   %r11d,%eax
  mov    %ebp,0x18(%rsp)
  mov    %r15d,%ebp
  imul   %ecx,%r10d
  mov    %eax,0x8(%rsp)
  mov    %r10d,%r15d
  nopw   0x0(%rax,%rax,1)
L3:
  mov    0xc(%rsp),%eax
  mov    0x10(%rsp),%edi
  mov    %ebp,%ecx
  mov    %r12d,%edx
  mov    %r13d,%esi
  add    $0x1,%ebp
  add    %ebx,%eax
  mov    %eax,(%rsp)
  call   <Draw_4_Ellipse_Points@plt>
  add    (%rsp),%r14d
  add    0x4(%rsp),%ebx
  lea    (%r15,%r14,2),%eax
  test   %eax,%eax
  jle    L38
  mov    0x18(%rsp),%edi
  sub    %edi,0x8(%rsp)
  sub    %r14d,%eax
  sub    $0x1,%r12d
  mov    %eax,%r14d
  add    %edi,%r15d
L38:
  cmp    %ebx,0x8(%rsp)
  jge    L3
  mov    0x1c(%rsp),%ebx
  mov    $0x1,%r10d
  xor    %r14d,%r14d
  xor    %r15d,%r15d
  mov    0x18(%rsp),%ebp
  mov    0x10(%rsp),%r12d
  mov    %r13d,0x10(%rsp)
  lea    (%rbx,%rbx,1),%eax
  sub    %eax,%r10d
  mov    0x4(%rsp),%eax
  imul   0xc(%rsp),%r10d
  mov    %ebp,0x8(%rsp)
  mov    %r12d,0xc(%rsp)
  mov    %ebx,%r12d
  imul   %ebx,%eax
  mov    %r14d,%ebx
  mov    %r10d,%r13d
  mov    %eax,(%rsp)
  xchg   %ax,%ax
L6:
  mov    0x14(%rsp),%eax
  mov    0x10(%rsp),%esi
  mov    %r14d,%edx
  mov    %r12d,%ecx
  mov    0xc(%rsp),%edi
  add    $0x1,%r14d
  lea    (%r15,%rax,1),%ebp
  call   <Draw_4_Ellipse_Points@plt>
  add    %ebp,%ebx
  add    0x8(%rsp),%r15d
  lea    0x0(%r13,%rbx,2),%eax
  test   %eax,%eax
  jle    L25
  mov    0x4(%rsp),%ecx
  sub    %ecx,(%rsp)
  sub    %ebx,%eax
  sub    $0x1,%r12d
  mov    %eax,%ebx
  add    %ecx,%r13d
L25:
  cmp    (%rsp),%r15d
  jle    L6
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void drawEllipse(int CX, int CY, int XRadius, int YRadius) {
  int X, Y, XChange, YChange, EllipseError, TwoASquare, TwoBSquare, StoppingX, StoppingY;
  if (XRadius < 0)
    XRadius = -XRadius;
  if (YRadius < 0)
    YRadius = -YRadius;
  TwoASquare = 2 * XRadius * XRadius;
  TwoBSquare = 2 * YRadius * YRadius;
  X = XRadius;
  Y = 0;
  XChange = YRadius * YRadius * (1 - 2 * XRadius);
  YChange = XRadius * XRadius;
  EllipseError = 0;
  StoppingX = TwoBSquare * XRadius;
  StoppingY = 0;
  while (StoppingX >= StoppingY) {
    Draw_4_Ellipse_Points(CX, CY, X, Y);
    Y++;
    StoppingY = StoppingY + TwoASquare;
    EllipseError = EllipseError + YChange;
    YChange = YChange + TwoASquare;
    if ((2 * EllipseError + XChange) > 0) {
      X--;
      StoppingX = StoppingX - TwoBSquare;
      EllipseError = EllipseError + XChange;
      XChange = XChange + TwoBSquare;
    }
  }
  X = 0;
  Y = YRadius;
  XChange = YRadius * YRadius;
  YChange = XRadius * XRadius * (1 - 2 * YRadius);
  EllipseError = 0;
  StoppingX = 0;
  StoppingY = TwoASquare * YRadius;
  while (StoppingX <= StoppingY) {
    Draw_4_Ellipse_Points(CX, CY, X, Y);
    X++;
    StoppingX = StoppingX + TwoBSquare;
    EllipseError = EllipseError + XChange;
    XChange = XChange + TwoBSquare;
    if ((2 * EllipseError + YChange) > 0) {
      Y--;
      StoppingY = StoppingY - TwoASquare;
      EllipseError = EllipseError + YChange;
      YChange = YChange + TwoASquare;
    }
  }
  return;
}
```