La siguiente ecuación matricial en Σ - para matrices B y C dadas - aparece en mi trabajo como una caracterización de una matriz de covarianza. He aprendido que esta ecuación es conocida, en particular en la teoría del control de tiempo continuo, como la ecuación de Lyapunov , y que existen varios algoritmos bien conocidos para resolverla que explotan la naturaleza especial de esta ecuación lineal.
Al buscar en Google también aprendí que existen implementaciones de Matlab y Fortran. He encontrado SLICOT y RECSY. Sin embargo, debido a problemas de licencia, se detuvo el acceso a la fuente SLICOT.
La mayor parte de mi trabajo se implementa en R, y como no he podido encontrar una interfaz R para un solucionador, considero escribir una. Mi pregunta es si SLICOT es la mejor biblioteca Fortran (o C) disponible con una implementación de un solucionador de la ecuación de Lyapunov. También estoy interesado en implementaciones que puedan manejar grandes matrices dispersas .