Why doesn't the sparse command accept uintx as inputs for the location vectors in the long form of the command?

Question

Nicholas Dinsmore 2013-5-14

2
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75776-why-doesn-t-the-sparse-command-accept-uintx-as-inputs-for-the-location-vectors-in-the-long-form-of

When I use the sparse command in the S = sparse(i,j,s,m,n) form, with vectors of type uintx for i & j, I get an error saying that type isn't supported. But since those are indexes would the sparse command have to convert to an integer format to use them? Now to be clear I am not trying to specify the contents ie the type of the sparse martix as an integer, the documentation clearly says only double and bool

As a note of the application. I have found that my function that generates a large sparse matrix, 25-50% faster if I change a single line of code to make the index arrays integers, which is understandable since integers=faster and smaller. I would have thought this would be supported. I could typecast the arrays but that seems like a redundant step.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Matt J 2013-5-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75776-why-doesn-t-the-sparse-command-accept-uintx-as-inputs-for-the-location-vectors-in-the-long-form-of#answer_85448

编辑：Matt J 2013-5-14

SPARSE uses doubles to store index data. Why it doesn't pre-cast uintx input to doubles for you is a fair enough question, but the index data will eventually need to be converted to doubles one way or another. It is also a fair question why after 10 years only double and logical matrix entries are supported. Probably because SPARSE originated as 3rd party code, and now it's hard to change.

7 个评论
显示 5更早的评论隐藏 5更早的评论

James Tursa 2013-5-15

编辑：James Tursa 2013-5-15

在 MATLAB Online 中打开

Formula for storage (not including overhead) of a sparse double matrix S assuming k bytes per integer for the index data is:

(8+k)*nnz(S) + k*(size(S,2)+1)

That is one 8-byte double and one k-byte row integer for each non-zero element, and then (number_of_columns+1) k-byte integers for the column data. So you should be able to deduce from this formula when compared to an actual sparse matrix what size integers are used for the sparse index data. After checking on Win32 and Win64 systems here are my conclusions:

Win32: 4-byte integer for index storage, 2^31-1 is largest index allowed

Win64: 8-byte integer for index storage, 2^48-1 is largest index allowed

As for the signed vs unsigned question, the mex header files show it as either an int (signed) or size_t (unsigned) depending on some settings. I don't know what all influences the settings. This only really matters if MATLAB allows indexes near the upper range into the "sign" bit (doesn't appear to based on tests) and in mex programming if one is subtracting indexes for some reason and possibly generating a negative result (which would get wrapped around to a large positive number if size_t is the underlying type).

Nicholas Dinsmore 2013-5-20

Guys, Thank you for your response. I do think that it is important for sparse to accept the int and uint as input. For a large sparse matrix it was surprising to me how much faster I could make the generation code just by using int instead of doubles of the indexing arithmetic(to see the benefit you need to give the JIT an opportunity to optimize). And because of the strength of the type system in matlab, it was very few lines of code to make things ints.

Matt J 2013-5-20

For a large sparse matrix it was surprising to me how much faster I could make the generation code just by using int instead of doubles of the indexing arithmetic(to see the benefit you need to give the JIT an opportunity to optimize)

You mentioned this in your original posting as well, but to me at least it's not clear what you mean by the "generation code". If the matrix is generated by calling sparse() and if sparse() responds to integer inputs with an error, then how can it be "faster".

请先登录，再进行评论。

Why doesn't the sparse command accept uintx as inputs for the location vectors in the long form of the command?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论
显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Why doesn't the sparse command accept uintx as inputs for the location vectors in the long form of the command?

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论 显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

7 个评论
显示 5更早的评论隐藏 5更早的评论