core
/
eccodes
mirror of https://github.com/ecmwf/eccodes.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

ECC-865: Implement roundf in grib_iterator_class_latlon.c

This commit is contained in:
Shahram Najm 2019-01-10 12:08:27 +00:00
parent 64173b17c4
commit b0a51bf272
1 changed files with 12 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
src/grib_iterator_class_latlon.c
View File

@ -88,6 +88,16 @@ static void init_class(grib_iterator_class* c)
#define RAD2DEG 57.29577951308232087684 /* 180 over pi */
#define DEG2RAD 0.01745329251994329576 /* pi over 180 */
/* Note: roundf implementation missing on MSVC! */
static float round_float(float x)
{
#ifdef ECCODES_ON_WINDOWS
return x >= 0.0f ? floorf(x + 0.5f) : ceilf(x - 0.5f);
#else
return roundf(x);
#endif
}
void unrotate(grib_handle* h,
const double inlat, const double inlon,
const double angleOfRot, const double southPoleLat, const double southPoleLon,
@ -127,11 +137,8 @@ void unrotate(grib_handle* h,
ret_lon = atan2(y, x) * RAD2DEG;
/* Still get a very small rounding error, round to 6 decimal places */
/* TODO: roundf implementation missing on MSVC! */
#ifndef ECCODES_ON_WINDOWS
ret_lat = roundf( ret_lat * 1000000.0 )/1000000.0;
ret_lon = roundf( ret_lon * 1000000.0 )/1000000.0;
#endif
ret_lat = round_float( ret_lat * 1000000.0 )/1000000.0;
ret_lon = round_float( ret_lon * 1000000.0 )/1000000.0;
ret_lon -= angleOfRot;