fix

Author: ValentinKaisermayer

test/Electrical/dc_motor.jl

@@ -66,11 +66,11 @@ D = Differential(t)
         # test steady-state values
         dc_gain = [f/(k^2 + f * R) k/(k^2 + f * R); k/(k^2 + f * R) -R/(k^2 + f * R)]
         idx_t = findfirst(sol.t .> 2.5)
-        @test sol[inertia.w][idx_t] ≈ (dc_gain * [V_step; 0])[2]
-        @test sol[emf.i][idx_t] ≈ (dc_gain * [V_step; 0])[1]
+        @test sol[inertia.w][idx_t]≈(dc_gain * [V_step; 0])[2] rtol=1e-3
+        @test sol[emf.i][idx_t]≈(dc_gain * [V_step; 0])[1] rtol=1e-3
         idx_t = findfirst(sol.t .> 5.5)
-        @test sol[inertia.w][idx_t] ≈ (dc_gain * [V_step; -tau_L_step])[2]
-        @test sol[emf.i][idx_t] ≈ (dc_gain * [V_step; -tau_L_step])[1]
+        @test sol[inertia.w][idx_t]≈(dc_gain * [V_step; -tau_L_step])[2] rtol=1e-3
+        @test sol[emf.i][idx_t]≈(dc_gain * [V_step; -tau_L_step])[1] rtol=1e-3
     end
 
     prob = DAEProblem(sys, D.(states(sys)) .=> 0.0, Pair[], (0, 6.0))
@@ -82,11 +82,11 @@ D = Differential(t)
     # test steady-state values
     dc_gain = [f/(k^2 + f * R) k/(k^2 + f * R); k/(k^2 + f * R) -R/(k^2 + f * R)]
     idx_t = findfirst(sol.t .> 2.5)
-    @test sol[inertia.w][idx_t] ≈ (dc_gain * [V_step; 0])[2]
-    @test sol[emf.i][idx_t] ≈ (dc_gain * [V_step; 0])[1]
+    @test sol[inertia.w][idx_t]≈(dc_gain * [V_step; 0])[2] rtol=1e-3
+    @test sol[emf.i][idx_t]≈(dc_gain * [V_step; 0])[1] rtol=1e-3
     idx_t = findfirst(sol.t .> 5.5)
-    @test sol[inertia.w][idx_t] ≈ (dc_gain * [V_step; -tau_L_step])[2]
-    @test sol[emf.i][idx_t] ≈ (dc_gain * [V_step; -tau_L_step])[1]
+    @test sol[inertia.w][idx_t]≈(dc_gain * [V_step; -tau_L_step])[2] rtol=1e-3
+    @test sol[emf.i][idx_t]≈(dc_gain * [V_step; -tau_L_step])[1] rtol=1e-3
 
     # p1 = Plots.plot(sol, vars=[inertia.w], ylabel="Angular Vel. in rad/s", label="")
     # p2 = Plots.plot(sol, vars=[emf.i], ylabel="Current in A", label="")
@@ -155,14 +155,14 @@ end
         # test steady-state values
         dc_gain = [f/(k^2 + f * R) k/(k^2 + f * R); k/(k^2 + f * R) -R/(k^2 + f * R)]
         idx_t = findfirst(sol.t .> 2.5)
-        @test sol[inertia.w][idx_t] ≈ (dc_gain * [V_step; 0])[2]
-        @test sol[emf.i][idx_t] ≈ (dc_gain * [V_step; 0])[1]
+        @test sol[inertia.w][idx_t]≈(dc_gain * [V_step; 0])[2] rtol=1e-3
+        @test sol[emf.i][idx_t]≈(dc_gain * [V_step; 0])[1] rtol=1e-3
         idx_t = findfirst(sol.t .> 5.5)
-        @test sol[inertia.w][idx_t] ≈ (dc_gain * [V_step; -tau_L_step])[2]
-        @test sol[emf.i][idx_t] ≈ (dc_gain * [V_step; -tau_L_step])[1]
+        @test sol[inertia.w][idx_t]≈(dc_gain * [V_step; -tau_L_step])[2] rtol=1e-3
+        @test sol[emf.i][idx_t]≈(dc_gain * [V_step; -tau_L_step])[1] rtol=1e-3
 
         # 
-        @test sol[inertia.w] == sol[speed_sensor.w]
+        @test all(sol[inertia.w] .== sol[speed_sensor.w.u])
     end
 
     prob = DAEProblem(sys, D.(states(sys)) .=> 0.0, Pair[], (0, 6.0))
@@ -175,13 +175,13 @@ end
     # test steady-state values
     dc_gain = [f/(k^2 + f * R) k/(k^2 + f * R); k/(k^2 + f * R) -R/(k^2 + f * R)]
     idx_t = findfirst(sol.t .> 2.5)
-    @test sol[inertia.w][idx_t] ≈ (dc_gain * [V_step; 0])[2]
-    @test sol[emf.i][idx_t] ≈ (dc_gain * [V_step; 0])[1]
+    @test sol[inertia.w][idx_t]≈(dc_gain * [V_step; 0])[2] rtol=1e-3
+    @test sol[emf.i][idx_t]≈(dc_gain * [V_step; 0])[1] rtol=1e-3
     idx_t = findfirst(sol.t .> 5.5)
-    @test sol[inertia.w][idx_t] ≈ (dc_gain * [V_step; -tau_L_step])[2]
-    @test sol[emf.i][idx_t] ≈ (dc_gain * [V_step; -tau_L_step])[1]
+    @test sol[inertia.w][idx_t]≈(dc_gain * [V_step; -tau_L_step])[2] rtol=1e-3
+    @test sol[emf.i][idx_t]≈(dc_gain * [V_step; -tau_L_step])[1] rtol=1e-3
     # 
-    @test sol[inertia.w] == sol[speed_sensor.w.u]
+    @test all(sol[inertia.w] .== sol[speed_sensor.w.u])
 end
 
 @testset "El. Heating Circuit" begin