[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Help-gsl] problem in using LU decomp

From: Serle Shuman
Subject: [Help-gsl] problem in using LU decomp
Date: Wed, 20 May 2009 22:38:01 +0200



Can anyone see what I am doing wrong. I seem to get and infinite v[0] value
when it should be -0.0000015531.


            int perm_sign = 0;

            gsl_matrix* C = gsl_matrix_calloc(2, 2);

            gsl_vector* b = gsl_vector_calloc(2);

            gsl_vector* v = gsl_vector_calloc(2);

            gsl_permutation* p = gsl_permutation_calloc(2);


            gsl_matrix_set(C, 0, 0,  1.0);

            gsl_vector_set(b, 0, -0.0000015531);

            gsl_linalg_LU_decomp(C, p, &perm_sign);

            gsl_linalg_LU_solve(C, p, b, v);

            for (int l = 0; l < n; l++) {

                  if (!gsl_finite(gsl_vector_get(v, l))) {

                        // we have a problem with v





Also can one use LU decomposition solver on a non-square matrix or does one
need to pad it up with zeros.


Thanks in advance


reply via email to

[Prev in Thread] Current Thread [Next in Thread]