debug: SimResults doctest

Author: franckgaga

src/estimator/execute.jl

@@ -7,7 +7,7 @@ Also store current inputs without operating points `u0` in `estim.lastu0`. This
 used for [`PredictiveController`](@ref) computations.
 """
 function remove_op!(estim::StateEstimator, u, ym, d)
-    u0, d0 = estim.model.buffer.u, estim.model.buffer.d
+    u0, d0 = estim.buffer.u, estim.buffer.d
     y0m = estim.buffer.ym
     u0  .= u  .- estim.model.uop
     y0m .= @views ym .- estim.model.yop[estim.i_ym]

src/plot_sim.jl

@@ -41,9 +41,11 @@ Simply call `plot` on them.
 
 # Examples
 ```jldoctest
-julia> model = LinModel(tf(1, [1, 1]), 1.0); N = 5; U_data = fill(1.0, 1, N);
+julia> model = LinModel(tf(1, [1, 1]), 1.0);
 
-julia> Y_data = reduce(hcat, (updatestate!(model, U_data[:, i]); model()) for i=1:N)
+julia> N = 5; U_data = fill(1.0, 1, N); Y_data = zeros(1, N);
+
+julia> for i=1:N; updatestate!(model, U_data[:, i]); Y_data[:, i] = model(); end; Y_data
 1×5 Matrix{Float64}:
  0.632121  0.864665  0.950213  0.981684  0.993262
 

src/sim_model.jl

@@ -234,10 +234,9 @@ julia> x = updatestate!(model, [1])
 """
 function updatestate!(model::SimModel{NT}, u, d=model.buffer.empty) where NT <: Real
     validate_args(model::SimModel, d, u)
-    u0, d0 = model.buffer.u, model.buffer.d
+    u0, d0, xnext0 = model.buffer.u, model.buffer.d, model.buffer.x
     u0 .= u .- model.uop
     d0 .= d .- model.dop
-    xnext0 = model.buffer.x
     f!(xnext0, model, model.x0, u0, d0)
     xnext0  .+= model.fop .- model.xop
     model.x0 .= xnext0
@@ -264,9 +263,8 @@ julia> y = evaloutput(model)
 """
 function evaloutput(model::SimModel{NT}, d=model.buffer.empty) where NT <: Real
     validate_args(model, d)
-    d0  = model.buffer.d
+    d0, y0  = model.buffer.d, model.buffer.y
     d0 .= d .- model.dop
-    y0 = model.buffer.y
     h!(y0, model, model.x0, d0)
     y   = y0
     y .+= model.yop