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

Динамическое выделение памяти - C++

Восстановить пароль Регистрация
 
Alberto_Timakov
1 / 1 / 0
Регистрация: 21.11.2011
Сообщений: 183
21.11.2011, 19:14     Динамическое выделение памяти #1
Почему распечатка матрицы Matr1 в f2 дает матрицу из ед., а в f1 матрицу из нулей? Ведь адрес первого эл. массива Matr2 сохранен в переменной Matr1, и динамическая память после завершения функции не очищается?

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
f2(Type** Matr1, int LineMatr, int ColumnMatr)
{
 
...
 
Type** Matr2 = new Type* [LineMatr];
for(int j = 0; j < ColumnMatr; j++)
     {
          Matr2[j] = new Type [ColumnMatr];
     }
 
for(int i = 0; i < LineMatr; i++)
    for(int j = 0; j < ColumnMatr; j++)
        Matr2[i][j] = 1;
 
...
 
Matr1 = Matr2;
 
printf(Matr1);
 
...
 
}
 
f1()
{
 
...
 
Type** Matr1 = new Type* [LineMatr];
for(int j = 0; j < ColumnMatr; j++)
     {
          Matr1[j] = new Type [ColumnMatr];
     }
 
for(int i = 0; i < LineMatr; i++)
    for(int j = 0; j < ColumnMatr; j++)
        Matr1[i][j] = 0;
 
...
 
f2(Matr1, LineMatr, ColumnMatr);
 
printf(Matr1);
 
...
 
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
21.11.2011, 19:14     Динамическое выделение памяти
Посмотрите здесь:

C++ динамическое выделение памяти
Динамическое выделение памяти C++
C++ динамическое выделение памяти
C++ Динамическое выделение памяти
Динамическое выделение памяти C++
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
mc.Duck
Заблокирован
21.11.2011, 19:18     Динамическое выделение памяти #2
Что за тип такой "Type" ? Вроде нигде не объявляли новый тип)
St1ck_man
9 / 9 / 1
Регистрация: 21.11.2011
Сообщений: 36
21.11.2011, 19:23     Динамическое выделение памяти #3
Внутри функций f1 и f2 идет присвоение всем элементам матрицы нулей или единиц
C
1
2
3
for(int i = 0; i < LineMatr; i++)
        for(int j = 0; j < ColumnMatr; j++)
                Matr2[i][j] = 1;
вот тут например
go
Эксперт C++
3582 / 1362 / 128
Регистрация: 16.04.2009
Сообщений: 4,528
21.11.2011, 19:24     Динамическое выделение памяти #4
Цитата Сообщение от Alberto_Timakov Посмотреть сообщение
printf(Matr1);
Что это?
St1ck_man
9 / 9 / 1
Регистрация: 21.11.2011
Сообщений: 36
21.11.2011, 19:27     Динамическое выделение памяти #5
если функция printf() написана вручную то дело наверняка в ней .
xAtom
 Аватар для xAtom
910 / 735 / 60
Регистрация: 09.12.2010
Сообщений: 1,346
Записей в блоге: 1
21.11.2011, 19:31     Динамическое выделение памяти #6
Alberto_Timakov, для полного объяснения нужен полный код.
go
Эксперт C++
3582 / 1362 / 128
Регистрация: 16.04.2009
Сообщений: 4,528
21.11.2011, 19:34     Динамическое выделение памяти #7
C
1
2
3
4
5
6
7
8
9
10
11
void f1 (type ***matr1,...)
{ 
...
 
 
*matr1=matr2;
...
 
 
}
f1 (&matr1,...);
mc.Duck
Заблокирован
21.11.2011, 19:49     Динамическое выделение памяти #8
Alberto_Timakov, какой-то дырявый у вас код, ничего в нем не понятно!
Alberto_Timakov
1 / 1 / 0
Регистрация: 21.11.2011
Сообщений: 183
21.11.2011, 20:33  [ТС]     Динамическое выделение памяти #9
это был псевдокод или что-то вроде этого)))
mc.Duck: Type - некоторый тип данных, к примеру, double или инт, не суть важно.
go : printf(Matr1); - так для краткости, я обозначил вывод матрицы:
C++
1
2
3
4
5
6
7
8
            for(int i = 0; i < LineMatr; i++)
            {
                for(int j = 0; j < ColumnMatr; j++)
                {
                    printf("%4.1lf",Matr1[i][j]);
                }
                printf("\n");
            }
St1ck_man : нет.

xAtom: закоментил все лишнее.

