mirror/hdf5
2
0
Fork 0
mirror of https://github.com/HDFGroup/hdf5.git synced 2025-02-23 16:20:57 +08:00
Code Issues Packages Projects Releases Wiki Activity

[svn-r3346] ./hdf5-1.4/src/H5Tconv.c

Browse Source 
2001-02-03 00:19:25 Robb Matzke  <matzke@llnl.gov>
	* H5T_conv_f_f: Fixed a bug where I had forgotten to increment the floating-point
	   exponent if rounding the significand resulted in a
	   carry. Thanks to Guillaume Colin de Verdiere for
	   finding this one!
This commit is contained in:
Robb Matzke 2001-02-02 23:38:09 -05:00
parent 500e15e77e
commit 5fa7f93263
1 changed files with 36 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

36
src/H5Tconv.c
View File

@ -2523,6 +2523,11 @@ H5T_conv_i_i (hid_t src_id, hid_t dst_id, H5T_cdata_t *cdata, hsize_t nelmts,
* then convert one value at each memory location advancing
* BUF_STRIDE bytes each time; otherwise assume both source and
* destination values are packed.
*
* Robb Matzke, 2001-02-02
* Oops, forgot to increment the exponent when rounding the
* significand resulted in a carry. Thanks to Guillaume Colin
* de Verdiere for finding this one!
*-------------------------------------------------------------------------
*/
herr_t
@ -2852,6 +2857,37 @@ H5T_conv_f_f (hid_t src_id, hid_t dst_id, H5T_cdata_t *cdata, hsize_t nelmts,
}
/* Write the exponent */
if (carry) {
expo++;
if (expo>=expo_max) {
/*
* The exponent is too large to fit in the available
* region or it results in the maximum possible value.
* Use positive or negative infinity instead unless the
* application specifies something else. Before
* calling the overflow handler make sure the source
* buffer we hand it is in the original byte order.
*/
if (H5T_overflow_g) {
uint8_t over_src[256];
assert(src_p->size<=sizeof over_src);
if (H5T_ORDER_BE==src.order) {
for (i=0; i<src_p->size; i++) {
over_src[src_p->size-(i+1)] = s[i];
}
} else {
for (i=0; i<src_p->size; i++) {
over_src[i] = s[i];
}
}
if ((H5T_overflow_g)(src_id, dst_id, over_src, d)>=0) {
goto next;
}
}
expo = expo_max;
H5T_bit_set(d, dst.u.f.mpos, dst.u.f.msize, FALSE);
}
}
H5_CHECK_OVERFLOW(expo,hssize_t,hsize_t);
H5T_bit_set_d(d, dst.u.f.epos, dst.u.f.esize, (hsize_t)expo);