[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,
Brijesh

**
**`test_lsolver.c`

*Description:* Text Data

**[Help-gsl] finite element assembly and calling UMFPACK for linear system solution**,
*Brijesh Upadhaya* **<=**