如何使用gsl打印cholesky中的下矩阵？ 我用GSL和C++来解决Cholesky方法的线性系统，所有的工作都很好，但是我想获得和打印下矩阵L。我在官方文档中找不到任何关于这个问题的东西。这是我的密码： #include <iostream> #include <gsl/gsl_matrix.h> #include <gsl/gsl_linalg.h> using namespace std; int main(void) { double a_data[] = { 4, -1, 1, -1, 4.25, 2.75, 1, 2.75, 3.5}; double b_data[] = {1, 2, 3}; gsl_matrix_view A = gsl_matrix_view_array (a_data, 3, 3); gsl_vector_view b = gsl_vector_view_array (b_data, 3); gsl_permutation * p = gsl_permutation_alloc (3); gsl_vector *x = gsl_vector_alloc (3); gsl_matrix *L = gsl_matrix_alloc (3, 3); gsl_linalg_pcholesky_decomp(& A.matrix, p); gsl_linalg_pcholesky_solve(& A.matrix, p, &b.vector, x); cout<<"x = "<<endl; gsl_vector_fprintf (stdout, x, "%g"); cout<<"p = "<<endl; gsl_permutation_fprintf (stdout, p, "%d"); cout<<endl; cout<<"L = "<<endl; gsl_matrix_fprintf (stdout, L, "%g"); gsl_permutation_free (p); return 0; } #包括 #包括 #包括 使用名称空间std； 内部主（空）{ 双a_数据[]={4，-1,1， -1, 4.25, 2.75, 1, 2.75, 3.5}; 双b_数据[]={1,2,3}； gsl_矩阵视图A=gsl_矩阵视图数组（A_数据，3，3）； gsl_矢量_视图b=gsl_矢量_视图_数组（b_数据，3）； gsl_置换*p=gsl_置换alloc（3）； gsl_向量*x=gsl_向量_alloc（3）； gsl_矩阵*L=gsl_矩阵_alloc（3,3）； gsl_linalg_pcholesky_decomp（&A.matrix，p）； gsl_linalg_pcholesky_solve（&A.矩阵，p，&b.向量，x）； cout

下矩阵L在原始矩阵A中被覆盖。矩阵A的下三角部分存储L