Как напечатать нижнюю матрицу в холески с помощью gsl? - PullRequest
Новая
       8

Как напечатать нижнюю матрицу в холески с помощью gsl?

0 голосов
/ 17 мая 2018

Я использую gsl с c ++ для решения линейной системы с помощью метода Холецкого, все работает нормально, но я хочу получить и напечатать нижнюю матрицу 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;
}

1 Ответ

0 голосов
/ 17 августа 2018

Нижняя матрица L перезаписывается в исходной матрице A. Нижняя треугольная часть матрицы A, хранит L.

...