Estoy buscando una biblioteca C
para resolver la ecuación de matriz lineal y, si es posible, no lineal de la forma Ax = b
. Para mí es importante que los paquetes no sean demasiado grandes y gratuitos. La velocidad no importa tanto como la simplicidad y la característica dispersa del almacenamiento matricial. Además, debería ser capaz de paralelizar los cálculos. Como soy bastante nuevo en el campo blas/lapack/... del álgebra lineal numérica, sería genial que contuviera una buena documentación, posiblemente con ejemplos. ¿Hay algún paquete que puedas recomendar? Estoy un poco abrumado por los resultados de Google y tampoco sé los nombres de las rutinas que debería buscar.Resolvador de ecuaciones lineales y no lineales de matriz dispersa
Acabo de encontrar SuperLu que parece tener una documentación muy buena, incluso con fragmentos de código. Está escrito en C y requiere CBLAS.
¿Cómo es Ax = b a * ecuación matricial no lineal *? – talonmies
No es lineal si A es una función de x. Se expresa mejor como A (x) * dx = db y x (i + 1) = x (i) + dx para capturar la carga incremental no lineal. Resuelve para dx y actualiza. – duffymo
SuperLU y Umfpack son solucionadores directos para problemas lineales. Si, como dices, tienes un problema no lineal, necesitarás algo más. – talonmies