Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.73/11: Рейтинг темы: голосов - 11, средняя оценка - 4.73
0 / 0 / 0
Регистрация: 18.05.2019
Сообщений: 14

Оператор скалярного произведения

04.12.2019, 18:05. Показов 2177. Ответов 12

Студворк — интернет-сервис помощи студентам
Имеется такой оператор скалярного произведения комплексных векторов:
C++
1
2
3
4
5
6
7
CComplex operator* (const CComplexVector& cv)
{
CComplex res = CComplex(0, 0);
for (size_t i = 0; i < MIN(this->length, cv.length); i++)
res += this->v[i] * CComplex(cv.v[i].re, -cv.v[i].im);
return res;
}
Компилятор выдает такую ошибку:

Ошибка: implicitly-declared "constexpr CComplex::CComplex(const CComplex&)" is deprecated
return res; (res подчеркнуто)

Что тут сделано не так?
Код полностью:

Добавлено через 2 минуты
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
#include <iostream>
 using namespace std;
class CComplex
{
public:
    double re, im;
    CComplex(double r = 0.0, double i = 0.0) { re = r; im = i; }
 
    CComplex& operator=(const CComplex& cc)
    {
        if (this == &cc)
            return *this;
        this->re = cc.re;
        this->im = cc.im;
        return *this;
    }
 
    CComplex& operator+=(const CComplex& cc)
    {
        return *this = *this + cc;
    }
 
    CComplex& operator-=(const CComplex& cc)
    {
        return *this = *this - cc;
    }
 
    friend const CComplex& operator+(const CComplex& left, const CComplex& right);
    friend const CComplex& operator-(const CComplex& left, const CComplex& right);
    friend const CComplex& operator*(const CComplex& left, const CComplex& right);
    friend ostream& operator<< (ostream& out, const CComplex& c);
};
 
const CComplex& operator+(const CComplex& left, const CComplex& right)
{
    return *new CComplex(left.re + right.re, left.im + right.im);
}
 
const CComplex& operator-(const CComplex& left, const CComplex& right)
{
    return *new CComplex(left.re - right.re, left.im - right.im);
}
 
const CComplex& operator*(const CComplex& left, const CComplex& right)
{
    return *new CComplex(
        left.re * right.re - left.im * right.im,
        left.im * right.re + left.re * right.im);
}
 
ostream& operator<< (ostream& out, const CComplex& c)
{
    return out << '(' << c.re << ',' << c.im << ')';
}
 
//Первый файл .h
//#pragma once
//#include "CComplex10.h"
//#include <iostream>
#define MAX(a, b) ((a)>(b)?(a):(b))
#define MIN(a, b) ((a)<(b)?(a):(b))
//using namespace std;
 
class CComplexVector
{
    CComplex* v = nullptr;
    size_t length = 0;
public:
 
    CComplexVector(const CComplexVector& cv)
    {
        if (length == cv.length)
        {
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cv.v[i];
        }
    }
 
    CComplexVector(size_t size = 0, CComplex cd = CComplex(0, 0))
    {
        if (size)
        {
            length = size;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cd;
        }
    }
 
    CComplexVector(size_t size, CComplex cd[])
    {
        if (size && cd)
        {
            length = size;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cd[i];
        }
    }
 
 
    ~CComplexVector() { Clear(); }
   size_t Length() { return length; }
void Clear()
    {
        if (v != nullptr)
            delete[]v;
        v = nullptr;
        length = 0;
    }
 
    CComplexVector& operator= (const CComplex& cd)
    {
        for (size_t i = 0; i < length; i++)
            this->v[i] = cd;
        return *this;
    }
 
    CComplexVector& operator= (const CComplexVector& cv)
    {
        if (this == &cv)
            return *this;
        Clear();
        if (cv.length)
        {   
            length = cv.length;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cv.v[i];
        }
        return *this;
    }
 
    CComplexVector& operator+= (const CComplexVector& cv)
    {
        return *this = *this + cv;
    }
 
    CComplexVector& operator-= (const CComplexVector& cv)
    {
        return *this = *this - cv;
    }
 
    CComplex& operator[] (size_t i)
    {
        if (i < this->length)
            return this->v[i];
        cout << "Array index out of beyond" << endl;
        exit(-1);
    }
 
    CComplexVector& operator+ (const CComplexVector& right)
    {
        CComplexVector* res = new CComplexVector(MAX(this->length, right.length));
        if (this->length > right.length)
            *res = *this;
        else
            *res = right;
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res->v[i] = this->v[i] + right.v[i];
        return *res;
    }
 
