Текст подпрограммы и версий ( Фортран ) aee4r.zip , aee4d.zip |
Тексты тестовых примеров ( Фортран ) taee4r.zip , taee4d.zip |
Текст подпрограммы и версий ( Си ) aee4r_c.zip , aee4d_c.zip |
Тексты тестовых примеров ( Си ) taee4r_c.zip , taee4d_c.zip |
Текст подпрограммы и версий ( Паскаль ) aee4r_p.zip , aee4e_p.zip |
Тексты тестовых примеров ( Паскаль ) taee4r_p.zip , taee4e_p.zip |
Вычисление К максимальных собственных значений симметричной трехдиагональной матрицы.
Подпрограма АЕЕ4R вычисляет матодом бисекций К максимальных собственных значений симметричной трехдиагональной матрицы.
Дж.Х.Уилкинсон, Алгебраическая проблема собственных значений, "Наука", М., 1970.
SUBROUTINE AEE4R ( A, B, N, RU, W, K)
Параметры
A - | вещественный вектор длины N, содержащий диагональные элементы симметричной трехдиагональной матрицы; |
B - | вещественный вектор длины N, содержащий в последних N - 1 компонентах поддиагональные элементы симметричной трехдиагональной матрицы; |
N - | заданный порядок симметричной трехдиагональной матрицы (тип: целый); |
RU - | вещественный вектор длины К, используемый как рабочий; |
W - | вещественный вектор длины К, содержащий вычисленные собственные значения в неубывающем порядке; |
K - | заданое число искомых собственных значений (тип: целый). |
Версии
AEE4D - | вычисление к максимальных собственных значений симметричной трехдиагональной матрицы, заданной с двойной точностью. |
Вызываемые подпрограммы : нет
Замечания по использованию
1. |
Подпрограмма АЕЕ4R не сохраняет исходные векторы А и В. | |
2. |
В подпрограмме АЕЕ4D параметры А, В, RU, W имеют тип DОUВLЕ РRЕСISIОN. | |
3. | Подпрограмма АЕЕ4R может быть использована и для вычисления К максимальных собственных значений комплексной эрмитовой трехдиагональной матрицы. Для этого нужно предварительно привести комплексную эрмитовую трехдиагональную матрицу к вещественной трехдиагональной симметричной матрице преобразованиями подобия (например, обратиться к подпрограмме АFE0С). | |
4. |
Подпрограмма АЕE4R использует служебную подпрограмму AMA1R. Подпрограмма АЕE4D использует служебную подпрограмму AMA1D. |
DIMENSION A(50), B(50), RU(5), S(5) DATA A /1., 49., 48*0.0/ DATA B /0., 7., 48*0.0/ N = 50 K = 5 CALL AEE4R (A, B, N, RU, S, K) Результаты: S(1) = S(2) = S(3) = S(4) = -2.55*10-12, S(5) = 50 .