Матрицы - перегрузка оператора *

@Fedor N · Регистрация: 22.06.2012

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

Здравствуйте, помогите пожалуйста. Имеется библиотека Matrixclass.h, необходимо в ней перегрузить оператор умножения. Как только я не пытался его реализовать, но выводит не понять что. Посоветуйте что-нибудь по этому поводу.

Matrixclass.h

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
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
#include <iostream>
#include <stdlib.h>
#include <time.h> // Srand, Randomize
#include <math.h>
using namespace std;
 
enum ElemStatus {minimal = -5000, maximal};
enum Action {sum = -6000, prod};
enum ConvertType {row, column};
 
class Matrix {
    private:
        // Размеры массива
        unsigned xsize, ysize;
 
        // Указатель на начало массива в динамической памяти
        double *Array;
 
        // Флаг временности
        int temporary;
    public:
 
        // Возвратить размеры матрицы
        unsigned getxsize(void) {
            return xsize;
        }
 
        unsigned getysize(void) {
            return ysize;
        }
 
        // Функция вычисления определителя
        double determinant(Matrix& M, unsigned size);
 
        // конструктор, заполняет матрицу
        Matrix(unsigned y, unsigned x);
        Matrix(unsigned y, unsigned x, double v);
 
        // конструктор копирования
        Matrix(Matrix &);
 
        // Конструктор "скалярной матрицы" xsize=ysize=1
        Matrix(double v);
 
        // деструктор, вывод "!" для отслеживания вызовов
        ~Matrix() {
//          delete [] Array;
//          cout << "!";
        }
 
        void add(unsigned x, unsigned y, double v);
 
        // "Примитивный" вывод
        void show(void);
 
        // Последовательный доступ к элементу
        double & operator()(unsigned n) const;
        Matrix & operator*(Matrix& M);
        Matrix & operator()(unsigned x1, unsigned x2, unsigned y1, unsigned y2) const;
};
 
double Matrix::determinant(Matrix& M, unsigned size) {
    double det = 0, sum1, sum2;
 
    // Если размерость 1, то первый элемент есть определитель
    if(size == 1)
        det = M(0);
    // Иначе если два элемента
    else if(size == 2)
        det = M(0) * M(3) - M(1) * M(2);
    else {
        // находим детерминант
        for(unsigned i = 0; i < size; i++) {
            sum2 = sum1 = 1;
            int l = 2 * size - 1 - i;
 
            for(unsigned j = 0; j < size; j++) {
                sum2 *= M(j * size + l % size);
                l--;
                sum1 *= M(j * size + ((j + i) % size));
            }
 
            det += sum1 - sum2;
        }
    }
 
    return det;
}
 
Matrix::Matrix(unsigned x, unsigned y) {
    Array = new double[(xsize = x)*(ysize = y)];
 
    temporary = 0;
    for(unsigned i = 0; i < xsize; i++)
        for(unsigned j = 0; j < ysize; j++)
//      cout << "Enter " << i + 1 << " elemet ";
//      cout << (v = rand() % 100) << endl;
            Array[i*x+j] = rand() % 100;
 
}
 
void Matrix::add(unsigned x, unsigned y, double v) {
    Array[x*xsize+y] = v;
}
 
void Matrix::show(void) {
    cout << endl;
    //cout << "{";
    for(int i = 0; i < xsize; i++) {
        cout << "{";
        for(int j = 0; j < ysize; j++)
            if(j != ysize - 1)
                cout << Array[i*xsize+j] << ", ";
            else if(j == 0)
                cout << Array[i*xsize+j];
            else
                cout << Array[i*xsize+j];
 
        cout << "}";
        if(i != xsize - 1)
            cout << "," << endl;
    }
 
    cout << " " << endl;
 
    if(temporary)
        delete this;
}
 
Matrix& Matrix::operator()(unsigned x1, unsigned x2, unsigned y1, unsigned y2) const {
    if((y2 < y1) || (x2 < x1) || (y2 > ysize) || (x2 > xsize)) {
        cout << "error";
        exit(1);
    }
 
    Matrix*A = new Matrix(y2 - y1 + 1, x2 - x1 + 1);
    A->temporary = 1;
    int i = 0;
 
    for(int y = y1; y <= y2; y++)
        for(int x = x1; x <= x2; x++)
            A->Array[i++] = Array[(y-1)*xsize+(x-1)];
 
    if(temporary)
        delete this;
 
    return *A;
}
 
