Estoy resolviendo un problema físico usando un esquema numérico implícito. Esto me lleva a resolver una ecuación lineal con matriz tridiagonal. He codificado este algoritmo de Wikipedia. Me pregunto si hay una biblioteca eficiente que permita resolver este tipo de ecuaciones de manera optimizada. Una nota importante es que la matriz en sí misma solo cambia cuando los parámetros del sistema están cambiando, por lo que tuve la oportunidad de calcular previamente algunos pasos del algoritmo para obtener una buena bonificación de rendimiento. Estoy usando C ++.