    CComplexVector& operator- (const CComplexVector& right)
    {
        CComplexVector* res = new CComplexVector(MAX(this->length, right.length));
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res->v[i] = this->v[i] - right.v[i];
        for (size_t i = MIN(this->length, right.length); i < MAX(this->length, right.length); i++)
            if (this->length > right.length)
                res->v[i] = this->v[i];
            else
                res->v[i] -= right.v[i];
        return *res;
    }
 
    CComplex operator* (const CComplexVector& cv)
    {
        CComplex res = CComplex(0, 0);
        for (size_t i = 0; i < MIN(this->length, cv.length); i++)
            res += this->v[i] * CComplex(cv.v[i].re, -cv.v[i].im);
        return res;
    }
    friend ostream& operator<< (ostream& out, const CComplexVector& cv);
};
ostream& operator<< (ostream& out, const CComplexVector& cv)
{
    if (cv.v)
    {
        size_t last = cv.length - 1;
        for (size_t i = 0; i < last; i++)
            out << cv.v[i] << ',';
        out << cv.v[last];
    }
    return out;
}
 
void print_vectors(CComplexVector& a, CComplexVector& b, CComplexVector& c)
{
    cout << "A size (" << a.Length() << "):\n\t" << a << endl;
    cout << "B size (" << b.Length() << "):\n\t" << b << endl;
    cout << "C size (" << c.Length() << "):\n\t" << c << endl;
}
 
int main()
{
    cout << "Create vectors" << endl;
    CComplexVector a = CComplexVector(3,CComplex(1,1));
    CComplexVector b = CComplexVector(3,CComplex(2,2));
    CComplex arr[3] = { (1,2), (1,3), (2,4) };
    CComplexVector c = CComplexVector(3,arr);
    print_vectors(a, b, c);
 
    b = CComplexVector(a);
 
    cout << "Create vectors" << endl;
    a = CComplexVector(3, CComplex(1,1));
    b = CComplexVector(3, CComplex(2, 2));
    c = CComplexVector(3, CComplex(3, 3));
    print_vectors(a, b, c);
 
    cout << "Demo '+' or '-'" << endl;
    a = a + a;
    b = b + a;
    c = c - a;
    print_vectors(a, b, c);
 
    cout << "Demo '+=' or '-='" << endl;
    a += c;
    b -= a;
    c += a;
    print_vectors(a, b, c);
 
    cout << "Print 'scalar product' of A and B using '*'" << endl;
    cout << '\t' << a * b << endl;
 
    return 0;
}
0
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
04.12.2019, 18:05
Ответы с готовыми решениями:

Обобщенное вычисление скалярного произведения
Здравствуйте! С помощью шаблонов (их рекурсивных вызовов) должен был посчитать произведение на этапе компиляции. Только ругается на...

Нахождение скалярного произведения векторов
Известно два вектора Х(х1,х2...х10) и Y(y1,y2...y10). Нужно ввести свои значения и найти скалярное произведение векторов X*Y

Вычисления скалярного произведения векторов
Напишите подпрограмму вычисления скалярного произведения векторов произвольной размерности. Первый раз с такой темой сталкиваюсь. Я...

12
Мозгоправ
 Аватар для L0M
1745 / 1039 / 468
Регистрация: 01.10.2018
Сообщений: 2,138
Записей в блоге: 2
04.12.2019, 19:04
Ouou12, вы наступаете на те же грабли с утечкой памяти, что и в теме Почему нельзя вызывать конструктор без аргументов?. Я обратил на это ваше внимание, но, видимо, это прошло мимо вас.

Цитата Сообщение от Ouou12 Посмотреть сообщение
Компилятор выдает такую ошибку:
Ошибка: implicitly-declared "constexpr CComplex::CComplex(const CComplex&)" is deprecated
Определите явно конструктор копии для CComplex.

Вообще-то есть стандартный заголовк <complex> для комплексных чисел.

Для MIN(a, b)/MAX(a, b) есть шаблонные функции.
0
0 / 0 / 0
Регистрация: 18.05.2019
Сообщений: 14
04.12.2019, 19:25  [ТС]
Исправил. Утечки памяти больше нет? Конструктор копирования вроде правилно написан.


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
#include <iostream>
 
using namespace std;
class CComplex
{
public:
    double re, im;
    CComplex(double r = 0.0, double i = 0.0) { re = r; im = i; }
    CComplex(const CComplex& cc)
    { 
        this->re = cc.re;
        this->im = cc.im; }
    CComplex& operator=(const CComplex& cc)
    {
        if (this == &cc)
            return *this;
        this->re = cc.re;
        this->im = cc.im;
        return *this;
    }
 
    CComplex& operator+=(const CComplex& cc)
    {
        return *this = *this + cc;
    }
 
    CComplex& operator-=(const CComplex& cc)
    {
        return *this = *this - cc;
    }
 
