cwpearson/astaroth
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Determine endianness.

Browse Source
This commit is contained in:
Miikka Vaisala
2019-09-24 12:52:52 +08:00
parent 88c3cdf0a5
commit d2e13d0139
2 changed files with 15 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
analysis/python/astar/data/read.py
View File

@@ -33,8 +33,8 @@ def read_bin(fname, fdir, fnum, minfo, numtype=np.longdouble):
timestamp = array[0] timestamp = array[0]
array = np.reshape(array[1:], (minfo.contents['AC_mx'], array = np.reshape(array[1:], (minfo.contents['AC_mx'],
minfo.contents['AC_my'], minfo.contents['AC_my'],
minfo.contents['AC_mz']), order='F') minfo.contents['AC_mz']), order='F')
else: else:
array = None array = None
timestamp = None timestamp = None

17
src/standalone/simulation.cc
View File

@@ -60,8 +60,17 @@ write_mesh_info(const AcMeshInfo* config)
infotxt = fopen("mesh_info.list", "w"); infotxt = fopen("mesh_info.list", "w");
// Determine endianness
unsigned int EE = 1;
char *CC = (char*) ⅇ
const int endianness = (int) *C;
// endianness = 0 -> big endian
// endianness = 1 -> little endian
fprintf(infotxt, "size_t %s %lu \n", "AcRealSize", sizeof(AcReal)); fprintf(infotxt, "size_t %s %lu \n", "AcRealSize", sizeof(AcReal));
fprintf(infotxt, "int %s %i \n", "endian", endianness);
// JP: this could be done shorter and with smaller chance for errors with the following // JP: this could be done shorter and with smaller chance for errors with the following
// (modified from acPrintMeshInfo() in astaroth.cu) // (modified from acPrintMeshInfo() in astaroth.cu)
// MV: Now adapted into working condition. E.g. removed useless / harmful formatting. // MV: Now adapted into working condition. E.g. removed useless / harmful formatting.
@@ -115,13 +124,13 @@ save_mesh(const AcMesh& save_mesh, const int step, const AcReal t_step)
save_ptr = fopen(bin_filename, "wb"); save_ptr = fopen(bin_filename, "wb");
// Start file with time stamp // Start file with time stamp
long double write_long_buf = (long double)t_step; AcReal write_long_buf = (AcReal)t_step;
fwrite(&write_long_buf, sizeof(long double), 1, save_ptr); fwrite(&write_long_buf, sizeof(AcReal), 1, save_ptr);
// Grid data // Grid data
for (size_t i = 0; i < n; ++i) { for (size_t i = 0; i < n; ++i) {
const AcReal point_val = save_mesh.vertex_buffer[VertexBufferHandle(w)][i]; const AcReal point_val = save_mesh.vertex_buffer[VertexBufferHandle(w)][i];
long double write_long_buf = (long double)point_val; AcReal write_long_buf = (AcReal)point_val;
fwrite(&write_long_buf, sizeof(long double), 1, save_ptr); fwrite(&write_long_buf, sizeof(AcReal), 1, save_ptr);
} }
fclose(save_ptr); fclose(save_ptr);
} }