double &Matrix::operator()(unsigned n) const {
    if(n > xsize * ysize) {
//  cout << n << " > " << xsize << " * " << ysize << endl;
        cerr << "Index exceeds matrix dimensions!";
        exit(1);
    }
 
    if(temporary) {
        cerr << "Operation is not allowed!";
        exit(1);
    }
 
//  cout << n << " > " << xsize << " * " << ysize << endl;
 
    return Array[n];
}
 
Matrix &Matrix::operator*(Matrix &M) const //умножение матриц
 
{
    Matrix result(*this);
 
    for (int i = 0; i < 3; i++)
     {
        for (int j = 0; j < 3; j++)
         {
            result.A[i][j] *= n.A[i][j];
         }
     }
    return result;
}
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
{
   /** Matrix result(*this);
    //Matrix T(ysize, M.xsize);
     double temp=0;
    if(xsize == M.ysize) {
        //реализация алгоритма перемножения матриц
         for(int j=0;j<M.ysize;j++)
        {
            for(int i = 0; i < M.xsize; i++)
            {
                temp+=Array[i][j]*M(i*M.xsize+j);
                //result.Array[i][j] *= M.Array[i][j];
            }
          Array[j]=temp;
          temp=0;
        }
                    /**for(int j = 0; j < M.ysize; j++)
                        for(int k = 0; k < xsize; k++)
                            Array[i *j] += Array[j * k] * M(k * i);**/
 
 
    /** for(unsigned a = 0; a < ysize; a++)
            for(unsigned b = 0; b < M.xsize; b++) {
                T(a * M.xsize + b) = 0;
 
                for(unsigned k = 0; k < M.xsize; k++) {
                    cout << T(a * M.ysize + b) << " = " << Array[b * ysize + k] << " * " << M(k + b) << endl;
                    T(a* M.ysize + b) += Array[b * ysize + k] * M(k + b);
 
//                  M(a * M.xsize + b) += Array[b * ysize + k] * M(k * M.ysize + a);
                }
            }**/
    } else
        cout << "Wrong matrix size" << endl;
 
    //return T;
    return result;
}

main

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
#include <iostream>
#include <time.h> // Srand, Randomize
#include "Matrixclass_.h"
#include <conio.h>
using namespace std;
 
int main(int argc, char **argv) {
    cout << "program begin" << endl;
    srand((unsigned) time(NULL)); // Randomize
 
    unsigned n;
 
    cout << "Enter size matrix: ";
    cout << (n = rand() % 1 + 3) << endl;
 
    Matrix A(n, n);
    A.show();
 
    cout << endl << "Output determinant of matrix " << endl;
    cout << A.determinant(A, n) << endl << endl;
 
 
    cout << "Output matrix free elements" << endl;
    Matrix B(n, 1);
 
    B.show();
 
    B*A;
 
    B.show();
 
    cout << "program end" << endl;
    getch();
}

@Jtalk · 26.06.2012, 16:24

Вы возвращаете из оператора ссылку на статический объект (который удалится после возврата из функции).

silent_1991 · 26.06.2012, 18:05

Сообщение от Jtalk

ссылку на статический объект

На локальный. Статический как раз продолжит жить до завершения программы.

IGPIGP · 26.06.2012, 18:41

Сообщение от Fedor N

C++
1
2
3
4
5
6
7
8
Matrix &Matrix::operator*(Matrix &M) const //умножение матриц
 { Matrix result(*this);
 for (int i = 0; i < 3; i++) 
{ for (int j = 0; j < 3; j++)
 { result.A[i][j] *= n.A[i][j];
 } 
} return result; 
}

