Use muladd in complex evalpoly (#35562)

JuliaLang · Apr 23, 2020 · eaa2c58 · eaa2c58
1 parent 407caee
commit eaa2c58
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/base/math.jl b/base/math.jl
@@ -153,7 +153,7 @@ function evalpoly(z::Complex, p::Tuple)
  ai = Symbol("a", i)
  push!(as, :($ai = $a))
  a = :(muladd(r, $ai, $b))
- b = :(p[$i] - s * $ai)
+ b = :(muladd(-s, $ai, p[$i]))
  end
  ai = :a0
  push!(as, :($ai = $a))
@@ -186,7 +186,7 @@ function _evalpoly(z::Complex, p)
  for i in N-2:-1:1
  ai = a
  a = muladd(r, ai, b)
- b = p[i] - s * ai
+ b = muladd(-s, ai, p[i])
  end
  ai = a
  muladd(ai, z, b)