Use ca^2 instead of ca*ca; important for Intervals

JuliaDiff · Apr 18, 2018 · 4d5ac34 · 4d5ac34
1 parent 2a646b8
commit 4d5ac34
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/power.jl b/src/power.jl
@@ -348,7 +348,7 @@ Return `c = a*a` with no allocation; all parameters are `HomogeneousPolynomial`.
         iszero(ca) && continue
         inda = idxTb[na]
         pos = posTb[2*inda]
-        c[pos] += ca * ca
+        c[pos] += ca^2
         @inbounds for nb = na+1:num_coeffs_a
             cb = a[nb]
             iszero(cb) && continue