Pourquoi ai-je une erreur de segmentation (core dumped) lorsque j'essaie d'allouer une matrice 4D?

Question

J'essaye d'allouer une grande matrice 4D mais je veux le faire dynamiquement. Quand je crée simplement la matrice statique, tout fonctionne bien, donc je sais que j'ai suffisamment de mémoire pour l'instant. Cependant, quand j'essaie de mettre en œuvre dynamiquement la même chose, ça casse chaque fois que j'entre dans la troisième dimension et je dois en arriver à une quatrième! Quelqu'un peut-il me dire pourquoi ce code ne fonctionne pas?

#include <iostream> 

using namespace std; 

static const int time1 = 7; 
static const int tlat = 15; 
static const int tlon = 17; 
static const int outlev = 3; 

int main(void) 
{ 
    //allocate four dimensional dataIn 
    int ****dataIn; 
    dataIn = new int ***[time1]; 
    if (dataIn == NULL) { return 1; } 

    for(int i = 0 ; i < time1 ; i++) { 
     dataIn[i] = new int **[tlat]; 
     if (dataIn[i] == NULL) { return 1; } 

     for(int j = 0 ; j < tlat ; j++) { 
      dataIn[i][j] = new int *[tlon]; 
      if (dataIn[i][j] == NULL) { return 1; } 

      for(int k = 0 ; k < tlon ; k++) { 
       dataIn[i][j][k] = new int[outlev]; 
       if (dataIn[i][j][k] == NULL) { return 1; } 
      } 
     } 
    } 
    //there is more code that happens here to add values to dataIn 
    //and eventually output it but I know all of that works   
    return 0; 
} 

J'ai essayé beaucoup de différentes variations sur ce code, et même utilisé malloc au lieu de nouveau, mais je ne peux pas le faire fonctionner. Toute aide est la bienvenue.

Answer 1

Avez-vous exécuté ceci dans un débogueur? À mon avis, le code semble bien, mais un débogueur vous dira où il se bloque, ce qui peut suffire à vous aider à le résoudre.

Answer 2

Vous devriez probablement allouer toute la mémoire dans un tableau plat, puis calculer les indices toi même. Envelopper le tout dans un objet pour l'encapsulation.

class Matrix { 
private: 
     int* data;   
     int[] sizes; 
     int nDimensions; 

public: 
     // allocates the data pointer and copies the other parameters 
     Matrix(int[] sizes, int nDimensions); 

     // frees the data and sizes arrays 
     ~Matrix(); 

     // calculates the cell position and returns it 
     int getCell(int[] coordinates); 

     // calcultes the cell position and sets its value 
     void setCell(int[] coordinates, int value); 

private: 
     // used by getCell and setCell, calculates the cell's 
     // location in the data array 
     size_t calculateCellPosition(int[] coordinates); 
}; 

Answer 3

Il compile et fonctionne correctement sur ma machine Linux. BTW, avec C++ standard, lancerait une exception std :: bad_alloc (au lieu de renvoyer NULL) quand il est incapable d'allouer la mémoire. Il peut donc être intéressant d'attraper cette exception au lieu de tester un pointeur NULL.

Answer 4

Comme souligné par cmeerw, un test sur NULL ou 0 n'est pas nécessaire dans std-C++.

Il serait utile si vous pouviez ajouter un commentaire à l'endroit exact où vous obtenez un segv.

En outre: Quel compilateur utilisez-vous et sur quel système d'exploitation?

Answer 5

Check out boost::multiarray, il fait exactement ce que vous voulez, sauf plus efficacement en ne faisant qu'une seule allocation de tas.