C++
1
2
3
4
5
6
7
8
9
10
Matrix &Matrix::operator*(Matrix &M) const //зачем const ?
 { Matrix result(*this);//доступ к членам объекта слева от * прямо по A[i][j] т.е. this уже здесь:)
//а временный объект бы пригодился:  Matrix tmpObj;
 
 for (int i = 0; i < 3; i++) 
{ for (int j = 0; j < 3; j++)
 { result.A[i][j] *= n.A[i][j];//что за n? передано же Matrix &M, тогда =M.A[i][j];
 } 
} return result; 
}

а если:

C++
1
2
3
4
5
6
7
8
9
10
Matrix &Matrix::operator*(Matrix &M) //зачем const ?
 { 
Matrix result;
 
 for (int i = 0; i < 3; i++) 
{ for (int j = 0; j < 3; j++)
 {result.A[i][j]= A[i][j]*M.A[i][j];
 } 
} return result; 
}

И еще, без перегрузки оператора "=" как это будет работать?
Память под массив выделяете динамически, а определения конструктора копирования не увидел.
Обратите внимание. Извините, за краткость. Что успел, - заметил.

@Fedor N · 28.06.2012, 15:01 **[ТС]**

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
Matrix &Matrix::operator*(Matrix &M)
{
    Matrix*res=new Matrix( this->ysize, M.xsize);
    for (int i=0; i<ysize; i++)
{
    for (int j=0; j<M.xsize;i++)
{
    res->Array(i*xsize+j)=0;
    for(int k=0;k<xsize; k++)
{
    res->Array(i*xsize+j)+=Array(i*xsize+k)*M.Array(k*M.xsize+j);
}
}
}
  return res;
}

Добавлено через 3 минуты
А если так. Компилятор ругается.

@Ksan · 28.06.2012, 15:06

Fedor N, а зачем вы используете одномерный массив, если это матрица?

@Fedor N · 28.06.2012, 15:13 **[ТС]**

В том то и дело, что нужен одномерный, иначе не передашь его функциям.

silent_1991 · 28.06.2012, 15:31

Fedor N, двумерный динамический массив нельзя передать в функцию?

C++
1
void func(type **array, size_t rows, size_t cols);

Нэ?

@Fedor N · 28.06.2012, 15:38 **[ТС]**

А разве одномерный массив вообще нельзя использовать? Мне преподаватель на этот счет ничего не говорил.

silent_1991 · 28.06.2012, 15:44

Fedor N, можно использовать всё, что позволяет выполнить задание. Человек удивился, почему не используется двумерный массив, в то время как это было бы логичнее с точки зрения человека (матрица - таблица, двумерный массив - тоже). Мой же комментарий был просто ответом на то, что двумерный динамический массив не получится передать в функцию.

@DU · 28.06.2012, 15:50

стоит обратить внимание на более существенные вещи, которые сейчас ограничивают использование этой матрицы в правильном коде.
Все функции, которые не изменяют матрицу, должны быть константными.
Если в функцию передается матрица и она в ней не изменяется, то надо передавать в нее матрицу по константной ссылке, а не просто по ссылке.
Те операторы, которые есть сейчас должны возвращать матрицу по значению (объект), а не по ссылке или указателю.

@Fedor N · 29.06.2012, 14:55 **[ТС]**

Я понял, что можно по ссылке передать, просто я с двумерным массивом не умею работать.

