Segmentation fault from dcsradd

Segmentation fault from dcsradd

Аватар пользователя Duke_K
Duke_K
чт, 28/06/2012 - 07:42

Hello. I get segmentation fault after this calling

call mkl_dcsradd('T', 1, sort, size(imat1)-1, maxval(jmat1), mat1, jmat1, imat1, bb, mat3, jmat3, imat3, Wv, jWv, iWv, size(Wv), info).

sort isn't set, bb=-1.d0. All arrays are allocatable, jWv and Wv are not allocated. Matrices mat1 and mat3 are in the attached files*. MKL version is 12.0.3. All parameters and matrices seem fine to me, and I can't figure what's it all about. Moreover, on MKL version 10.2.7 all work properly.

Can you, please, help me with some advice?

*I attached them to my next comment
notes about matrices storage:
order of arrays - ja, ia, a
matrices are 15x15. mat1 has 15 nonzero elements, mat3 has 225 nonzero elements, most of them are zeros.

RSS Вверх
6 сообщений / 0 новое
Последнее сообщение
Пожалуйста, обратитесь к странице Уведомление об оптимизации для более подробной информации относительно производительности и оптимизации в программных продуктах компании Intel.
Аватар пользователя Chao Y (Intel)
Chao Y (Intel)
чт, 28/06/2012 - 19:56

Hello,

Do you have any code that could show this problem? That will help to look at this problem.

Thanks,
Chao

Вверх
Аватар пользователя Duke_K
Duke_K
вт, 03/07/2012 - 01:17

Thank you for answer.

The thing is that I get this error just like that - read matrices from files and run the routine. I've made small separate program, and attached it to post along with updated matrices. I compile this program by "ifort -mkl segfault.f90", maybe it's not exactly the right way, but before that case I have no problems.

Also I'd like to add, that version of MKL, which crushes on the case, is 10.3.3, not 12.0.3.

Вложения: 

ВложениеРазмер
Скачать mat1.dat743 байта
Скачать mat3.dat8.21 КБ
Скачать segfault.f90777 байтов
Вверх
Аватар пользователя Chao Y (Intel)
Chao Y (Intel)
пн, 09/07/2012 - 20:35

Hello,

For the attached code, it computed withthe request=1:

According to the manual:
If request=1, the routine computes only values of the array ic of length m + 1, the memory for this array must be allocated beforehand.

In thecode, resmat, jresmat do not allocate the memory. Also, the sort is undefined, it should be a integer value.

Thanks,
Chao

Вверх
Аватар пользователя Duke_K
Duke_K
пн, 03/09/2012 - 01:05

Thank you for your answer, Chao.
But I didn't understand your suggestions quite well, will you, please, explain it to me one more time?
If request=1, array ic must be allocated, so in my code. After this calling, we get number of nonzero elements in matrix (nnz), so we can allocate arrays c and jc. But we may not allocate them, am I right? For example, if we just want to know nnz. And anyway, thing is that segmentation fault prevents procedure with request=1 from finishing, and this is my problem.
Thanks in advance for your answer and your time.

Вверх
Аватар пользователя Chao Y (Intel)
Chao Y (Intel)
вт, 11/09/2012 - 19:54

Hello,
>But we may not allocate them, am I right? For example, if we just want to know nnz. And anyway, thing is that segmentation fault prevents procedure with request=1 from finishing,

yes, those arrays are not allocated in your program, and it needs to allocated before calling the function.

Thanks,
Chao

Вверх

Зарегистрируйтесь, чтобы оставить комментарий.