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

@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;
}

Новые блоги и статьи Все статьи Все блоги /
Реалии 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 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.

@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