- Timestamp:
- Mar 1, 2016, 12:16:46 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ohana.20160226/src/libfits/table/F_uncompress_T.c
r39409 r39428 13 13 int gfits_dump_raw_table (FTable *table, char *message) { 14 14 # if (VERBOSE_DUMP) 15 int i;16 15 fprintf (stderr, "%s data:\n", message); 17 for (i = 0; i < table->header->Naxis[0]*table->header->Naxis[1]; i++) {16 for (int i = 0; i < table->header->Naxis[0]*table->header->Naxis[1]; i++) { 18 17 // for (i = 0; i < 16; i++) { 19 18 fprintf (stderr, "%02hhx", table->buffer[i]); … … 31 30 int gfits_dump_cmp_table (FTable *table, char *message) { 32 31 # if (VERBOSE_DUMP) 33 int i;34 32 fprintf (stderr, "%s pntr:\n", message); 35 for (i = 0; i < table->header->Naxis[0]*table->header->Naxis[1]; i++) {33 for (int i = 0; i < table->header->Naxis[0]*table->header->Naxis[1]; i++) { 36 34 fprintf (stderr, "%02hhx", table->buffer[i]); 37 35 if (i % 2) fprintf (stderr, " "); … … 41 39 42 40 fprintf (stderr, "%s data:\n", message); 43 for (i = 0; i < table->header->pcount; i++) {41 for (int i = 0; i < table->header->pcount; i++) { 44 42 fprintf (stderr, "%02hhx", table->buffer[table->heap_start + i]); 45 43 if (i % 2) fprintf (stderr, " "); … … 83 81 int gfits_uncompress_table (FTable *srctable, FTable *tgttable) { 84 82 85 int i;86 87 83 char keyword[256]; 88 84 … … 120 116 121 117 ALLOCATE (fields, TableField, Nfields); 122 for (i = 0; i < Nfields; i++) {118 for (int i = 0; i < Nfields; i++) { 123 119 snprintf (keyword, 80, "TTYPE%d", i+1); 124 120 if (!gfits_scan (srcheader, keyword, "%s", 1, fields[i].ttype)) ESCAPE; … … 176 172 // calling with an empty data array, data values are not actually copied (this is done 177 173 // below after decompression) 178 for (i = 0; i < Nfields; i++) {174 for (int i = 0; i < Nfields; i++) { 179 175 if (!gfits_set_bintable_column (tgtheader, tgttable, fields[i].ttype, NULL, Ny)) ESCAPE; 180 176 } … … 206 202 gettimeofday (&startTimer, (void *) NULL); 207 203 208 off_t row; 209 for (row = 0; row < srcheader->Naxis[1]; row++) { 210 211 for (i = 0; i < Nfields; i++) { 204 for (off_t row = 0; row < srcheader->Naxis[1]; row++) { 205 206 for (int i = 0; i < Nfields; i++) { 212 207 213 208 // gfits_uncompress_data can take values specific to the compression mode … … 224 219 225 220 if (VERBOSE_DUMP) { 226 int k;227 221 fprintf (stderr, "u4: "); 228 for ( k = 0; k < Nzdata; k++) {222 for (int k = 0; k < Nzdata; k++) { 229 223 fprintf (stderr, "%02hhx", zdata[k]); 230 224 if (k % 2) fprintf (stderr, " "); … … 256 250 257 251 if (VERBOSE_DUMP) { 258 int k;259 252 fprintf (stderr, "u3: "); 260 for ( k = 0; k < Nzdata; k++) {253 for (int k = 0; k < Nzdata; k++) { 261 254 fprintf (stderr, "%02hhx", zdata[k]); 262 255 if (k % 2) fprintf (stderr, " "); … … 277 270 278 271 if (VERBOSE_DUMP) { 279 unsigned long int k;280 272 fprintf (stderr, "u2: "); 281 for ( k = 0; k < Nraw*fields[i].pixsize; k++) {273 for (unsigned long int k = 0; k < Nraw*fields[i].pixsize; k++) { 282 274 fprintf (stderr, "%02hhx", raw[k]); 283 275 if (k % 2) fprintf (stderr, " "); … … 305 297 306 298 if (VERBOSE_DUMP) { 307 unsigned long int k;308 299 fprintf (stderr, "u1: "); 309 for ( k = 0; k < Nraw*fields[i].pixsize; k++) {300 for (unsigned long int k = 0; k < Nraw*fields[i].pixsize; k++) { 310 301 fprintf (stderr, "%02hhx", raw[k]); 311 302 if (k % 2) fprintf (stderr, " "); … … 362 353 int gfits_distribute_table_data (FTable *table, TableField *field, char *raw, int row_start, int Nrows) { 363 354 364 off_t i;365 366 355 // we are copying NN rows into the column which starts at XX and has MM bytes per row 367 356 … … 373 362 374 363 if (VERBOSE) fprintf (stderr, "distribute: "); 375 for ( i = 0; i < Nrows; i++, tblbuffer += Nx) {364 for (off_t i = 0; i < Nrows; i++, tblbuffer += Nx) { 376 365 memcpy (tblbuffer, &raw[i*rowsize], rowsize); 377 366 # if (VERBOSE) 378 int j; for (j = 0; j < field->rowsize; j++) fprintf (stderr, "0x%02hhx ", table->buffer[Nx*row + field->offset + j]);367 for (int j = 0; j < field->rowsize; j++) fprintf (stderr, "0x%02hhx ", table->buffer[Nx*row + field->offset + j]); 379 368 # endif 380 369 } … … 385 374 386 375 int gfits_distribute_table_gzp2 (FTable *table, TableField *field, char *raw, int row_start, int Nrows) { 387 388 off_t i, j, k;389 376 390 377 // we are copying NN rows into the column which starts at XX and has MM bytes per row … … 397 384 char *rawptr = raw; 398 385 399 for ( k = 0; k < field->pixsize; k++) {386 for (off_t k = 0; k < field->pixsize; k++) { 400 387 # ifdef BYTE_SWAP 401 388 char *tblptr_start = &table->buffer[Nx*row_start + offset + (pixsize - k - 1)]; … … 403 390 char *tblptr_start = &table->buffer[Nx*row_start + offset + k]; 404 391 # endif 405 for ( i = 0; i < Nrows; i++, tblptr_start += Nx) {392 for (off_t i = 0; i < Nrows; i++, tblptr_start += Nx) { 406 393 char *tblptr = tblptr_start; 407 for ( j = 0; j < Nvalues; j++, tblptr += pixsize, rawptr++) {394 for (off_t j = 0; j < Nvalues; j++, tblptr += pixsize, rawptr++) { 408 395 *tblptr = *rawptr; 409 396 } … … 415 402 int gfits_distribute_table_gzp2_alt (FTable *table, TableField *field, char *raw, int row_start, int Nrows) { 416 403 417 off_t i, j, k;418 419 404 // we are copying NN rows into the column which starts at XX and has MM bytes per row 420 405 421 406 off_t Nx = table->header->Naxis[0]; 422 407 423 for ( k = 0; k < field->pixsize; k++) {424 for ( i = 0; i < Nrows; i++) {408 for (off_t k = 0; k < field->pixsize; k++) { 409 for (off_t i = 0; i < Nrows; i++) { 425 410 int row = row_start + i; 426 411 char *rawptr = &raw[i*field->Nvalues + k*Nrows*field->Nvalues]; … … 430 415 char *tblptr = &table->buffer[Nx*row + field->offset + k]; 431 416 # endif 432 for ( j = 0; j < field->Nvalues; j++, tblptr += field->pixsize, rawptr++) {417 for (off_t j = 0; j < field->Nvalues; j++, tblptr += field->pixsize, rawptr++) { 433 418 *tblptr = *rawptr; 434 419 }
Note:
See TracChangeset
for help on using the changeset viewer.