код:

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
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
файл LUP-decomposition.h:
 
 
template <class Type>
    Type** Decomposition
        (
            /*Type** Matr,*///Матрица
            const int LineMatr,//количество ее строк, 
            const int ColumnMatr,//количество ее столбцов.
        /*  Type** MatrUp,*///Верхнетреугольная матрица, получившиеся в ходе LUP-разложения.
            Type** MatrLow//Нижнетреугольная матрица, получившиеся в ходе LUP-разложения.
        )
 
        {
 
            //Объявление перестановочной матрицы.
 
            Type** MatrPermutation;
            Type** _MatrPermutation = new Type* [LineMatr];
            for(int j = 0; j < ColumnMatr; j++)
                {
                    _MatrPermutation[j] = new Type [ColumnMatr];
                }
 
            Type** MatrMk = new Type* [LineMatr];
            for(int j = 0; j < ColumnMatr; j++)
                {
                    MatrMk[j] = new Type [ColumnMatr];
                }
 
            for(int i = 0; i < LineMatr; i++)
                {
                for(int j = 0; j < ColumnMatr; j++)
                    {
                        if(i == j)
                        {
                            _MatrPermutation[i][j] = 1;
                            MatrMk[i][j] = 1;
                        }
                        else
                        {
                            MatrMk[i][j] = 0;
                            _MatrPermutation[i][j] = 0;
                        }
 
                    }
                }
 
            //for(int k = 0; k < ColumnMatr; k++)
            //{
 
                //Определение перестановочной матрицы и перестановка соотв. строк
                //исходной матрицы.
                //int MaxElColumn = SearchMaxElColumnStartWithDiag(Matr,LineMatr - 1,k);
                //MatrPermutation = PermutationRows
                //(
                //  Matr,
                //  ColumnMatr,
                //  k,
                //  MaxElColumn//SearchMaxElColumnStartWithDiag(Matr,LineMatr - 1,k)
                //);
 
                // _MatrPermutation = MatrMul
                //( 
                //  MatrPermutation,
                //  LineMatr,
                //  ColumnMatr,
                //  _MatrPermutation,
                //  LineMatr,
                //  ColumnMatr
                //);
                //На главной диагонали у нижнетреугольной матрицы, должны стоять ед.
 
            //  MatrLow[k][k] =  1;
 
                //Уничтожение поддиагональных элементов текущего столбца
 
                //for(int i = k + 1; i < LineMatr; i++)
                //{
                //  MatrMk[i][k] =  Matr[i][k]/Matr[k][k];
 
                //  for(int j = k; j < ColumnMatr; j++)
                //  {
                //      Matr[i][j] = Matr[i][j] - MatrMk[i][k] * Matr[k][j];
                //  }
                //}
 
                //if(k > 0)
                //{;
                    //MatrLow = MatrMul
                    //( 
                    //  MatrMul
                    //  ( 
                    //      MatrMul
                    //      ( 
                    //          MatrPermutation,
                    //          LineMatr,
                    //          ColumnMatr,
                    //          MatrLow,
                    //          LineMatr,
                    //          ColumnMatr
                    //      ),
                    //      LineMatr,
                    //      ColumnMatr,
                    //      MatrPermutation,
                    //      LineMatr,
                    //      ColumnMatr
                    //  ),
                    //  LineMatr,
                    //  ColumnMatr,
                    //  MatrMk,
                    //  LineMatr,
                    //  ColumnMatr
                    //);
                    //Type temp = MatrLow[k][k - 1];
                    //MatrLow[k][k - 1] = MatrLow[MaxElColumn][k - 1];
                    //MatrLow[MaxElColumn][k - 1] = temp;
                }
                //else
                    MatrLow = MatrMk;
            //}
 
 
            //Определение верхнетреугольной матрицы 
            //по получившийся в результате прямого хода Гауса.
 
            //for(int i = 0; i < LineMatr; i++)
            //  {
            //      for(int j = 0; j < ColumnMatr; j++)
            //      {
            //          MatrUp[i][j] = Matr[i][j];
            //      }
            //  }
 
            // Возврат перестановочной матрицы в вызывающую функцию.
            printf("\nf1\n");
            for(int i = 0; i < LineMatr; i++)
            {
                for(int j = 0; j < ColumnMatr; j++)
                {
                    printf("%4.1lf",MatrLow[i][j]);
                }
                printf("\n");
            }
            return _MatrPermutation;
        }
 
 
файл использующий предыдущий:
 
 
#include "LUP-decomposition.h"
 
