Why the sparse pattern of a sparse matrix changed

Question

0 个投票

I'm building a sparse matrix through the mex function. Both the Ir and Jc properties are appointed in the mex function A. Later, this sparse matrix is returned to matlab, and is further transferred to another mex function B together with other parameters to assemble a large matrix.

For one test case, part of the Jc property of the created sparse matrix in the mex function A is as follows:

368

392

416

440

464

488

512

But when this sparse matrix is further transferred to mex function B, this part is altered and is as follows:

368

392

416

440

463

487

511

Actually, after this sparse matrix is established, I did no alternation about it. I want to know what may case this.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Bruno Luong 2021-4-11

You MUST NOT store 0 elements internally. Use spok I indicate bellow to check if your mex-build result is valid.

Rylan 2021-4-11

(9,22) -0.0000

(10,22) 0

(11,22) 0.0000

(12,22) -0.0000

(13,22) -0.2500

(14,22) 0.5000

It seems problem was caused by (10,22), and use of spok indicate that "0 jumbled row indices, 3 explicit zeros". To remove this problem, may be a vary small number, for example, 1e-16, should be added to all the data member.

Thanks for your suggestion!

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question