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

Hello,

I don't get why in the attached code, if line 36 I set transa to 't', then I get a segmentation fault. The matrix A is symmetric, so there shouldn't be any other changes to do to make this work ? In its Fortran version, dcsr_multiplication.f, there is no problem when I switch trans line 82.

Thanks in advance for any help.

Link Copied

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

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

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

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

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

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

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

*- 1; j < ic[i + 1] - 1; ++j) std::cout << i + 1 << " " << jc* << " " << c << std::endl;
}[/cpp] You will see that when transa='n', the results are correct and the corresponding CSR is upper triangular, while when transa='t', the CSR is not lower triangular, but since A is upper triangular, A' is lower triangular, hence A'*A' should be lower triangular.
Do you agree with what I am saying ? Do you get incorrect results too ?
Thanks in advance for your help.

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

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

qweasd q. wrote:Thanks for all your help !Another quick question, is it planed to support dcsrmultcsr with arbitrary 'matdescr' (for example symmetric input, unsymmetric output). I'm using this routine to compute A*D where A is symmetric (and so only its lower part is stored) and D is diagonal (but might as well be a general matrix). What I'm doing right now is csradd(csrmultcsr(A, D), csrmultcsr(A', D)) (I don't forget to remove the diagonal terms of csrmultcsr(A', D)), but not only it is making the MKL segfaults (not anymore with your fix though), I don't think this is very efficient.

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

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

hello, pls check the issue with the latest 11.0 update3 and let us know the result.

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

Hello,

It seems Update 3 breaks ARPACK. Because I also use it in my application, I won't update until Update 4. Thanks anyway, I hope everything will be alright in the next update.

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