Форум программистов, компьютерный форум CyberForum.ru

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

Восстановить пароль Регистрация
 
 
Рейтинг: Рейтинг темы: голосов - 17, средняя оценка - 4.88
Fedor N
 Аватар для Fedor N
0 / 0 / 0
Регистрация: 22.06.2012
Сообщений: 17
26.06.2012, 15:38     Матрицы - перегрузка оператора * #1
Здравствуйте, помогите пожалуйста. Имеется библиотека 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();
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
26.06.2012, 15:38     Матрицы - перегрузка оператора *
Посмотрите здесь:

Перегрузка оператора C++
Перегрузка оператора == C++
Перегрузка оператора - C++
Перегрузка оператора *= C++
Перегрузка оператора + для динамической матрицы C++
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Jtalk
93 / 79 / 4
Регистрация: 13.05.2011
Сообщений: 279
26.06.2012, 16:24     Матрицы - перегрузка оператора * #2
Вы возвращаете из оператора ссылку на статический объект (который удалится после возврата из функции).
silent_1991
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
26.06.2012, 18:05     Матрицы - перегрузка оператора * #3
Цитата Сообщение от Jtalk Посмотреть сообщение
ссылку на статический объект
На локальный. Статический как раз продолжит жить до завершения программы.
IGPIGP
Комп_Оратор)
 Аватар для IGPIGP
6172 / 2901 / 284
Регистрация: 04.12.2011
Сообщений: 7,721
Записей в блоге: 3
26.06.2012, 18:41     Матрицы - перегрузка оператора * #4
Цитата Сообщение от 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
 Аватар для Fedor N
0 / 0 / 0
Регистрация: 22.06.2012
Сообщений: 17
28.06.2012, 15:01  [ТС]     Матрицы - перегрузка оператора * #5
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
26 / 26 / 0
Регистрация: 02.11.2010
Сообщений: 370
28.06.2012, 15:06     Матрицы - перегрузка оператора * #6
Fedor N, а зачем вы используете одномерный массив, если это матрица?
Fedor N
 Аватар для Fedor N
0 / 0 / 0
Регистрация: 22.06.2012
Сообщений: 17
28.06.2012, 15:13  [ТС]     Матрицы - перегрузка оператора * #7
В том то и дело, что нужен одномерный, иначе не передашь его функциям.
silent_1991
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
28.06.2012, 15:31     Матрицы - перегрузка оператора * #8
Fedor N, двумерный динамический массив нельзя передать в функцию?
C++
1
void func(type **array, size_t rows, size_t cols);
Нэ?
Fedor N
 Аватар для Fedor N
0 / 0 / 0
Регистрация: 22.06.2012
Сообщений: 17
28.06.2012, 15:38  [ТС]     Матрицы - перегрузка оператора * #9
А разве одномерный массив вообще нельзя использовать? Мне преподаватель на этот счет ничего не говорил.
silent_1991
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
28.06.2012, 15:44     Матрицы - перегрузка оператора * #10
Fedor N, можно использовать всё, что позволяет выполнить задание. Человек удивился, почему не используется двумерный массив, в то время как это было бы логичнее с точки зрения человека (матрица - таблица, двумерный массив - тоже). Мой же комментарий был просто ответом на то, что двумерный динамический массив не получится передать в функцию.
DU
1477 / 1053 / 45
Регистрация: 05.12.2011
Сообщений: 2,279
28.06.2012, 15:50     Матрицы - перегрузка оператора * #11
стоит обратить внимание на более существенные вещи, которые сейчас ограничивают использование этой матрицы в правильном коде.
Все функции, которые не изменяют матрицу, должны быть константными.
Если в функцию передается матрица и она в ней не изменяется, то надо передавать в нее матрицу по константной ссылке, а не просто по ссылке.
Те операторы, которые есть сейчас должны возвращать матрицу по значению (объект), а не по ссылке или указателю.
Fedor N
 Аватар для Fedor N
0 / 0 / 0
Регистрация: 22.06.2012
Сообщений: 17
29.06.2012, 14:55  [ТС]     Матрицы - перегрузка оператора * #12
Я понял, что можно по ссылке передать, просто я с двумерным массивом не умею работать.

Добавлено через 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
 Аватар для Fedor N
0 / 0 / 0
Регистрация: 22.06.2012
Сообщений: 17
01.07.2012, 19:54  [ТС]     Матрицы - перегрузка оператора * #13
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\mingw32\4.4.1\..\..\..\..\include\Matrixclass.h|180|error: 'class Matrix' has no member named 'tmp'|
silent_1991
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
01.07.2012, 20:04     Матрицы - перегрузка оператора * #14
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
 Аватар для Fedor N
0 / 0 / 0
Регистрация: 22.06.2012
Сообщений: 17
01.07.2012, 20:24  [ТС]     Матрицы - перегрузка оператора * #15
Цитата Сообщение от 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
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
01.07.2012, 20:36     Матрицы - перегрузка оператора * #16
Цитата Сообщение от Fedor N Посмотреть сообщение
Вот эти варианты работают
Аа, ну да, не заметил, что Array - одномерный. Тогда можно вместо первого варианта:
C++
1
tmp(i,j) +=Array[i * xsize + k]*M(k,j);
Fedor N
 Аватар для Fedor N
0 / 0 / 0
Регистрация: 22.06.2012
Сообщений: 17
01.07.2012, 20:44  [ТС]     Матрицы - перегрузка оператора * #17
Цитата Сообщение от 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
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
01.07.2012, 20:52     Матрицы - перегрузка оператора * #18
Ошибок в коде много, и очень серьёзных, но я не заметил причины неверной работы. Покажите весь код (в том числе клиентскую часть).
Fedor N
 Аватар для Fedor N
0 / 0 / 0
Регистрация: 22.06.2012
Сообщений: 17
01.07.2012, 21:02  [ТС]     Матрицы - перегрузка оператора * #19
Библиотека 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();
}
Да и в принципе вот и все с чем я сейчас работаю=). Потом другие функции появятся (они вроде норм работают)
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
01.07.2012, 21:37     Матрицы - перегрузка оператора *
Еще ссылки по теме:

C++ Перегрузка оператора =
C++ Перегрузка оператора ++
Перегрузка оператора () для разреженной матрицы в CSR-формате C++

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

Или воспользуйтесь поиском по форуму:
silent_1991
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
01.07.2012, 21:37     Матрицы - перегрузка оператора * #20
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;
}
Yandex
Объявления
01.07.2012, 21:37     Матрицы - перегрузка оператора *
Ответ Создать тему
Опции темы

Текущее время: 00:07. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru