gaussian elimination without pivoting

Open discussion regarding features, bugs, issues, vendors, etc.

gaussian elimination without pivoting

Postby steven82 » Mon Aug 22, 2005 8:40 am

Hello,

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?

Thx,

Steven
steven82
 
Posts: 1
Joined: Mon Aug 22, 2005 8:33 am

Postby sven » Mon Aug 22, 2005 8:56 am

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.
sven
 
Posts: 144
Joined: Wed Dec 22, 2004 4:28 am


Return to User Discussion

Who is online

Users browsing this forum: Yahoo [Bot] and 4 guests