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

ECC-1403: WIP

This commit is contained in:
Shahram Najm 2022-06-19 14:00:39 +01:00
parent ebddfc7205
commit e7b5750113
1 changed files with 24 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
src/grib_accessor_class_data_apply_boustrophedonic.c
View File

@ -273,7 +273,6 @@ static int unpack_double_element(grib_accessor* a, size_t idx, double* val)
int err = 0; int err = 0;
/* GRIB-564: The index idx relates to codedValues NOT values! */ /* GRIB-564: The index idx relates to codedValues NOT values! */
err = grib_get_size(a->parent->h, "codedValues", &size); err = grib_get_size(a->parent->h, "codedValues", &size);
if (err) if (err)
return err; return err;
@ -293,10 +292,31 @@ static int unpack_double_element(grib_accessor* a, size_t idx, double* val)
static int unpack_double_element_set(grib_accessor* a, const size_t* index_array, size_t len, double* val_array) static int unpack_double_element_set(grib_accessor* a, const size_t* index_array, size_t len, double* val_array)
{ {
Assert(!"unpack_double_element_set: apply boustrophedonic packing "); size_t size = 0, i = 0;
return GRIB_DECODING_ERROR; double* values;
} int err = 0;
/* GRIB-564: The indexes in index_array relate to codedValues NOT values! */
err = grib_get_size(grib_handle_of_accessor(a), "codedValues", &size);
if (err)
return err;
for (i = 0; i < len; i++) {
if (index_array[i] > size) return GRIB_INVALID_ARGUMENT;
}
values = (double*)grib_context_malloc_clear(a->context, size * sizeof(double));
err = grib_get_double_array(grib_handle_of_accessor(a), "codedValues", values, &size);
if (err) {
grib_context_free(a->context, values);
return err;
}
for (i = 0; i < len; i++) {
val_array[i] = values[index_array[i]];
}
grib_context_free(a->context, values);
return GRIB_SUCCESS;
}
static int pack_double(grib_accessor* a, const double* val, size_t* len) static int pack_double(grib_accessor* a, const double* val, size_t* len)
{ {