IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Mar 11, 2016, 10:23:42 PM (10 years ago)
Author:
eugene
Message:

modify to pass with extremely pedantic build; force consistency for signed vs unsigned and int sizes; various relastro updates

Location:
trunk/Ohana
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/Ohana

  • trunk/Ohana/src/libfits/table/F_uncompress_T.c

    r39358 r39457  
    1313int gfits_dump_raw_table (FTable *table, char *message) {
    1414# if (VERBOSE_DUMP)
    15   int i;
    1615  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++) {
    1817    // for (i = 0; i < 16; i++) {
    1918    fprintf (stderr, "%02hhx", table->buffer[i]);
     
    2221  }
    2322  fprintf (stderr, "\n");
     23# else
     24  OHANA_UNUSED_PARAM(table);
     25  OHANA_UNUSED_PARAM(message);
    2426# endif
    2527  return TRUE;
     
    2830int gfits_dump_cmp_table (FTable *table, char *message) {
    2931# if (VERBOSE_DUMP)
    30   int i;
    3132  fprintf (stderr, "%s pntr:\n", message);
    32   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++) {
    3334    fprintf (stderr, "%02hhx", table->buffer[i]);
    3435    if (i % 2) fprintf (stderr, " ");
     
    3839
    3940  fprintf (stderr, "%s data:\n", message);
    40   for (i = 0; i < table->header->pcount; i++) {
     41  for (int i = 0; i < table->header->pcount; i++) {
    4142    fprintf (stderr, "%02hhx", table->buffer[table->heap_start + i]);
    4243    if (i % 2) fprintf (stderr, " ");
     
    4445  }
    4546  fprintf (stderr, "\n");
     47# else
     48  OHANA_UNUSED_PARAM(table);
     49  OHANA_UNUSED_PARAM(message);
    4650# endif
    4751  return TRUE;
     
    7781int gfits_uncompress_table (FTable *srctable, FTable *tgttable) {
    7882
    79   int i;
    80 
    8183  char keyword[256];
    8284
     
    114116
    115117  ALLOCATE (fields, TableField, Nfields);
    116   for (i = 0; i < Nfields; i++) {
     118  for (int i = 0; i < Nfields; i++) {
    117119    snprintf (keyword, 80, "TTYPE%d", i+1);
    118120    if (!gfits_scan (srcheader, keyword, "%s", 1, fields[i].ttype)) ESCAPE;
     
    170172  // calling with an empty data array, data values are not actually copied (this is done
    171173  // below after decompression)
    172   for (i = 0; i < Nfields; i++) {
     174  for (int i = 0; i < Nfields; i++) {
    173175    if (!gfits_set_bintable_column (tgtheader, tgttable, fields[i].ttype, NULL, Ny)) ESCAPE;
    174176  }
     
    200202  gettimeofday (&startTimer, (void *) NULL);
    201203
    202   off_t row;
    203   for (row = 0; row < srcheader->Naxis[1]; row++) {
    204 
    205     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++) {
    206207
    207208      // gfits_uncompress_data can take values specific to the compression mode
     
    218219   
    219220      if (VERBOSE_DUMP) {
    220         int k;
    221221        fprintf (stderr, "u4: ");
    222         for (k = 0; k < Nzdata; k++) {
     222        for (int k = 0; k < Nzdata; k++) {
    223223          fprintf (stderr, "%02hhx", zdata[k]);
    224224          if (k % 2) fprintf (stderr, " ");
     
    250250 
    251251      if (VERBOSE_DUMP) {
    252         int k;
    253252        fprintf (stderr, "u3: ");
    254         for (k = 0; k < Nzdata; k++) {
     253        for (int k = 0; k < Nzdata; k++) {
    255254          fprintf (stderr, "%02hhx", zdata[k]);
    256255          if (k % 2) fprintf (stderr, " ");
     
    266265      gettimeofday (&startTimer, (void *) NULL);
    267266 
    268       int Nrows = (row == Ntile - 1) ? ztilelast : ztilelen;
    269       int Nraw = Nrows*fields[i].Nvalues*fields[i].pixsize; // expected number of bytes
     267      unsigned long int Nrows = (row == Ntile - 1) ? ztilelast : ztilelen;
     268      unsigned long int Nraw = Nrows*fields[i].Nvalues*fields[i].pixsize; // expected number of bytes
    270269      if (!gfits_uncompress_data (zdata, Nzdata, fields[i].zctype, NULL, NULL, 0, raw, &Nraw, Nraw_alloc, fields[i].pixsize)) ESCAPE;
    271270
    272271      if (VERBOSE_DUMP) {
    273         int k;
    274272        fprintf (stderr, "u2: ");
    275         for (k = 0; k < Nraw*fields[i].pixsize; k++) {
     273        for (unsigned long int k = 0; k < Nraw*fields[i].pixsize; k++) {
    276274          fprintf (stderr, "%02hhx", raw[k]);
    277275          if (k % 2) fprintf (stderr, " ");
     
    299297 
    300298      if (VERBOSE_DUMP) {
    301         int k;
    302299        fprintf (stderr, "u1: ");
    303         for (k = 0; k < Nraw*fields[i].pixsize; k++) {
     300        for (unsigned long int k = 0; k < Nraw*fields[i].pixsize; k++) {
    304301          fprintf (stderr, "%02hhx", raw[k]);
    305302          if (k % 2) fprintf (stderr, " ");
     
    356353int gfits_distribute_table_data (FTable *table, TableField *field, char *raw, int row_start, int Nrows) {
    357354
    358   off_t i;
    359 
    360355  // we are copying NN rows into the column which starts at XX and has MM bytes per row
    361356
     
    367362
    368363  if (VERBOSE) fprintf (stderr, "distribute: ");
    369   for (i = 0; i < Nrows; i++, tblbuffer += Nx) {
     364  for (off_t i = 0; i < Nrows; i++, tblbuffer += Nx) {
    370365    memcpy (tblbuffer, &raw[i*rowsize], rowsize);
    371366# if (VERBOSE)
    372     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]);
    373368# endif
    374369  }
     
    379374
    380375int gfits_distribute_table_gzp2 (FTable *table, TableField *field, char *raw, int row_start, int Nrows) {
    381 
    382   off_t i, j, k;
    383376
    384377  // we are copying NN rows into the column which starts at XX and has MM bytes per row
     
    391384  char *rawptr = raw;
    392385
    393   for (k = 0; k < field->pixsize; k++) {
     386  for (off_t k = 0; k < field->pixsize; k++) {
    394387# ifdef BYTE_SWAP     
    395388    char *tblptr_start = &table->buffer[Nx*row_start + offset + (pixsize - k - 1)];
     
    397390    char *tblptr_start = &table->buffer[Nx*row_start + offset + k];
    398391# endif
    399     for (i = 0; i < Nrows; i++, tblptr_start += Nx) {
     392    for (off_t i = 0; i < Nrows; i++, tblptr_start += Nx) {
    400393      char *tblptr = tblptr_start;
    401       for (j = 0; j < Nvalues; j++, tblptr += pixsize, rawptr++) {
     394      for (off_t j = 0; j < Nvalues; j++, tblptr += pixsize, rawptr++) {
    402395        *tblptr = *rawptr;
    403396      }
     
    409402int gfits_distribute_table_gzp2_alt (FTable *table, TableField *field, char *raw, int row_start, int Nrows) {
    410403 
    411   off_t i, j, k;
    412  
    413404  // we are copying NN rows into the column which starts at XX and has MM bytes per row
    414405 
    415406  off_t Nx = table->header->Naxis[0];
    416407 
    417   for (k = 0; k < field->pixsize; k++) {
    418     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++) {
    419410      int row = row_start + i;
    420411      char *rawptr = &raw[i*field->Nvalues + k*Nrows*field->Nvalues];
     
    424415      char *tblptr = &table->buffer[Nx*row + field->offset + k];
    425416# endif
    426       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++) {
    427418        *tblptr = *rawptr;
    428419      }
Note: See TracChangeset for help on using the changeset viewer.