fixed bug in euclidean_div

This commit is contained in:
Andreas Fabri 2003-04-04 12:42:50 +00:00
parent 74c73c4815
commit e009ffb7b7
1 changed files with 1 additions and 1 deletions

View File

@ -1930,7 +1930,7 @@ void Polynomial<NT>::euclidean_div(
qd = rd - gd + 1;
q = Polynomial<NT>(std::size_t(qd));
}
while ( rd >= gd ) {
while ( rd >= gd && !(r.is_zero())) {
NT S = r[rd] / g[gd];
qd = rd-gd;
q.coeff(qd) += S;