cholmod_factorize ошибка сегментации. Как локализовать ошибку?

Извините, если этот вопрос слишком двусмысленный или наивный

Я реализую алгоритм деформации сетки, основанный на оптимизации. Последний шаг оптимизации - решение большой разреженной матрицы с помощью разреженного решателя.

Я использую cholmod для этой задачи. Но в некоторых случаях я получаю segfault на cholmod_factorize

cholmod_sparse *A = composeASparseCholmodMatrix(........);//Creating my matrix. 
cholmod_factor *L = cholmod_analyze(A, commonCholmodPtr);
cholmod_factorize(A, L, commonCholmodPtr);//receive Segmentation fault on some cases
cholmod_free_sparse(&A, commonCholmodPtr);
//L are used for solving 

Так что мой вопрос в каких случаях cholmod_factorize может вызвать segfault? И как я могу найти источник проблемы?

Я не могу предоставить матрицу здесь, потому что эта матрица очень большая (9000+ строк и столбцов, 45000 ненулевых элементов). Состав матрицы A подходит во многих случаях (8000+ строк / столбцов), алгоритм работает отлично. Я пытаюсь использовать error_handler из cholmod_commmon, но ничего не получил

PS извините за мой английский

UPD Кажется, что программа выходит из строя, тогда L->is_super true

0 ответов

Другие вопросы по тегам