Как проверить dpotrf

Question

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

Должна ли эта матрица иметь специальный шаблон? Могу ли я просто сгенерировать простую двойную матрицу NxN?

Answer 1

из dpotrf.f , DPOTRF вычисляет факторизацию Холецкого реальной симметричной положительно определенной матрицы A.

Следовательно, ожидается, что входное значение будет положительно определенным и симметричным.

Кроме того, dpotrf ожидает, что входные данные будут массивом, т. Е. Ваш входной сигнал равен (N * N) x1. Я обычно изменяю вывод в NxN.