Текст подпрограммы и версий ( Фортран )
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 

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

Назначение

Обращение положительно определенной симметричной ленточной матрицы с компактной формой представления с использованием метода квадратного корня (метода Холецкого).

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

Подпрограмма А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