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

tridiagonal matrix solution

LRaim
New Contributor I
‎04-15-2016 02:53 AM
1,963 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

Link Copied

3 Replies
mecej4
Honored Contributor III
‎04-15-2016 12:42 PM
1,963 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
Copy link
Gennady_F_Intel
Moderator
‎04-15-2016 11:52 PM
1,963 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
Copy link
Royi
Novice
‎08-11-2018 02:37 PM
1,963 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
Copy link
Reply

Community support is provided Monday to Friday. Other contact methods are available here.

Intel does not verify all solutions, including but not limited to any file transfers that may appear in this community. Accordingly, Intel disclaims all express and implied warranties, including without limitation, the implied warranties of merchantability, fitness for a particular purpose, and non-infringement, as well as any warranty arising from course of performance, course of dealing, or usage in trade.

For more complete information about compiler optimizations, see our Optimization Notice.