Текст подпрограммы и версий ( Фортран )
astar.zip , astad.zip , astac.zip
Тексты тестовых примеров ( Фортран )
tastar.zip , tastad.zip , tastac.zip
Текст подпрограммы и версий ( Си )
astar_c.zip , astad_c.zip , astac_c.zip
Тексты тестовых примеров ( Си )
tastar_c.zip , tastad_c.zip , tastac_c.zip
Текст подпрограммы и версий ( Паскаль )
astar_p.zip , astae_p.zip , astac_p.zip
Тексты тестовых примеров ( Паскаль )
tastar_p.zip , tastae_p.zip , tastac_p.zip

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

Назначение

Решение системы линейных алгебраических уравнений с теплицевой матрицей.

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

Решается система Аx = b, где А - теплицева матрица размера М на М. Эффективно используется теплицевость матрицы А.

С.Н.Воеводина. Решение системы уравнений с клеточно - теплицевыми матрицами. В сб. "Вычислительные методы и программирование", вып. 24, Изд-во МГУ, 1975.

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

    SUBROUTINE  ASTAR (A, X, R, M) 

Параметры

A - вещественный одномерный массив длины 2М - 1, в котором задается матрица системы следующим образом: в массиве А записывается ее первая строка, за которой следуют элементы первого столбца начиная со второго элемента;
X - вещественный одномерный массив длины М, в котором задается правая часть; на выходе из подпрограммы в Х получается решение;
R - вещественный одномерный рабочий массив длины 2М;
M - порядок матрицы (тип: целый).

Версии

ASTAD - решение системы линейных алгебраических уравнений с теплицевой матрицей, заданной с двойной точностью;
ASTAC - решение системы линейных алгебраических уравнений с комплексной теплицевой матрицей.

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

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

  1. 

При обращении к подпрограмме АSКАD параметры А, Х, R имеют тип DОUВLЕ РRЕСISIОN.

  2. 

При обращении к подпрограмме АSКАС параметры А, Х, R имеют тип СОМРLЕХ.

  3.  В подпрограммах АSТАR, АSТАD, АSТАС используются соответствующие версии служебной подпрограммы АSТ0R.

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

      REAL  A(7), X(4), R(8)
      DATA  A /1., 2., 3., 4., 4., 3., 2./
      DATA  X /4*10./
      M = 4
      CALL  ASTAR (A, X, R, M)

Результат:

      X  =   ( 1., 1., 1., 1. )