diff options
Diffstat (limited to 'ripser.cpp')
-rw-r--r-- | ripser.cpp | 33 |
1 files changed, 25 insertions, 8 deletions
@@ -830,10 +830,27 @@ enum file_format { BINARY }; -template <typename T> T read(std::istream& input_stream) { +template <typename T> inline void reverse_endianness(T& x) { + uint8_t* p = reinterpret_cast<uint8_t *>(&x); + std::reverse(p, p + sizeof(T)); +} + +template <typename T> T read_le(std::istream& input_stream) { + T result; + input_stream.read(reinterpret_cast<char*>(&result), sizeof(T)); + #ifdef BIGENDIAN + reverse_endianness(result); + #endif + return result; +} + +template <typename T> T read_be(std::istream& input_stream) { T result; input_stream.read(reinterpret_cast<char*>(&result), sizeof(T)); - return result; // on little endian: boost::endian::little_to_native(result); + #ifndef BIGENDIAN + reverse_endianness(result); + #endif + return result; } compressed_lower_distance_matrix read_point_cloud(std::istream& input_stream) { @@ -929,33 +946,33 @@ compressed_lower_distance_matrix read_distance_matrix(std::istream& input_stream } compressed_lower_distance_matrix read_dipha(std::istream& input_stream) { - if (read<int64_t>(input_stream) != 8067171840) { + if (read_le<int64_t>(input_stream) != 8067171840) { std::cerr << "input is not a Dipha file (magic number: 8067171840)" << std::endl; exit(-1); } - if (read<int64_t>(input_stream) != 7) { + if (read_le<int64_t>(input_stream) != 7) { std::cerr << "input is not a Dipha distance matrix (file type: 7)" << std::endl; exit(-1); } - index_t n = read<int64_t>(input_stream); + index_t n = read_le<int64_t>(input_stream); std::vector<value_t> distances; for (int i = 0; i < n; ++i) for (int j = 0; j < n; ++j) if (i > j) - distances.push_back(read<double>(input_stream)); + distances.push_back(read_le<double>(input_stream)); else - read<double>(input_stream); + read_le<double>(input_stream); return compressed_lower_distance_matrix(std::move(distances)); } compressed_lower_distance_matrix read_binary(std::istream& input_stream) { std::vector<value_t> distances; - while (!input_stream.eof()) distances.push_back(read<value_t>(input_stream)); + while (!input_stream.eof()) distances.push_back(read_le<value_t>(input_stream)); return compressed_lower_distance_matrix(std::move(distances)); } |