Optimize poly division (#985)

* reproduce issue

* implement suggestion from issue

* refactor

---------

Co-authored-by: jotabulacios <[email protected]>

Author: nicole-graus

crates/math/src/polynomial/mod.rs

@@ -248,8 +248,12 @@ impl<F: IsField> Polynomial<FieldElement<F>> {
             Polynomial::new(&[FieldElement::zero()])
         } else {
             for i in 0..=factor.degree() {
-                for j in 0..=self.degree() {
-                    coefficients[i + j] += &factor.coefficients[i] * &self.coefficients[j];
+                if factor.coefficients[i] != FieldElement::zero() {
+                    for j in 0..=self.degree() {
+                        if self.coefficients[j] != FieldElement::zero() {
+                            coefficients[i + j] += &factor.coefficients[i] * &self.coefficients[j];
+                        }
+                    }
                 }
             }
             Polynomial::new(&coefficients)