Умножение треугольных матриц«Методы обработки разреженных матриц»

@8VAid8 · Регистрация: 21.04.2012

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

Нужно перемножить треугольные матрицы в обычном виде и в свёрнутом. С обычным проблем нет.
Доступ к элементам свёрнутой матрицы осуществляется по формуле : k= j*(j-1)/2+i. Подстановкой этой формулы в алгоритм умножения обычных матриц пытаюсь добится результата, но он не совпадает. Возможно потому что обычные матрицы также ниже главной диагонали заполнены нулями...
Каким еще способом можно перемножить свёрнутые разряженные матрицы? Или где у меня ошибки?

C++

#include "stdafx.h"
#include "targetver.h"
#include <stdlib.h>
#include <iostream>
#include <iomanip>
#include <time.h>
#define n 4
using namespace std;
int main()
{ 
        setlocale(LC_ALL,"Russian"); //русская кодировка
        int a[n][n]; //начальная матрица (1)
        int b[n][n]; //начальная матрица (2)
        const int m=n*(n+1)/2;//размерность свёрнутой матрицы
        //int a1[m];   //начальнай свёрнутая матрица (1)
        //int b1[m];   //начальнай свёрнутая матрица (2)
        //int c[n][n]; //произведение (1) и (2)
        //int c1[m];   //произведение свёрнутых (1) и (2)
        clock_t t1;  //переменная таймера
        clock_t t2;
 
        //Матрица 1 
        for(int i=0;i<n;i++)
        {
             for(int j=0;j<n;j++)
             {
                if(j<i)a[i][j]=0;
                else
                {
                    a[i][j]=rand()%10;
                }
             }
        }
        //-------------Свёрнутая матрица 1---------
        int a1[m];   //начальнай свёрнутая матрица (1)
        int k1=0;
        for(int j=0;j<n;j++)
        {
            for(int i=0;i<n;i++)
            {
                if(i<=j)
                {
                  a1[k1]=a[i][j];k1++;
                }
            
            }
        }
 
            //Вывод матрицы 1
        //--------------------------------
        cout<<"Матрица 1\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                cout<<a[i][j]<<"  ";
            }
            cout<<endl;
        }
        cout<<"\n---------------------------------\n";
        //------------------------------------------
        //Вывод свёрнутой матрицы 1
        cout<<"\nСвёрнутая матрица 1\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<m;i++)
        {
            cout<<a1[i]<<" ";
        }
            cout<<"\n---------------------------------\n";
            cout<<"\n\n\n";
        //------------------------------------------
            //Матрица 2
        for(int i=0;i<n;i++)
        {
             for(int j=0;j<n;j++)
             {
                if(j<i) b[i][j]=0;
                else
                {
                    b[i][j]=rand()%10;
                }
             }
        }
        //-------------Свёрнутая матрица 2---------
        int b1[m];   //начальнай свёрнутая матрица (2)
        int k2=0;
        for(int j=0;j<n;j++)
        {
            for(int i=0;i<n;i++)
            {
                if(i<=j)
                {
                  b1[k2]=b[i][j];k2++;
                }
            
            }
        }
            //Вывод матрицы 2   
        //-------------------------------------------
        cout<<"\nМатрица 2\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                cout<<b[i][j]<<"  ";
            }
            cout<<endl;
        }
        cout<<"\n---------------------------------\n";
        //-------------------------------------------
        //Вывод свёрнутой матрицы 2
        cout<<"\nСвёрнутая матрица 2\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<m;i++)
        {
            cout<<b1[i]<<" ";
        }
            cout<<"\n---------------------------------\n";
            cout<<"\n\n\n";
        //-----------------------------------------
 
 
///---------------------умножение простых матриц----------------------
 
            int c[n][n]; //произведение (1) и (2)
            for(int q = 0;q<n;q++)
            {
                    for(int w = 0; w<n;w++)
                    {
                        c[q][w]=0;
                    }
            }
    
 
   for (int i = 0; i < n; i++) 
    {
        for (int j = 0; j < n; j++)
        {
            // Multiply the row of A by the column of B to get the row, column of product.
            for (int inner = 0; inner < n; inner++)
            {
                c[i][j] += a[i][inner] * b[inner][j];
            }
        }
     }
 
 
 
  //--------------------Вывод результата умножения простых матриц------
 
 //-------------------------------------------
        cout<<"\nРезультат умножения простых матриц\n";
        cout<<"\n---------------------------------\n";
        for(int i = 0; i < n; i++)
                            {
                                for(int j = 0; j < n; j++)
                                    {
                                        cout << setw(3) << c[i][j];
                                    }
 
                                cout << "\n";
                            }
        cout<<"\n---------------------------------\n";
        //-------------------------------------------
 
        ///---------------------умножение свёрнутых матриц----------------------
 
        //-------------------------------------------k=j*(j-1)/2 +i
 
