merge common code for linear cpost

JuliaReach · Jan 10, 2025 · 8f41b33 · 8f41b33
1 parent 6cdf144
commit 8f41b33
Show file tree

Hide file tree

Showing 9 changed files with 30 additions and 140 deletions.
diff --git a/src/Algorithms/A20/post.jl b/src/Algorithms/A20/post.jl
@@ -1,23 +1,4 @@
-# continuous post
-function post(alg::A20, ivp::IVP{<:AbstractContinuousSystem}, tspan;
-              Δt0::TimeInterval=zeroI, kwargs...)
-    δ = alg.δ
-
-    NSTEPS = get(kwargs, :NSTEPS, compute_nsteps(δ, tspan))
-
-    # normalize system to canonical form
-    ivp_norm = _normalize(ivp)
-
-    # homogenize the initial-value problem
-    if get(kwargs, :homogenize, false)
-        ivp_norm = homogenize(ivp_norm)
-    end
-
-    # discretize system
-    ivp_discr = discretize(ivp_norm, δ, alg.approx_model)
-
-    return post(alg, ivp_discr, NSTEPS; Δt0=Δt0, kwargs...)
-end
+# default continuous post
 
 # discrete post
 function post(alg::A20{N}, ivp::IVP{<:AbstractDiscreteSystem}, NSTEPS=nothing;

diff --git a/src/Algorithms/BFFPSV18/post.jl b/src/Algorithms/BFFPSV18/post.jl
@@ -1,23 +1,4 @@
-# continuous post
-function post(alg::BFFPSV18, ivp::IVP{<:AbstractContinuousSystem}, tspan;
-              Δt0::TimeInterval=zeroI, kwargs...)
-    δ = alg.δ
-
-    NSTEPS = get(kwargs, :NSTEPS, compute_nsteps(δ, tspan))
-
-    # normalize system to canonical form
-    ivp_norm = _normalize(ivp)
-
-    # homogenize the initial-value problem
-    if get(kwargs, :homogenize, false)
-        ivp_norm = homogenize(ivp_norm)
-    end
-
-    # discretize system
-    ivp_discr = discretize(ivp_norm, δ, alg.approx_model)
-
-    return post(alg, ivp_discr, NSTEPS; Δt0=Δt0, kwargs...)
-end
+# default continuous post
 
 # discrete post
 function post(alg::BFFPSV18{N,ST}, ivp::IVP{<:AbstractDiscreteSystem}, NSTEPS=nothing;

diff --git a/src/Algorithms/BOX/post.jl b/src/Algorithms/BOX/post.jl
@@ -1,23 +1,4 @@
-# continuous post
-function post(alg::BOX, ivp::IVP{<:AbstractContinuousSystem}, tspan;
-              Δt0::TimeInterval=zeroI, kwargs...)
-    δ = alg.δ
-
-    NSTEPS = get(kwargs, :NSTEPS, compute_nsteps(δ, tspan))
-
-    # normalize system to canonical form
-    ivp_norm = _normalize(ivp)
-
-    # homogenize the initial-value problem
-    if get(kwargs, :homogenize, false)
-        ivp_norm = homogenize(ivp_norm)
-    end
-
-    # discretize system
-    ivp_discr = discretize(ivp_norm, δ, alg.approx_model)
-
-    return post(alg, ivp_discr, NSTEPS; Δt0=Δt0, kwargs...)
-end
+# default continuous post
 
 # discrete post
 function post(alg::BOX{N}, ivp::IVP{<:AbstractDiscreteSystem}, NSTEPS=nothing;

diff --git a/src/Algorithms/GLGM06/post.jl b/src/Algorithms/GLGM06/post.jl
@@ -1,23 +1,4 @@
-# continuous post
-function post(alg::GLGM06, ivp::IVP{<:AbstractContinuousSystem}, tspan;
-              Δt0::TimeInterval=zeroI, kwargs...)
-    δ = alg.δ
-
-    NSTEPS = get(kwargs, :NSTEPS, compute_nsteps(δ, tspan))
-
-    # normalize system to canonical form
-    ivp_norm = _normalize(ivp)
-
-    # homogenize the initial-value problem
-    if get(kwargs, :homogenize, false)
-        ivp_norm = homogenize(ivp_norm)
-    end
-
-    # discretize system
-    ivp_discr = discretize(ivp_norm, δ, alg.approx_model)
-
-    return post(alg, ivp_discr, NSTEPS; Δt0=Δt0, kwargs...)
-end
+# default continuous post
 
 # discrete post
 function post(alg::GLGM06{N}, ivp::IVP{<:AbstractDiscreteSystem}, NSTEPS=nothing;

diff --git a/src/Algorithms/INT/post.jl b/src/Algorithms/INT/post.jl
@@ -1,23 +1,4 @@
-# continuous post
-function post(alg::INT, ivp::IVP{<:AbstractContinuousSystem}, tspan;
-              Δt0::TimeInterval=zeroI, kwargs...)
-    δ = alg.δ
-
-    NSTEPS = get(kwargs, :NSTEPS, compute_nsteps(δ, tspan))
-
-    # normalize system to canonical form
-    ivp_norm = _normalize(ivp)
-
-    # homogenize the initial-value problem
-    if get(kwargs, :homogenize, false)
-        ivp_norm = homogenize(ivp_norm)
-    end
-
-    # discretize system
-    ivp_discr = discretize(ivp_norm, δ, alg.approx_model)
-
-    return post(alg, ivp_discr, NSTEPS; Δt0=Δt0, kwargs...)
-end
+# default continuous post
 
 # discrete post
 function post(alg::INT{N}, ivp::IVP{<:AbstractDiscreteSystem}, NSTEPS=nothing;

diff --git a/src/Algorithms/LGG09/post.jl b/src/Algorithms/LGG09/post.jl
@@ -1,23 +1,4 @@
-# continuous post
-function post(alg::LGG09, ivp::IVP{<:AbstractContinuousSystem}, tspan;
-              Δt0::TimeInterval=zeroI, kwargs...)
-    δ = alg.δ
-
-    NSTEPS = get(kwargs, :NSTEPS, compute_nsteps(δ, tspan))
-
-    # normalize system to canonical form
-    ivp_norm = _normalize(ivp)
-
-    # homogenize the initial-value problem
-    if get(kwargs, :homogenize, false)
-        ivp_norm = homogenize(ivp_norm)
-    end
-
-    # discretize system
-    ivp_discr = discretize(ivp_norm, δ, alg.approx_model)
-
-    return post(alg, ivp_discr, NSTEPS; Δt0=Δt0, kwargs...)
-end
+# default continuous post
 
 # discrete post
 function post(alg::LGG09{N,AM,VN,TN}, ivp::IVP{<:AbstractDiscreteSystem}, NSTEPS=nothing;

diff --git a/src/Algorithms/VREP/post.jl b/src/Algorithms/VREP/post.jl
@@ -1,23 +1,4 @@
-# continuous post
-function post(alg::VREP, ivp::IVP{<:AbstractContinuousSystem}, tspan;
-              Δt0::TimeInterval=zeroI, kwargs...)
-    δ = alg.δ
-
-    NSTEPS = get(kwargs, :NSTEPS, compute_nsteps(δ, tspan))
-
-    # normalize system to canonical form
-    ivp_norm = _normalize(ivp)
-
-    # homogenize the initial-value problem
-    if get(kwargs, :homogenize, false)
-        ivp_norm = homogenize(ivp_norm)
-    end
-
-    # discretize system
-    ivp_discr = discretize(ivp_norm, δ, alg.approx_model)
-
-    return post(alg, ivp_discr, NSTEPS; Δt0=Δt0, kwargs...)
-end
+# default continuous post
 
 # discrete post
 function post(alg::VREP{N}, ivp::IVP{<:AbstractDiscreteSystem}, NSTEPS=nothing;

diff --git a/src/Algorithms/linear_post.jl b/src/Algorithms/linear_post.jl
@@ -0,0 +1,22 @@
+# default continuous post for linear systems
+function post(alg::Union{A20,BFFPSV18,BOX,GLGM06,INT,LGG09,VREP},
+              ivp::IVP{<:AbstractContinuousSystem}, tspan;
+              Δt0::TimeInterval=zeroI, kwargs...)
+    δ = alg.δ
+
+    NSTEPS = get(kwargs, :NSTEPS, compute_nsteps(δ, tspan))
+
+    # normalize system to canonical form
+    ivp_norm = _normalize(ivp)
+
+    # homogenize system
+    if get(kwargs, :homogenize, false)
+        ivp_norm = homogenize(ivp_norm)
+    end
+
+    # discretize system
+    ivp_discr = discretize(ivp_norm, δ, alg.approx_model)
+
+    # apply discrete post to discretized system
+    return post(alg, ivp_discr, NSTEPS; Δt0=Δt0, kwargs...)
+end
diff --git a/src/ReachabilityAnalysis.jl b/src/ReachabilityAnalysis.jl
@@ -45,6 +45,7 @@ include("Algorithms/ASB07/ASB07.jl")
 include("Algorithms/ORBIT/ORBIT.jl")
 include("Algorithms/VREP/VREP.jl")
 include("Algorithms/FLOWSTAR/FLOWSTAR.jl")
+include("Algorithms/linear_post.jl")
 
 # Nonlinear
 include("Algorithms/CARLIN/CARLIN.jl")