Добавлено через 22 часа 59 минут

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
Matrix &Matrix::operator*(Matrix &M) {
 
    Matrix c(M.ysize, xsize);
 
    if (ysize != M.xsize) {
        cerr << "Размерности матриц не верны" << endl;
        exit(1);
    } else {
 
        int i, j, k;
 
        for (i = 0; i < M.ysize; ++i) {
            for (j = 0; j < ysize; ++j) {
                c(i * c.xsize + j) = 0;
                for (k = 0; k < M.xsize; ++k) {
                                       c(i * c.xsize + j) += M(i * xsize + k) * Array[k * M.ysize + j];
                    //cout << i * M.xsize + k;
                }
            }
        }
    }
 
    show(c);

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

@Fedor N · 01.07.2012, 19:54 **[ТС]**

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
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
#include <iostream>
#include <stdlib.h>
#include <conio.h>
using namespace std;
enum ElemStatus {minimal=-5000, maximal};
enum Action {sum=-6000, prod};
enum ConvertType{row, column};
class Matrix
{
    private:
        // Размеры массива
        unsigned xsize, ysize;
        // Указатель на начало массива в динамической памяти
        double *Array;
      // Флаг временности
      int temporary;
    public:
        // Возвратить размеры матрицы
        unsigned getxsize(void) {return xsize;}
        unsigned getysize(void) {return ysize;}
 
        // конструктор, заполняет матрицу значением v
        Matrix(unsigned y = 1, unsigned x = 1, double v = 0);
        // конструктор копирования
        Matrix(Matrix &);
      // Конструктор "скалярной матрицы" xsize=ysize=1
        Matrix(double v);
        // деструктор, вывод "!" для отслеживания вызовов
        ~Matrix(){delete []Array; cout << "!";}
        // "Примитивный" вывод
        void show(void);
      // Последовательный доступ к элементу
      double & operator()(unsigned n) const;
        Matrix operator *(const Matrix&);
        // ЗАДАНИЕ1
      // Параллельный доступ
      double & operator()(unsigned y, unsigned x) const;
    /** // ЗАДАНИЕ2
      // Присваивание матрицы
        Matrix & operator=(Matrix &M);
        // Присваивание скаляра
        Matrix & operator=(double v);
      // Присваивание массива
      Matrix & operator=(double v[]);
        // ЗАДАНИЕ3
      // Аддитивные операции
        Matrix & operator+(const Matrix& M) const;
        Matrix & operator-(const Matrix& M) const;
        Matrix & operator+(double v) const;
        Matrix & operator-(double v) const;
        // ЗАДАНИЕ4
      // Аддитивные операции
      friend Matrix & operator+(double v, const Matrix& M);
        friend Matrix & operator-(double v, const Matrix& M);
        // ЗАДАНИЕ5
      // Умножение матрицы на матрицу
      Matrix & Matrix::operator*(const Matrix& M) const;
        // ЗАДАНИЕ6
      // Умножение матрицы на скаляр
      Matrix & Matrix::operator*(double v) const;
      // Умножение скаляра на матрицу
      friend Matrix & operator*(double v, const Matrix& M);
      // ЗАДАНИЕ7
      // Сравнение матриц
        Matrix & operator==(const Matrix& M) const;
        Matrix & operator>(const Matrix& M) const;
        Matrix & operator<(const Matrix& M) const;
      // ЗАДАНИЕ8
      // Сравнение матрицы со скаляром
      Matrix & operator==(double v) const;
        Matrix & operator>(double v) const;
        Matrix & operator<(double v) const;
        // ЗАДАНИЕ9
      // Сравнение скаляра с матрицей
      friend Matrix & operator==(double v, const Matrix& M);
      friend Matrix & operator>(double v, const Matrix& M);
      friend Matrix & operator<(double v, const Matrix& M);
        // ЗАДАНИЕ10
      // Перегрузка унарный операций
      Matrix & operator!() const;
        Matrix & operator-() const;
        // ЗАДАНИЕ11
      // Перегрузка потокового вывода
      friend ostream & operator << (ostream & out, const Matrix & M);
      // ЗАДАНИЕ12
      // Перегрузка потокового ввода
      friend istream & operator >> (istream & in, const Matrix & M);
      // ЗАДАНИЕ13
      // Транспонирование матрицы
      Matrix & operator~() const;
      // ЗАДАНИЕ14
      // Вертикальное объединение матриц
        Matrix & operator % (Matrix &M) const;
      // ЗАДАНИЕ15
      // Горизонтальное объединение матриц
        Matrix & Matrix::operator || (Matrix &M) const;
        // ЗАДАНИЕ16
      // Нахождение минимума, максимума
        double & operator [](ElemStatus es)  const;
        // ЗАДАНИЕ17
      // Нахождение суммы, произведения
        double operator [](Action act)  const;
      // ЗАДАНИЕ18
      // Преобразование в строку, столбец
        Matrix & operator [](ConvertType ct)  const;
      // ЗАДАНИЕ19
      // нахождение ближайшего по значению элемента
        double&  operator [](double v)  const;
        // ЗАДАНИЕ 20
      // Вырезка прямоугольного фрагмента
      Matrix&  operator ()(unsigned y1, unsigned y2,
                           unsigned x1, unsigned x2)  const;
**/
};
 
Matrix::Matrix(unsigned y, unsigned x, double v)
{
    unsigned asize = x*y;
    Array = new double[asize];
    xsize = x;
    ysize = y;
   temporary = 0;
    for(int i = 0; i<asize; i++) Array[i] = v;
}
 
Matrix::Matrix(Matrix & M)
{
    xsize = M.xsize;
    ysize = M.ysize;
    unsigned asize = xsize*ysize;
    Array = new double[asize];
    for(int i = 0; i<asize; i++) Array[i] = M.Array[i];
   temporary = 0;
    if (M.temporary) delete &M;
}
 
Matrix::Matrix(double v)
{
    xsize = ysize = 1;
   Array = new double;
    *Array = v;
   temporary = 0;
}
 
void Matrix::show(void)
{
    cout << '\n';
    for (int i = 0; i<ysize; i++)
   {
        for (int j = 0; j<xsize; j++)
                cout << Array[i*xsize+j]<< '\t';
        cout << '\n';
   }
    if (temporary) delete this;
}
double & Matrix::operator()(unsigned y, unsigned x) const
{
   if (y>ysize-1 || x>xsize-1)         // Проверка индексов массива. Они не должны быть больше максимального значения
   {
    cerr << "Index exceeds matrix dimensions!";
    exit(1);
   }
   if (temporary)                     // Если объект временный, то операция не выполняется
   {
    cerr << "Missing operator!";
    exit(1);
   }
   double *elem;                       // Создание указателя типа double
   elem = &Array[x*xsize+y];           // Настройка указателя на нужный элемент
   return *elem;
}
 
 
Matrix Matrix::operator *(const Matrix& M)
{       Matrix tmp(xsize,ysize);
        for(int i = 0; i < xsize; i++)
            for(int j = 0; j < ysize; j++)
        {    tmp(i,j) = 0;
                for(int k = 0; k < xsize; k++)
        tmp(i,j) +=M(i,k)*M(k,j);
        }
 
        return tmp;
}

Здравствуйте, привожу код библиотеки Matrixclass.h но хочу заострить внимание на операторе умножения, если компилировать в таком варианте, то в результате матрицы В перемножаться. Я же хочу перемножить А и В матрицы. И тут возникает вопрос. Как?) пробую повсякому использовать указатель this,

