... В подпрограмме amazd_c параметры a, b, и c должны иметь тип double. int main(void) { /* Initialized data */ static float a[16] /* was [4][4] */ = { 1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f, 1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f }; static float b[10] = { 1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f }; /* Local variables */ static float c__[16] /* was [4][4] */; static int n, i__; extern int amazr_c(float *, int *, float *, float *); n = 4; amazr_c(a, , b, c__); for (i__ = 0; i__ ...
... Нахождение максимального по модулю элемента строки (или столбца) симметричной матрицы, заданной в компактной форме. real amy6r_c (real *a, integer *n, integer *i, integer *j) . ... вещественный вектор длины n на (n + 1)/2, в котором задается симметричная матрица в компактной форме; . ... целая переменная, содержащая второй (или первый) индекс максимального по модулю элемента исходной матрицы в предположении, что она задана в виде квадратной матрицы, а не в компактной форме. amy6d_c - . ...
... Вычисление суммы абсолютных значений элементов строки (или столбца) симметричной матрицы, заданной в компактной форме. real amy5r_c (real *a, integer *n, integer *i) . ... вещественный вектор длины n на (n + 1)/2, в котором задается симметричная матрица в компактной форме; . ... порядок матрицы (тип: целый); . ... вычисление суммы абсолютных значений элементов строки (или столбца) симметричной матрицы, заданной в компактной форме с удвоенной точностью. amy5c_c - . ...
... Вычисление первой нормы симметричной матрицы, заданной в компактной форме. АМY3R вычисляет первую норму симметричной матрицы А порядка N, заданной в компактной форме N max | ... вычисление первой нормы симметричной матрицы, заданной в компактной форме с удвоенной точностью. ... вычисление суммы абсолютных значений элементов строки (или столбца) симметричной матрицы, заданной в компактной форме с обычной и удвоенной точностью, и эрмитовой матрицы, заданной в компактной форме, соответственно. ...
... Вычисление первой нормы симметричной матрицы, заданной в компактной форме. amy3r_c вычисляет первую норму симметричной матрицы А порядка N, заданной в компактной форме N max | ... вычисление первой нормы эрмитовой матрицы, заданной в компактной форме. amy5r_c - . ... вычисление суммы абсолютных значений элементов строки (или столбца) симметричной матрицы, заданной в компактной форме с обычной и удвоенной точностью, и эрмитовой матрицы, заданной в компактной форме, соответственно. ...
... amy2r.zip , amy2d.zip , amy2c.zip . ... amy2r_p.zip , amy2e_p.zip , amy2c_p.zip . ... Вычисление первой нормы матрицы. АМY2R вычисляет первую норму матрицы А порядка N: N max | ... j i =1 REAL FUNCTION AMY2R (A, N) . ... вещественный двумерный массив размера N на N, в котором задается исходная матрица; . ... заданный порядок матрицы (тип: целый). ... вычисление первой нормы матрицы, заданной с удвоенной точностью. ... вычисление первой нормы комплексной матрицы. ...
... Вычисление первой нормы матрицы. amy2r_c вычисляет первую норму матрицы А порядка N: N max | ... В функции amy2c_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 amy2r_c(float *, int *); static int n; static float y; n = 4; y = (float)amy2r_c(a, ); printf("\n %16.7e \n", y); return 0; } /* main */ Результат: y = 11. ...
... amy1r.zip , amy1d.zip , amy1c.zip . ... amy1r_p.zip , amy1e_p.zip , amy1c_p.zip . ... Вычисление - нормы матрицы. АМY1R вычисляет - норму матрицы А порядка N: N max | ... i j =1 REAL FUNCTION AMY1R (A, N) . ... вещественный двумерный массив размера N на N, в котором задается исходная матрица; . ... порядок матрицы (тип: целый). ... вычисление - нормы матрицы, заданной с удвоенной точностью. ... вычисление - нормы комплексной матрицы. ... В функции АМY1D массив А должен иметь тип DОUВLЕ РRЕСISIОN. ...
... Вычисление - нормы матрицы. amy1r_c вычисляет - норму матрицы А порядка N: N max | ... В функции 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, ); printf("\n %16.7e \n", y); return 0; } /* main */ Результат: y = 11. ...
... amtpr.zip , amtpd.zip . ... Транспонирование прямоугольной разреженной матрицы, заданной в формате RR (С) U, с одновременной перестановкой строк. ... Технология разреженных матриц. М.: Мир, 1988 SUBROUTINE AMTPR ( IA, JA, AN, N, M, K, IAT, JAT, ANT) . ... IAT, JAT, - . ANT . ... транспонирование прямоугольной разреженной матрицы, заданной в формате RR (С) U, с одновременной перестановкой строк в режиме удвоенной точности; при этом параметры AN и ANT должны быть описаны как DOUBLE PRECISION. ...
... Вычисление произведения вещественной теплицевой или ханкелевой матрицы на вещественный вектор на основе алгоритма Быстрого Преобразования Фурье. ... a 0 a -1 a -2 ... a -(N-2) a -(N-1) | ... Дополним теплицеву матрицу А порядка N до циркулянтной матрицы А 0 порядка М, которая определяется М элементами своего первого столбца: a -(N-1) , a -(N-2) , ... , a -1 , a 0 , a 1 , ... , a N-1 , 0, 0, ... , 0 , . ... Аналогичный метод умножения теплицевой матрицы на вектор изложен в работе [3]. ...
... Вычисление произведения вещественной теплицевой или ханкелевой матрицы на вещественный вектор на основе алгоритма Быстрого Преобразования Фурье. ... Аналогичный метод умножения теплицевой матрицы на вектор изложен в работе [3]. ... подпрограммы умножения вещественной циркулянтной матрицы на вещественный вектор и компексной циркулянтной матрицы на комплексный вектор на основе алгоритма Быстрого Преобразования Фурье (используются в ams1r_c и ams1c_c соответственно). ...
... Умножение прямоугольной разреженной матрицы, заданной в формате RR (C) U, на диагональную матрицу . ... умножение прямоугольной разреженной матрицы , заданной в формате RR (C) U, на диагональную матрицу в режиме удвоенной точности; при этом параметры an и d должны иметь тип double int main(void) { /* Initialized data */ static int ia[5] = { 1,3,4,7,10 }; static int ja[9] = { 4,1,3,4,2,3,4,3,1 }; static float an[9] = { 1.f , 1.f , 1.f , 1.f , 1.f , 1.f , 1.f , 1.f , ...
... Технология разреженных матриц. ... численное умножение обратной для нижней треугольной разреженной матрицы с единичной диагональю, заданной в формате RR (L) U, на прямоугольную разреженную матрицу, заданную в формате RR (C) U, в режиме удвоенной точности; при этом параметры unt, bn, xn и p должны иметь тип double int main(void) { /* Initialized data */ static int iut[5] = { 1,1,1,2,4 }; static int jut[3] = { 2,3,1 }; static int ib[5] = { 1,3,4,5,6 }; static ...
Текст подпрограммы и версий ( Фортран ) . ... Численное умножение прямоугольных разреженных матриц, заданных в формате RR (C) U . ... Пусть в формате RR (C) U заданы матрица A размеров p на q и матрица B размеров q на r. Пусть известен также портрет результирующей матрицы C = A * B в формате RR (C) U, полученный, например, предварительно при помощи подпрограммы AMMIR . Подпрограмма AMMJR вычисляет ненулевые элементы матрицы C посредством алгоритма, изложенного в описании подпрограммы AMMIR . ...
... заданное число столбцов матрицы A и строк матрицы B (тип: целый); . ... численное умножение прямоугольных разреженных матриц , заданных в формате RR (C) U, в режиме удвоенной точности; при этом параметры an, bn, cn и x должны иметь тип double int main(void) { /* Initialized data */ static int ia[5] = { 1,4,4,6,7 }; static int ja[6] = { 1,5,4,4,2,1 }; static int ib[6] = { 1,2,4,6,6,7 }; static int jb[6] = { 2,2,1,2,1,2 }; static int ic[5] = { 1,2,2,4,5 }; static ...
... Умножение матрицы, обратной к заданной в компактной форме вещественной ленточной матрице A, на прямоугольную матрицу и оценка числа обусловленности матрицы A. Для заданной в компактной форме ленточной вещественной . ... М.: Мир, 1980. int ammhr_c (real *a, integer *ma, integer *n, integer *ml, integer *mu, integer *nlead, real *b, integer *mb, integer *nn, real *rcond, real *z, integer *ierr) . ... подпрограмма треугольной факторизации и оценки числа обусловленности ленточной матрицы A. utafsi_c - ....
... Умножение матрицы, обратной к заданной вещественной матрице общего вида, на прямоугольную матрицу. ... Изд-во "Мир", М: 1980. int ammgr_c (real *a, integer *m, integer *n, integer *nlead, real *b, integer *mm, integer *nn, integer *ierr) . ... порядок матрицы A и число строк матриц B и C (тип: целый); . ... ierr - . ... умножение матрицы, обратной к вещественной матрице, заданной с удвоенной точностью, на вещественную прямоугольную матрицу, заданную с удвоенной точностью. ammgc_c - . ...