Извините, если этот вопрос слишком неоднозначный или наивный
Я реализую некоторый алгоритм деформации сетки на основе оптимизации
Последний шаг оптимизации — решение большой разреженной матрицы с помощью разреженного решателя.
Я использую 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 верно
Задача ещё не решена.