Home My Page Projects PaStiX
Summary Activity Forums Lists Docs News Files

Forum: open-discussion

Monitor Forum | Start New Thread Start New Thread
RE: How to get the D in LDL' after factorization [ Reply ]
By: Nobody on 2015-05-25 00:54
[forum:149583]
Thank you!

RE: How to get the D in LDL' after factorization [ Reply ]
By: Théophile Terraz on 2015-05-20 12:54
[forum:149575]

pastix_release_bce84c2.tar.bz2 (8) downloads
Hello,

Here is a new release of PaStiX. You can now get the D of the LDLt factorization by setting iparm[IPARM_INERTIA] to API_YES, then run only the factorization step of PaStiX. The diagonal will be returned in the vector X, don't forget to allocate it before the PaStiX call.

The step 4 of the examples "step-by-step.c" and "step-by-step_dist.c" show you how to get the diagonal of the factorized matrix. It also works with LLt and LU factorizations.

We would be very interested in your feedback.

Regards,

Théophile

RE: How to get the D in LDL' after factorization [ Reply ]
By: Shaojun Fu on 2015-05-13 01:26
[forum:149558]
Thanks!

Looking forward to your patch.

Shaojun

RE: How to get the D in LDL' after factorization [ Reply ]
By: Pierre Ramet on 2015-05-12 06:07
[forum:149556]
Dear Shaojun,
it is not yet possible to get the D matrix of the LDLt factorization through the current interface. But it should be easy to add this functionality with a new iparm option, at the end of the factorization step, if the vector X is allocated by the user. Of course, we will take care of reordering and data distribution (in the same time, we can also retrieve the diagonal of L from LLt and LU factorization).
We will provide a patch that you can try in few days, then it could be committed in the develop branch of PaStiX.
Regard,
Pierre.

How to get the D in LDL' after factorization [ Reply ]
By: Shaojun Fu on 2015-05-11 11:06
[forum:149555]
Hello!

I want to get the D in LDL' after factorization.
Is there an interface? Or where is it?
It seems to be inside pastix_data.

Your help would be very appriciated!

Shaojun