reducción de filas de lapack

Question

Estoy intentando escribir una función que produce una solución única a un sistema de ecuaciones subrepresentados (por ejemplo, la matriz que describe el sistema es más ancho que alto). Con el fin de hacer esto, he estado buscando en la documentación LAPACK una manera de reducir fila de una matriz a su forma reducida escalonada, similar a la rref () en ambas calculadoras Mathematica y TI. El más cercano que encontré fue http: / /software.intel.com/en-us/forums/intel-math-kernel-library/topic/53107/ este diminuto hilo. Este hilo, sin embargo, parece implicar que simplemente tomando la "U" matriz triangular superior (y dividiendo cada fila por la diagonal) es la misma que la forma escalonada reducida de una matriz, que no creo ser el caso. Podría codificar hasta rref () yo mismo, pero no creo que pudiera lograr el rendimiento LAPACK es famoso.

1) ¿Hay una mejor manera de conseguir simplemente cualquier solución específica a un sistema subrepresentados?
2) Si no, hay una manera para LAPACK a la fila-a reducir una matriz?

Gracias!

Answer 1

Un método frecuentemente utilizado para esta es la solución de mínimos cuadrados, véase sgelsx de lapack