Текст подпрограммы и версий ( Фортран )
aee3r.zip , aee3d.zip
Тексты тестовых примеров ( Фортран )
taee3r.zip , taee3d.zip
Текст подпрограммы и версий ( Си )
aee3r_c.zip , aee3d_c.zip
Тексты тестовых примеров ( Си )
taee3r_c.zip , taee3d_c.zip
Текст подпрограммы и версий ( Паскаль )
aee3r_p.zip , aee3e_p.zip
Тексты тестовых примеров ( Паскаль )
taee3r_p.zip , taee3e_p.zip

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

Назначение

Вычисление К минимальных собственных значений симметричной трехдиагональной матрицы.

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

Подпрограма АЕЕ3R вычисляет К минимальных собственных значений симметричной трехдиагональной матрицы методом бисекций.

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

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

    SUBROUTINE  AEE3R ( A, B, N, RU, W, K) 

Параметры

A - вещественный вектор длины N, содержащий диагональные элементы симметричной трехдиагональной матрицы;
B - вещественный вектор длины N, содержащий в последних N - 1 компонентах поддиагональные элементы симметричной трехдиагональной матрицы;
N - заданный порядок симметричной трехдиагональной матрицы (тип: целый);
RU - вещественный вектор длины К, используемый как рабочий;
W - вещественный вектор длины К, содержащий вычисленные собственные значения в неубывающем порядке;
K - заданное число искомых минимальных собственных значений (тип: целый).

Версии

AEE3D - вычисление К минимальных собственных значений симметричной трехдиагональной матрицы, заданной с двойной точностью.

Вызываемые подпрограммы : нет

Замечания по использованию

  1. 

Подпрограмма АЕЕ3R не сохраняет исходные векторы А и В.

  2. 

В подпрограмме АЕЕ3D параметры А, В, RU, W должны иметь тип DОUВLЕ РRЕСISIОN.

  3.  Подпрограмма АЕЕ3R может быть использована и для вычисления К минимальных собственных значений комплексной эрмитовой трехдиагональной матрицы. Для этого нужно предварительно привести комлексную эрмитову трехдиагональную матрицу к вещественной симметричной трехдиагональной матрице преобразованиями подобия (например, обратиться к подпрограмме АЕF0С).
  4.  Подпрограмма АЕE3R использует служебную подпрограмму AMA1R.
Подпрограмма АЕE3D использует служебную подпрограмму AMA1D.

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

      DIMENSION  A(50), B(50), RU(15), W(15)
      DATA  A /1., 49., 48*0.0/
      DATA  B /0., 7., 48*0./
      N = 50
      K = 15
      CALL  AEE3R (A, B, N, RU, W, K)

Результаты:

      W(1)  =  W(2)  =  ...  =  W(15)  =  -2.55*10-12