    friend  CComplex& operator+(const CComplex& left, const CComplex& right);
    friend  CComplex& operator-(const CComplex& left, const CComplex& right);
    friend  CComplex& operator*(const CComplex& left, const CComplex& right);
    friend ostream& operator<< (ostream& out, const CComplex& c);
};
 
CComplex& operator+(const CComplex& left, const CComplex& right)
{
    return *new CComplex(left.re + right.re, left.im + right.im);
}
 
CComplex& operator-(const CComplex& left, const CComplex& right)
{
    return *new CComplex(left.re - right.re, left.im - right.im);
}
 
 CComplex& operator*(const CComplex& left, const CComplex& right)
{
    return *new CComplex(
        left.re * right.re - left.im * right.im,
        left.im * right.re + left.re * right.im);
}
 
ostream& operator<< (ostream& out, const CComplex& c)
{
    return out << '(' << c.re << ',' << c.im << ')';
}
 
 
#define MAX(a, b) ((a)>(b)?(a):(b))
#define MIN(a, b) ((a)<(b)?(a):(b))
 
 
class CComplexVector
{
    CComplex* v = nullptr;
    size_t length = 0;
public:
 
    CComplexVector(const CComplexVector& cv)
    {
        if (length == cv.length)
        {
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cv.v[i];
        }
    }
 
    CComplexVector(size_t size = 0, CComplex cd = CComplex(0, 0))
    {
        if (size)
        {
            length = size;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cd;
        }
    }
 
 
 
    CComplexVector(size_t size, CComplex cd[])
    {
        if (size && cd)
        {
            length = size;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cd[i];
        }
    }
 
 
    ~CComplexVector() { Clear(); }
   size_t Length() { return length; }
void Clear()
    {
        if (v != nullptr)
            delete[]v;
        v = nullptr;
        length = 0;
    }
 
    CComplexVector& operator= (const CComplex& cd)
    {
        for (size_t i = 0; i < length; i++)
            this->v[i] = cd;
        return *this;
    }
 
    CComplexVector& operator= (const CComplexVector& cv)
    {
        if (this == &cv)
            return *this;
        Clear();
        if (cv.length)
        {length = cv.length
           v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cv.v[i];
        }
        return *this;
    }
 
    CComplexVector& operator+= (const CComplexVector& cv)
    {
        return *this = *this + cv;
    }
 
    CComplexVector& operator-= (const CComplexVector& cv)
    {
        return *this = *this - cv;
    }
 
    CComplex& operator[] (size_t i)
    {
        if (i < this->length)
            return this->v[i];
        cout << "Array index out of beyond" << endl;
        exit(-1);
    }
 
    CComplexVector& operator+ (const CComplexVector& right)
    {
        CComplexVector* res = new CComplexVector(MAX(this->length, right.length));
        if (this->length > right.length)
            *res = *this;
        else
            *res = right;
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res->v[i] = this->v[i] + right.v[i];
        return *res;
    }
 
    CComplexVector& operator- (const CComplexVector& right)
    {
        CComplexVector* res = new CComplexVector(MAX(this->length, right.length));
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res->v[i] = this->v[i] - right.v[i];
        for (size_t i = MIN(this->length, right.length); i < MAX(this->length, right.length); i++)
            if (this->length > right.length)
                res->v[i] = this->v[i];
            else
               res->v[i] -= right.v[i];
        return *res;
    }
 
    CComplex operator* (const CComplexVector& cv)
    {
        CComplex res = CComplex(0, 0);
        for (size_t i = 0; i < MIN(this->length, cv.length); i++)
            res += this->v[i] * CComplex(cv.v[i].re, -cv.v[i].im);
        return res;
    }
    friend ostream& operator<< (ostream& out, const CComplexVector& cv);
};
ostream& operator<< (ostream& out, const CComplexVector& cv)
{
    if (cv.v)
    {
        size_t last = cv.length - 1;
        for (size_t i = 0; i < last; i++)
            out << cv.v[i] << ',';
        out << cv.v[last];
    }
    return out;
}
 
void print_vectors(CComplexVector& a, CComplexVector& b, CComplexVector& c)
{
    cout << "A size (" << a.Length() << "):\n\t" << a << endl;
    cout << "B size (" << b.Length() << "):\n\t" << b << endl;
    cout << "C size (" << c.Length() << "):\n\t" << c << endl;
}
 
