diff --git a/Project.toml b/Project.toml
index 12b89dc..49e9429 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,6 +1,6 @@
 name = "FastTransforms"
 uuid = "057dd010-8810-581a-b7be-e3fc3b93f78c"
-version = "0.17.1"
+version = "0.17.2"
 
 
 [deps]
diff --git a/src/libfasttransforms.jl b/src/libfasttransforms.jl
index 320b9ee..fc3fdef 100644
--- a/src/libfasttransforms.jl
+++ b/src/libfasttransforms.jl
@@ -1293,3 +1293,5 @@ function modified_jacobi_matrix(R, XP)
     end
     return XQ
 end
+
+inv(P::FTPlan{T,2,DISKANALYSIS}) where T = plan_disk_synthesis(T, P.n, P.m)
\ No newline at end of file
diff --git a/test/libfasttransformstests.jl b/test/libfasttransformstests.jl
index 545da7d..7828ab0 100644
--- a/test/libfasttransformstests.jl
+++ b/test/libfasttransformstests.jl
@@ -171,6 +171,9 @@ FastTransforms.ft_set_num_threads(ceil(Int, Base.Sys.CPU_THREADS/2))
     ps = plan_disk_synthesis(A)
     pa = plan_disk_analysis(A)
     test_nd_plans(p, ps, pa, A)
+    @test inv(pa) isa typeof(ps)
+    @test inv(pa)*(p*A) ≈ ps*(p*A)
+    @test p\(pa*(inv(pa)*(p*A))) ≈ A
     A = diskones(Float64, n, 4n-3) + im*diskones(Float64, n, 4n-3)
     p = plan_disk2cxf(A, α, β)
     ps = plan_disk_synthesis(A)