[EricVanWyk]

You only half integrated my changes to divide.

q and r are coming in as references. You do not need to declare them again.

Also,
divide(base, q[j-1], &q[j], &r[j]);
should read as
divide(base, q[j-1], q[j], r[j]);

I'm going to be offline for a while. Good luck!