Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.84/19: Рейтинг темы: голосов - 19, средняя оценка - 4.84
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313

Перемножение матриц.

21.07.2011, 11:08. Показов 3585. Ответов 5
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Нужен код для перемножения матрицы и столбца за минимально возможное время.
Порядок матрицы ( и столбца ) огромен - около 100000. Помогите найти оптимальный код!
Ищу уже 2 день... что-то безрезультатно как то(((
0
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
21.07.2011, 11:08
Ответы с готовыми решениями:

Возведение матриц в степени, перемножение и сложение матриц
Даны 3-ех элементные вещественные векторы x и y и квадратные матрицы A и B. Вычислить: (A 2 * x, y) + (B 2 * y, x) + (A * B * x, y) ...

Перемножение матриц
Доброго времени суток. Подскажите как разобраться с перемножением матриц, если есть 2 матрицы размерностью МхМ (М = 2000-5000),...

Перемножение 2-ух матриц
Помогите пожалуйста с перемножением 2-ух матриц? в долгу не останусь) Выделение памяти написал, а вот как использовать с умножением, и...

5
 Аватар для Mиxaил
542 / 447 / 162
Регистрация: 10.12.2009
Сообщений: 1,857
21.07.2011, 11:16
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
int multMatrixBlock(double *a, const int n, double *b, const int m, double *res, double *c, double *d, double*e){
    int n_block=n/m;
    int m_small;
    int i,j,k,o;
    m_small=n-n_block*m;
    clearMatrix(res, n*n);
    for ( i = 0; i <= n_block; i++)
    {
        for ( j = 0; j <= n_block; j++)
        {
            for ( k = 0; k <=n_block; k++)
            {
                if (i<n_block && j<n_block && k<n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixSq(c, m, d, e);
                    setSubMatrixAdded(res, n, e, m, i, j);
                }
                else if (i<n_block && j<n_block && k==n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m, d,m_small, m, e);
                    setSubMatrixAdded(res, n, e, m, i, j);  
                }
                else if (i==n_block && j<n_block && k<n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m_small, d,m, m, e);
                    setSubMatrixAdded(res, n, e, m, i, j);
                }
                else if (i==n_block && j<n_block && k==n_block) {   
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m_small, d,m_small, m, e);
                    setSubMatrixAdded(res, n, e, m, i, j);  
                }
                else if (i<n_block && j==n_block && k<n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m, d,m, m_small, e);
                    setSubMatrixAdded(res, n, e, m, i, j);
                }
                else if (i<n_block && j==n_block && k==n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m, d,m_small, m_small, e);
                    setSubMatrixAdded(res, n, e, m, i, j);
                }
                else if (i==n_block && j==n_block && k<n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m_small, d,m, m_small, e);
                    setSubMatrixAdded(res, n, e, m, i, j);  
                }
                else if (i==n_block && j==n_block && k==n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixSq(c, m_small, d, e);
                    setSubMatrixAdded(res, n, e, m, i, j);  
                }
            }
        }
    }
    return 0;
}
нашел в интернете

Добавлено через 1 минуту
матрица хранится в одномерном массиве.
0
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
21.07.2011, 11:21  [ТС]
Вот только там нет определния кучи ф-ций...
0
 Аватар для Mиxaил
