diff --git a/Project.toml b/Project.toml index 4da61866..7cb0ab7b 100644 --- a/Project.toml +++ b/Project.toml @@ -1,7 +1,7 @@ name = "DynamicalSystemsBase" uuid = "6e36e845-645a-534a-86f2-f5d4aa5a06b4" repo = "https://github.com/JuliaDynamics/DynamicalSystemsBase.jl.git" -version = "3.2.1" +version = "3.2.2" [deps] ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210" diff --git a/src/derived_systems/poincare/hyperplane.jl b/src/derived_systems/poincare/hyperplane.jl index d8405a57..e1d21cf2 100644 --- a/src/derived_systems/poincare/hyperplane.jl +++ b/src/derived_systems/poincare/hyperplane.jl @@ -120,6 +120,7 @@ function poincaresos(A::StateSpaceSet, planecrossing::PlaneCrossing, j) return data end +using LinearAlgebra function interpolate_crossing(A, B, pc::PlaneCrossing{<:AbstractVector}) # https://en.wikipedia.org/wiki/Line%E2%80%93plane_intersection t = LinearAlgebra.dot(pc.n, (pc.p₀ .- A))/LinearAlgebra.dot((B .- A), pc.n) diff --git a/test/poincare.jl b/test/poincare.jl index b42e0c7a..7828bf74 100644 --- a/test/poincare.jl +++ b/test/poincare.jl @@ -96,4 +96,9 @@ end A = poincaresos(X, plane1) @test dimension(A) == 3 @test all(x -> abs(x) < 1e-12, A[:, 1]) + vec_plane = [1.0,0.0,0.0,0.0] + B = poincaresos(X, vec_plane) + @test dimension(B) == 3 + @test all(x -> abs(x) < 1e-12, B[:, 1]) + @test vec(A) == vec(B) end