New LR-ADI Lyapunov solver.

Author: andreasvarga

Project.toml

@@ -1,14 +1,25 @@
 name = "MatrixEquations"
 uuid = "99c1a7ee-ab34-5fd5-8076-27c950a045f4"
 authors = ["Andreas Varga <varga.andreas@gmail.com>"]
-version = "2.4.5"
+version = "2.5.0"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 LinearMaps = "7a12625a-238d-50fd-b39a-03d52299707e"
 
+[extras]
+BandedMatrices = "aae01518-5342-5314-be14-df237901396f"
+DoubleFloats = "497a8b3b-efae-58df-a0af-a86822472b78"
+GenericLinearAlgebra = "14197337-ba66-59df-a3e3-ca00e7dcff7a"
+GenericSchur = "c145ed77-6b09-5dd9-b285-bf645a82121e"
+IterativeSolvers = "42fd0dbc-a981-5370-80f2-aaf504508153"
+Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
+SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
+
 [compat]
-DoubleFloats = "1.2"
+BandedMatrices = "1.9"
+DoubleFloats = "1.3"
 GenericLinearAlgebra = "0.3"
 GenericSchur = "0.5"
 IterativeSolvers = "0.9"
@@ -17,16 +28,8 @@ LinearMaps = "3.3"
 Random = "1"
 SparseArrays = "1"
 Test = "1"
-julia = "1.8"
+julia = "1.10"
 
-[extras]
-DoubleFloats = "497a8b3b-efae-58df-a0af-a86822472b78"
-GenericLinearAlgebra = "14197337-ba66-59df-a3e3-ca00e7dcff7a"
-GenericSchur = "c145ed77-6b09-5dd9-b285-bf645a82121e"
-IterativeSolvers = "42fd0dbc-a981-5370-80f2-aaf504508153"
-Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
-SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
-Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["Random", "Test", "GenericSchur", "GenericLinearAlgebra", "DoubleFloats", "IterativeSolvers", "SparseArrays"]
+test = ["Random", "Test", "GenericSchur", "GenericLinearAlgebra", "DoubleFloats", "IterativeSolvers", "SparseArrays","BandedMatrices"]

ReleaseNotes.md

@@ -1,13 +1,17 @@
 # Release Notes
 
+## Version 2.5.0
+Version bump to enforce compatibility with Julia's LTS version.
+New iterative solvers for positive continuous-time Lyapunov equations based on low-rank ADI methods. 
+
 ## Version 2.4.5
 Enhanced Lyapunov solvers `lyapc` and `lyapd` and Sylvester solvers `sylvc` and `sylvd` to efficiently handle symmetric/Hermitian/diagonal inputs.  
 
 ## Version 2.4.4
-Back to using the standard Sylvester solver sylvc with the old wrappers (the *trsyl family of solvers), until issue [#150](https://github.com/JuliaLinearAlgebra/libblastrampoline/issues/150) will be fixed.
+Back to using the standard Sylvester solver `sylvc` with the old wrappers (the `*trsyl` family of solvers), until issue [#150](https://github.com/JuliaLinearAlgebra/libblastrampoline/issues/150) will be fixed.
 
 ## Version 2.4.3
-New wrappers for BLAS Level 3 based LAPACK family *trsylv3 of solvers for Sylvester equation. Updating the standard Sylvester solver sylvc to use the new wrappers (instead the *trsyl family of solvers).    
+New wrappers for BLAS Level 3 based LAPACK family `*trsylv3` of solvers for Sylvester equation. Updating the standard Sylvester solver `sylvc` to use the new wrappers (instead the `*trsyl` family of solvers).    
 
 ## Version 2.4.2
 Version bump to fix type piracy detected by Aqua.  

docs/src/index.md

@@ -61,6 +61,7 @@ The current version of the package includes the following functions:
 | Function | Description |
 | :--- | :--- |
 | **[`lyapci`](@ref)** |Iterative solution of continuous Lyapunov equations|
+| **[`plyapci`](@ref)** |Iterative solution of low-rank non-negative continuous Lyapunov equations|
 | **[`lyapdi`](@ref)** |Iterative solution of discrete Lyapunov equations|
 | **[`tlyapci`](@ref)** |Iterative solution of the continuous T-Lyapunov equations|
 | **[`hlyapci`](@ref)** |Iterative solution of the continuous H-Lyapunov equations|

docs/src/iterative.md

@@ -10,6 +10,7 @@ cgls
 
 ```@docs
 lyapci
+plyapci
 lyapdi
 tlyapci
 hlyapci

src/MatrixEquations.jl

@@ -28,7 +28,7 @@ export sylvckr, sylvdkr, gsylvkr, sylvsyskr, dsylvsyskr,
 export opnorm1, opnorm1est, oprcondest, opsepest
 export lyapop, invlyapop, sylvop, invsylvop, sylvsysop, invsylvsysop, trmatop, eliminationop, duplicationop
 export tulyaplikeop, hulyaplikeop, lyaplikeop, gsylvop
-#export plyapci
+export plyapci
 
 include("meutil.jl")
 include("sylvester.jl")
@@ -39,6 +39,5 @@ include("plyapunov.jl")
 include("meoperators.jl")
 include("condest.jl")
 include("iterative_methods.jl")
-#include("iterative_methods_new.jl")
 
 end