С Новым годом! Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
 
Рейтинг: Рейтинг темы: голосов - 12, средняя оценка - 4.67
Mike_Texnik
4 / 4 / 0
Регистрация: 25.03.2011
Сообщений: 28
#1

Резервирование памяти/освобождение памяти для трехмерного массива - C++

17.10.2012, 13:37. Просмотров 1605. Ответов 17
Метки нет (Все метки)

Необходимо создать трехмерный массив (A[Nx][Ny][Nz]), в котором элементы вдоль направления Z выли бы выровнены по 16 байт.

Есть две проблемы:

1. Как скрестить строку вроде
C++
1
float *Uz __attribute__ ((aligned(16))) = (float*)_mm_malloc(Nz * sizeof(float), 16);
С кодом типа (то есть как корректно подставить выше написанную строчку вместо строки 8):

Кликните здесь для просмотра всего текста
C++
1
2
3
4
5
6
7
8
9
float ***U = (float ***)malloc(Nx * sizeof(float **));
     for(i = 0; i<=Nx; i++)
     {
         U[i] = (float **)malloc(Ny * sizeof(float *));
         for(j = 0; j <=Ny; j++)
         {
                      U[i][j] = (float *) malloc(Nz * sizeof(float));
         }
     }


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
#include <iostream>
#include <malloc.h>
 
const int Nx=108,Ny=103,Nz=250;
float ***U = (float ***)malloc(Nx * sizeof(float **));
float ***U0 = (float ***)malloc(Nx * sizeof(float **));
int k,i,j;
 
int main()
{
 
    //выделяем память
     for(i = 0; i<=Nx; i++)
     {
         U[i] = (float **)malloc(Ny * sizeof(float *));
         U0[i] = (float **)malloc(Ny * sizeof(float *));
         for(j = 0; j <=Ny; j++)
         {
               U[i][j] = (float *) malloc(Nz * sizeof(float));
               U0[i][j] = (float *) malloc(Nz * sizeof(float));
         }
     }
    
 
    //очищаем память
     for(i = 0; i<=Nx; i++)
     {
         for(j = 0; j <=Ny; j++)
         {
               free(U[i][j]);
               free(U0[i][j]);
         }
     }
 
    free (U); free(U0);
    return 0;
}


Вылетает сообщение:
0
Миниатюры
Резервирование памяти/освобождение памяти для трехмерного массива  
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
17.10.2012, 13:37
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Резервирование памяти/освобождение памяти для трехмерного массива (C++):

Динамическая загрузка и освобождение памяти для двумерного массива - C++
Всем добрый день:victory: У меня возник небольшой вопрос , и надеюсь что отзывчивые форумчане 8-) помогут мне с ним разобраться. ...

Как правильно записать в виде функции выделение памяти для двумерного массива и ее освобождение - C++
Здравствуйте! Подскажите, как правильно записать в виде функции выделение памяти для двумерного массива и ее освобождение. ...

Освобождение памяти динамического массива структур - C++
Поставлена следующая задача: сделать любую структуру, сгененировать в неё рандомные данные, записать в бинарник. Потом считать из...

Освобождение памяти после динамического массива - C++
Есть динамический массив, но почему-то на строке delete FirstPoint; студия (2010) пишет что-то про попытку записи в память после конца кучи...

Освобождение памяти динамического массива. Деструктор - C++
Почему выдает ошибку при написании деструктора? Если его убрать, то все работает. #pragma once #ifndef MATRIX_H #define...

