- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hello
In the attached example (creating dense matrix sv - > transforming into sparse matrix v -> multiplying with vector u) there is a problem with the mkl_dcsrmv routine. I.e. the value stored in the product vector y is wrong.
I would be very happy if someone could tell me what I do wrong.
Heres the code:
test.cpp
Kind regards
Wolfram
In the attached example (creating dense matrix sv - > transforming into sparse matrix v -> multiplying with vector u) there is a problem with the mkl_dcsrmv routine. I.e. the value stored in the product vector y is wrong.
I would be very happy if someone could tell me what I do wrong.
Heres the code:
test.cpp
Kind regards
Wolfram
test.cpp
(Virus scan in progress ...)
Link Copied
1 Reply
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
I found the reason: as the solution which is solved is y = alpha*A*x + beta*y (and not just y=A*x) it is necessary to initialize y. In the given case all y values hasve to be set to 0 before calling to mkl_dcsrmv.
Kind regards
Wolfram
Kind regards
Wolfram

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