|
Текст подпрограммы и версий amy1r_c.zip , amy1d_c.zip , amy1c_c.zip |
Тексты тестовых примеров tamy1r_c.zip , tamy1d_c.zip , tamy1c_c.zip |
Вычисление ∞ - нормы матрицы.
amy1r_c вычисляет ∞ - норму матрицы А порядка N:
N
max ∑ | a i j | , i = 1, 2, ..., N .
i j =1
real amy1r_c (real *a, integer *n)
Параметры
| a - | вещественный двумерный массив размера n на n, в котором задается исходная матрица; |
| n - | порядок матрицы (тип: целый). |
Версии
| amy1d_c - | вычисление ∞ - нормы матрицы, заданной с удвоенной точностью. |
| amy1c_c - | вычисление ∞ - нормы комплексной матрицы. |
Вызываемые подпрограммы
|
avz1r_c - avz1d_c avz1c_c | вычисление суммы абсолютных значений компонент вектора или заданного подмножества компонент этого вектора для случаев, когда вектор задан с обычной или удвоенной точностью и комплексного вектора соответственно. Используются в amy1r_c, amy1d_c и amy1c_c соответственно. |
Замечания по использованию
|
В функции amy1d_c массив a должен иметь тип double. В функции amy1c_c массив a должен иметь тип complex. |
int main(void)
{
/* Initialized data */
static float a[16] /* was [4][4] */ = { 1.f,2.f,3.f,5.f,2.f,3.f,4.f,1.f,
3.f,4.f,1.f,2.f,4.f,1.f,2.f,3.f };
/* Local variables */
extern float amy1r_c(float *, int *);
static int n;
static float y;
n = 4;
y = (float)amy1r_c(a, &n);
printf("\n %16.7e \n", y);
return 0;
} /* main */
Результат: y = 11.