Polynomial long division: Difference between revisions

Content added Content deleted
(Updated D entry)
(Updated D entry)
Line 705: Line 705:
throw new Error("ZeroDivisionError");
throw new Error("ZeroDivisionError");
if (dN >= dD) {
if (dN >= dD) {
//q = [0.0].replicate(dN);
q = [0.0].replicate(dN);
q = std.array.replicate([0.0], dN);
while (dN >= dD) {
while (dN >= dD) {
auto d = std.array.replicate([0.0], dN - dD) ~ D;
auto d = [0.0].replicate(dN - dD) ~ D;
immutable mult = q[dN - dD] = N[$ - 1] / d[$ - 1];
immutable mult = q[dN - dD] = N[$ - 1] / d[$ - 1];
d[] *= mult;
d[] *= mult;