Open discussion regarding features, bugs, issues, vendors, etc.
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?
- Posts: 1
- Joined: Mon Aug 22, 2005 8:33 am
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!)
- Posts: 144
- Joined: Wed Dec 22, 2004 4:28 am
Return to User Discussion
Who is online
Users browsing this forum: Baidu [Spider], Bing [Bot] and 1 guest