Fix bug in ultra2ultra with integer difference

JuliaApproximation · Nov 16, 2023 · d4fd97a · d4fd97a
1 parent a1024f2
commit d4fd97a
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,6 +1,6 @@
 name = "FastTransforms"
 uuid = "057dd010-8810-581a-b7be-e3fc3b93f78c"
-version = "0.15.14"
+version = "0.15.15"
 
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"

diff --git a/src/toeplitzhankel.jl b/src/toeplitzhankel.jl
@@ -299,7 +299,7 @@ function plan_th_ultra2ultra!(::Type{S}, mn, λ₁, λ₂, dims) where {S}
 
     if isapproxinteger(λ₂ - λ₁)
         # TODO: don't make extra plan
-        plans = typeof(_good_plan_th_ultra2ultra!(S, mn, λ₁+0.1, λ₂, dims))[]
+        plans = typeof(_good_plan_th_ultra2ultra!(S, mn, λ₂+0.1, λ₂, dims))[]
     else
         plans = [_good_plan_th_ultra2ultra!(S, mn, λ₁, c, dims)]
     end

diff --git a/test/toeplitzhankeltests.jl b/test/toeplitzhankeltests.jl
@@ -15,6 +15,7 @@ Random.seed!(0)
         @test th_ultra2ultra(x,1, 2) ≈ lib_ultra2ultra(x, 1, 2)
         @test th_ultra2ultra(x,0.1, 2.2) ≈ lib_ultra2ultra(x, 0.1, 2.2)
         @test th_ultra2ultra(x, 2.2, 0.1) ≈ lib_ultra2ultra(x, 2.2, 0.1)
+        @test th_ultra2ultra(x, 1, 3) ≈ lib_ultra2ultra(x, 1, 3)
         @test @inferred(th_jac2jac(x,0.1, 0.2,0.1,0.4)) ≈ lib_jac2jac(x, 0.1, 0.2,0.1,0.4)
         @test th_jac2jac(x,0.1, 0.2,0.3,0.2) ≈ lib_jac2jac(x, 0.1, 0.2,0.3,0.2)
         @test th_jac2jac(x,0.1, 0.2,0.3,0.4) ≈ lib_jac2jac(x, 0.1, 0.2,0.3,0.4)