C++
1
 tmp(i,j) +=this->M(i,k)*M(k,j);

- думаю, что нужно сделать это так, но в итоге: c:\program files (x86)\codeblocks\mingw\bin\..\lib\gcc\mi ngw32\4.4.1\..\..\..\..\include\Matrixcl ass.h|180|error: 'class Matrix' has no member named 'tmp'|

silent_1991 · 01.07.2012, 20:04

C++
1
tmp(i,j) +=M(i,k)*M(k,j);

-->

C++
1
tmp(i,j) +=Array[i][k]*M(k,j);

либо

C++
1
tmp(i,j) +=(*this)(i,k)*M(k,j);

либо

C++
1
tmp(i,j) +=this->operator()(i,k)*M(k,j);

@Fedor N · 01.07.2012, 20:24 **[ТС]**

Сообщение от silent_1991

Код C++
1 tmp(i,j) +=(*this)(i,k)*M(k,j);
либоКод C++
1 tmp(i,j) +=this->operator()(i,k)*M(k,j);

Вот эти варианты работают,теперь надо разобраться пчм некоторые на 1 больше , а другие на 1 меньше элементы матрицы)

silent_1991 · 01.07.2012, 20:36

Сообщение от Fedor N

Вот эти варианты работают

Аа, ну да, не заметил, что Array - одномерный. Тогда можно вместо первого варианта:

C++
1
tmp(i,j) +=Array[i * xsize + k]*M(k,j);

@Fedor N · 01.07.2012, 20:44 **[ТС]**

Сообщение от silent_1991

Код C++
1 tmp(i,j) +=Array[i * xsize + k]*M(k,j);

