[Top][All Lists]

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

[Help-gsl] finite element assembly and calling UMFPACK for linear system

From: Brijesh Upadhaya
Subject: [Help-gsl] finite element assembly and calling UMFPACK for linear system solution
Date: Wed, 19 Jun 2019 18:15:24 +0300

Hi everyone,

I am working on a finite element problem and wanted to ask if anyone of you
have tried using UMFPACK for the linear algebra. I have following findings
when using sparse matrix module from gsl-2.5.

1. In CCS the row indices are not sorted so UMFPACK gives error (see
attached example)
2. Adding boundary conditions were also difficult once compressed, so I
couldn't utilize gsl_spmatrix_add() to add mass matrix with stiffness
matrix. Instead I had to do it with triplet and calling the get/set
routines. Could only add boundary condition at triplet level and compressed
at the end and passed the pointers to UMFPACK routines.
3. Another thing (not related to gsl) was that UMFPACK don't use (size_t
*), instead it has (int *) and (long *) for column pointers and row
indices. I tried to cast (size_t *) to (long *) to utilize the pointers of
gsl_spmatrix object.

I am bit curious to know if anyone have a better experience. I wanted to
use direct solver instead of iterative solver.

King Regards,

Attachment: test_lsolver.c
Description: Text Data

reply via email to

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