Текст подпрограммы и версий ( Фортран )
aeh2c.zip , aeh2p.zip 		Тексты тестовых примеров ( Фортран )
taeh2c.zip , taeh2p.zip
Текст подпрограммы и версий ( Си )
aeh2c_c.zip , aeh2p_c.zip 		Тексты тестовых примеров ( Си )
taeh2c_c.zip , taeh2p_c.zip
Текст подпрограммы и версий ( Паскаль )
aeh2c_p.zip , aeh2z_p.zip 		Тексты тестовых примеров ( Паскаль )
taeh2c_p.zip , taeh2z_p.zip

Подпрограмма:  AEH2C

Назначение

Вычисление всех собственных значений комплексной эрмитовой матрицы QL - алгоpитмом со сдвигом.

Математическое описание

Заданная матрица приводится унитарными преобразованиями к симметрической вещественной трехдиагональной форме, которая затем используется для определения собственных значений исходной матрицы.

Дж.Х.Уилкинсон, Алгебраическая проблема собственных значений, "Наука", М., 1970.

Использование 

    SUBROUTINE  AEH2C (N, AR, AI, EV, RAB, IERR) 
   Параметры
N - порядок исходной матрицы (тип: целый);
AR, AI - вещественные двумерные массивы размерности N на N, содержащие соответственно действительную и мнимую части исходной матрицы;
EV - вещественный одномерный массив размерности  N, содержащий вычисленные собственные значения, расположенные в возрастающем порядке;
RAB - вещественный двумерный массив размерности 4 на N, используемый как рабочий;
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы. Значение IЕRR полагается равным номеру собственного значения, для вычисления которого потребовалось более 30 итераций, при этом собственные значения с индексами 1, 2, ..., IЕRR - 1 вычислены правильно и расположены в возрастающем порядке, но они не обязательно являются самыми меньшими из всех  N собственных значений. 
   Версии
AEH2P - вычисление всех собственных значений комплексной эрмитовой матрицы с повышенной точностью. Массивы АR, АI, ЕV, RАВ имеют тип DОUВLЕ РRЕСISIОN. 
   Вызываемые подпрограммы
UTAE10 - подпрограмма выдачи диагностических сообщений при работе подпрограмм AEH2C и AEH2P. 
   Замечания по использованию
  Подпрограммы AEH2C(P) сохраняют элементы исходной матрицы.
  Подпрограмма АЕH2P использует служебную подпрограмму AA01P.

Пример использования 

       DIMENSION  AR(4, 4), AI(4, 4), EV(4), RAB(16)
       DATA  AR /1., 0.42, 0.54, 0.66, 0.42, 1., 0.32, 0.44, 0.54, 0.32, 
      *                  1., 0.22, 0.66, 0.44, 0.22, 1./
       DATA  AI /16*0./
       N = 4
       CALL  AEH2C (N, AR, AI, EV, RAB, IERR)

Результаты:

       EV(1)  =  0.2422607 , 
       EV(2)  =  0.6382838 , 
       EV(3)  =  0.7967067 , 
       EV(4)  =  2.3227488 ;

       IERR  =  0