Освобождение памяти у динамического массива char - C++
Добрый день. Имеем код: class _ArrayChar { typedef int Ivalue; typedef char* Pchar; typedef const char Cchar;

17
alsav22
5428 / 4823 / 442
Регистрация: 04.06.2011
Сообщений: 13,587
17.10.2012, 13:47 #2
Что в глаза бросается : for(i = 0; i<=Nx; i++). Нужно: for(i = 0; i<Nx; i++)
0
Mike_Texnik
4 / 4 / 0
Регистрация: 25.03.2011
Сообщений: 28
17.10.2012, 13:49  [ТС] #3
если не секрет, то почему так? планируется использовать элементы с 0 по Nx
0
Croessmah
Ушел
Эксперт CЭксперт С++
13558 / 7708 / 872
Регистрация: 27.09.2012
Сообщений: 18,996
Записей в блоге: 3
Завершенные тесты: 1
17.10.2012, 13:53 #4
Память освобождать в порядке, обратном выделению.
и sizeof(float **) == sizeof(float *)

Добавлено через 1 минуту
Цитата Сообщение от Mike_Texnik Посмотреть сообщение
если не секрет, то почему так? планируется использовать элементы с 0 по Nx
Вы выделяете Nx элементов, а индексация в массиве начинается с 0, поэтому последний индекс будет (Nx-1)
1
Mike_Texnik
4 / 4 / 0
Регистрация: 25.03.2011
Сообщений: 28
17.10.2012, 13:55  [ТС] #5
Croessmah, большое спасибо! а сможете подсказать как быть с выравниванием?
0
Croessmah
Ушел
Эксперт CЭксперт С++
13558 / 7708 / 872
Регистрация: 27.09.2012
Сообщений: 18,996
Записей в блоге: 3
Завершенные тесты: 1
17.10.2012, 14:00 #6
Здесь у Вас выход за пределы массива из-за i<=Nx
C++
1
2
3
4
5
6
7
8
9
     for(i = 0; i<=Nx; i++)
     {
         for(j = 0; j <=Ny; j++)
         {
               free(U[i][j]);
               free(U0[i][j]);
         }
     } 
    free (U); free(U0);
и утечка памяти.

C++
1
2
3
4
5
6
7
8
9
10
11
     for(i = 0; i<Nx; i++)
     {
         for(j = 0; j <Ny; j++)
         {
               free(U[i][j]);
               free(U0[i][j]);
         }
               free(U[i]);
               free(U0[i]);
     } 
    free (U); free(U0);
Добавлено через 3 минуты
Цитата Сообщение от Mike_Texnik Посмотреть сообщение
а сможете подсказать как быть с выравниванием?
что-то я самого смысла задания не догоняю
0
Mike_Texnik
4 / 4 / 0
Регистрация: 25.03.2011
Сообщений: 28
17.10.2012, 14:02  [ТС] #7
смысл в том, чтобы воспользоваться радостью SIMD (типа SSE) для действий над элементами массива по направлению Z
0
castaway
Эксперт С++
4916 / 3024 / 370
Регистрация: 10.11.2010
Сообщений: 11,081
Записей в блоге: 10
Завершенные тесты: 1
17.10.2012, 14:03 #8
C
1
2
3
4
...
void *ptr = malloc( size + 15 );
void *ptr16 = ((char *)ptr + 15) & ~0x0f;
...
0
Mike_Texnik
4 / 4 / 0
Регистрация: 25.03.2011
Сообщений: 28
17.10.2012, 14:09  [ТС] #9
lazybiz, извиняюсь, но не очень понял Вашу мысль. как это можно прикрутить к строкам массива?

Цитата Сообщение от lazybiz Посмотреть сообщение
C
1
2
3
4
...
void *ptr = malloc( size + 15 );
void *ptr16 = ((char *)ptr + 15) & ~0x0f;
...
0
castaway
Эксперт С++
4916 / 3024 / 370
Регистрация: 10.11.2010
Сообщений: 11,081
Записей в блоге: 10
Завершенные тесты: 1
17.10.2012, 14:25 #10
C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
#include <stdio.h>
#include <malloc.h>
 
void *aligned_malloc( size_t size, unsigned alignment, void **ptr )
{
    *ptr = malloc( size + alignment - 1 );
    return (void *)(((size_t)*ptr + (alignment - 1)) & ~(alignment - 1));
}
 
int main()
{
    void *  p;
    void *  p16;
 
    p16 = aligned_malloc( 1000, 16, &p );
 
    printf( "%08x\n", (unsigned int)p );
    printf( "%08x\n", (unsigned int)p16 );
 
    free( p );
 
    return 0;
}
0
Mike_Texnik
4 / 4 / 0
Регистрация: 25.03.2011
Сообщений: 28
17.10.2012, 15:31  [ТС] #11
lazybiz, большое спасибо! я, наверное, нифига не понял... но верно ли, что применительно к моей задаче это будет выглядеть как:
Кликните здесь для просмотра всего текста
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
#include <stdio.h>
#include <malloc.h>
 
const int Nx=108,Ny=103,Nz=254;
float ***U = (float ***)malloc((Nx+1) * sizeof(float **));
float ***U0 = (float ***)malloc((Nx+1) * sizeof(float **));
float ***U0p16 = (float ***)malloc((Nx+1) * sizeof(float **));
 
int k,i,j;
 
void *aligned_malloc( size_t size, unsigned alignment, void **ptr )
{
    *ptr = malloc( size + alignment - 1 );
    return (void *)(((size_t)*ptr + (alignment - 1)) & ~(alignment - 1));
}
 
int main()
{
 
    //выделяем память
     for(i = 0; i<=Nx; i++)
     {
         U[i] = (float **)malloc((Ny+1) * sizeof(float *));
         U0[i] = (float **)malloc((Ny+1) * sizeof(float *));
         U0p16[i] = (float **)malloc((Ny+1) * sizeof(float *));
         for(j = 0; j <=Ny; j++)
         {
               U[i][j] = (float *) malloc((Nz+1) * sizeof(float));
               U0[i][j] = (float *) malloc((Nz+1) * sizeof(float));
               U0p16[i][j] = (float *)malloc((Nz+1) * sizeof(float *));
               U0p16[i][j] = aligned_malloc( (Nz+1), 16, &U0[i][j] );
         }
     }
    //очищаем память
     for(i = 0; i<=Nx; i++)
     {
         for(j = 0; j <=Ny; j++)
         {
               free(U[i][j]);
               free(U0[i][j]);
               free(U0p16[i][j]);
         }
     }
 
    free (U); free(U0); free(U0p16);
    return 0;
}
 
    return 0;
}
0
castaway
Эксперт С++
4916 / 3024 / 370
Регистрация: 10.11.2010
Сообщений: 11,081
Записей в блоге: 10
Завершенные тесты: 1
17.10.2012, 15:54 #12
Ты вообще какую-то глупость написал. Тебе всего то надо заменить в твоём же первом коде последний malloc на aligned_malloc.
0
Mike_Texnik
4 / 4 / 0
Регистрация: 25.03.2011
Сообщений: 28
17.10.2012, 16:07  [ТС] #13
то есть просто вместо строки:
U0[i][j] = (float *) malloc((Nz+1) * sizeof(float));
поставить:
U0[i][j] = aligned_malloc( Nz+1, 16, &U0[i][j] );

