From fc11d80c7539af0c4bcc714fef0dfe4f8b6fc86d Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Sun, 29 Sep 2024 06:55:35 -0400
Subject: [PATCH 01/13] removed defaults from SISO component

---
 src/Blocks/utils.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/Blocks/utils.jl b/src/Blocks/utils.jl
index 10ad0261f..7b845fb14 100644
--- a/src/Blocks/utils.jl
+++ b/src/Blocks/utils.jl
@@ -130,8 +130,8 @@ Single input single output (SISO) continuous system block.
         y_start = 0.0
     end
     @variables begin
-        u(t) = u_start, [description = "Input of SISO system"]
-        y(t) = y_start, [description = "Output of SISO system"]
+        u(t), [guess = u_start, description = "Input of SISO system"]
+        y(t), [guess = y_start, description = "Output of SISO system"]
     end
     @components begin
         input = RealInput(guess = u_start)

From 5baefd5794cbd356e7e3e096f20e9322f734ff01 Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Sun, 29 Sep 2024 07:25:52 -0400
Subject: [PATCH 02/13] Added a test for SISO initial value set

---
 test/Blocks/utils.jl | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)

diff --git a/test/Blocks/utils.jl b/test/Blocks/utils.jl
index 6bde3e1d4..a88d0a637 100644
--- a/test/Blocks/utils.jl
+++ b/test/Blocks/utils.jl
@@ -1,5 +1,6 @@
 using ModelingToolkitStandardLibrary.Blocks
 using ModelingToolkit
