Текст подпрограммы и версий ( Фортран ) aih2r.zip aih2d.zip |
Тексты тестовых примеров ( Фортран ) taih2r.zip taih2d.zip |
Текст подпрограммы и версий ( Си ) aih2r_c.zip aih2d_c.zip |
Тексты тестовых примеров ( Си ) taih2r_c.zip taih2d_c.zip |
Текст подпрограммы и версий ( Паскаль ) aih2r_p.zip aih2e_p.zip |
Тексты тестовых примеров ( Паскаль ) taih2r_p.zip taih2e_p.zip |
Обращение положительно определенной симметричной ленточной матрицы с компактной формой представления с использованием метода квадратного корня (метода Холецкого).
Подпрограмма АIН2R вычисляет матрицу В (симметричная, в компактной форме), такую, что А*В = Е, где А - исходная симметричная ленточная положительно определенная заданная в компактной форме, а Е - единичная матрицы. Используется метод квадратного корня (метод Холецкого).
В.В.Воеводин, Численные методы алгебры (теория и алгорифмы), Наука, М., 1966 г.
SUBROUTINE AIH2R (A, AINV, N, NC, IERR)
Параметры
A - | вещественный двумерный массив размера N на (NС+1), в котором содержится заданная матрица, записанная в компактной форме; |
AINV - | вещественный вектор размерности N на (N+1)/2, в котором запоминается вычисленная обратная матрица в компактной форме; |
N - | порядок матрицы (тип: целый); |
NC - | число ненулевых диагоналей заданной матрицы по одну сторону от главной диагонали (тип: целый); |
IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы; при этом: |
IЕRR=65 - | если заданная матрица не является положительно определенной. |
Версии
AIH2D - | обращение положительно определенной симметричной ленточной матрицы с компактной формой представления с использованием метода квадратного корня (метод Холецкого) с повышенной точностью. |
Вызываемые подпрограммы
UTAI10 - | подпрограмма печати диагностических сообщений при работе подпрограммы АIН2R; |
UTAI11 - | подпрограмма печати диагностических сообщений при работе подпрограммы АIН2D; |
Замечания по использованию
При работе подпрограммы исходная матрица не сохраняется. На место вектора А помещается ее треугольное разложение, полученное методом квадратного корня (методом Холецкого). Для подпрограммы АIН2D параметры А, АINV должны иметь тип DОUВLЕ РRЕСISIОN. |
DIMENSION A(3, 2), AINV(6) DATA A /0., -1., -1., 2., 2., 2./ N = 3 NC = 1 CALL AIH2R (A, AINV, N, NC, IERR) Результаты: АINV = (0.75, 0.5, 1., 0.25, 0.5, 0.75), IERR = 0