http://www.gnu.org/software/gsl/manual/html_node/
ю, что есть несколько оптимизированных алгоритмов для всех видов матричных разложений (QR-разложение, SVD, ...), умножения и тому подобное. Тем не менее, я не мог найти хороший обзор. Для C ++ есть довольно полезная информацияэтот вопрос, но я ищу эти вещи в C.