Superl LU MT Cant получить решение X
Я пытался реализовать библиотеку SuperLU MT в моем проекте, но я не понимаю, как она работает.
В документации сказано, что драйвер возвращает решение в плотной матрице "B", которую вы дали через аргументы, но когда я попытался напечатать эту матрицу B, я получил вектор единиц.
Я печатаю матрицу B в файле pzlinsol.c, из примеров.
pzgssv(nprocs, &A, perm_c, perm_r, &L, &U, &B, &info);
zPrint_Dense_Matrix(&B);
Я решил эту проблему с помощью scipy, и результат не является вектором единиц, я также попробовал еще несколько примеров, но ничего.
Должно быть что-то, что я делаю неправильно!
Спасибо!