From cbd818849b494d43860568c94bd8accf2b8dc9bc Mon Sep 17 00:00:00 2001 From: Uziel Linares Date: Sun, 28 Jun 2020 17:37:06 -0500 Subject: [PATCH 1/2] modify TaylorModel1s product --- src/arithmetic.jl | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/src/arithmetic.jl b/src/arithmetic.jl index 060b9816..e80001f8 100644 --- a/src/arithmetic.jl +++ b/src/arithmetic.jl @@ -126,9 +126,12 @@ end function remainder_product(a, b, aux, Δnegl) Δa = a.pol(aux) Δb = b.pol(aux) - Δ = Δnegl + Δb * a.rem + Δa * b.rem + a.rem * b.rem + Δ₁ = Δa * b.rem + a.rem * (Δb + b.rem) + Δ₂ = Δb * a.rem + b.rem * (Δa + a.rem) + Δ = Δnegl + (Δ₁ ∩ Δ₂) return Δ end + function remainder_product(a::TaylorModel1{TaylorModelN{N,T,S},S}, b::TaylorModel1{TaylorModelN{N,T,S},S}, aux, Δnegl) where {N,T,S} Δa = a.pol(aux) From e4ae2244cd5f661a41e088190a03a33647786c8d Mon Sep 17 00:00:00 2001 From: Uziel Linares Date: Mon, 6 Jul 2020 15:38:45 -0500 Subject: [PATCH 2/2] fix remainder_product function --- src/arithmetic.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/arithmetic.jl b/src/arithmetic.jl index e80001f8..aa1d3339 100644 --- a/src/arithmetic.jl +++ b/src/arithmetic.jl @@ -128,7 +128,7 @@ function remainder_product(a, b, aux, Δnegl) Δb = b.pol(aux) Δ₁ = Δa * b.rem + a.rem * (Δb + b.rem) Δ₂ = Δb * a.rem + b.rem * (Δa + a.rem) - Δ = Δnegl + (Δ₁ ∩ Δ₂) + Δ = (Δnegl + Δ₁) ∩ (Δnegl + Δ₂) return Δ end