Использование научной библиотеки GNU для поиска ядра матрицы
Учитывая gsl_matrix * Объект, который является M-от-N матрица, какой самый простой способ найти ядро А? Я пытался использовать разложение по сингулярным значениям (в частности, gsl_linalg_SV_decomp метод), но это не работает, когда M<N, что иногда имеет место для моих матриц.
Решение
Самый простой способ (если вы ограничены этой библиотекой) состоит в том, чтобы заполнить вашу матрицу рядами нулей до N x N, которые затем вы можете передать своей функции. Заполненная матрица будет иметь то же пустое пространство.
Другие решения