и добавить функцию aligned_malloc в программу

Верно?

P.S. заранее очень извиняюсь если глупости спрашиваю
0
castaway
Эксперт С++
4916 / 3024 / 370
Регистрация: 10.11.2010
Сообщений: 11,081
Записей в блоге: 10
Завершенные тесты: 1
17.10.2012, 16:27 #14
Так. aligned_malloc возвращает два указателя, тот который передается как 3-й параметр нужен для того чтобы по нему освобождать память, т.е. тебе для этого нужен еще один массив указателей размером Nz.
И еще, зачем ты прибавляешь единицу к размеру?
0
Mike_Texnik
4 / 4 / 0
Регистрация: 25.03.2011
Сообщений: 28
17.10.2012, 16:46  [ТС] #15
Цитата Сообщение от lazybiz Посмотреть сообщение
Так. aligned_malloc возвращает два указателя, тот который передается как 3-й параметр нужен для того чтобы по нему освобождать память, т.е. тебе для этого нужен еще один массив указателей размером Nz.
то есть сперва нужно сделать массив:

float *U0zp16 = (float *)malloc((Nz+1) * sizeof(float));

а потом затолкать его в aligned_malloc:

U0[i][j] = (float *) aligned_malloc( Nz+1, 16, &U0zp16 );

Цитата Сообщение от lazybiz Посмотреть сообщение
И еще, зачем ты прибавляешь единицу к размеру?
Вся эта радость будет встраиваться в уже готовый и рабочий код в котором обработка массивов идет от 0 до Nz, то есть нужно Nz+1 элементов, при этом само Nz очень часто поминается в ходе работы программы, поэтому его задаю как константу,
0
17.10.2012, 16:46
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
17.10.2012, 16:46
Привет! Вот еще темы с ответами:

Освобождение памяти для члена класса - C++
Добрый вечер :victory: Проблема, чувствую, тривиальная, но что-то я не могу сообразить :coffee2: В моём самописном классе есть член...

Освобождение памяти для объектов производного класса - C++
Всем привет. Не могу найти инфу в интернете, всё перерыл, но понял что случай специфичный у меня. Придумал себе конструкцию, но не уверен,...

Освобождение памяти (удаление массива char) и raised exception class EAccessViolation - C++
Подскажите плиз, есть программа, вот кусок AnsiString inputText; int inputTextLength; inputText= Form1-&gt;Edit1-&gt;Text; ...

Освобождение памяти для динамического масива типа char - C++
Здравствуйте, помогите пожалуйста разобраться с одним моментом. При освобождении памяти двумерного динамического массива типа char, в конце...


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

Или воспользуйтесь поиском по форуму:
15
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.