topic Re: dss routines - sequence in Intel® oneAPI Math Kernel Library

dss routines - sequence

Harald31415 — Wed, 19 Nov 2025 12:06:08 GMT

Hi everyone,
I am testing the dss routines. Everything seems fine: I can call dss_create, dss_define_structure, dss_reorder, dss_factor_real and dss_slove_real in that sequence and get the expected result. I can also repeat the sequence.

However, if I again call define_structure after factor_real and before solve_real I get error -13 ( MKL_DSS_TERM_LVL_ERR )

Did I miss something ? Is this not allowed ? Is there a way to undo the results of dss_reorder and dss_factor_real ?
Any help appreciated

Harald

Re: dss routines - sequence

Harald31415 — Wed, 19 Nov 2025 18:53:10 GMT

I solved it pragmatically: I start from scratch before each define_structure. Yet I am interested in the concept.

Re: dss routines - sequence

mecej4O — Sun, 23 Nov 2025 13:33:39 GMT

Please read the documentation, which lists the order in which the individual DSS routines should be called. Usually, the structure of the matrix does not change during a calculation. Often, even the matrix is constant but the right hand side vector changes. If so, the matrix may be factorized just once, and only the solution plase needs to be repeated for each right hand side. The factorization is the most time-consuming part, so it makes sense to do it as few times as possible -- just once if that is correct to do.

Re: dss routines - sequence

Harald31415 — Sun, 23 Nov 2025 14:19:46 GMT

thank you mecey40. I did read the documentation; my question was, if the same handle can be used for different matrices, even before finishing the cycle.

Re: dss routines - sequence

Fengrui — Mon, 24 Nov 2025 21:40:48 GMT

No, it can not. Internal results are stored in the handle. Using it for different matrices before finishing the cycle leads to undefined behavior.

Re: dss routines - sequence

Harald31415 — Mon, 24 Nov 2025 21:48:58 GMT

Thank you. That is a clear answer. I thought, there is some clear or free or release command which I missed.