static int unpack_long (grib_accessor* a, long* val, size_t *len) { grib_accessor_signed_bits* self = (grib_accessor_signed_bits*)a; int i; int ret=0; long pos = a->offset*8; long rlen = value_count(a); long numberOfBits = 0; if(*len < rlen) { grib_context_log(a->parent->h->context, GRIB_LOG_ERROR, " wrong size (%ld) for %s it contains %d values ",*len, a->name , rlen); *len = 0; return GRIB_ARRAY_TOO_SMALL; } ret=grib_get_long(a->parent->h,self->numberOfBits,&numberOfBits); if (ret) return ret; if (numberOfBits==0) { int i; for (i=0;i<rlen;i++) val[i]=0; return GRIB_SUCCESS; } for (i=0;i<rlen;i++) val[i] = grib_decode_signed_longb(a->parent->h->buffer->data, &pos, numberOfBits); *len = rlen; return GRIB_SUCCESS; }
static int unpack_double(grib_accessor* a, double* val, size_t *len) { grib_accessor_data_2order_packing* self = (grib_accessor_data_2order_packing*)a; size_t i = 0; size_t j = 0; size_t n_vals = 0; int err=0; long vcount = 0, bias = 0; double reference_value; long binary_scale_factor; long bits_per_value, decimal_scale_factor; long n1 = 0, n2 = 0, extraValues = 0, p1 = 0, p2 = 0; long offsetsection = 0, snd_bitmap = 0, snd_ordr_wdiff = 0; long matrix_values =0; long general_ext =0; long boustrophedonic =0; long two_ordr_spd =0; long plus1_spd =0; long nbits_per_width =0; long nbits_per_group_size =0; long octet_start_group =0; long width_spd_sp_desc =0; grib_handle* gh = grib_handle_of_accessor(a); unsigned char* buf_size_of_groups = (unsigned char*)gh->buffer->data; unsigned char* buf_width_of_group = (unsigned char*)gh->buffer->data; unsigned char* bufrefs = (unsigned char*)gh->buffer->data; unsigned char* bufvals = (unsigned char*)gh->buffer->data; double s = 0; double d = 0; double max = 0; double min = 0; unsigned long* sec_val = NULL; long ref_vals = 0; short n_sp_diff = 0; short f_size_of_group = 0; short f_width_of_group = 0; long bitp = 0; long pointer_of_group_size = 0; long pointer_of_group_width = 0; long refsp = 0; long nap = 0; long nn=0; unsigned char* bitmap=NULL; grib_accessor* abitmap=NULL; size_t bitmap_len=0; err=grib_value_count(a,&nn); n_vals=nn; if((err = grib_get_long_internal(gh,self->offsetsection,&offsetsection)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->bits_per_value,&bits_per_value)) != GRIB_SUCCESS) return err; if((err = grib_get_double_internal(gh,self->reference_value, &reference_value)) != GRIB_SUCCESS)return err; if((err = grib_get_long_internal(gh,self->binary_scale_factor, &binary_scale_factor)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->decimal_scale_factor, &decimal_scale_factor)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->n1,&n1)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->n2, &n2)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->p1, &p1)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->extraValues, &extraValues)) != GRIB_SUCCESS) return err; p1=p1+65536*extraValues; if((err = grib_get_long_internal(gh,self->p2, &p2)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->matrix_values, &matrix_values)) != GRIB_SUCCESS)return err; if((err = grib_get_long_internal(gh,self->snd_bitmap, &snd_bitmap)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->snd_ordr_wdiff, &snd_ordr_wdiff)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->general_ext, &general_ext)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->boustrophedonic, &boustrophedonic)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->two_ordr_spd, &two_ordr_spd)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->plus1_spd, &plus1_spd)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->width_widths, &nbits_per_width)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->width_lengths, &nbits_per_group_size)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->octet_start_group, &octet_start_group)) != GRIB_SUCCESS) return err; if((err = grib_get_long_internal(gh,self->width_spd_sp_desc, &width_spd_sp_desc)) != GRIB_SUCCESS) width_spd_sp_desc=-1; if((err = grib_get_long_internal(gh,self->nap, &nap)) != GRIB_SUCCESS) return err; self->dirty=0; n_sp_diff = two_ordr_spd*2+plus1_spd; Assert(bits_per_value < (sizeof(unsigned long)*8)-1); if ((abitmap=grib_find_accessor(gh,self->bitmap))!=NULL) { bitmap_len=grib_byte_count(abitmap); bitmap=(unsigned char*)grib_context_malloc_clear(a->context,sizeof(char)*bitmap_len); err=grib_unpack_bytes(abitmap,bitmap,&bitmap_len); if (err) {grib_context_free(a->context,bitmap); return err;} } if(bits_per_value == 0) return GRIB_NOT_IMPLEMENTED; /* I have everything now start decoding */ /* fprintf(stdout,"\n****************************************\n"); fprintf(stdout," bits_per_value = %ld\n", bits_per_value); fprintf(stdout," reference_value = %g\n", reference_value); fprintf(stdout," binary_scale_factor = %ld\n", binary_scale_factor); fprintf(stdout," decimal_scale_factor = %ld\n", decimal_scale_factor); fprintf(stdout," n1 = %ld\n", n1); fprintf(stdout," n2 = %ld\n", n2); fprintf(stdout," p1 = %ld\n", p1); fprintf(stdout," p2 = %ld\n", p2); fprintf(stdout," matrix_values = %ld\n", matrix_values); fprintf(stdout," snd_bitmap = %ld\n", snd_bitmap); fprintf(stdout," snd_ordr_wdiff = %ld\n", snd_ordr_wdiff); fprintf(stdout," general_ext = %ld\n", general_ext); fprintf(stdout," boustrophedonic = %ld\n", boustrophedonic); fprintf(stdout," two_ordr_spd = %ld \n", two_ordr_spd); fprintf(stdout," plus1_spd = %ld\n", plus1_spd); fprintf(stdout," n_sp_diff = %d\n", n_sp_diff); fprintf(stdout," width_widths = %ld\n", nbits_per_group_size); fprintf(stdout," width_lengths = %ld\n", nbits_per_width); fprintf(stdout," octet_start_group = %ld\n", octet_start_group); fprintf(stdout," width_spd_sp_desc = %ld\n", width_spd_sp_desc); fprintf(stdout," offsetsection = %ld\n", offsetsection); fprintf(stdout," offset w = %ld\n", octet_start_group + offsetsection); fprintf(stdout,"\n****************************************\n"); */ if(snd_bitmap || matrix_values) return GRIB_NOT_IMPLEMENTED; sec_val = (unsigned long*)grib_context_malloc(a->context,(n_vals)*sizeof(unsigned long)); buf_width_of_group += a->offset; buf_size_of_groups += offsetsection+(octet_start_group-1); /* -1 because of documented starting at 1(and not 0)*/ bufrefs += offsetsection+n1-1; /* -1 because of documented starting at 1(and not 0)*/ pointer_of_group_size = 0; pointer_of_group_width = 0; refsp = 0; for(i=0;i < n_sp_diff;i++) sec_val[i] = grib_decode_unsigned_long(buf_width_of_group, &pointer_of_group_width, width_spd_sp_desc); bias = grib_decode_signed_longb(buf_width_of_group, &pointer_of_group_width, width_spd_sp_desc); bufvals += offsetsection+n2-1; bitp = 0; vcount = n_sp_diff; if(pointer_of_group_width%8) pointer_of_group_width = 8+(pointer_of_group_width-(pointer_of_group_width%8)); #if KEEP_OLD == 1 if(sd == NULL){ sd = grib_context_malloc_clear(a->context,sizeof(second_order_packed)); sd->packed_byte_count = 0; sd->nbits_per_group_size = nbits_per_group_size; sd->nbits_per_widths = nbits_per_width; sd->size_of_group_array = p1; sd->array_of_group_size = grib_context_malloc_clear(a->context,sizeof(unsigned long)*sd->size_of_group_array); sd->array_of_group_width = grib_context_malloc_clear(a->context,sizeof(unsigned long)*sd->size_of_group_array); sd->array_of_group_refs = grib_context_malloc_clear(a->context,sizeof( long)*sd->size_of_group_array); } #endif for(i=0;i < p1;i++){ f_width_of_group = (short) grib_decode_unsigned_long(buf_width_of_group, &pointer_of_group_width, nbits_per_width); f_size_of_group = (short) grib_decode_unsigned_long(buf_size_of_groups, &pointer_of_group_size, nbits_per_group_size); ref_vals = grib_decode_unsigned_long(bufrefs, &refsp, bits_per_value); #if KEEP_OLD == 1 if(sd->packed_byte_count == 0){ sd->array_of_group_width[i] = f_width_of_group; sd->array_of_group_size[i] = f_size_of_group; sd->array_of_group_refs[i] = ref_vals; } #endif for(j=0; j < f_size_of_group;j++){ sec_val[vcount+j] = ref_vals + grib_decode_unsigned_long(bufvals, &bitp, f_width_of_group); } vcount += f_size_of_group; } #if KEEP_OLD == 1 if(sd->packed_byte_count == 0) sd->packed_byte_count = (bitp+7)/8; #endif Assert (n_vals == vcount); /*for(i=0;i < 10;i++) printf("readvalue [%d] %ld %ld bias %ld <<\n", i,sec_val[i],binary_scale_factor,bias );*/ if(snd_ordr_wdiff) de_spatial_difference(a->context,sec_val, n_vals, n_sp_diff, bias); if(boustrophedonic) reverse_rows(sec_val,n_vals,nap,bitmap,bitmap_len); s = grib_power(binary_scale_factor,2); d = grib_power(-decimal_scale_factor,10) ; for(i=0;i < n_vals;i++) val[i] = (double) ((((double)sec_val[i])*s)+reference_value)*d; max = val[0]; min = max; for(i=0;i< n_vals;i++) { if ( val[i] > max ) max = val[i]; else if ( val[i] < min ) min = val[i]; } min *= d; max *= d; grib_context_free(a->context,sec_val); if (bitmap!=NULL) grib_context_free(a->context,bitmap); return err; }