# Plans to port SLICOT Toeplitz solvers to Julia?

**URL:** https://discourse.julialang.org/t/plans-to-port-slicot-toeplitz-solvers-to-julia/55996
**Category:** General Usage
**Tags:** question
**Created:** [February 25, 2021, 9:21am UTC](https://discourse.julialang.org/t/plans-to-port-slicot-toeplitz-solvers-to-julia/55996 "2021-02-25T09:21:52Z")
**Posts on this page:** 6
**Page:** 1

<div class="post-metadata">

### Author: ![BambOoxX](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/bambooxx/32/22179_2.png) [@BambOoxX](https://discourse.julialang.org/u/BambOoxX)
#### Post date: [February 25, 2021, 9:21am UTC](https://discourse.julialang.org/t/plans-to-port-slicot-toeplitz-solvers-to-julia/55996/1 "2021-02-25T09:21:52Z")

</div>

I saw that the MatrixEquations.jl package was released to reproduce some solvers related to Riccati and Lyapunov equations, linked with the SLICOT library.

Are there any plans to link the Toeplitz specialized solvers from SLICOT to a julia package ?

---

<div class="post-metadata">

### Author: ![baggepinnen](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/baggepinnen/32/693_2.png) [@baggepinnen](https://discourse.julialang.org/u/baggepinnen)
#### Post date: [February 25, 2021, 10:05am UTC](https://discourse.julialang.org/t/plans-to-port-slicot-toeplitz-solvers-to-julia/55996/2 "2021-02-25T10:05:09Z")

</div>

Does this package offer what you need?

> **[GitHub - JuliaLinearAlgebra/ToeplitzMatrices.jl: Fast matrix multiplication...](https://github.com/JuliaLinearAlgebra/ToeplitzMatrices.jl)**
>
> Fast matrix multiplication and division for Toeplitz matrices in Julia - GitHub - JuliaLinearAlgebra/ToeplitzMatrices.jl: Fast matrix multiplication and division for Toeplitz matrices in Julia

---

<div class="post-metadata">

### Author: ![BambOoxX](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/bambooxx/32/22179_2.png) [@BambOoxX](https://discourse.julialang.org/u/BambOoxX)
#### Post date: [February 25, 2021, 10:20am UTC](https://discourse.julialang.org/t/plans-to-port-slicot-toeplitz-solvers-to-julia/55996/3 "2021-02-25T10:20:41Z")

</div>

I am currently using this package to form some toeplitz matrices, but as far as I understand the code, there are also some solvers shipped, yet they do not seem to call SLICOT.  
I am particulrly interested in the solution of block Toeplitz matrices problems as with e.g [MB02ID](http://slicot.org/objects/software/shared/doc/MB02ID.html) from SLICOT, but I did not see anything about block Toeplitz matrices in ToeplitzMatrices.jl

---

<div class="post-metadata">

### Author: ![baggepinnen](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/baggepinnen/32/693_2.png) [@baggepinnen](https://discourse.julialang.org/u/baggepinnen)
#### Post date: [February 25, 2021, 11:19am UTC](https://discourse.julialang.org/t/plans-to-port-slicot-toeplitz-solvers-to-julia/55996/4 "2021-02-25T11:19:25Z")

</div>

I’m not aware of any package targeting this particular matrix structure, but

> **[GitHub - JuliaLinearAlgebra/BlockBandedMatrices.jl: A Julia package for...](https://github.com/JuliaLinearAlgebra/BlockBandedMatrices.jl)**
>
> A Julia package for representing block-banded matrices and banded-block-banded matrices - GitHub - JuliaLinearAlgebra/BlockBandedMatrices.jl: A Julia package for representing block-banded matrices ...

and the other packages in JuliaMatrices might be able to solve this problem efficiently, maybe by composing a toeplitz matrix with a block array etc.

---

<div class="post-metadata">

### Author: ![BambOoxX](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/bambooxx/32/22179_2.png) [@BambOoxX](https://discourse.julialang.org/u/BambOoxX)
#### Post date: [February 26, 2021, 2:33pm UTC](https://discourse.julialang.org/t/plans-to-port-slicot-toeplitz-solvers-to-julia/55996/5 "2021-02-26T14:33:50Z")

</div>

Do you mean that with julia the composition of matrix types may lead to performance improvements ? As in e.g. `BlockDiagonal(Toeplitz(Array)` will optimize output computation using both blokc matrix and Toepitz matrices ?

---

<div class="post-metadata">

### Author: ![baggepinnen](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/baggepinnen/32/693_2.png) [@baggepinnen](https://discourse.julialang.org/u/baggepinnen)
#### Post date: [February 26, 2021, 6:35pm UTC](https://discourse.julialang.org/t/plans-to-port-slicot-toeplitz-solvers-to-julia/55996/6 "2021-02-26T18:35:53Z")

</div>

I am not familiar with the internals of BlockArrays, but if their backsolve is implemented in terms of backsolve on the blocks (I don’t know if this is feasible) then yes, the composition will be performant.  
They do have a bandedBlockBandedMatrix, so they might have some custom solvers for that.