int main()
{  
    cout << "Create vectors" << endl;
    CComplexVector a = CComplexVector(3,CComplex(1,1));
    CComplexVector b = CComplexVector(3,CComplex(2,2));
    CComplex arr[3] = { CComplex(1,2), CComplex(1,3), CComplex(2,4) };
    CComplexVector c = CComplexVector(3,arr);
    print_vectors(a, b, c);
 
    b = CComplexVector(a);
 
    cout << "Create vectors" << endl;
    a = CComplexVector(3);
    b = CComplexVector(3, CComplex(2, 2));
    c = CComplexVector(3, CComplex(3, 3));
    print_vectors(a, b, c);
 
    cout << "Demo '+' or '-'" << endl;
    a = a + a;
    b = b + a;
    c = c - a;
    print_vectors(a, b, c);
 
    cout << "Demo '+=' or '-='" << endl;
    a += c;
    b -= a;
    c += a;
    print_vectors(a, b, c);
 
    cout << "Print 'scalar product' of A and B using '*'" << endl;
    cout << '\t' << a * b << endl;
 
    return 0;
}
0
Мозгоправ
 Аватар для L0M
1745 / 1039 / 468
Регистрация: 01.10.2018
Сообщений: 2,138
Записей в блоге: 2
04.12.2019, 21:05
Цитата Сообщение от Ouou12 Посмотреть сообщение
Утечки памяти больше нет?
А что вы у меня спрашиваете? Проверяйте. Есть она в вашем коде или нет.

Спойлер: есть.

И синтаксическая ошибка в вашем коде есть. Вы перед тем как показать код его компилируете? запускаете? тестируете?

И ещё раз:
C++
1
2
3
    friend  CComplex operator+(const CComplex& left, const CComplex& right);
    friend  CComplex operator-(const CComplex& left, const CComplex& right);
    friend  CComplex operator*(const CComplex& left, const CComplex& right);
