Текст подпрограммы и версий
am09r_c.zip , am09d_c.zip , am09c_c.zip 		Тексты тестовых примеров
tam09r_c.zip , tam09d_c.zip , tam09c_c.zip

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

Назначение

Левосторонние преобразования вектора последовательностью матриц вращения, упакованных под нижней двухдиагональной частью прямоугольной матрицы размера N * M (N ≤ M).

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

Вектор B длины N последовательно умножается слева на матрицы вращения  R3 2, R4 3, ..., RN N - 1,  упакованные под нижней двухдиагональной частью заданной прямоугольной матрицы А размера N * М (N ≤ М) таким образом, что если  t = А (i, j-1)  для i > j , то матрица Ri j размера N*N является марицей плоского вращения в плоскости (i, j), у которой элементы  Ri i =  Rj j = С,   Ri j = - S,   Rj i = S, а остальные элементы главной диагонали равны 1,
 где  С = (1 - | t |2) / (1 + | t |2) ,  S = 2t / (1 + | t |2).

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

    int am09r_c (real *a, integer *n, integer *m, real *b)

   Параметры
a - двумерный массив размера n * m, содержащий прямоугольную матрицу, под нижней двухдиагональной частью которой хранится в упакованном виде последовательность матриц вращения (тип: вещественный);
n, m - число строк и столбцов матрицы  A, причем n ≤ m (тип: целый);
b - одномерный массив длины  n, в котором задается исходный вектор; в результате работы подпрограммы на месте  b запоминается преобразованный вектор (тип: вещественный). 
   Версии
am09d_c - левосторонние преобразования вектора последовательностью матриц вращения с повышенной точностью.
am09c_c - левосторонние преобразования комплексного вектора последовательностью комплексных матриц вращения. 
   Вызываемые подпрограммы:  нет. 
   Замечания по использованию
  1. 

Заданная в необходимой форме последовательность матриц вращения, используемых данной подпрограммой, получается, например, в результате работы подпрограммы, реализующей приведение прямоугольной матрицы размера N на M (N ≤ M) к нижнему треугольному виду (afp7r_c).

  2. 

В подпрограмме am09d_c пapaмeтpы  a и  b имeют тип double.

  3.  В подпрограмме am09c_c пapaмeтpы  a и  b имeют тип complex.

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

int main(void)
{
    /* Initialized data */
    static float a[20] /* was [4][5] */ = { 0.f,0.f,.5f,1.f,0.f,0.f,0.f,1.f,
                            0.f,1.f,1.f,1.f,.5f,1.f,1.f,1.f,1.f,1.f,1.f,.5f };
    static float b[4] = { 1.f,1.f,1.f,1.f };

    /* Local variables */
    extern int am09r_c(float *, int *, int *, float *);
    static int i__;

    am09r_c(a, &c__4, &c__5, b);

    for (i__ = 1; i__ <= 4; ++i__) {
         printf("\n  %16.7e \n", b[i__-1]);
    }
    for (i__ = 0; i__ <= 15; i__+=3) {
         printf("\n  %16.7e %16.7e %16.7e \n",
                a[i__], a[i__+1], a[i__+2]);
    }
    printf("\n  %16.7e %16.7e \n", a[18], a[19]);
    return 0;
} /* main */


Результат:  b  =  ( 1.0, - 1.0, 0.2, 1.4 )