Текст подпрограммы и версий
av12r_c.zip  av12d_c.zip  av12c_c.zip  		Тексты тестовых примеров
tav12r_c.zip  tav12d_c.zip  tav12c_c.zip 

Функция:  av12r_c

Назначение

Вычисление евклидовой нормы вещественного вектора.

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

Функция av12r_c вычисляет евклидову норму (длину)  n - мерного вектора по формуле: 

     || X ||2 = ( | x1 |2 + | x2 |2 + ... + | xN |2 )1/2 ,

     где  X = ( x1, x2, ... , xN ) .

Вектор X не обязательно должен занимать последовательные слова памяти машины и может располагаться в памяти с некоторыми постоянными интервалами между компонентами, которые называются параметром приращения.

Этот параметр может быть меньше нуля, если обработка компонент вектора должна проводится в обратном порядке.

В некоторых случаях, когда компоненты вектора имеют достаточно большое значение, при нахождении суммы квадратов компонент может возникнуть переполнение центрального процессора, хотя конечный результат меньше максимального числа, представимого на ЭВМ.

Для того, чтобы переполнение не возникало, проводится масштабирование по формуле: 

     evnorm_c = rnorm_c * ( (x12 + x22 +...+ xN2) / rnorm2 )1/2
  где    rnorm_c = 2lmax + 1 ,        lmax = log2AMAX ,
            AMAX = max ( | x1 |, | x2 |, ... ,| xn | ) ,

Масштабирование проводится в случае АМАХ > 1*106, когда возникает возможность переполнения.

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

    real av12r_c (integer *n, real *sx, integer *incx, integer *ierr)

   Параметры
n - размерность вектора (тип: целый);
sx - вещественный одномерный массив длины n * | incx |, содержащий заданный вектор;
incx - параметр приращения, задающий расположение вектора в массиве sx (тип: целый);
ierr - целая переменная, указывающая характер ошибки; при этом:
ierr=65 - когда n ≤ 0, 
   Версии
av12d_c - вычисление евклидовой нормы вектора, заданного с удвоенной точностью;
av12c_c - вычисление евклидовой нормы комплексного вектоpа. 
   Вызываемые подпрограммы
utav10_c - подпрограмма выдачи диагностических сообщений при работе подпрограмм av12r_c и av12c_c.
utav11_c - подпрограмма выдачи диагностических сообщений при работе подпрограммы av12d_c. 
   Замечания по использованию
 

В функции av12d_c одномерный массив, содержащий вектор, обозначен dx и имеет тип double.

В функции av12c_c одномерный массив, содержащий вектор, обозначен cx и имеет тип complex.

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

int main(void)
{
    /* Initialized data */
    static float sx[11] = { 2.f,0.f,2.f,0.f,2.f,0.f,2.f,0.f,2.f,0.f,2.f };

    /* Local variables */
    extern float av12r_c(int *, float *, int *, int *);
    static int incx, ierr, n;
    static float r12;

    n = 6;
    incx = 2;
    ierr = 0;
    r12 = (float)av12r_c(&n, sx, &incx, &ierr);

    printf("\n  %12.3e \n", r12);
    printf("\n  %5i \n", ierr);
    return 0;
} /* main */


Результаты:    r12  =  4.898979 ,     ierr  =  0