вычислительный детерминант комплексной матрицы fortran90

Мне нужно вычислить детерминант комплексной матрицы, которая является симметричной. Размер матрицы варьируется от 500 * 500 до 2000 * 2000. Есть ли у меня подпрограмма? Кстати, я использую ifort для компиляции.вычислительный детерминант комплексной матрицы fortran90

источник

2013-09-26 someone

Вы уверены, что вам действительно нужно вычислить определитель? Скорее всего, вы хотите получить соответствующее количество, такое как номер условия. В LAPACK существует несколько соответствующих подпрограмм (это входит в библиотеку Intel MKL, которую вы можете иметь, если у вас есть ifort). Согласно частому варианту LAPACK, существуют проблемы стабильности с использованием детерминант http://www.netlib.org/lapack/faq.html#_are_there_routines_in_lapack_to_compute_determinants – DaveP

Самый простой способ - сделать LU-разложение as described here. Я предложил бы использовать LAPACK для решения этой задачи ...

Это article имеет некоторый код в C делает, что для вещественной симметричной матрицы, поэтому вам необходимо обмениваться dspsv по zspsv для обработки двойной точности комплексных матриц.

источник

2013-09-26 05:22:35

вычислительный детерминант комплексной матрицы fortran90

ответ

Смежные вопросы