|
Текст подпрограммы и версий aee2r_c.zip , aee2d_c.zip |
Тексты тестовых примеров taee2r_c.zip , taee2d_c.zip |
Вычисление всех собственных значений вещественной симметрической трехдиагональной матрицы при помощи неявного Q*L - алгоритма.
Подпрограмма aee2r_c определяет собственные значения симметрической трехдиагональной матрицы, располагая их в возрастающем порядке.
Дж.Х.Уилкинсон, Алгебраическая проблема собственных значений, "Наука", М., 1970.
int aee2r_c (integer *n, real *a, real *ev, real *rab,
integer *ierr)
Параметры
| n - | порядок исходной матрицы (тип: целый); |
| a - | вещественный двумерный массив размерности n на 2, содержащий в последних n - 1 компонентах первого столбца элементы нижней диагонали, а во втором столбце - элементы главной диагонали; |
| ev - | вещественный одномерный массив размерности n, содержащий вычисленные собственные значения в возрастающем порядке; |
| rab - | вещественный одномерный массив размерности n, используемый как рабочий; |
| ierr - | целая переменая, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы; значение ierr полагается равным номеру собственного значения, для вычисления которого потребовалось более 30 итераций, при этом собственные значения с индексами 1, 2, ..., ierr - 1 вычислены правильно и расположены в возрастающем порядке, но они не обязательно являются самыми меньшими из всех n собственных значений. |
Версии
| aee2d_c - | вычисление всех собственных значений вещественной симметрической трехдиагональной матрицы с повышенной точностью. Массивы a, ev, rab имеют тип double. |
Вызываемые подпрограммы
| utae10_c - | подпрограмма выдачи диагностических сообщений при работе подпрограмм aee2r_c и aee2d_c. |
Замечания по использованию
| Подпрограмма aee2r_c сохраняет исходную информацию. |
int main(void)
{
/* Initialized data */
static float a[10] /* was [5][2] */ = { 0.f,1.f,1.f,1.f,1.f,0.f,0.f,0.f,
0.f,0.f };
/* Local variables */
static int ierr;
extern int aee2r_c(int *, float *, float *, float *, int *);
static int i__;
static float ev[5], rab[5];
#define a_ref(a_1,a_2) a[(a_2)*5 + a_1 - 6]
for (i__ = 1; i__ <= 5; ++i__) {
printf("\n %16.7e %16.7e \n",
a_ref(i__, 1), a_ref(i__, 2));
}
aee2r_c(&c__5, a, ev, rab, &ierr);
for (i__ = 1; i__ <= 5; ++i__) {
printf("\n %16.7e %16.7e \n",
a_ref(i__, 1), a_ref(i__, 2));
}
for (i__ = 1; i__ <= 5; ++i__) {
printf("\n %16.7e \n", ev[i__-1]);
}
printf("\n %5i \n", ierr);
return 0;
} /* main */
Результаты:
ev(1) = -1.7320508 ,
ev(2) = -1.0 ,
ev(3) = 6.28*10-12 ,
ev(4) = 1.0 ,
ev(5) = 1.7320508 ;
ierr = 0