how to convert a matrix to a diagonally dominant matrix using pivoting in Matlab
3 views (last 30 days)
Show older comments
Hi I am trying to solve a linear system of the following type: A*x=b, where A is the coefficient matrix, x is the vectors of unknowns and b is the vector of solution. The coefficient matrix (A) is a n-by-n sparse matrix, with even zeros in the diagonal. In order to solve this system in an accurate way I am using an iterative method in Matlab called bicgstab (Biconjugate gradients stabilized method ). This coefficient matrix (A) has a det(A)=-4.1548e-05 and a rcond(A)= 1.1331e-04. Therefore the matrix is ill-conditioned. I first try to perform a scaling and the results where: det(A)= -1.2612e+135 but the rcond(A)=5.0808e-07... Therefore the matrix is still ill-conditioned... I verify and the sum of all absolute value of the non-diagonal elements where 163.60 and the sum of all absolute value of the diagonal elements where 32.49... Therefore the matrix of coefficient is not diagonally dominant and will not converge using my function bicgstab...
I am looking for someone that can help me with performing a pivoting to the coefficient matrix (A) so it can be diagonally dominant. Or any advice to solve this problem....
Thanks for the help.
0 Comments
Answers (0)
See Also
Categories
Find more on Linear Algebra in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!