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

Showing results for

- Intel Community
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- unit diagnal elements in a upper triangular matrix

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

anjuta

Beginner

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

03-22-2010
02:51 PM

79 Views

only unit diagnal elements. I found that I still need to include the diagnal elements in the matrix storage vectors passed to mkl_dcscsv, otherwise the result will be wrong. I'm wondering if there's any specific reason for this? I mean isn't it intuitive not to include the unit diagnal elements?

Thanks

Link Copied

Accepted Solutions

Gennady_F_Intel

Moderator

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

03-25-2010
06:14 AM

79 Views

Thanks for the example. all right with this example.

The unit diagonal elements can be omitted accordingly the description of mkl_dcscsv and therefore that's the bug.This issue has been submitted to our internal development tracking database for further investigation, we will inform you once a new update becomes available.

--Gennady Fedorov

7 Replies

ArturGuzik

Valued Contributor I

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

03-22-2010
08:40 PM

79 Views

I believe you might have overlooked the description of the input parameters. Take a look at

anjuta

Beginner

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

03-23-2010
11:42 AM

79 Views

from my testing, I think it's not.

ArturGuzik

Valued Contributor I

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

03-23-2010
05:56 PM

79 Views

A.

anjuta

Beginner

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

03-24-2010
11:53 AM

79 Views

#include

#include "mkl.h"

#include "mkl_types.h"

#include "mkl_spblas.h"

int main(int argc, char **argv)

{

/*double values[] = {1.0, 1.0, -0.000999, 1.0, 0.2, -0.2, 1.0, -1.0, 1.0};

MKL_INT rows[] = {0, 1, 0, 2, 0, 2, 3, 1, 4};

MKL_INT colIndex[] = {0, 1, 2, 4, 7, 9};*/

double values[] = {-0.000999, 0.2, -0.2, -1.0};

MKL_INT rows[] = {0, 0, 2, 1};

MKL_INT colIndex[] = {0, 0, 0, 1, 3, 4};

double rhs[5] = {0, 0, 0, 0, 0.001998};

double sol[5] = {1.0, 1.0, 1.0, 1.0, 1.0};

char transa = 'n';

char matdescra[6];

MKL_INT m = 5;

double alpha = 1.0;

matdescra[0] = 't';

matdescra[1] = 'u';

matdescra[2] = 'u';

matdescra[3] = 'c';

mkl_dcscsv(&transa, &m, α, matdescra, values, rows, colIndex, &(colIndex[1]), rhs, sol);

printf("sol = %.7e %.7e %.7e %.7e %.7e\n", sol[0], sol[1], sol[2], sol[3], sol[4]);

return 0;

}

The correct results would be:

0, 0.001998, 0, 0, 0.001998

but with the above code, I'm getting

0, 0, 0, 0, 0.001998

once I put the unit diagonal elements in, the result is correct.

Thanks

Gennady_F_Intel

Moderator

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

03-25-2010
06:14 AM

80 Views

Thanks for the example. all right with this example.

The unit diagonal elements can be omitted accordingly the description of mkl_dcscsv and therefore that's the bug.This issue has been submitted to our internal development tracking database for further investigation, we will inform you once a new update becomes available.

--Gennady Fedorov

Gennady_F_Intel

Moderator

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

04-01-2010
05:18 AM

79 Views

Hello,

just for your informatin:

This issue has been submitted to our internal development tracking database for further implementation,

we will inform you once a new update becomes available.

Here is a bug tracking number for your reference: 200092519

Regards, Gennady

Gennady_F_Intel

Moderator

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

09-01-2010
04:43 AM

79 Views

Dear Customer,

The problem you reported has been fixed
into the latest MKL **10.2 Update 6** version released recently.

You can Download this version from Intel Registration Center: https://registrationcenter.intel.com/regcenter/register.aspx

See announcement about thathere.

Please let us know if the problem is still exist.

--Gennady

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