... Технология разреженных матриц. ... умножение прямоугольной разреженной матрицы, заданной в формате RR (C) U, на заполненный вектор - столбец в режиме удвоенной точности; при этом параметры an, b и c должны быть описаны в режиме double int main(void) { /* Initialized data */ static int ia[5] = { 1,3,7,9,11 }; static int ja[10] = { 5,3,4,3,1,5,1,6,4,2 }; static float an[10] = { - 1.f ,2.f,3.f,3.f,4.f,7.f,-2.f,- 1.f , 1.f , 1.f }; static float b[6] = { 1.f , 1.f , ...
... В подпрограмме am18c_c пapaмeтpы a и b имeют тип complex. int main(void) { /* Initialized data */ static float a[30] /* was [5][6] */ = { 1.f,3.f,0.f,1.f,3.f,3.f,-2.f,0.f, 1.f,5.f,2.f,0.f,1.f,1.f,1.f,1.f,1.f,2.f,0.f,1.f,-1.f,1.f,1.f,2.f, 0.f,0.f,1.f,1.f,3.f,2.f }; static float x[6] = { 1.f,0.f,2.f,-1.f,1.f,0.f }; /* Local variables */ extern int am18r_c(float *, int *, int *, float *), afp3r_c(float *, int *, int *, float *, int *); static float t[6]; static int is[5], i__; for (i__ = 0; i__ ...
... Умножение симметричной матрицы, заданной в компактной форме, на вектор. ... int am16r_c (real *a, real *b, real *c, integer *n) . ... умножение с повышенной точностью вещественной симметричной матрицы, заданной в компактной форме, на вектор; . ... В подпрограмме am16c_c пapaмeтpы a, b, c имeют тип complex. int main(void) { /* Local variables */ extern int am16r_c(float *, float *, float *, int *); static float a[10], b[4], c__[4]; static int k, n, i__; static float r__; n = 4; for (k = 1; k ...
... Преобразование вектора последовательностью матриц отражения, упакованных над главной диагональю прямоугольной матрицы размера N * M (N M). ... Q M - 1 - матрицы отражения размера M * M, B - заданный вектор длины M. Матрицы отражения Q i имеют вид Q i = Е - W i W i H , . ... Ненулевые компоненты векторов W i , порождающих матрицы Q i , задаются в строках над главной диагональю прямоугольной матрицы А размера N * М (N М), т.е. W i = ( 0, 0, ... , 0, a i, i+1 , a i, i+2 , ... , a i, M ) H . ...
... Умножение вектора на прямоугольную матрицу. Вычисляется призведение С = А В вектора А длины N и прямоугольной матрицы B размера N * М. Результатом этого произведения является вектор длины М. int am06r_c (real *a, real *b, real *c, integer *n, integer *m) . ... одномерный массив длины n, в котором задается исходный вектор A (тип: вещественный); . ... умножение вектора на прямоугольную матрицу с повышенной точностью; . ... умножение комплексного вектора на комплексную прямоугольную матрицу; . ...
... Умножение прямоугольной матрицы на вектор. Вычисляется произведение С = А В прямоугольной матрицы А размера N * М и вектора B длины М. Результатом этого произведения является вектор длины N. int am05r_c (real *a, real *b, real *c, integer *n, integer *m) . ... двумерный массив размера n * m, в котором задается исходная матрица A (тип: вещественный); . ... умножение прямоугольной матрицы на вектор с повышенной точностью; . ... умножение комплексной прямоугольной матрицы на комплексный вектор; . ...
... Вычисление евклидовой нормы симметричной матрицы, заданной в компактной форме. Для симметричной матрицы А порядка N, заданной в компактной форме, вычисляется евклидова норма N N || ... 2 ) 1/2 k =1 m =1 real am01r_c (real *a, integer *n) . ... вектор длины n * (n + 1)/2, в котором задается симметричная матрица A в компактной форме (тип: вещественный); . ... вычисление евклидовой нормы вещественной симметричной матрицы, заданной в компактной форме с удвоенной точностью; . ...
... ama8r.zip , ama8d.zip . ... Численное сложение двух прямоугольных разреженных матриц, заданных в формате RR (C) U. Описание формата RR (C) U приведено в описании подпрограммы AMTSR . Пусть заданы прямоугольные разреженные матрицы A и B размеров N на M в формате RR (C) U и портрет матрицы C, равной A + B, также в формате RR (C) U, который может быть получен подпрограммой AMA7R на символическом этапе сложения. Подпрограмма AMA8R по этим исходным данным вычисляет ненулевые элементы матрицы C . ...
... Сложение симметричной матрицы и ленточной симметричной, заданных в компактной форме. ama6r_c вычисляет сумму симметричной матрицы А порядка N, заданной в компактной форме, и ленточной симметричной матрицы B того же порядка, заданной также в компактной форме. Результатом этого сложения является результирующая симметричная матрица C порядка N. Матрица C запоминается в компактной форме. int ama6r_c (real *a, integer *n, real *b, integer *nc, real *c) . ...
... ama5r.zip , ama5d.zip . ... Сложение симметричной матрицы и ленточной, заданных в компактной форме. ... Результатом этого сложения является квадратная матрица C порядка N. SUBROUTINE AMA5R (A, N, B, NLB, NUB, C) . ... сложение с повышенной точностью симметричной матрицы и ленточной, заданных в компактной форме. ... DIMENSION A(15), B(5, 4), C(5, 5) DATA A /15*1.0/ DATA B /2*0.0, 3*1.0, 0.0, 13*1.0, 0.0/ N = 5 NLB = 2 NUB = 1 CALL AMA5R (A, N, B, NLB, NUB, C) Результаты: | 2.0 2.0 1.0 1.0 1.0 | ...
... Сложение ленточной и ленточной симметричной матриц, заданных в компактной форме. ama4r_c вычисляет сумму ленточной матрицы А порядка N, заданной в компактной форме, и ленточной симметричной матрицы B того же порядка, заданной также в компактной форме. ... Матрица C запоминается в компактной форме. int ama4r_c (real *a, integer *n, integer *nla, integer *nua, real *b, integer *nc, real *c) . ... заданный порядок матриц A и B (тип: целый); . ... заданное число ко - диагоналей матрицы B (тип: целый); ....
... В подпрограмме ama2d_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[8] /* was [4][2] */ = { 0.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f }; /* Local variables */ extern int ama2r_c(float *, int *, float *, int *, float *); static float c__[16] /* was [4][4] */; static int n, nc, i__; n = 4; nc = 1; ama2r_c(a, , b, , c__); for (i__ = 0; i__ ...