Умножение матриц. SSE2

@fedosenk0 · Регистрация: 25.10.2017

Студворк — интернет-сервис помощи студентам

Есть 2 матрицы 1000х1000, элементами которых являются матрицы float 12х12. Как с помощью команд SSE2 перемножить матрицы и записать результат в 3-ю матрицу? Транспонировать матрицы запрещено.
Вот как я это делаю, но что-то не правильно.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
__m128 aaa;
    __m128 bbb;
    __m128 ccc;
    for (int i = 0; i < 500; i++) {
        for (int k = 0; k < 500; k++) {
            for (int j = 0; j < 500; j++) {
                temp = q[i][j];
                temp1 = w[j][k];
                temp2 = e[i][k];
 
                for (int m = 0; m < 12; m++) {
                    for (int l = 0; l < 12; l += 4) {
                        aaa = _mm_load_ps(&temp[m][l]);
                        bbb = _mm_load_ps(&temp1[m][l]);
                        ccc = _mm_mul_ps(aaa, bbb);
                        temp2[m][l] = ccc.m128_f32[l];
 
                    }
                }
                e[i][k] = temp2;
            }
        }
    }

@nonedark2008 · 12.02.2019, 23:31

fedosenk0, посмотри как умножаются блочные матрицы. Ты явно где-то растерял все свои сложения.
Раз транспонировать нельзя, то стоит глянуть на _mm_load1_pd и _mm_add_ps, да и вообще изучи гайд.

Новые блоги и статьи Все статьи Все блоги /
Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .
Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.

@nonedark2008 1394 / 1023 / 325 Регистрация: 28.07.2012 Сообщений: 2,813
	12.02.2019, 23:31
	fedosenk0, посмотри как умножаются блочные матрицы. Ты явно где-то растерял все свои сложения. Раз транспонировать нельзя, то стоит глянуть на `_mm_load1_pd` и `_mm_add_ps`, да и вообще изучи гайд. 0