template <class Type>
    void SolutionSystemLUP
        (
            Type** Matr,
            const int LineMatr,
            const int ColumnMatr,
            Type* VectorX,
            Type* VectorB,
            FILE* fout
        )
 
        {
 
            Type** MatrUp=new Type* [LineMatr];for(int i=0;i<LineMatr;i++){MatrUp[i]=new Type [ColumnMatr];}
            Type** MatrLow=new Type* [LineMatr];for(int i=0;i<LineMatr;i++){MatrLow[i]=new Type [ColumnMatr];}
            //Type** MatrPermutation;
            //Type* _VectorB=new Type [LineMatr];for(int i=0;i<LineMatr;i++){_VectorB[i] = 0;}
            for(int i=0;i<LineMatr;i++)
            {
                for(int j=0;j<ColumnMatr;j++)
                {
                    MatrUp[i][j] = 0;
                    MatrLow[i][j] = 0;
                }
            }
 
 
            /*MatrPermutation = */Decomposition(/*Matr,*/LineMatr,ColumnMatr,/*MatrUp,*/MatrLow);
 
            printf("\nf2\n");
            for(int i = 0; i < ColumnMatr; i++)
            {
                for(int j = 0; j < ColumnMatr; j++)
                {
                    printf("%4.1lf",MatrLow[i][j]);
                }
                printf("\n");
            }
 
            //Type** M = MatrMul
            //  ( 
            //      MatrLow,
            //      LineMatr,
            //      ColumnMatr,
            //      MatrUp,
            //      LineMatr,
            //      ColumnMatr
            //  );
 
            //fprintf(fout,"\n\n");
            //for(int i = 0; i < ColumnMatr; i++)
            //{
            //  for(int j = 0; j < ColumnMatr; j++)
            //  {
            //      fprintf(fout,"%4.1lf",M[i][j]);
            //  }
            //  fprintf(fout,"\n");
            //}
 
            /*printf("\nssss2\n");
            for(int i = 0; i < ColumnMatr; i++)
            {
                for(int j = 0; j < ColumnMatr; j++)
                {
                    printf("%4.1lf",MatrLow[i][j]);
                }
                printf("\n");
            }   */
            //printf("\n\n\n");
            //for(int i = 0; i < LineMatr; i++)
            //{
            //  for(int j = 0; j< ColumnMatr; j++)
            //  {
            //      printf("%4.5lf",MatrLow[i][j]);
            //  }
            //  printf("\n");
            //}
 
            //for(int i = 0; i < LineMatr; i++)
            //{
            //  for(int j = 0; j < ColumnMatr; j++)
            //  {
            //      _VectorB[i] += MatrPermutation[i][j]*VectorB[j];
            //  }
            //}
 
            //VectorB = _VectorB;
            // Опеределение вектора решения системы соотв. нижнетреугольной матр. 
            //с правой частью VectorB.
 
            //Type* VectorZ = new Type[ColumnMatr];
            //for(int i = 0; i < ColumnMatr; i++)
            //  VectorZ[i] = 0;
 
            // Решение системы соотв. нижнетреугольной матр. с правой частью VectorB.
 
            //Type Sum = 0;
            //for(int i = 0; i < LineMatr; i++, Sum = 0)
            //{
            //  for(int j = 0; j < i; j++)
            //      Sum += MatrLow[i][j]*VectorZ[j];
 
            //  VectorZ[i] = VectorB[i] - Sum;
            //}
 
            // Решение системы соотв. верхнетреугольной матр. с правой частью VectorZ.
            
            //Sum = 0;
            //for(int i = LineMatr - 1; i >= 0; i--, Sum = 0)
            //{
            //  for(int j = i + 1; j < ColumnMatr; j++)
            //      Sum += MatrUp[i][j]*VectorX[j];
 
            //  VectorX[i] = (VectorZ[i] - Sum)/MatrUp[i][i];
            //}
 
            //Вернуть вектор решения в вызывающую программу.
 
            //return VectorX;
        }
Добавлено через 11 минут
go :

C++
1
2
3
4
5
6
7
8
9
10
11
void f1 (type ***matr1,...)
{ 
...
 
 
*matr1=matr2;
...
 
 
}
f1 (&matr1,...);
- работает
, но я не могу понять, почему мое не работает
Сыроежка
Заблокирован
21.11.2011, 20:37     Динамическое выделение памяти #10
Цитата Сообщение от Alberto_Timakov Посмотреть сообщение
Почему распечатка матрицы Matr1 в f2 дает матрицу из ед., а в f1 матрицу из нулей? Ведь адрес первого эл. массива Matr2 сохранен в переменной Matr1, и динамическая память после завершения функции не очищается?

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
f2(Type** Matr1, int LineMatr, int ColumnMatr)
{
 
...
 
Type** Matr2 = new Type* [LineMatr];
for(int j = 0; j < ColumnMatr; j++)
     {
          Matr2[j] = new Type [ColumnMatr];
     }
 
for(int i = 0; i < LineMatr; i++)
    for(int j = 0; j < ColumnMatr; j++)
        Matr2[i][j] = 1;
 
...
 
Matr1 = Matr2;
 
printf(Matr1);
 
...
 
}
 
