Page 1 of 1

gaussian elimination without pivoting

PostPosted: Mon Aug 22, 2005 8:40 am
by steven82

I've already used the DGESV function to perform gaussian elimination of a general linear system. This function uses LU decomposition with partial pivoting and row interchanges. Now I want to solve my system without any pivoting to compare it to the solution with pivoting. Now my question is which function can I use to solve my system using gaussian elimination without pivoting?



PostPosted: Mon Aug 22, 2005 8:56 am
by sven
Dear Steven,

There is no routine in LAPACK to perform Gaussian elimination without pivoting on a general matrix. It is not available because, as I am sure you are aware, that would not be numerically stable unless the matrix had some special property.

If you wanted to implement such a method for illustrative purposes it would not be hard to take DGETRF and remove the pivoting. (But please change the name of the routine if you do that!)

Best wishes,

Sven Hammarling.