Текст подпрограммы и версий ( Фортран )
afg0r.zip , afg0d.zip , afg0c.zip
Тексты тестовых примеров ( Фортран )
tafg0r.zip , tafg0d.zip , tafg0c.zip
Текст подпрограммы и версий ( Си )
afg0r_c.zip , afg0d_c.zip , afg0c_c.zip
Тексты тестовых примеров ( Си )
tafg0r_c.zip , tafg0d_c.zip , tafg0c_c.zip
Текст подпрограммы и версий ( Паскаль )
afg0r_p.zip , afg0e_p.zip , afg0c_p.zip
Тексты тестовых примеров ( Паскаль )
tafg0r_p.zip , tafg0e_p.zip , tafg0c_p.zip

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

Назначение

Треугольное разложение вещественной матрицы методом Гаусса без выбора ведущего элемента.

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

Выполняется факторизация матрицы А вида А = LU , где L - нижняя, U - верхняя треугольные матрицы. Матрица L имеет единичные диагональные элементы.

В.В.Воеводин, Р.В.Петрина, Комплекс алгоритмов, основанных на преобразованиях Гаусса, в пакете линейной алгебры, Сб. "Численный анализ на ФОРТРАНе", вып.3, Изд-во МГУ, 1973.

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

    SUBROUTINE  AFG0R (A, N) 

Параметры

A - двумерный N*N массив, в котором задается исходная матрица (тип: RЕАL); по окончании работы подпрограммы в массиве А на соответствующих местах запоминаются элементы матрицы U и поддиагональные элементы матрицы L;
N - заданный порядок исходной матрицы (тип: целый).

Версии

AFG0D - треугольное разложение матрицы, заданной с удвоенной точностью, методом Гаусса без выбора ведущего элемента.
AFG0C - треугольное разложение комплексной матрицы методом Гаусса без выбора ведущего элемента.

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

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

  1. 

В подпрограмме АFG0D массив А имеет тип DОUВLЕ РRЕСISIОN.

  2.  В подпрограмме АFG0С массив А имеет тип СОМРLЕХ.

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

      DIMENSION   A(4, 4)
      DATA  A /7.9, 8.5, 4.3, 3.2, 5.6, -4.8, 4.2, -1.4, 5.7, .8, -3.2, 
   *                -8.9, -7.2, 3.5, 9.3, 3.3/
      CALL  AFG0R (A, 4)

Результат:

                 | 7.900     5.600    5.700   -7.200 |
                 | 1.076  -10.825  -5.333   11.247 |
      A  =    | 0.544    -0.106  -6.870   14.415 |   , 
                 | 0.405      0.339   1.369  -17.323 |

это означает, что

                 | 1           0          0            0       |
                 | 1.076    1          0            0       |
      L  =    | 0.544   -0.106   1            0       |   , 
                 | 0.405    0.339   1.369     1       |
                 | 7.900    5.600   5.700   -7.200 |

                 | 0      -10.825  -5.333    11.247 |
      U  =    | 0        0           -6.870    14.415 |  .
                 | 0        0            0          -17.323 |