0
0 / 0 / 0
Регистрация: 18.05.2019
Сообщений: 14
06.12.2019, 12:34  [ТС]
Ошибка исчезла, спасибо! Но память все равно течет(
Миниатюры
Оператор скалярного произведения  
0
Мозгоправ
 Аватар для L0M
1745 / 1039 / 468
Регистрация: 01.10.2018
Сообщений: 2,138
Записей в блоге: 2
06.12.2019, 13:38
Ouou12, в таких операторах не надо выделять память:
C++
1
2
3
4
CComplex operator+(const CComplex& left, const CComplex& right)
{
    return CComplex(left.re + right.re, left.im + right.im);
}
Проверьте и исправьте все подобные операторы в CComplex и в CComplexVector.
0
0 / 0 / 0
Регистрация: 18.05.2019
Сообщений: 14
08.12.2019, 17:24  [ТС]
Все равно есть небольшая утечка
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
#include <iostream>
 
using namespace std;
class CComplex
{
public:
    double re, im;
    CComplex(double r = 0.0, double i = 0.0) { re = r; im = i; }
    CComplex(const CComplex& cc)
    { 
        this->re = cc.re;
        this->im = cc.im; }
    CComplex& operator=(const CComplex& cc)
    {
        if (this == &cc)
            return *this;
        this->re = cc.re;
        this->im = cc.im;
        return *this;
    }
 
    CComplex& operator+=(const CComplex& cc)
    {
        return *this = *this + cc;
    }
 
    CComplex& operator-=(const CComplex& cc)
    {
        return *this = *this - cc;
    }
 
    friend  CComplex operator+(const CComplex& left, const CComplex& right);
    friend  CComplex operator-(const CComplex& left, const CComplex& right);
    friend  CComplex operator*(const CComplex& left, const CComplex& right);
    friend ostream& operator<< (ostream& out, const CComplex& c);
};
 
CComplex operator+(const CComplex& left, const CComplex& right)
{
    return CComplex(left.re + right.re, left.im + right.im);
}
 
CComplex operator-(const CComplex& left, const CComplex& right)
{
    return CComplex(left.re - right.re, left.im - right.im);
}
 
 CComplex operator*(const CComplex& left, const CComplex& right)
{
    return CComplex(
        left.re * right.re - left.im * right.im,
        left.im * right.re + left.re * right.im);
}
 
ostream& operator<< (ostream& out, const CComplex& c)
{
    return out << '(' << c.re << ',' << c.im << ')';
}
 
 
#define MAX(a, b) ((a)>(b)?(a):(b))
#define MIN(a, b) ((a)<(b)?(a):(b))
 
 
class CComplexVector
{
    CComplex* v = nullptr;
    size_t length = 0;
public:
 
    CComplexVector(const CComplexVector& cv)
    {
        if (length == cv.length)
        {
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cv.v[i];
        }
    }
 
    CComplexVector(size_t size = 0, CComplex cd = CComplex(0, 0))
    {
        if (size)
        {
            length = size;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cd;
        }
    }
 
 
 
    CComplexVector(size_t size, CComplex cd[])
    {
        if (size && cd)
        {
            length = size;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cd[i];
        }
    }
 
 
    ~CComplexVector() { Clear(); }
   size_t Length() { return length; }
void Clear()
    {
        if (v != nullptr)
            delete[]v;
        v = nullptr;
        length = 0;
    }
 
    CComplexVector& operator= (const CComplex& cd)
    {
        for (size_t i = 0; i < length; i++)
            this->v[i] = cd;
        return *this;
    }
 
    CComplexVector& operator= (const CComplexVector& cv)
    {
        if (this == &cv)
            return *this;
        Clear();
        if (cv.length)
        {length = cv.length;
           v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cv.v[i];
        }
        return *this;
    }
 
    CComplexVector& operator+= (const CComplexVector& cv)
    {
        return *this = *this + cv;
    }
 
    CComplexVector& operator-= (const CComplexVector& cv)
    {
        return *this = *this - cv;
    }
 
    CComplex& operator[] (size_t i)
    {
        if (i < this->length)
            return this->v[i];
        cout << "Array index out of beyond" << endl;
        exit(-1);
    }
 
    CComplexVector& operator+ (const CComplexVector& right)
    {
        CComplexVector* res = new CComplexVector(MAX(this->length, right.length));
        if (this->length > right.length)
            *res = *this;
        else
            *res = right;
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res->v[i] = this->v[i] + right.v[i];
        return *res;
    }
 
    CComplexVector& operator- (const CComplexVector& right)
    {
        CComplexVector* res = new CComplexVector(MAX(this->length, right.length));
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res->v[i] = this->v[i] - right.v[i];
        for (size_t i = MIN(this->length, right.length); i < MAX(this->length, right.length); i++)
            if (this->length > right.length)
                res->v[i] = this->v[i];
            else
               res->v[i] -= right.v[i];
        return *res;
    }
 
    CComplex operator* (const CComplexVector& cv)
    {
        CComplex res = CComplex(0, 0);
        for (size_t i = 0; i < MIN(this->length, cv.length); i++)
            res += this->v[i] * CComplex(cv.v[i].re, -cv.v[i].im);
        return res;
    }
    friend ostream& operator<< (ostream& out, const CComplexVector& cv);
};
ostream& operator<< (ostream& out, const CComplexVector& cv)
{
    if (cv.v)
    {
        size_t last = cv.length - 1;
        for (size_t i = 0; i < last; i++)
            out << cv.v[i] << ',';
        out << cv.v[last];
    }
    return out;
}
 
void print_vectors(CComplexVector& a, CComplexVector& b, CComplexVector& c)
{
    cout << "A size (" << a.Length() << "):\n\t" << a << endl;
    cout << "B size (" << b.Length() << "):\n\t" << b << endl;
    cout << "C size (" << c.Length() << "):\n\t" << c << endl;
}
 
int main()
{  
    cout << "Create vectors" << endl;
    CComplexVector a = CComplexVector(3,CComplex(1,1));
    CComplexVector b = CComplexVector(3,CComplex(2,2));
    CComplex arr[3] = { CComplex(1,2), CComplex(1,3), CComplex(2,4) };
    CComplexVector c = CComplexVector(3,arr);
    print_vectors(a, b, c);
 
    b = CComplexVector(a);
 
    cout << "Create vectors" << endl;
    a = CComplexVector(3);
    b = CComplexVector(3, CComplex(2, 2));
    c = CComplexVector(3, CComplex(3, 3));
    print_vectors(a, b, c);
 
    cout << "Demo '+' or '-'" << endl;
    a = a + a;
    b = b + a;
    c = c - a;
    print_vectors(a, b, c);
 
    cout << "Demo '+=' or '-='" << endl;
    a += c;
    b -= a;
    c += a;
    print_vectors(a, b, c);
 
    cout << "Print 'scalar product' of A and B using '*'" << endl;
    cout << '\t' << a * b << endl;
 
    return 0;
}
0
Мозгоправ
 Аватар для L0M
1745 / 1039 / 468
Регистрация: 01.10.2018
Сообщений: 2,138
Записей в блоге: 2
08.12.2019, 19:42
Ouou12, неправильно написан конструктор копии CComplexVector.
Память течёт из-за
C++
1
2
CComplexVector& operator+ (const CComplexVector& right);
CComplexVector& operator- (const CComplexVector& right);
которые должны возвращать не ссылку, а значение.
Я вам об этом говорил:
Цитата Сообщение от L0M Посмотреть сообщение
Проверьте и исправьте все подобные операторы в CComplex и в CComplexVector.
0
0 / 0 / 0
Регистрация: 18.05.2019
Сообщений: 14
09.12.2019, 16:19  [ТС]
Почему неправильно? Эта строчка лишняя:
C++
1
if (length == cv.length)
?
0
Мозгоправ
 Аватар для L0M
1745 / 1039 / 468
Регистрация: 01.10.2018
Сообщений: 2,138
Записей в блоге: 2
09.12.2019, 18:12
Цитата Сообщение от Ouou12 Посмотреть сообщение
Эта строчка лишняя:
Да. Если у вас длины не равны, а они не равны, то конструктор копии ничего не делает.
0
0 / 0 / 0
Регистрация: 18.05.2019
Сообщений: 14
09.12.2019, 18:48  [ТС]
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
CComplexVector operator+ (const CComplexVector& right)
    {
        CComplexVector res = CComplexVector(MAX(this->length, right.length));
        if (this->length > right.length)
            res = *this;
        else
            res = right;
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res.v[i] = this->v[i] + right.v[i];
        return res;
    }
 
    CComplexVector operator- (const CComplexVector& right)
    {
        CComplexVector res = CComplexVector(MAX(this->length, right.length));
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res.v[i] = this->v[i] - right.v[i];
        for (size_t i = MIN(this->length, right.length); i < MAX(this->length, right.length); i++)
            if (this->length > right.length)
                res.v[i] = this->v[i];
            else
                res.v[i] -= right.v[i];
        return res;
    }
Что-то не то. После такого исправления, векторы почему-то очищаются.
0
 Аватар для zayats80888
6352 / 3523 / 1428
Регистрация: 07.02.2019
Сообщений: 8,995
09.12.2019, 19:17
Ouou12, исправлял не очень внимательно и не всё, мог накосячить
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
class CComplexVector
{
    CComplex* v = nullptr;
    size_t length = 0;
public:
 
    CComplexVector(const CComplexVector& cv):v(new CComplex[cv.length]), length(cv.length)
    {
        for (size_t i = 0; i < length; i++)
        {
            this->v[i] = cv.v[i];
        }
    }
 
    CComplexVector(size_t size = 0, CComplex cd = CComplex(0, 0)) :v(new CComplex[size]), length(size)
    {
        for (size_t i = 0; i < length; i++)
        {
            this->v[i] = cd;
        }
    }
    
    template<size_t size> CComplexVector(CComplex (&cd)[size]):CComplexVector(size)
    {
        for (size_t i = 0; i < length; i++)
            this->v[i] = cd[i];
    }
    
 
    ~CComplexVector() { Clear(); }
    size_t Length() const { return length; }
    void Clear()
    {
        delete[] v;
        v = nullptr;
        length = 0;
    }
 
    CComplexVector& operator= (const CComplex& cd)
    {
        for (size_t i = 0; i < length; i++)
            this->v[i] = cd;
        return *this;
    }
 
    CComplexVector& operator=(const CComplexVector& cv)
    {
        if (this == &cv)
            return *this;
        Clear();
        if (cv.length)
        {
            length = cv.length;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cv.v[i];
        }
        return *this;
    }
 
    CComplexVector& operator+= (const CComplexVector& cv)
    {
        return *this = *this + cv;
    }
 
    CComplexVector& operator-= (const CComplexVector& cv)
    {
        return *this = *this - cv;
    }
 
    CComplex& operator[] (size_t i)
    {
        if (i < this->length)
            return this->v[i];
        cout << "Array index out of beyond" << endl;
        exit(-1);
    }
 
    CComplexVector operator+ (const CComplexVector& right) const
    {
        CComplexVector res;
        const CComplexVector* add = &right;
        if (this->length >= right.length)
        {
            res = *this;
        }
        else
        {
            res = right;
            add = this;
        }
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res.v[i] += add->v[i];
        return res;
    }
 
    CComplexVector operator- (const CComplexVector& right) const
    {
        CComplexVector res;
        const CComplexVector* add = &right;
        if (this->length >= right.length)
        {
            res = *this;
        }
        else
        {
            res = right;
            add = this;
        }
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res.v[i] -= add->v[i];
        return res;
    }
 
    CComplex operator* (const CComplexVector& cv) const
    {
        CComplex res = CComplex(0, 0);
        for (size_t i = 0; i < MIN(this->length, cv.length); i++)
            res += this->v[i] * CComplex(cv.v[i].re, -cv.v[i].im);
        return res;
    }
    friend ostream& operator<< (ostream& out, const CComplexVector& cv);
};
0
Мозгоправ
 Аватар для L0M
1745 / 1039 / 468
Регистрация: 01.10.2018
Сообщений: 2,138
Записей в блоге: 2
09.12.2019, 20:18
Цитата Сообщение от Ouou12 Посмотреть сообщение
Что-то не то. После такого исправления, векторы почему-то очищаются.
Видимо потому, что не исправили конструктор копии.

Ваш код без утечек памяти, который вроде работает:
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
243
244
245
246
#include <iostream>
 
using namespace std;
class CComplex
{
public:
    double re, im;
    CComplex(double r = 0.0, double i = 0.0) { re = r; im = i; }
    CComplex(const CComplex& cc)
    {
        this->re = cc.re;
        this->im = cc.im;
    }
    CComplex& operator=(const CComplex& cc)
    {
        if (this == &cc)
            return *this;
        this->re = cc.re;
        this->im = cc.im;
        return *this;
    }
 
    CComplex& operator+=(const CComplex& cc)
    {
        return *this = *this + cc;
    }
 
    CComplex& operator-=(const CComplex& cc)
    {
        return *this = *this - cc;
    }
 
    friend  CComplex operator+(const CComplex& left, const CComplex& right);
    friend  CComplex operator-(const CComplex& left, const CComplex& right);
    friend  CComplex operator*(const CComplex& left, const CComplex& right);
    friend ostream& operator<< (ostream& out, const CComplex& c);
};
 
CComplex operator+(const CComplex& left, const CComplex& right)
{
    return CComplex(left.re + right.re, left.im + right.im);
}
 
CComplex operator-(const CComplex& left, const CComplex& right)
{
    return CComplex(left.re - right.re, left.im - right.im);
}
 
CComplex operator*(const CComplex& left, const CComplex& right)
{
    return CComplex(
        left.re * right.re - left.im * right.im,
        left.im * right.re + left.re * right.im);
}
 
ostream& operator<< (ostream& out, const CComplex& c)
{
    return out << '(' << c.re << ',' << c.im << ')';
}
 
 
#define MAX(a, b) ((a)>(b)?(a):(b))
#define MIN(a, b) ((a)<(b)?(a):(b))
 
 
class CComplexVector
{
    CComplex* v = nullptr;
    size_t length = 0;
public:
 
    CComplexVector(const CComplexVector& cv)
    {
        length = cv.length;
        v = new CComplex[length];
        for (size_t i = 0; i < length; i++)
            this->v[i] = cv.v[i];
    }
 
    CComplexVector(size_t size = 0, CComplex cd = CComplex(0, 0))
    {
        if (size)
        {
            length = size;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cd;
        }
    }
 
 
 
    CComplexVector(size_t size, CComplex cd[])
    {
        if (size && cd)
        {
            length = size;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cd[i];
        }
    }
 
 
    ~CComplexVector() { Clear(); }
    size_t Length() { return length; }
    void Clear()
    {
        if (v != nullptr)
            delete[]v;
        v = nullptr;
        length = 0;
    }
 
    CComplexVector& operator= (const CComplex& cd)
    {
        for (size_t i = 0; i < length; i++)
            this->v[i] = cd;
        return *this;
    }
 
    CComplexVector& operator= (const CComplexVector& cv)
    {
        if (this == &cv)
            return *this;
        Clear();
        if (cv.length)
        {
            length = cv.length;
            v = new CComplex[length];
            for (size_t i = 0; i < length; i++)
                this->v[i] = cv.v[i];
        }
        return *this;
    }
 
    CComplexVector& operator+= (const CComplexVector& cv)
    {
        return *this = *this + cv;
    }
 
    CComplexVector& operator-= (const CComplexVector& cv)
    {
        return *this = *this - cv;
    }
 
    CComplex& operator[] (size_t i)
    {
        if (i < this->length)
            return this->v[i];
        cout << "Array index out of beyond" << endl;
        exit(-1);
    }
 
    CComplexVector operator+ (const CComplexVector& right)
    {
        CComplexVector res = CComplexVector(MAX(this->length, right.length));
        if (this->length > right.length)
            res = *this;
        else
            res = right;
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res.v[i] = this->v[i] + right.v[i];
        return res;
    }
 
    CComplexVector operator- (const CComplexVector& right)
    {
        CComplexVector res = CComplexVector(MAX(this->length, right.length));
        for (size_t i = 0; i < MIN(this->length, right.length); i++)
            res.v[i] = this->v[i] - right.v[i];
        for (size_t i = MIN(this->length, right.length); i < MAX(this->length, right.length); i++)
            if (this->length > right.length)
                res.v[i] = this->v[i];
            else
                res.v[i] -= right.v[i];
        return res;
    }
 
    CComplex operator* (const CComplexVector& cv)
    {
        CComplex res = CComplex(0, 0);
        for (size_t i = 0; i < MIN(this->length, cv.length); i++)
            res += this->v[i] * CComplex(cv.v[i].re, -cv.v[i].im);
        return res;
    }
    friend ostream& operator<< (ostream& out, const CComplexVector& cv);
};
ostream& operator<< (ostream& out, const CComplexVector& cv)
{
    if (cv.v)
    {
        size_t last = cv.length - 1;
        for (size_t i = 0; i < last; i++)
            out << cv.v[i] << ',';
        out << cv.v[last];
    }
    return out;
}
 
void print_vectors(CComplexVector& a, CComplexVector& b, CComplexVector& c)
{
    cout << "A size (" << a.Length() << "):\n\t" << a << endl;
    cout << "B size (" << b.Length() << "):\n\t" << b << endl;
    cout << "C size (" << c.Length() << "):\n\t" << c << endl;
}
 
int main()
{
    cout << "Create vectors" << endl;
    CComplexVector a = CComplexVector(3, CComplex(1, 1));
    CComplexVector b = CComplexVector(3, CComplex(2, 2));
    CComplex arr[3] = { CComplex(1,2), CComplex(1,3), CComplex(2,4) };
    CComplexVector c = CComplexVector(3, arr);
    print_vectors(a, b, c);
 
    b = CComplexVector(a);
 
    cout << "Create vectors" << endl;
    a = CComplexVector(3);
    b = CComplexVector(3, CComplex(2, 2));
    c = CComplexVector(3, CComplex(3, 3));
    print_vectors(a, b, c);
 
    cout << "Demo '+' or '-'" << endl;
    a = a + a;
    b = b + a;
    c = c - a;
    print_vectors(a, b, c);
 
    a = a + b;
    b = b + c;
    c = c - b;
    print_vectors(a, b, c);
 
    cout << "Demo '+=' or '-='" << endl;
    a += c;
    b -= a;
    c += a;
    print_vectors(a, b, c);
 
    cout << "Print 'scalar product' of A and B using '*'" << endl;
    cout << '\t' << a * b << endl;
 
    return 0;
}
Надо тестировать. Особенно для векторов разной длины.
1
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
09.12.2019, 20:18
Помогаю со студенческими работами здесь

Возникли трудности с поиском скалярного произведения
По условию дана целочисленная матрица X размера mxn. скалярное произведение двух векторов B и C = b1c1+b2c2+...bcn. Вывести на экран...

Функции скалярного произведения и длины вектора
Здравствуйте, помогите с лабой, у самого не выходит. С++ В программе создать функцию, грамотно определив ее параметры и тип,...

Перегрузка операции умножения, сложения и скалярного произведения
Помогите пожалуйста написать программу! Создать шаблон класса N-мерного вектора (с элементами неизвестного заранее типа) и перегрузить...

Вычисление суммы, разности и скалярного произведения векторов
Ввести два вектора. Вектора – динамические массивы. Вычислить сумму, разность и скалярное произведение. Ввод размерности 0 – выход из...

ScalarProduct для вычисления скалярного произведения двух векторов
Написать программу для проверки работы процедуры ScalarProduct для вычисления скалярного произведения двух векторов. Заранее спасибо...


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

Или воспользуйтесь поиском по форуму:
13
Ответ Создать тему
Новые блоги и статьи
SDL3 для Web (WebAssembly): сборка C/C++ проекта из консоли
8Observer8 30.01.2026
Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .
Установка Emscripten SDK (emsdk) и CMake на Windows для сборки C и C++ приложений в WebAssembly (Wasm)
8Observer8 30.01.2026
Чтобы скачать Emscripten SDK (emsdk) необходимо сначало скачать и уставить Git: Install for Windows. Следуйте стандартной процедуре установки Git через установщик. Система контроля версиями Git. . .
Подключение Box2D v3 к SDL3 для Android: физика и отрисовка коллайдеров
8Observer8 29.01.2026
Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами. Версия v3 была полностью переписана на Си, в. . .
Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT
bedvit 28.01.2026
Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .
Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL_LoadPNG (без SDL3_image)
8Observer8 28.01.2026
Содержание блога SDL3 имеет собственные средства для загрузки и отображения PNG-файлов с альфа-каналом и базовой работы с ними. В этой инструкции используется функция SDL_LoadPNG(), которая. . .
Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL3_image
8Observer8 27.01.2026
Содержание блога SDL3_image - это библиотека для загрузки и работы с изображениями. Эта пошаговая инструкция покажет, как загрузить и вывести на экран смартфона картинку с альфа-каналом, то есть с. . .
Влияние грибов на сукцессию
anaschu 26.01.2026
Бифуркационные изменения массы гриба происходят тогда, когда мы уменьшаем массу компоста в 10 раз, а скорость прироста биомассы уменьшаем в три раза. Скорость прироста биомассы может уменьшаться за. . .
Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android
8Observer8 26.01.2026
Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru