БЧА НИВЦ МГУ. AIG5R. Обращение матриц общего вида

Документ взят из кэша поисковой машины. Адрес оригинального документа : http://num-anal.srcc.msu.su/lib_na/cat/ai/aig5r.htm
Дата изменения: Thu Dec 3 12:35:58 2015
Дата индексирования: Sun Apr 10 00:34:50 2016
Кодировка: Windows-1251

Текст подпрограммы и версий ( Фортран ) aig5r.zip aig5d.zip	Тексты тестовых примеров ( Фортран ) taig5r.zip taig5d.zip
Текст подпрограммы и версий ( Си ) aig5r_c.zip aig5d_c.zip	Тексты тестовых примеров ( Си ) taig5r_c.zip taig5d_c.zip
Текст подпрограммы и версий ( Паскаль ) aig5r_p.zip aig5e_p.zip	Тексты тестовых примеров ( Паскаль ) taig5r_p.zip taig5e_p.zip

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

Назначение

Модификация обратной матрицы после изменения одного элемента исходной матрицы.

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

Пусть известна обратная матрица AI = A^- 1, где A - вещественная квадратная матрица порядка N. Пусть матрица AM получена из матрицы A увеличением элемента a_i j на величину D. Подпрограмма AIG5R выполняет такую модификацию матрицы AI, что в результате формируется матрица AIM = AM^- 1 без непосредственного обращения матрицы AM.

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

    SUBROUTINE  AIG5R (AI, N, I, J, D, AIM)

   Параметры

AI -	вещественный двумерный массив размеров N на N, в котором задается известная матрица, обратная к исходной матрице A;
N -	порядок матрицы AI (тип: целый);
I, J -	заданные номера строки и столбца элемента a_i j исходной матрицы A, значение которого увеличивается на величину D (тип: целый);
D -	заданное значение константы, на которую увеличивается элемент a_i j (тип: вещественный);
AIM -	вещественный двумерный массив размеров N на N, в котором помещается модифицированная обратная матрица.

   Версии

AIG5D -

модификация обратной матрицы после изменения одного элемента исходной матрицы в режиме удвоенной точности; при этом параметры AI, D и AIM должны иметь тип DOUBLE PRECISION.

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

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

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

      DIMENSION  AI(2,2), AIM(2,2) 
      DATA  AI(1,1) /-2.0/, AI(1,2) /1.0/, AI(2,1) /1.5/, AI(2,2) /-0.5/ 
      N = 2 
      I = 1 
      J = 2 
      D = 3.0 
      CALL  AIG5R (AI, N, I, J, D, AIM) 

Результаты: 

                       |   -0.363636     0.454545    |
    AIM    =     |
                       |    0.272727    -0.0909091  |