Теперь работает, но вот перемножают криво-все 3 варианта ( , а должно быть
ans =

32 33 34
56 60 64
56 60 64

silent_1991 · 01.07.2012, 20:52

Ошибок в коде много, и очень серьёзных, но я не заметил причины неверной работы. Покажите весь код (в том числе клиентскую часть).

@Fedor N · 01.07.2012, 21:02 **[ТС]**

Библиотека Matrixclass.h есть выше, преподаватель бы обиделся если бы ему сказали, что в его библиотеке много серьезных ошибок )) Я же пытаюсь оператор умножения по его библиотеке сделать. А так это вскоре будет целая программа для расчета СЛАУ, осталось только как раз это умножение доделать... Вот main.cpp

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
#include <Matrixclass.h>
#include <iostream>
#include <stdlib.h>
#include <conio.h>
 
using namespace std;
int main()
{
    Matrix A(3,3);
    A(0,0)=1;A(1,0)=2;A(2,0)=3;
    A(0,1)=4;A(1,1)=4;A(2,1)=4;
    A(0,2)=4; A(1,2)=4;A(2,2)=4;
    cout<<"A:";A.show();
    Matrix B(3,3);
    B(0,0)=2; B(1,0)=3; B(2,0)=4;
    B(0,1)=6;B(1,1)=6;B(2,1)=6;
    B(0,2)=6;B(1,2)=6;B(2,2)=6;
    cout<<"B:";B.show();
    cout<<"A*B:";cout<<endl;
    (A*B).show();
}

Да и в принципе вот и все с чем я сейчас работаю=). Потом другие функции появятся (они вроде норм работают)

silent_1991 · 01.07.2012, 21:37

Fedor N, в коде преподавателя большинство из этих ошибок отсутствует (хотя мне всё же не понятен глубинный смысл поля temporary и адовых инструкций delete &M; и delete this; (

)). Вы же сами себя запутали, в одних местах поменяв порядок индексов (y, x вместо x, y), в других нет. Таким образом, инициализируете матрицу вы в одном порядке следования индексов (неверном), перемножаете же в верном... Вот подправленный, отформатированный и вычищенный код:

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
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
#include <iostream>
#include <stdlib.h>
#include <conio.h>
using namespace std;
 
class Matrix
{
private:
    unsigned xsize;
    unsigned ysize;
    double *Array;
    
public:
    unsigned getxsize(void) const
    {
        return xsize;
    }
    
    unsigned getysize(void) const
    {
        return ysize;
    }
    
    Matrix(unsigned y = 1, unsigned x = 1, double v = 0);
    Matrix(const Matrix& M);
    Matrix(double v);
    
    ~Matrix()
    {
        delete [] Array;
        //cout << "!";
    }
    
    void show(void) const;
    
    double& operator()(unsigned y, unsigned x);
    double  operator()(unsigned y, unsigned x) const;
    
    Matrix operator*(const Matrix& M) const;
};
 
Matrix::Matrix(unsigned x, unsigned y, double v)
{
    unsigned asize = x * y;
    Array = new double[asize];
    
    xsize = x;
    ysize = y;
    
    for (unsigned i = 0; i < asize; i++)
        Array[i] = v;
}
 
Matrix::Matrix(const Matrix& M)
{
    xsize = M.xsize;
    ysize = M.ysize;
    
    unsigned asize = xsize * ysize;
    Array = new double[asize];
    
    for (unsigned i = 0; i < asize; i++)
        Array[i] = M.Array[i];
}
 
Matrix::Matrix(double v)
{
    xsize = ysize = 1;
    Array = new double;
    *Array = v;
}
 
void Matrix::show(void) const
{
    cout << '\n';
    
    for (unsigned i = 0; i < xsize; i++)
    {
        for (unsigned j = 0; j < ysize; j++)
            cout << (*this)(i, j) << '\t';
        
        cout << '\n';
    }
}
 
double& Matrix::operator()(unsigned x, unsigned y)
{
    if (y > ysize - 1 || x > xsize - 1)
    {
        cerr << "Index exceeds matrix dimensions!";
        exit(1);
    }
    
    return Array[x * xsize + y];
}
 
