Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- BUG: Buffer Overrun in mkl_dcsrcsc

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

bamberg

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

06-17-2010
02:36 AM

60 Views

BUG: Buffer Overrun in mkl_dcsrcsc

This was verified with mkl10.2.5.0.35 and mkl 10.3.0.025.

Best regards,

Brian

Brian

This is the output I get for the test program

Link Copied

1 Reply

Ying_H_Intel

Employee

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

06-17-2010
08:02 PM

60 Views

Hi Brian,

Thanks a lot for the small test case! We can reproduce the problem with it. I get same result as yours.

But the result may be expected if we consider the fact of the routinetakethe sparse matrix as square matrices by default.(Thisinfo seems be missed mkldocumentation. I will ask ourdoc developer to add such claim.).

For example, the parameter m INTEGER. Dimension of the matrix A.We give it 4 as input,

as a result, it thought the sparse matrix is 4x4.Thus, theresult matrix should be

11 ** *

*12**

1314 * *

*15 * *

The colIndex: 1 32 3 4 is explained as row in CSC format ( Valume, Row, ColumnIndex)

Value: 11, 13, 12, 14, 15

The out_colStart will be explained as column Index in CSC . It is array with m+1=5 elements. Thusit is 1, 3, 6, 6, 6. according to the below definition.

ia1 INTEGER. Array of length m + 1, containing indices of elements in the array acsc, such that ia(I) is

the index in the array acsc of the first non-zero element ia1 from thecolumn I. The value of the last element ia(m + 1)

is equal to the number of non-zeros plus one.

(It is ok to explainthem as ( Value,Column,RowIndex) inCSCifthe matrix is transposed in your case)

^{ThanksYing }

ia1 INTEGER. Array of length m + 1, containing indices of elements in the array acsc, such that ia(I) is

the index in the array acsc of the first non-zero element ia1 from thecolumn I. The value of the last element ia(m + 1)

is equal to the number of non-zeros plus one.

(It is ok to explainthem as ( Value,Column,RowIndex) inCSCifthe matrix is transposed in your case)

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

For more complete information about compiler optimizations, see our Optimization Notice.