Не совсем понятна постановка задачи: реализовать умножение матриц ленточным способом с распределением строк матрицы B. Если с распределением столбцов все легко-мы умножаем строку на столбец, то какой алгоритм тут хотят увидеть-я не представляю. Быть может, у кого-нибудь есть примеры кода?
1 ответ
Кода у меня нет (и не будет), но могу предложить идею.
Матрица B
разбивается построчно, а затем процедура скалярного умножения строк A
на столбцы B
реализуется параллельно через MPI.
Пусть у нас два процесса S
и T
. Тогда процесс S
получает первую половину всех строк B
, а процесс T
- вторую. Далее, мы должны перемножить каждую строку матрицы A
на каждый столбец матрицы B
При этом половина столбца находится в S
, а а вторая половина в T
. Что делать, если нам надо перемножить строку a
матрицы A
на столбец b
матрицы B
, который уже разбит на две части bS и bT?
Очень просто строку a
также разбиваем на две части aS и aT. aS идет в S
и умножается на bS, одновременно с этим вторая половина строки aT идёт в T
, и перемножается там на bT. Далее оба результата складываются, получается результирующий элемент итоговой матрицы. И так нужно повторить для каждой строки и каждого столбца исходных матриц.
Надеюсь, мне никогда не придётся кодировать подобную ерунду, это непрактично и неэффективно.
B
сначала транспонируется (есть такая форма оптимизации), либо вас просят реализовать достаточно сомнительный по эффективности и корявый по красоте алгоритм, просто чтобы проверить ваши способности.