double Matrix::operator()(unsigned x, unsigned y) const
{
    if (y > ysize - 1 || x > xsize - 1)
    {
        cerr << "Index exceeds matrix dimensions!";
        exit(1);
    }
    
    return Array[x * xsize + y];
}
 
Matrix Matrix::operator*(const Matrix& M) const
{
    Matrix tmp(xsize, ysize);
    
    for (unsigned i = 0; i < xsize; i++)
    {
        for (unsigned j = 0; j < ysize; j++)
        {
            tmp(i, j) = 0;
            
            for (unsigned k = 0; k < xsize; k++)
                tmp(i, j) += (*this)(i, k) * M(k, j);
        }
    }
 
    return tmp;
}
 
int main()
{
    Matrix A(3, 3);
    
    A(0, 0) = 1;    A(0, 1) = 2;    A(0, 2) = 3;
    A(1, 0) = 4;    A(1, 1) = 4;    A(1, 2) = 4;
    A(2, 0) = 4;    A(2, 1) = 4;    A(2, 2) = 4;
    
    cout << "A:";
    A.show();
    
    Matrix B(3, 3);
    
    B(0, 0) = 2;    B(0, 1) = 3;    B(0, 2) = 4;
    B(1, 0) = 6;    B(1, 1) = 6;    B(1, 2) = 6;
    B(2, 0) = 6;    B(2, 1) = 6;    B(2, 2) = 6;
    
    cout << "B:";
    B.show();
    
    cout << "A*B:" << endl;
    (A * B).show();
    
    return 0;
}

Новые блоги и статьи Все статьи Все блоги /
Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .	Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Учёным и волонтёрам проекта «Einstein@home» удалось обнаружить четыре гамма-лучевых пульсара в джете Млечного Пути Programma_Boinc 01.01.2026 Учёным и волонтёрам проекта «Einstein@home» удалось обнаружить четыре гамма-лучевых пульсара в джете Млечного Пути Сочетание глобально распределённой вычислительной мощности и инновационных. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .
Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .	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. . .

@Jtalk 95 / 81 / 3 Регистрация: 13.05.2011 Сообщений: 279
	26.06.2012, 16:24
	Вы возвращаете из оператора ссылку на статический объект (который удалится после возврата из функции). 1

@Ksan 27 / 27 / 4 Регистрация: 02.11.2010 Сообщений: 370
	28.06.2012, 15:06
	Fedor N, а зачем вы используете одномерный массив, если это матрица? 0

@Fedor N 0 / 0 / 0 Регистрация: 22.06.2012 Сообщений: 17
	28.06.2012, 15:13 [ТС]
	В том то и дело, что нужен одномерный, иначе не передашь его функциям. 0

@Fedor N 0 / 0 / 0 Регистрация: 22.06.2012 Сообщений: 17
	28.06.2012, 15:38 [ТС]
	А разве одномерный массив вообще нельзя использовать? Мне преподаватель на этот счет ничего не говорил. 0

silent_1991 5058 / 3118 / 271 Регистрация: 11.11.2009 Сообщений: 7,044
	28.06.2012, 15:44
	Fedor N, можно использовать всё, что позволяет выполнить задание. Человек удивился, почему не используется двумерный массив, в то время как это было бы логичнее с точки зрения человека (матрица - таблица, двумерный массив - тоже). Мой же комментарий был просто ответом на то, что двумерный динамический массив не получится передать в функцию. 0

@DU 1500 / 1146 / 165 Регистрация: 05.12.2011 Сообщений: 2,279
	28.06.2012, 15:50
	стоит обратить внимание на более существенные вещи, которые сейчас ограничивают использование этой матрицы в правильном коде. Все функции, которые не изменяют матрицу, должны быть константными. Если в функцию передается матрица и она в ней не изменяется, то надо передавать в нее матрицу по константной ссылке, а не просто по ссылке. Те операторы, которые есть сейчас должны возвращать матрицу по значению (объект), а не по ссылке или указателю. 1

silent_1991 5058 / 3118 / 271 Регистрация: 11.11.2009 Сообщений: 7,044
	01.07.2012, 20:52
	Ошибок в коде много, и очень серьёзных, но я не заметил причины неверной работы. Покажите весь код (в том числе клиентскую часть). 0