+using OrdinaryDiffEq
 
 @testset "Array Guesses" begin
     for (block, guess) in [
@@ -23,6 +24,22 @@ end
     end
 end
 
+@testset "SISO Check" begin
+    k, w, d = 1.0, 1.0, 0.5
+    @named c = Constant(; k = 1)
+    @named so = SecondOrder(; k = k, w = w, d = d, xd = 1)
+    @named iosys = ODESystem(connect(c.output, so.input), t, systems = [so, c])
+    sys = structural_simplify(iosys)
+
+    initsys = ModelingToolkit.generate_initializesystem(sys)
+    initsys = structural_simplify(initsys)
+    initprob = NonlinearProblem(initsys, [t=>0])
+    initsol = solve(initprob)
+
+    @test initsol[sys.so.xd] == 1.0 
+    @test initsol[sys.so.u] == 1.0
+end
+
 @test_deprecated RealInput(; name = :a, u_start = 1.0)
 @test_deprecated RealInput(; name = :a, nin = 2, u_start = ones(2))
 @test_deprecated RealOutput(; name = :a, u_start = 1.0)

From 0b8b50143eb4c68f560861ff1c8af7710a24ec09 Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Sun, 29 Sep 2024 07:34:11 -0400
Subject: [PATCH 03/13] format

---
 test/Blocks/utils.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/test/Blocks/utils.jl b/test/Blocks/utils.jl
index a88d0a637..c87177611 100644
--- a/test/Blocks/utils.jl
+++ b/test/Blocks/utils.jl
@@ -33,10 +33,10 @@ end
 
     initsys = ModelingToolkit.generate_initializesystem(sys)
     initsys = structural_simplify(initsys)
-    initprob = NonlinearProblem(initsys, [t=>0])
+    initprob = NonlinearProblem(initsys, [t => 0])
     initsol = solve(initprob)
 
-    @test initsol[sys.so.xd] == 1.0 
+    @test initsol[sys.so.xd] == 1.0
     @test initsol[sys.so.u] == 1.0
 end
 

From 45a9f3a5ee5a35d55fa49cde706b2a76f62be1b8 Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Sun, 29 Sep 2024 08:03:58 -0400
Subject: [PATCH 04/13] fixed missing t

---
 test/Blocks/utils.jl | 1 +
 1 file changed, 1 insertion(+)

diff --git a/test/Blocks/utils.jl b/test/Blocks/utils.jl
index c87177611..624a84acb 100644
--- a/test/Blocks/utils.jl
+++ b/test/Blocks/utils.jl
@@ -1,6 +1,7 @@
 using ModelingToolkitStandardLibrary.Blocks
 using ModelingToolkit
 using OrdinaryDiffEq
+using ModelingToolkit: t_nounits as t, D_nounits as D
 
 @testset "Array Guesses" begin
     for (block, guess) in [

From 6c3269202e488059a91cface4cd83787faea66b1 Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Sun, 29 Sep 2024 06:55:35 -0400
Subject: [PATCH 05/13] removed defaults from SISO component

---
 src/Blocks/utils.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/Blocks/utils.jl b/src/Blocks/utils.jl
index 10ad0261f..7b845fb14 100644
--- a/src/Blocks/utils.jl
+++ b/src/Blocks/utils.jl
@@ -130,8 +130,8 @@ Single input single output (SISO) continuous system block.
         y_start = 0.0
     end
     @variables begin
-        u(t) = u_start, [description = "Input of SISO system"]
-        y(t) = y_start, [description = "Output of SISO system"]
+        u(t), [guess = u_start, description = "Input of SISO system"]
+        y(t), [guess = y_start, description = "Output of SISO system"]
     end
     @components begin
         input = RealInput(guess = u_start)

From c37236153700cb55748334dc82a1d1b59f5a2757 Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Sun, 29 Sep 2024 07:25:52 -0400
Subject: [PATCH 06/13] Added a test for SISO initial value set

---
 test/Blocks/utils.jl | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)

diff --git a/test/Blocks/utils.jl b/test/Blocks/utils.jl
index 6bde3e1d4..a88d0a637 100644
--- a/test/Blocks/utils.jl
+++ b/test/Blocks/utils.jl
@@ -1,5 +1,6 @@
 using ModelingToolkitStandardLibrary.Blocks
 using ModelingToolkit
+using OrdinaryDiffEq
 
 @testset "Array Guesses" begin
     for (block, guess) in [
@@ -23,6 +24,22 @@ end
     end
 end
 
+@testset "SISO Check" begin
+    k, w, d = 1.0, 1.0, 0.5
+    @named c = Constant(; k = 1)
+    @named so = SecondOrder(; k = k, w = w, d = d, xd = 1)
+    @named iosys = ODESystem(connect(c.output, so.input), t, systems = [so, c])
+    sys = structural_simplify(iosys)
+
+    initsys = ModelingToolkit.generate_initializesystem(sys)
+    initsys = structural_simplify(initsys)
+    initprob = NonlinearProblem(initsys, [t=>0])
+    initsol = solve(initprob)
+
+    @test initsol[sys.so.xd] == 1.0 
+    @test initsol[sys.so.u] == 1.0
+end
+
 @test_deprecated RealInput(; name = :a, u_start = 1.0)
 @test_deprecated RealInput(; name = :a, nin = 2, u_start = ones(2))
 @test_deprecated RealOutput(; name = :a, u_start = 1.0)

From f76ce5730d017be9c5a741beb079922b3788daf9 Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Sun, 29 Sep 2024 07:34:11 -0400
Subject: [PATCH 07/13] format

---
 test/Blocks/utils.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/test/Blocks/utils.jl b/test/Blocks/utils.jl
index a88d0a637..c87177611 100644
--- a/test/Blocks/utils.jl
+++ b/test/Blocks/utils.jl
@@ -33,10 +33,10 @@ end
 
     initsys = ModelingToolkit.generate_initializesystem(sys)
     initsys = structural_simplify(initsys)
-    initprob = NonlinearProblem(initsys, [t=>0])
+    initprob = NonlinearProblem(initsys, [t => 0])
     initsol = solve(initprob)
 
-    @test initsol[sys.so.xd] == 1.0 
+    @test initsol[sys.so.xd] == 1.0
     @test initsol[sys.so.u] == 1.0
 end
 

From 3efe41184e664fcd520aab760d6fab9467086260 Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Sun, 29 Sep 2024 08:03:58 -0400
Subject: [PATCH 08/13] fixed missing t

---
 test/Blocks/utils.jl | 1 +
 1 file changed, 1 insertion(+)

diff --git a/test/Blocks/utils.jl b/test/Blocks/utils.jl
index c87177611..624a84acb 100644
--- a/test/Blocks/utils.jl
+++ b/test/Blocks/utils.jl
@@ -1,6 +1,7 @@
 using ModelingToolkitStandardLibrary.Blocks
 using ModelingToolkit
 using OrdinaryDiffEq
+using ModelingToolkit: t_nounits as t, D_nounits as D
 
 @testset "Array Guesses" begin
     for (block, guess) in [

From ef096b97ce6b13425cb1863a0d121d4c4a8647af Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Wed, 2 Oct 2024 10:55:58 -0400
Subject: [PATCH 09/13] fixed connections.md

---
 docs/src/connectors/connections.md         | 28 +++++++++-------
 src/Mechanical/Translational/components.jl | 38 ++++++++++------------
 2 files changed, 34 insertions(+), 32 deletions(-)

diff --git a/docs/src/connectors/connections.md b/docs/src/connectors/connections.md
index 01381db43..1111a7778 100644
--- a/docs/src/connectors/connections.md
+++ b/docs/src/connectors/connections.md
@@ -136,8 +136,8 @@ using ModelingToolkitStandardLibrary
 const TV = ModelingToolkitStandardLibrary.Mechanical.Translational
 
 systems = @named begin
-    damping = TV.Damper(d = 1, flange_a.v = 1)
-    body = TV.Mass(m = 1, v = 1)
+    damping = TV.Damper(d = 1)
+    body = TV.Mass(m = 1)
     ground = TV.Fixed()
 end
 
@@ -155,7 +155,7 @@ nothing # hide
 As expected, we have a similar solution…
 
 ```@example connections
-prob = ODEProblem(sys, [], (0, 10.0), [])
+prob = ODEProblem(sys, [], (0, 10.0), [];   initialization_eqs=[sys.body.s ~ 0, sys.body.v ~ 1])
 sol_v = solve(prob)
 
 p1 = plot(sol_v, idxs = [body.v])
@@ -228,7 +228,7 @@ In this problem, we have a mass, spring, and damper which are connected to a fix
 The damper will connect the flange/flange 1 (`flange_a`) to the mass, and flange/flange 2 (`flange_b`) to the fixed point.  For both position- and velocity-based domains, we set the damping constant `d=1` and `va=1` and leave the default for `v_b_0` at 0.  For the position domain, we also need to set the initial positions for `flange_a` and `flange_b`.
 
 ```@example connections
-@named dv = TV.Damper(d = 1, flange_a.v = 1)
+@named dv = TV.Damper(d = 1)
 @named dp = TP.Damper(d = 1, va = 1, vb = 0.0, flange_a__s = 3, flange_b__s = 1)
 nothing # hide
 ```
@@ -238,7 +238,7 @@ nothing # hide
 The spring will connect the flange/flange 1 (`flange_a`) to the mass, and flange/flange 2 (`flange_b`) to the fixed point.  For both position- and velocity-based domains, we set the spring constant `k=1`.  The velocity domain then requires the initial velocity `va` and initial spring stretch `delta_s`.  The position domain instead needs the initial positions for `flange_a` and `flange_b` and the natural spring length `l`.
 
 ```@example connections
-@named sv = TV.Spring(k = 1, flange_a__v = 1, delta_s = 1)
+@named sv = TV.Spring(k = 1)
 @named sp = TP.Spring(k = 1, flange_a__s = 3, flange_b__s = 1, l = 1)
 nothing # hide
 ```
@@ -248,7 +248,7 @@ nothing # hide
 For both position- and velocity-based domains, we set the mass `m=1` and initial velocity `v=1`. Like the damper, the position domain requires the position initial conditions set as well.
 
 ```@example connections
-@named bv = TV.Mass(m = 1, v = 1)
+@named bv = TV.Mass(m = 1)
 @named bp = TP.Mass(m = 1, v = 1, s = 3)
 nothing # hide
 ```
@@ -270,7 +270,7 @@ As can be seen, the position-based domain requires more initial condition inform
 Let's define a quick function to simplify and solve the 2 different systems. Note, we will solve with a fixed time step and a set tolerance to compare the numerical differences.
 
 ```@example connections
-function simplify_and_solve(damping, spring, body, ground)
+function simplify_and_solve(damping, spring, body, ground; initialization_eqs = Equation[])
     eqs = [connect(spring.flange_a, body.flange, damping.flange_a)
            connect(spring.flange_b, damping.flange_b, ground.flange)]
 
@@ -280,8 +280,8 @@ function simplify_and_solve(damping, spring, body, ground)
 
     println.(full_equations(sys))
 
-    prob = ODEProblem(sys, [], (0, 10.0), [])
-    sol = solve(prob; dt = 0.1, adaptive = false, reltol = 1e-9, abstol = 1e-9)
+    prob = ODEProblem(sys, [], (0, 10.0), []; initialization_eqs)
+    sol = solve(prob; abstol=1e-9, reltol=1e-9)
 
     return sol
 end
@@ -291,7 +291,12 @@ nothing # hide
 Now let's solve the velocity domain model
 
 ```@example connections
-solv = simplify_and_solve(dv, sv, bv, gv);
+initialization_eqs = [
+  bv.s ~ 3
+  bv.v ~ 1
+  sv.delta_s ~ 1
+]
+solv = simplify_and_solve(dv, sv, bv, gv; initialization_eqs);
 nothing # hide
 ```
 
@@ -346,12 +351,11 @@ By definition, the spring stretch is
 \Delta s = s - s_{b_0} - l
 ```
 
-Which means both systems are actually solving the same exact system.  We can plot the numerical difference between the 2 systems and see the result is negligible.
+Which means both systems are actually solving the same exact system.  We can plot the numerical difference between the 2 systems and see the result is negligible (much less than the tolerance of 1e-9).
 
 ```@example connections
 plot(title = "numerical difference: vel. vs. pos. domain", xlabel = "time [s]",
     ylabel = "solv[bv.v] .- solp[bp.v]")
 time = 0:0.1:10
 plot!(time, (solv(time)[bv.v] .- solp(time)[bp.v]), label = "")
-Plots.ylims!(-1e-15, 1e-15)
 ```
diff --git a/src/Mechanical/Translational/components.jl b/src/Mechanical/Translational/components.jl
index c29da00e0..b890b6711 100644
--- a/src/Mechanical/Translational/components.jl
+++ b/src/Mechanical/Translational/components.jl
@@ -65,9 +65,9 @@ Sliding mass with inertia
     @named flange = MechanicalPort()
 
     vars = @variables begin
-        s(t)
-        v(t)
-        f(t)
+        s(t), [guess=0]
+        v(t), [guess=0]
+        f(t), [guess=0]
     end
 
     eqs = [flange.v ~ v
@@ -98,22 +98,21 @@ Linear 1D translational spring
   - `flange_a`: 1-dim. translational flange on one side of spring
   - `flange_b`: 1-dim. translational flange on opposite side of spring
 """
-@component function Spring(; name, k, delta_s = 0.0, flange_a__v = 0.0, flange_b__v = 0.0)
-    Spring(REL; name, k, delta_s, flange_a__v, flange_b__v)
+@component function Spring(; name, k)
+    Spring(REL; name, k)
 end # default
 
-@component function Spring(::Val{:relative}; name, k, delta_s = 0.0, flange_a__v = 0.0,
-        flange_b__v = 0.0)
+@component function Spring(::Val{:relative}; name, k)
     pars = @parameters begin
         k = k
     end
     vars = @variables begin
-        delta_s(t) = delta_s
-        f(t)
+        delta_s(t), [guess = 0]
+        f(t), [guess=0]
     end
 
-    @named flange_a = MechanicalPort(; v = flange_a__v)
-    @named flange_b = MechanicalPort(; v = flange_b__v)
+    @named flange_a = MechanicalPort()
+    @named flange_b = MechanicalPort()
 
     eqs = [D(delta_s) ~ flange_a.v - flange_b.v
            f ~ k * delta_s
@@ -125,20 +124,19 @@ end # default
 end
 
 const ABS = Val(:absolute)
-@component function Spring(::Val{:absolute}; name, k, sa = 0, sb = 0, flange_a__v = 0.0,
-        flange_b__v = 0.0, l = 0)
+@component function Spring(::Val{:absolute}; name, k, l = 0)
     pars = @parameters begin
         k = k
         l = l
     end
     vars = @variables begin
-        sa(t) = sa
-        sb(t) = sb
-        f(t)
+        sa(t), [guess=0]
+        sb(t), [guess=0]
+        f(t), [guess=0]
     end
 
-    @named flange_a = MechanicalPort(; v = flange_a__v)
-    @named flange_b = MechanicalPort(; v = flange_b__v)
+    @named flange_a = MechanicalPort()
+    @named flange_b = MechanicalPort()
 
     eqs = [D(sa) ~ flange_a.v
            D(sb) ~ flange_b.v
@@ -169,8 +167,8 @@ Linear 1D translational damper
         d
     end
     @variables begin
-        v(t)
-        f(t)
+        v(t), [guess=0]
+        f(t), [guess=0]
     end
 
     @components begin

From b2f6c93e3848c0ebb617f8e01643f7d0ecc6858d Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Wed, 2 Oct 2024 11:05:15 -0400
Subject: [PATCH 10/13] format

---
 docs/src/connectors/connections.md         | 13 ++++++-------
 src/Mechanical/Translational/components.jl | 18 +++++++++---------
 2 files changed, 15 insertions(+), 16 deletions(-)

diff --git a/docs/src/connectors/connections.md b/docs/src/connectors/connections.md
index 1111a7778..1e5a3270f 100644
--- a/docs/src/connectors/connections.md
+++ b/docs/src/connectors/connections.md
@@ -155,7 +155,8 @@ nothing # hide
 As expected, we have a similar solution…
 
 ```@example connections
-prob = ODEProblem(sys, [], (0, 10.0), [];   initialization_eqs=[sys.body.s ~ 0, sys.body.v ~ 1])
+prob = ODEProblem(
+    sys, [], (0, 10.0), []; initialization_eqs = [sys.body.s ~ 0, sys.body.v ~ 1])
 sol_v = solve(prob)
 
 p1 = plot(sol_v, idxs = [body.v])
@@ -281,7 +282,7 @@ function simplify_and_solve(damping, spring, body, ground; initialization_eqs =
     println.(full_equations(sys))
 
     prob = ODEProblem(sys, [], (0, 10.0), []; initialization_eqs)
-    sol = solve(prob; abstol=1e-9, reltol=1e-9)
+    sol = solve(prob; abstol = 1e-9, reltol = 1e-9)
 
     return sol
 end
@@ -291,11 +292,9 @@ nothing # hide
 Now let's solve the velocity domain model
 
 ```@example connections
-initialization_eqs = [
-  bv.s ~ 3
-  bv.v ~ 1
-  sv.delta_s ~ 1
-]
+initialization_eqs = [bv.s ~ 3
+                      bv.v ~ 1
+                      sv.delta_s ~ 1]
 solv = simplify_and_solve(dv, sv, bv, gv; initialization_eqs);
 nothing # hide
 ```
diff --git a/src/Mechanical/Translational/components.jl b/src/Mechanical/Translational/components.jl
index b890b6711..2df7f8afb 100644
--- a/src/Mechanical/Translational/components.jl
+++ b/src/Mechanical/Translational/components.jl
@@ -65,9 +65,9 @@ Sliding mass with inertia
     @named flange = MechanicalPort()
 
     vars = @variables begin
-        s(t), [guess=0]
-        v(t), [guess=0]
-        f(t), [guess=0]
+        s(t), [guess = 0]
+        v(t), [guess = 0]
+        f(t), [guess = 0]
     end
 
     eqs = [flange.v ~ v
@@ -108,7 +108,7 @@ end # default
     end
     vars = @variables begin
         delta_s(t), [guess = 0]
-        f(t), [guess=0]
+        f(t), [guess = 0]
     end
 
     @named flange_a = MechanicalPort()
@@ -130,9 +130,9 @@ const ABS = Val(:absolute)
         l = l
     end
     vars = @variables begin
-        sa(t), [guess=0]
-        sb(t), [guess=0]
-        f(t), [guess=0]
+        sa(t), [guess = 0]
+        sb(t), [guess = 0]
+        f(t), [guess = 0]
     end
 
     @named flange_a = MechanicalPort()
@@ -167,8 +167,8 @@ Linear 1D translational damper
         d
     end
     @variables begin
-        v(t), [guess=0]
-        f(t), [guess=0]
+        v(t), [guess = 0]
+        f(t), [guess = 0]
     end
 
     @components begin

From d2667802c064e194481ddea35680c1fdcdf068f0 Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Thu, 3 Oct 2024 14:42:02 -0400
Subject: [PATCH 11/13] format

---
 src/Blocks/nonlinear.jl | 38 ++++++++++++++++++++++++--------------
 1 file changed, 24 insertions(+), 14 deletions(-)

diff --git a/src/Blocks/nonlinear.jl b/src/Blocks/nonlinear.jl
index 8f2466f72..8dd3eaf79 100644
--- a/src/Blocks/nonlinear.jl
+++ b/src/Blocks/nonlinear.jl
@@ -93,20 +93,30 @@ Initial value of state `Y` can be set with `int.y`
   - `input`
   - `output`
 """
-@mtkmodel SlewRateLimiter begin
-    @parameters begin
-        rising = 1.0, [description = "Maximum rising slew rate of SlewRateLimiter"]
-        falling = -rising, [description = "Derivative time constant of SlewRateLimiter"]
-        Td = 0.001, [description = "Derivative time constant"]
-    end
-    begin
-        getdefault(rising) ≥ getdefault(falling) ||
-            throw(ArgumentError("`rising` must be smaller than `falling`"))
-        getdefault(Td) > 0 ||
-            throw(ArgumentError("Time constant `Td` must be strictly positive"))
+@component function SlewRateLimiter(;
+        name, y_start = 0.0, rising = 1.0, falling = -rising, Td = 0.001)
+    pars = @parameters begin
+        rising = rising, [description = "Maximum rising slew rate of SlewRateLimiter"]
+        falling = falling, [description = "Derivative time constant of SlewRateLimiter"]
+        Td = Td, [description = "Derivative time constant"]
+        y_start = y_start
     end
-    @extend u, y = siso = SISO(; y_start)
-    @equations begin
+
+    getdefault(rising) ≥ getdefault(falling) ||
+        throw(ArgumentError("`rising` must be smaller than `falling`"))
+    getdefault(Td) > 0 ||
+        throw(ArgumentError("Time constant `Td` must be strictly positive"))
+
+    @named siso = SISO(; y_start)
+    @unpack y, u = siso
+
+    eqs = [
         D(y) ~ max(min((u - y) / Td, rising), falling)
-    end
+    ]
+
+    initialization_eqs = [
+        y ~ y_start
+    ]
+
+    return extend(ODESystem(eqs, t, [], pars; name, initialization_eqs), siso)
 end

From d881c14216f0d96a4f8d08b756a346ffc7c9d58a Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Thu, 3 Oct 2024 14:54:55 -0400
Subject: [PATCH 12/13] fixed translational test

---
 test/Mechanical/translational.jl | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/test/Mechanical/translational.jl b/test/Mechanical/translational.jl
index 2ed07eb05..0c72e2c53 100644
--- a/test/Mechanical/translational.jl
+++ b/test/Mechanical/translational.jl
@@ -34,7 +34,7 @@ end
     @named dv = TV.Damper(d = 1)
     @named dp = TP.Damper(d = 1, va = 1, vb = 0.0, flange_a.s = 3, flange_b.s = 1)
 
-    @named sv = TV.Spring(k = 1, flange_a__v = 1, delta_s = 1)
+    @named sv = TV.Spring(k = 1)
     @named sp = TP.Spring(k = 1, flange_a__s = 3, flange_b__s = 1, l = 1)
 
     @named bv = TV.Mass(m = 1)
@@ -43,7 +43,7 @@ end
     @named gv = TV.Fixed()
     @named gp = TP.Fixed(s_0 = 1)
 
-    function simplify_and_solve(damping, spring, body, ground)
+    function simplify_and_solve(damping, spring, body, ground; initialization_eqs=Equation[])
         eqs = [connect(spring.flange_a, body.flange, damping.flange_a)
                connect(spring.flange_b, damping.flange_b, ground.flange)]
 
@@ -51,13 +51,13 @@ end
 
         sys = structural_simplify(model)
 
-        prob = ODEProblem(sys, [body.s => 0], (0, 20.0), [])
-        sol = solve(prob, ImplicitMidpoint(), dt = 0.01)
+        prob = ODEProblem(sys, [], (0, 20.0), []; initialization_eqs)
+        sol = solve(prob; abstol = 1e-9, reltol = 1e-9)
 
         return sol
     end
 
-    solv = simplify_and_solve(dv, sv, bv, gv)
+    solv = simplify_and_solve(dv, sv, bv, gv; initialization_eqs=[bv.s ~ 3, bv.v ~ 1, sv.delta_s ~ 1])
     solp = simplify_and_solve(dp, sp, bp, gp)
 
     @test solv[bv.v][1] == 1.0

From 13b930827f7fdea489af86b2cf93204c2aa1917c Mon Sep 17 00:00:00 2001
From: Brad Carman <bradleygcarman@outlook.com>
Date: Thu, 3 Oct 2024 14:59:05 -0400
Subject: [PATCH 13/13] format

---
 test/Mechanical/translational.jl | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/test/Mechanical/translational.jl b/test/Mechanical/translational.jl
index 0c72e2c53..ec077de13 100644
--- a/test/Mechanical/translational.jl
+++ b/test/Mechanical/translational.jl
@@ -43,7 +43,8 @@ end
     @named gv = TV.Fixed()
     @named gp = TP.Fixed(s_0 = 1)
 
-    function simplify_and_solve(damping, spring, body, ground; initialization_eqs=Equation[])
+    function simplify_and_solve(
+            damping, spring, body, ground; initialization_eqs = Equation[])
         eqs = [connect(spring.flange_a, body.flange, damping.flange_a)
                connect(spring.flange_b, damping.flange_b, ground.flange)]
 
@@ -57,7 +58,8 @@ end
         return sol
     end
 
-    solv = simplify_and_solve(dv, sv, bv, gv; initialization_eqs=[bv.s ~ 3, bv.v ~ 1, sv.delta_s ~ 1])
+    solv = simplify_and_solve(
+        dv, sv, bv, gv; initialization_eqs = [bv.s ~ 3, bv.v ~ 1, sv.delta_s ~ 1])
     solp = simplify_and_solve(dp, sp, bp, gp)
 
     @test solv[bv.v][1] == 1.0