Документ взят из кэша поисковой машины. Адрес оригинального документа : 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. Обращение матриц общего вида
Текст подпрограммы и версий ( Фортран )
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 увеличением элемента  ai j на величину D. Подпрограмма AIG5R выполняет такую модификацию матрицы AI, что в результате формируется матрица AIM = AM - 1 без непосредственного обращения матрицы AM.

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

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

Параметры

AI - вещественный двумерный массив размеров N на N, в котором задается известная матрица, обратная к исходной матрице A;
N - порядок матрицы AI (тип: целый);
I, J - заданные номера строки и столбца элемента  ai j  исходной матрицы A, значение которого увеличивается на величину D (тип: целый);
D - заданное значение константы, на которую увеличивается элемент  ai 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  |