542 / 447 / 162
Регистрация: 10.12.2009
Сообщений: 1,857
21.07.2011, 11:38
Цитата Сообщение от hello19 Посмотреть сообщение
Вот только там нет определния кучи ф-ций...
исходые функции не удалось найти, сорри... Зато есть еще вот это:
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
int multMatrixSq(double *a, const int n, double *b, double *c){
    double s1=0, s2=0, s3=0, s4=0;
    double *f, *g, *h;
    int i, j, k;
    if (n%2==0){
        for (i = 0; i < n-1; i+=2)
        {
            f=a+i*n;
            for (j = 0; j < n-1; j+=2)
            {
                g=b+j;
                for (k=0;k<n;k++){
                    s1+=f[k]*g[n*k];    
                    s2+=f[n+k]*g[n*k];  
                    s3+=f[k]*g[n*k+1];  
                    s4+=f[n+k]*g[n*k+1];    
                }
                h=c+i*n+j;
                h[0]=s1; s1=0;
                h[n]=s2; s2=0;
                h[1]=s3; s3=0;
                h[n+1]=s4; s4=0;
            }
        }
    }
    else {
        for (i = 0; i < n-1; i+=2)
            {
                f=a+i*n;
                for (j = 0; j < n-1; j+=2)
                {
                    g=b+j;
                    for (k=0;k<n;k++){
                        s1+=f[k]*g[n*k];
                        s2+=f[n+k]*g[n*k];
                        s3+=f[k]*g[n*k+1];
                        s4+=f[n+k]*g[n*k+1];
                    }
                    h=c+i*n+j;
                    h[0]=s1; s1=0;
                    h[n]=s2; s2=0;
                    h[1]=s3; s3=0;
                    h[n+1]=s4; s4=0;
                }
                if  (j==n-1){
                    g=b+j;
                    for (k=0;k<n;k++){
                        s1+=f[k]*g[n*k];    
                        s2+=f[n+k]*g[n*k];  
                    }
                    h=c+i*n+j;
                    h[0]=s1; s1=0;
                    h[n]=s2; s2=0;
                }
            }
            if (i==n-1){
                f=a+i*n;
                for (j = 0; j < n-1; j+=2)
                {
                    g=b+j;
                    for (k=0;k<n;k++){
                        s1+=f[k]*g[n*k];    
                        s3+=f[k]*g[n*k+1];  
                    }
                    h=c+i*n+j;
                    h[0]=s1; s1=0;
                    h[1]=s3; s3=0;
                }
                if  (j==n-1){
                    g=b+j;
                    for (k=0;k<n;k++){
                        s1+=f[k]*g[n*k];    
                    }
                    h=c+i*n+j;
                    h[0]=s1; s1=0;
                }
            }
        }
    return 0;
}
тут, вроде, без сторонних функций
1
Каратель
Эксперт С++
6610 / 4029 / 401
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
21.07.2011, 11:42

Не по теме:

Цитата Сообщение от hello19 Посмотреть сообщение
около 100000
лол) я представляю сколько веремени будут перемножаться такие матрицы кодом приведеным выше, а ТС рекомендую смирится с этим заданием, к сожалению это не студенческая задача


а по теме: вам нужно ленточное перемножение матриц
1
Freelance
Эксперт С++
 Аватар для asics
2891 / 1826 / 356
Регистрация: 09.09.2010
Сообщений: 3,841
21.07.2011, 11:49
Вот подобная тема -- Распараллеливание
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
21.07.2011, 11:49
Помогаю со студенческими работами здесь

Перемножение матриц
Доброго времени суток! У меня есть две матрицы: первая k x m и вторая u x s размерностей. Нужно их перемножить. Подскажите, пожалуйста, в...

Перемножение матриц
Как сделать перегрузку *= для матриц ?! Был такой вариант: Matrix&amp; Matrix::operator *= (const Matrix &amp;rhs) { Matrix...

Перемножение матриц
Помогите решить проблему суть задачи проверить матрицу бинарного отношения на свойства. Все сделал кроме Транзитивности, суть того, что...

Перемножение матриц
Даны две квадратные матрицы.Напишите программу перемножающую матрицы по правилам перемножения матриц. У меня есть программа. но мне...

Перемножение матриц
Помогите исправить ошибки. //--------------------------------------------------------------------------- #include &lt;iostream&gt; ...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
6
Ответ Создать тему
Новые блоги и статьи
PhpStorm 2025.3: WSL Terminal всегда стартует в ~
and_y87 14.12.2025
PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .
Access
VikBal 11.12.2025
Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.
Новый ноутбук
volvo 07.12.2025
Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro
Музыка, написанная Искусственным Интеллектом
volvo 04.12.2025
Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .
От async/await к виртуальным потокам в Python
IndentationError 23.11.2025
Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов
Argus19 22.11.2025
Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .
Сколько Государство потратило денег на меня, обеспечивая инсулином.
Programma_Boinc 20.11.2025
Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .
Ломающие изменения в C#.NStar Alpha
Etyuhibosecyu 20.11.2025
Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .
Мысли в слух
kumehtar 18.11.2025
Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .
Создание Single Page Application на фреймах
krapotkin 16.11.2025
Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2025, CyberForum.ru