Блог

У меня есть очень простой фрагмент кода:

 #include <iostream>
#include <random>

using namespace std;

const int K=1;
const int N=1024*K;
const int M=1024;

class C {
  public:
  C(){
    minstd_rand0 g(1);
    for(int i=0; i<N;   i) for(int j=0; j<M;   j)
        a[i][j]=generate_canonical<float, numeric_limits<float>::digits>(g);
     }
  private:
  float a[N][M];
};

int main()
{
  C c;
  cout<<"size="<<N*M/(1<<20)<<endl;
  return 0;
}
 

При K = 1 код работает, но при K = 2 происходит немедленный сбой из-за ошибки сегментации.
Итак, как я понимаю, существует некоторое ограничение на объем данных, который может содержать произвольный класс (не более этого объема МБ). Как точно узнать эту сумму?

В своей работе я имею дело с массивами, подобными a, которые обозначают распределение функции вероятности на некотором отрезке линии. У меня есть десятки таких массивов. Я предполагал поместить все эти массивы в 1 класс менеджера, который будет обрабатывать их все, но, как показывает этот пример, это невозможно, возникает ошибка сегментации. Итак, у меня есть, например, 10 массивов (я хотел бы инкапсулировать их все внутри какого-то одного объекта) типа double с количеством элементов 2x1024x1024. Как правильно управлять ими в программе?

P. S. Если распределить их по куче:

 class C {
  public:
  C(){
    p = new float*[N];
    for(int i=0; i<N;   i) p[i]=new float[M];
    minstd_rand0 g(1);
    for(int i=0; i<N;   i) for(int j=0; j<M;   j)
        p[i][j]=generate_canonical<float, 
                numeric_limits<float>::digits>(g);
  }
  ~C(){
    for(int i=0; i<N;   i) delete [] p[i];
    delete [] p;
  }
  private:
  float ** p;
};
 

программа работает даже при K=100. Правильно ли это управлять ими?
В конце концов, есть ли способ выделить массив статически, не делая его статическим членом класса?

P. P. S. Этот фрагмент кода работает:

 const int K=10;
const int N=1024*K;
const int M=1024;
float a[N][M];

class C {
  public:
  C(){
    p=(float *)a;
    minstd_rand0 g(1);
    for(int i=0; i<N;   i) for(int j=0; j<M;   j)
      p[i*M j]=generate_canonical<float, 
               numeric_limits<float>::digits> 
               (g);
    }
  public:
  float *p;
};
 

Но я не знаю, правильно ли это обрабатывать такие большие массивы.