I'm with a problem when I try to write 3D int array using HDF.
The problem happens when I use malloc to dynamically allocate the 3D int
array, the hdf5 file shows wrong numbers.
When I allocate 3D int arrays with: *int data[NX][NY][NZ]* all goes fine.
I try to print on screen all data with printf and I got correct answers.
I would appreciate some advices
Here is the source code:
*/* *
* * This example writes data to the HDF5 file.*
* * Data conversion is performed during write operation. *
* */*
* *
*#include<hdf5.h>*
*#include <stdio.h>*
*#include <stdlib.h>*
*
*
*#define FILE "SDS16.h5"*
*#define DATASETNAME "IntArray" *
*#define NX 5 /* dataset dimensions */*
*#define NY 5*
*#define NZ 5*
*#define RANK 3*
*
*
*int*** allocate3D(int l,int m,int n);*
*void deallocate3D(int*** arr3D,int l,int m);*
*
*
*
*
*int main (void)*
*{*
* hid_t file, dataset; /* file and dataset handles */*
* hid_t datatype, dataspace; /* handles */*
* hsize_t dimsf[3]; /* dataset dimensions */*
* herr_t status; *
* int i,j,k;*
*// int data[NX][NY][NZ];*
* int*** data; /* data to write */*
* *
*
*
* //Data and output buffer initialization. *
* data = allocate3D(NX,NY,NZ);*
*
*
* /**
* * Create a new file using H5F_ACC_TRUNC access,*
* * default file creation properties, and default file*
* * access properties.*
* */*
* file = H5Fcreate(FILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);*
* /**
* * Describe the size of the array and create the data space for fixed*
* * size dataset. *
* */*
* dimsf[0] = NX;*
* dimsf[1] = NY;*
* dimsf[2] = NZ;*
* dataspace = H5Screate_simple(RANK, dimsf, NULL); *
* *
* //Define datatype for the data in the file. We will store little
endian INT numbers.*
* datatype = H5Tcopy(H5T_NATIVE_INT);*
* status = H5Tset_order(datatype, H5T_ORDER_LE);*
*
*
* /*Create a new dataset within the file using defined dataspace and*
* * datatype and default dataset creation properties.*/*
* *
* dataset = H5Dcreate(file, DATASETNAME, datatype, dataspace,
H5P_DEFAULT);*
* //Write the data to the dataset using default transfer properties.*
*
*
* *
* for (k = 0; k < NZ; k++) *
* {*
* for (j = 0; j < NY; j++) *
* {*
* for (i = 0; i < NX; i++)*
* {*
* data[i][j][k] = k;//campo[i*j*NX];*
* }*
* }*
* }*
* *
* *
* status = H5Dwrite(dataset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL,
H5P_DEFAULT, data);*
* *
* deallocate3D(data,NX,NY);*
* *
* *
* //Close/release resources.*
* *
* H5Sclose(dataspace);*
* H5Tclose(datatype);*
* H5Dclose(dataset);*
* H5Fclose(file);*
* *
* return 0;*
*}*
*
*
*
*
*//allocate a 3D array*
*int ***allocate3D(int l,int m,int n)*
*{*
*int ***arr3D;*
*int i,j,k;*
*
*
*arr3D = (int***)malloc(l * sizeof(int **));*
*
*
*for(i=0;i<l;i++)*
*{*
* arr3D[i] = (int**)malloc(m * sizeof(int*));*
* for(j=0;j<m;j++)*
* {*
* arr3D[i][j] = (int*)malloc(n*sizeof(int));*
* }*
*}*
*
*
*return arr3D;*
*}*
*
*
*//deallocate a 3D array*
*void deallocate3D(int*** arr3D,int l,int m)*
*{*
* int i,j;*
*
*
* for(i=0;i<l;i++)*
* {*
* for(j=0;j<m;j++)*
* {*
* free(arr3D[i][j]);*
* }*
* free(arr3D[i]);*
* }*
* free(arr3D);*
*} *
_______________________________________________
Hdf-forum is for HDF software users discussion.
[email protected]
http://mail.hdfgroup.org/mailman/listinfo/hdf-forum_hdfgroup.org
