Add steady state matrices.

[samuelsonric]

Example use.

using AlgebraicDynamics
using AlgebraicDynamics.Steady

n = 100

# disrete dynamical systems
M1 = Machine(rand(1:n, n, n), rand(1:n, n))
M2 = Machine(rand(1:n, n, n), rand(1:n, n))
M3 = Machine(rand(1:n, n, n), rand(1:n, n))

# steady-state matrices
S1 = ss(M1)
S2 = ss(M2)
S3 = ss(M3)

# compose steady-state matrices
S3 * S2 * S1

We can also compose the steady-state matrices using WiringDiagrams.jl.

using WiringDiagrams, OMEinsum

# composition pattern
d = WiringDiagram([[1, 2], [2, 3], [3, 4]], [1, 4], Dict(1 => n, 2 => n, 3 => n, 4 => n))

# wiring diagram algebra
a = ArrayAlgebra{Array{Float64}}()

# compose steady-state matrices
a(d)(S3, S2, S1)