Текст подпрограммы и версий ( Фортран )
askar.zip , askad.zip , askac.zip 		Тексты тестовых примеров ( Фортран )
taskar.zip , taskad.zip , taskac.zip
Текст подпрограммы и версий ( Си )
askar_c.zip , askad_c.zip , askac_c.zip 		Тексты тестовых примеров ( Си )
taskar_c.zip , taskad_c.zip , taskac_c.zip
Текст подпрограммы и версий ( Паскаль )
askar_p.zip , askae_p.zip , askac_p.zip 		Тексты тестовых примеров ( Паскаль )
taskar_p.zip , taskae_p.zip , taskac_p.zip

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

Назначение

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

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

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

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

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

    SUBROUTINE  ASKAR ( A, X, R, M, L) 
   Параметры
A - вещественный одномерный массив длины М2 на (2L - 1), в котором задается матрица системы следующим образом: массив А делится на 2L - 1 кусков длины М2. В I - ом куске,  I = 1, ..., L, записывается I - ая клетка первой клеточной строки матрицы системы. В (L + I) - ом куске,  I = 1, ..., L - 1, записывается (I + 1) - ая клетка первого клеточного столбца. Каждая клетка в своем куске записывается по столбцам;
X - вещественный одномерный массив длины М на L, в котором задается правая часть; на выходе из подпрограммы в Х получается решение;
R - вещественный одномерный рабочий массив длины 2М2L + 4М2 + М;
M - порядок клетки (тип: целый);
L - клеточный порядок исходной матрицы (тип: целый). 
   Версии
ASKAD - решение системы линейных алгебраических уравнений с комплексной клеточно - теплицевой матрицей, заданной с двойной точностью.
ASKAC - решение системы линейных алгебраических уравнений с комплексной клеточно - теплицевой матрицей. 
   Вызываемые подпрограммы
 AIG1R -
 AIG1D  
 AIG1C   		обращение матрицы методом Жордана с выбором по столбцу. 
   Замечания по использованию
  1. 

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

  2. 

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

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

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

      REAL  A(12), X(4), R(34)
      DATA  A /1., 2., 2., 1., 2., 4., 4., 2., 2., 4., 4., 2./
      DATA  X /4*9./
      M = 2
      L = 2
      CALL  ASKAR (A, X, R, M, L)

Результат:    X  =   (1., 1., 1., 1.)