通过M<N的GSL计算矩阵的零空间会产生错误

遵循gsl中奇异值分解(SVD)函数的文档,例如A=USVt

对于秩亏矩阵,a的零空间由对应于零奇异值的V列给出

int gsl_linalg_SV_decomp(gsl_matrix * A, gsl_matrix * V, gsl_vector * S, gsl_vector * work)

语法如下:

#include <gsl/gsl_linalg.h>
#include <gsl/gsl_blas.h>
#include <stdio.h>
#include <stdlib.h>
...

int main(int argc, char **argv){
    ...
    gsl_linalg_SV_decomp(A,V,S,work);

return 0;
}

gsl: svd.c:60: ERROR: svd of MxN matrix, M<N, is not implemented
Default GSL error handler invoked

事实上,A有M<N

你知道另一个库会计算M<矩阵的零空间吗;Ngsl有解决方法吗?