Описание тега matrix-multiplication

Вопросы, связанные с умножением матриц, особенно реализация. Математические вопросы должны учитывать тег [linear-algebra].
5 ответов

Умножение матриц с потоками Java

Я пытаюсь создать программу на Java с потоками для умножения матриц. Это исходный код: import java.util.Random; public class MatrixTest { //Creating the matrix static int[][] mat = new int[3][3]; static int[][] mat2 = new int[3][3]; static int[][] r…
2 ответа

Матрица эффективного умножения матрицы (матрица грамм)

Я хочу размножаться B = A @ A.T в клочья. Очевидно, что ответом будет симметричная матрица (т.е. B[i, j] == B[j, i]). Тем не менее, мне не ясно, как использовать это легко, чтобы сократить время вычислений вдвое (путем вычисления только нижнего треу…
07 июн '18 в 04:50
1 ответ

R несоответствующие аргументы

Мои извинения за то, что задали вопрос, который задавался тысячу раз, но я все еще не могу решить свою проблему. Я пытаюсь запустить этот двойной цикл: rho=matrix(0,12,1) for (i in 1:k){ for (j in 1:n-k){ temp=((1/(sigma^2))*(1/(n-i)))*((as.matrix(v…
25 апр '16 в 16:49
1 ответ

Упрощенное умножение матриц - возможно ли это?

Есть ли способ упростить a=np.dot(a,b) так же, как вы пишете a=a+b как a+=b? (a,b оба np.array)
24 май '16 в 21:09
0 ответов

Как разбросать столбцы матрицы с использованием C и MPI при выполнении умножения матрицы на матрицу NxN

У меня есть проект, чтобы написать программу, которая выполняет умножение матриц на две матрицы NxN A и B. Я могу использовать любой алгоритм, но мы узнали об алгоритме Кэннона в классе, так что, вероятно, я буду его использовать. У меня проблемы с …
1 ответ

Ошибка переполнения стека при получении целого размера?

Как часть проекта кода класса, часть назначения состоит в создании алгоритма, который будет умножать две матрицы, используя метод Штрассена. Я создал алгоритмы, необходимые для этого, однако я получаю ошибку переполнения стека при вызове getSize(). …
24 фев '14 в 00:26
4 ответа

Как оптимизировать этот продукт из трех матриц в C++ для x86?

У меня есть ключевой алгоритм, в котором большая часть времени выполнения тратится на вычисление продукта с плотной матрицей: A*A'*Y, where: A is an m-by-n matrix, A' is its conjugate transpose, Y is an m-by-k matrix Typical characteristics: - k is …
1 ответ

Невозможно сохранить значения массива в ocaml

Я пишу простую программу для расчета умножения двух матриц. (* matrix multiplication *) let multiply_row_per_col (m1 : int array array) (m2 : int array array) = (* create matrix *) let nrows = get_num_rows m1 in let ncols = get_num_rows m2 in let ma…
17 авг '16 в 13:18
4 ответа

Как использовать этот код C для умножения двух матриц с использованием алгоритма Штрассена?

Я искал реализацию алгоритма Штрассена в C, и я нашел этот код в конце. Чтобы использовать multiply функция: void multiply(int n, matrix a, matrix b, matrix c, matrix d); который умножает две матрицы a, b и помещает результат в c (d является промежу…
02 мар '12 в 19:19
2 ответа

Умножение матрицы на PHP

Я написал алгоритм Matrix Chained Multiplication в C около 3 лет назад, и я только начал изучать PHP. вот алгоритм в C: #include<stdio.h> #include<limits.h> int MatrixChainOrder(int p[], int n) { int m[n][n]; int i, j, k, L, q; for (i = …
05 июн '14 в 12:55
1 ответ

Умножение матриц с использованием одномерных массивов в C

У меня были некоторые проблемы с умножением матрицы (в виде двух одномерных массивов) в C. Никакие переменные не могут быть созданы или удалены. Правильный результат должен быть: http://s32.postimg.org/nrkigh2l1/IMG_20160516_183316.jpg Циклы печати …
16 май '16 в 16:47
1 ответ

Эффективное умножение матрицы SSE NxN

Я пытаюсь реализовать SSE версия большой матрицы с умножением матрицы. Я ищу эффективный алгоритм, основанный на SIMD Реализации. Мой желаемый метод выглядит так: A(n x m) * B(m x k) = C(n x k) И все матрицы считаются выровненными 16-байтовыми масси…
01 фев '14 в 21:27
1 ответ

Почему умножение матриц быстрее с Repa, чем с hmatrix?

Интересно, что Data.Array.Repa на самом деле быстрее, чем hmatrix, что неожиданно, поскольку hmatrix реализован с использованием LAPACK. Это потому, что Repa использует распакованный тип? import Data.Array.Repa import Data.Array.Repa.Algorithms.Matr…
01 ноя '13 в 18:49
1 ответ

Матричное умножение со случайными числами больше 1000

Я пытаюсь получить две матрицы со случайно сгенерированными числами для умножения вместе, но я продолжаю получать ошибку сегментации и знаю, как правильно распределить память. любая помощь с этим была бы велика // вот код #include<iostream> #i…
29 фев '16 в 10:48
1 ответ

Обобщение векторного внешнего произведения: примените его к каждому столбцу матрицы

У меня есть матрица A = [x1, x2, ..., xm] где каждый xi является вектором столбца размера [n, 1], Так что А имеет форму [n, m], Я пытаюсь найти ковариационную матрицу каждого вектора столбца, чтобы, если в результате получилась другая матрица C, C и…
0 ответов

Умножение матрицы на чистый кватернион

У меня есть уравнение, в котором мне нужно умножить матрицу 3 x 3 на вектор V 3 x 1, который сохраняется как чистый кватернион Q = [0 V]. Простой способ сделать это - просто получить мой вектор V из Q и затем умножить M на V. Это наиболее логичный с…
01 янв '16 в 13:43
2 ответа

Как оптимизировать код умножения матриц (matmul) для быстрой работы на одном ядре процессора

Я работаю над концепциями параллельного программирования и пытаюсь оптимизировать пример умножения матриц на одном ядре. Самая быстрая реализация, которую я придумал, заключается в следующем: /* This routine performs a dgemm operation * C := C + A *…
1 ответ

Умножение подмножества матрицы на другую матрицу, используя lapack/blas

Я хочу умножить подмножество матрицы A на другую матрицу, используя dgemm или любую другую функцию lapack/blas. Я думаю, что, поскольку элементы подматрицы не могут быть смежными, я не могу использовать dgemm напрямую, не копируя подматрицу в другое…
26 сен '17 в 07:10
1 ответ

Производительность умножения матрицы Eigen3 зависит от процессора?

В последние месяцы я работал с умножением компьютерных матриц и провел несколько тестов с использованием openMP и eigen3. Испытания проводились на следующих машинах: Компьютер 1: Процессор Intel Core i7-3610QM @ 2,30 ГГц / 6 ГБ ddr3 Компьютер 2: Шес…
2 ответа

OpenGL glMatrixMode(GL_PROJECTION) против glMatrixMode(GL_MODELVIEW)

В чем разница между размещением glRotatef() после glMatrixMode(GL_PROJECTION); glMatrixMode(GL_PROJECTION); glLoadIdentity(); glRotatef(red, green, blue); и размещение glRotatef() после glMatrixMode(GL_MODELVIEW); glMatrixMode(GL_MODELVIEW); glLoadI…