f1()
{
 
...
 
Type** Matr1 = new Type* [LineMatr];
for(int j = 0; j < ColumnMatr; j++)
     {
          Matr1[j] = new Type [ColumnMatr];
     }
 
for(int i = 0; i < LineMatr; i++)
    for(int j = 0; j < ColumnMatr; j++)
        Matr1[i][j] = 0;
 
...
 
f2(Matr1, LineMatr, ColumnMatr);
 
printf(Matr1);
 
...
 
}
Я думаю, что вы некорректно представили свой псевдокод. То есть он не соответсвует тому коду, на основе которого вы привели здесь свой псевдокод.
Чтобы убедиться в этом. вы сами переведите свой псевдокод в нормальный код, и у вас все должно печаться, как надо.
go
Эксперт C++
3582 / 1362 / 128
Регистрация: 16.04.2009
Сообщений: 4,528
21.11.2011, 20:39     Динамическое выделение памяти #11
Цитата Сообщение от Alberto_Timakov Посмотреть сообщение
но я не могу понять, почему мое не работает
вы присваивали копии указателя адрес на первый элемент матрицы, а я указателю
Alberto_Timakov
1 / 1 / 0
Регистрация: 21.11.2011
Сообщений: 183
21.11.2011, 20:56  [ТС]     Динамическое выделение памяти #12
Возможно) А я уже скинул действующий код, который не верно работает.

Добавлено через 13 минут
Непонял, почему копия? Объясните поподробней, пожалуйста)
go
Эксперт C++
3582 / 1362 / 128
Регистрация: 16.04.2009
Сообщений: 4,528
21.11.2011, 21:01     Динамическое выделение памяти #13
Цитата Сообщение от Alberto_Timakov Посмотреть сообщение
Непонял, почему копия? Объясните поподробней, пожалуйста)
ну в функцию передается копия переменной, и изменять мы можем только копию, находясь в функции (или нужно использовать глобальные переменные, что в вашем случае не лучший вариант)
Alberto_Timakov
1 / 1 / 0
Регистрация: 21.11.2011
Сообщений: 183
21.11.2011, 21:03  [ТС]     Динамическое выделение памяти #14
Matr1 не по значению же передается, а по ссылке...
Это же матрица...
go
Эксперт C++
3582 / 1362 / 128
Регистрация: 16.04.2009
Сообщений: 4,528
21.11.2011, 21:05     Динамическое выделение памяти #15
Alberto_Timakov, а я про значение ничего не говорил
Передается Копия ссылки
Alberto_Timakov
1 / 1 / 0
Регистрация: 21.11.2011
Сообщений: 183
22.11.2011, 01:13  [ТС]     Динамическое выделение памяти #16
Ну, если я тебя правильно понял, то следующая замена, не должна была бы изменить результат распечатки в f1.

Если я заменю строку
C++
1
Matr1 = Matr2;
на

C++
1
2
3
4
5
6
7
    for(int i = 0; i < ColumnMatr; i++)
    {
          for(int j = 0; j < ColumnMatr; j++)
          {
        MatrLow[i][j] = 1;
          }
    }
то результат распечатки в f1 - матрица из ед. Если, я что-то недогоняю, то объясни поподробней, пожалуйста))) я хочу понять, но никак не могу)

Добавлено через 44 минуты
Просто для меня "Передается Копия ссылки" звучит, как работа не с исходным массивом, а с его копией, чего быть не может...
go
Эксперт C++
3582 / 1362 / 128
Регистрация: 16.04.2009
Сообщений: 4,528
22.11.2011, 13:46     Динамическое выделение памяти #17
Цитата Сообщение от Alberto_Timakov Посмотреть сообщение
Просто для меня "Передается Копия ссылки" звучит, как работа не с исходным массивом, а с его копией, чего быть не может...
у вас есть исходный указатель на массив, когда вы его передаете в функцию, то в функции создается другой указатель, но на этот же массив
Вывод: вы можете менять значения массива, но не можете менять значение исходного указателя.
вообще если это не поймете, функциями лучше не пользоваться
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
23.11.2011, 00:12     Динамическое выделение памяти
Еще ссылки по теме:

Динамическое выделение памяти C++
C++ Динамическое выделение памяти
Динамическое выделение памяти C++

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

Или воспользуйтесь поиском по форуму:
Alberto_Timakov
1 / 1 / 0
Регистрация: 21.11.2011
Сообщений: 183
23.11.2011, 00:12  [ТС]     Динамическое выделение памяти #18
и то и другое понял)
Yandex
Объявления
23.11.2011, 00:12     Динамическое выделение памяти
Ответ Создать тему
Опции темы

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