Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.

tridiagonal matrix solution

LRaim
New Contributor I
874 Views

I have tried to understand if MKL provides a faster method for the solution of simple linear tri-diagonal systems which can be easily solved with simple Gaussian elimination.
​I am using PARDISO in different cases.

 

 

 

 

0 Kudos
3 Replies
mecej4
Honored Contributor III
874 Views

Pardiso is not the best choice for solving banded linear equations. Please consider using the MKL/Lapack banded equations solvers ?gbslv with bandwidths kl = ku = 1, see https://software.intel.com/en-us/node/468882 .

0 Kudos
Gennady_F_Intel
Moderator
874 Views

you may also try to evaluate ?gtsv ( Computes the solution to the system of linear equations with a tridiagonal coefficient matrix A and multiple right-hand sides )

0 Kudos
Royi
Novice
874 Views

@Gennady,

Documentation of ?gtsv states:

The routine uses Gaussian elimination with partial pivoting.

How come it doesn't use Thomas Algorithm?
How does it compare to vanilla Thomas Algorithm in performance?

0 Kudos
Reply