Erreur de segmentation ...

Bonjour,
j'écris un petit programme en C pour lire un fichier binaire très volumineux, qui contient un vecteur très long, et le transformer en matrice 3 D.

Tou se passe bien jusqu'à ce que je m'attaque à des vraiment gros fichiers, typiquement, quand je depasse environ 500 000 lignes (quand le paramètre nslice depasse 10, et donc le nombre total de ligne pass au dessus de 456020).

Est ce qu'il y a une limite de taille quelque part, que je ne vois pas ?
je ne vois pas où ça bloque...

Merci d'avance

marie

le programme:
#include <iostream>
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <cstdio>
using namespace std;

#define ncar 151
#define nslice 10

int main()
{

// Indices pour trier data
long int i,j,xi,yi,si,s;
long int Ndatatot;

// Données sur le format du fichier

int NdataR,NdataC;
NdataR=NdataC=ncar*ncar;
Ndatatot=ncar*ncar*nslice*2;

cout << "Ndatatot : " << Ndatatot << endl;

double data[ncar*ncar*nslice*2]={0};

double dataout[ncar*ncar*nslice*2]={0};

// Données sur le format de sortie
double FieldR[ncar][ncar][nslice]={0};
double FieldC[ncar][ncar][nslice]={0};


// Ouverture du fichier binaire à traiter
FILE* fichier = NULL;
fichier=fopen("templateTD.out.fld","rb");

// templateTD.out.fld
printf("opened \n");

// Lecture du fichier
if (fichier != NULL)
{
for(i=0;i<Ndatatot;i++)
{
// fscanf(fichier,"%lf", &data[i]);
fread(&data[i],sizeof(double),1,fichier);

}
printf("Data %e, %e, %e \n", data[0], data[1],data[2]);
fclose(fichier);
}


// Lecture des données du fichier et écriture en matrices
for(s=0;s<nslice;s++)
{
si=s;
for(i=0;i<NdataR;i++)
{
xi=floor(i/ncar);
yi=i-ncar*(floor(i/ncar));
// si=0;
FieldR[xi][yi][si]=data[i+s*ncar*ncar*2];
}
cout << "fieldR : " << FieldR[0][0][s] << endl;
// cout << "fieldR : " << FieldR[0][1][0] << endl;
// cout << "fieldR : " << FieldR[0][2][0] << endl;

for(j=NdataR;j<NdataR+NdataC;j++)
{
xi=floor((j-NdataR)/ncar);
yi=j-ncar*(floor((j)/ncar));
// si=0;
FieldC[xi][yi][si]=data[j+s*ncar*ncar*2];
}
// cout << "fieldC : " << FieldC[0][0][0] << endl;
// cout << "fieldC : " << data[NdataR] << endl;
}

// Ecriture de la matrice en une seule colonne
for(s=0;s<nslice;s++)
{
si=s;
for(i=0;i<NdataR;i++)
{
xi=floor(i/ncar);
yi=i-ncar*(floor(i/ncar));
// si=0;
dataout[i+s*ncar*ncar*2]=FieldR[xi][yi][si];
// cout << "datai " << dataout[i] << endl;

// cout << "fieldR : " << FieldR[xi][yi][si] << endl;

// cout << "dataout : " << dataout[0] << endl;
}


for(j=NdataR;j<NdataR+NdataC;j++)
{
xi=floor((j-NdataR)/ncar);
yi=j-ncar*(floor(j/ncar));
// si=0;
dataout[j+s*ncar*ncar*2]=FieldC[xi][yi][si];
}
}

// Ecriture de la colonne dans un fichier
cout << "Ecriture" << endl;

FILE* fichier2 = NULL;

fichier2=fopen("dataout.txt","w");
if (fichier2 != NULL)
{
for(i=0;i<Ndatatot;i++)
{
fprintf(fichier2,"%e \n", dataout[i]);
}
fprintf(fichier2,"\n");

fclose(fichier2);
}
printf("Data %e, %e, %e \n", dataout[0], dataout[1],dataout[2]);

return 0;
}