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

Detecting ill-condition matrix in PARDISO (instability in FE model)

Alemdar__Bulent
Beginner
431 Views

Hi

We are planning to use PARDISO in out finite element analysis package. Usually the matrix to be solved is positive-deifinite symmetric matrix. But sometimes, wemight haveill-conditioned matrix (this might be due to mesh discretization error or any unsupressed rigid body mode in our finite element model). If this is the case, we try to extract node (or the position of diagonal term that leads to negative pivoting) and then we check our FE model at this node to make any correction in our model.

I tried a few examples in PARDISO, and once it fails to solve the system but I coudn't find any more information returned by the solver.

Is there any way to get more info from the solver like node or diagonal number if PARDISO fails to solve an ill-conditioned matrix?

best regards

bulent

0 Kudos
1 Reply
Sergey_K_Intel1
Employee
431 Views

Dear Bulent,

At the present PARDISO doesn't returm the position of zero diagonalelement for ill-conditioned or singular matrices. We have plans to implement the feature in one of future releases.

All the best

Sergey

0 Kudos
Reply