Implementation
int aaellipseRgba(Pointer<SdlRenderer> renderer, int x, int y, int rx, int ry,
int r, int g, int b, int a) {
int result;
int i;
int a2, b2, ds, dt, dxt, t, s, d;
int xp, yp, xs, ys, dyt, od, xx, yy, xc2, yc2;
double cp;
double sab;
int weight, iweight;
/*
* Sanity check radii
*/
if ((rx < 0) || (ry < 0)) {
return (-1);
}
/*
* Special cases for rx=0 and/or ry=0: draw a hline/vline/pixel
*/
if (rx == 0) {
if (ry == 0) {
return (pixelRgba(renderer, x, y, r, g, b, a));
} else {
return (vlineRgba(renderer, x, y - ry, y + ry, r, g, b, a));
}
} else {
if (ry == 0) {
return (hlineRgba(renderer, x - rx, x + rx, y, r, g, b, a));
}
}
/* Variable setup */
a2 = rx * rx;
b2 = ry * ry;
ds = 2 * a2;
dt = 2 * b2;
xc2 = 2 * x;
yc2 = 2 * y;
sab = sqrt(a2 + b2);
od = (sab * 0.01).toInt() + 1; /* introduce some overdraw */
dxt = (a2 ~/ sab) + od;
t = 0;
s = -2 * a2 * ry;
d = 0;
xp = x;
yp = y - ry;
/* Draw */
result = 0;
result |= sdlSetRenderDrawBlendMode(
renderer, (a == 255) ? SDL_BLENDMODE_NONE : SDL_BLENDMODE_BLEND);
/* "End points" */
result |= pixelRgba(renderer, xp, yp, r, g, b, a);
result |= pixelRgba(renderer, xc2 - xp, yp, r, g, b, a);
result |= pixelRgba(renderer, xp, yc2 - yp, r, g, b, a);
result |= pixelRgba(renderer, xc2 - xp, yc2 - yp, r, g, b, a);
for (i = 1; i <= dxt; i++) {
xp--;
d += t - b2;
if (d >= 0) {
ys = yp - 1;
} else if ((d - s - a2) > 0) {
if ((2 * d - s - a2) >= 0) {
ys = yp + 1;
} else {
ys = yp;
yp++;
d -= s + a2;
s += ds;
}
} else {
yp++;
ys = yp + 1;
d -= s + a2;
s += ds;
}
t -= dt;
/* Calculate alpha */
if (s != 0) {
cp = d.abs() / s.abs();
if (cp > 1.0) {
cp = 1.0;
}
} else {
cp = 1.0;
}
/* Calculate weights */
weight = (cp * 255).toInt();
iweight = 255 - weight;
/* Upper half */
xx = xc2 - xp;
result |= pixelRgbaWeight(renderer, xp, yp, r, g, b, a, iweight);
result |= pixelRgbaWeight(renderer, xx, yp, r, g, b, a, iweight);
result |= pixelRgbaWeight(renderer, xp, ys, r, g, b, a, weight);
result |= pixelRgbaWeight(renderer, xx, ys, r, g, b, a, weight);
/* Lower half */
yy = yc2 - yp;
result |= pixelRgbaWeight(renderer, xp, yy, r, g, b, a, iweight);
result |= pixelRgbaWeight(renderer, xx, yy, r, g, b, a, iweight);
yy = yc2 - ys;
result |= pixelRgbaWeight(renderer, xp, yy, r, g, b, a, weight);
result |= pixelRgbaWeight(renderer, xx, yy, r, g, b, a, weight);
}
/* Replaces original approximation code dyt = abs(yp - yc); */
dyt = b2 ~/ sab + od;
for (i = 1; i <= dyt; i++) {
yp++;
d -= s + a2;
if (d <= 0) {
xs = xp + 1;
} else if ((d + t - b2) < 0) {
if ((2 * d + t - b2) <= 0) {
xs = xp - 1;
} else {
xs = xp;
xp--;
d += t - b2;
t -= dt;
}
} else {
xp--;
xs = xp - 1;
d += t - b2;
t -= dt;
}
s += ds;
/* Calculate alpha */
if (t != 0) {
cp = d.abs() / t.abs();
if (cp > 1.0) {
cp = 1.0;
}
} else {
cp = 1.0;
}
/* Calculate weight */
weight = (cp * 255).toInt();
iweight = 255 - weight;
/* Left half */
xx = xc2 - xp;
yy = yc2 - yp;
result |= pixelRgbaWeight(renderer, xp, yp, r, g, b, a, iweight);
result |= pixelRgbaWeight(renderer, xx, yp, r, g, b, a, iweight);
result |= pixelRgbaWeight(renderer, xp, yy, r, g, b, a, iweight);
result |= pixelRgbaWeight(renderer, xx, yy, r, g, b, a, iweight);
/* Right half */
xx = xc2 - xs;
result |= pixelRgbaWeight(renderer, xs, yp, r, g, b, a, weight);
result |= pixelRgbaWeight(renderer, xx, yp, r, g, b, a, weight);
result |= pixelRgbaWeight(renderer, xs, yy, r, g, b, a, weight);
result |= pixelRgbaWeight(renderer, xx, yy, r, g, b, a, weight);
}
return (result);
}