//     inn*(inn-1)/2 +i
 
//     j*(j-1)/2 +inn
int c1[m];   //произведение свёрнутых (1) и (2)
for(int q = 0;q<m;q++){c1[q]=0;}
   
 
  for (int i = 0; i < n; i++) 
   {
        for (int j = 0; j < n; j++) 
        {
            // Multiply the row of A by the column of B to get the row, column of product.
            for (int inn = 0; inn < n; inn++)
            {
                //if(i<=j)
                    c1[j*(j-1)/2 +i] += a1[inn*(inn-1)/2 +i] * b1[j*(j-1)/2 +inn];
            }
 
        }
        
    }
 
    //-------------------------------------------
//--------------------Вывод результата умножения свёрнутых матриц------
 
        cout<<"\nРезультат умножения свёрнутых матриц\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<m;i++)
        {
            cout<<c1[i]<<" ";
        }
            cout<<"\n---------------------------------\n";
        //-------------------------------------------
 
        
 
        system("pause");
}

@8VAid8 · 27.11.2012, 22:46 **[ТС]**

Имеет ли смысл с точки зрения производительности преобразовывать разряженную матрицу в одномерный массив, также заполненный и нулями?

Новые блоги и статьи
Это работает. Скорость асинхронной логики велика. Вопрос видимо останется в стабильности. Плата - огонь! Hrethgir 13.01.2025 По прошлому проекту в Logisim Evolution прилагаю файл архива проекта Gowin Eda и снимок. Восьмибитный счётчик из сумматора+ генератор сигнала согласования+ стартер fast регистров. Файлы прилагаю. . . .	UserScript для подсветки кнопок языков программирования в зависимости от текущего раздела volvo 13.01.2025 В результате работы этого скрипта подсвечиваются нужные кнопки не только в форме быстрого ответа, но и при редактировании сообщения: / / ==UserScript== / / @name CF_DefaultLangSelect / / . . .	Введение в модели и алгоритмы машинного обучения InfoMaster 12.01.2025 Машинное обучение представляет собой одну из наиболее динамично развивающихся областей искусственного интеллекта, которая фокусируется на разработке алгоритмов и методов, позволяющих компьютерам. . .	Как на Python создать нейросеть для решения задач InfoMaster 12.01.2025 В контексте стремительного развития современных технологий особое внимание уделяется таким инструментам, как нейросети. Эти структуры, вдохновленные биологическими нейронными сетями, используются для. . .	Как создать нейросеть для генерации картинок на Python InfoMaster 12.01.2025 Генерация изображений с помощью искусственных нейронных сетей стала одним из наиболее захватывающих направлений в области компьютерного зрения и машинного обучения. В этой статье мы рассмотрим. . .	Создание нейросети для генерации текста на Python InfoMaster 12.01.2025 Нейросети, или искусственные нейронные сети, представляют собой модели машинного обучения, вдохновленные работой человеческого мозга. Они состоят из множества взаимосвязанных узлов, или "нейронов",. . .
Как создать нейросеть распознавания изображений на Python InfoMaster 12.01.2025 Введение в распознавание изображений с помощью нейросетей Распознавание изображений с помощью нейронных сетей стало одним из самых впечатляющих достижений в области искусственного интеллекта. Эта. . .	Основы искуственного интеллекта InfoMaster 12.01.2025 Искусственный интеллект (ИИ) представляет собой одну из наиболее динамично развивающихся областей современной науки и технологий. В широком смысле под искусственным интеллектом понимается способность. . .	Python и нейросети InfoMaster 12.01.2025 Искусственные нейронные сети стали неотъемлемой частью современных технологий, революционизировав множество областей - от медицинской диагностики до автономных транспортных средств. Python, благодаря. . .	Python в машинном обучении InfoMaster 12.01.2025 Python стал неотъемлемой частью современного машинного обучения, завоевав позицию ведущего языка программирования в этой области. Его популярность обусловлена несколькими ключевыми факторами, которые. . .	Создание UI на Python с TKinter InfoMaster 12.01.2025 TKinter — это одна из наиболее популярных библиотек для создания графических интерфейсов пользователей (GUI) в языке программирования Python. TKinter входит в стандартную библиотеку Python, что. . .	HTML5 в разработке мобильных приложений InfoMaster 12.01.2025 Введение: Обзор роли HTML5 в мобильной разработке В современном мире мобильных технологий HTML5 стал ключевым инструментом для разработки кроссплатформенных приложений. Эта технология произвела. . .

@8VAid8 14 / 15 / 2 Регистрация: 21.04.2012 Сообщений: 192
	27.11.2012, 22:46 [ТС]	2
	Имеет ли смысл с точки зрения производительности преобразовывать разряженную матрицу в одномерный массив, также заполненный и нулями? 0