caiyu/openEMS
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

ProcessFields fix: hdf5 dumps need dynamic memory allocation

Browse Source
This commit is contained in:
Thorsten Liebig 2011-01-20 14:37:35 +01:00
parent 3f8f0156d9
commit ae3e15770e
1 changed files with 6 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Common/processfields.cpp
View File

@ -449,7 +449,8 @@ bool ProcessFields::DumpVectorArray2HDF5(string filename, string groupName, stri
// I have not the slightest idea why this array-copy action is necessary... but it's the only way hdf5 does what it is supposed to do anyway!!
// at least it is save in case FDTD_FLOAT was defined as double...
// why does hdf5 write the dimensions backwards??? or matlab???
float hdf5array[3][numLines[2]][numLines[1]][numLines[0]];
unsigned int reverse_numLines[] = {numLines[2], numLines[1], numLines[0]};
float ****hdf5array = Create_N_3DArray<float>(reverse_numLines);
for (int n=0; n<3; ++n)
{
for (unsigned int i=0; i<numLines[0]; ++i)
@ -464,6 +465,7 @@ bool ProcessFields::DumpVectorArray2HDF5(string filename, string groupName, stri
}
}
dataset.write( hdf5array, H5::PredType::NATIVE_FLOAT );
Delete_N_3DArray(hdf5array,reverse_numLines);
return true;
}
@ -496,7 +498,8 @@ bool ProcessFields::DumpVectorArray2HDF5(string filename, string groupName, stri
// I have not the slightest idea why this array-copy action is necessary... but it's the only way hdf5 does what it is supposed to do anyway!!
// at least it is save in case FDTD_FLOAT was defined as double...
// why does hdf5 write the dimensions backwards??? or matlab???
float hdf5array[3][numLines[2]][numLines[1]][numLines[0]];
unsigned int reverse_numLines[] = {numLines[2], numLines[1], numLines[0]};
float ****hdf5array = Create_N_3DArray<float>(reverse_numLines);
for (int n=0; n<3; ++n)
{
for (unsigned int i=0; i<numLines[0]; ++i)
@ -533,6 +536,7 @@ bool ProcessFields::DumpVectorArray2HDF5(string filename, string groupName, stri
}
}
dataset.write( hdf5array, H5::PredType::NATIVE_FLOAT );
Delete_N_3DArray(hdf5array,reverse_numLines);
return true;
}