IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
May 3, 2006, 12:01:09 PM (20 years ago)
Author:
eugene
Message:

converted fits_ to gfits_

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r7039 r7054  
    33
    44/*********************** fits read table ***********************************/
    5 int fits_read_ftable (char *filename, FTable *table, char *extname) {
     5int gfits_read_ftable (char *filename, FTable *table, char *extname) {
    66
    77  int status;
    … …  
    1111  if (f == NULL) return (FALSE);
    1212
    13   status = fits_fread_ftable (f, table, extname);
     13  status = gfits_fread_ftable (f, table, extname);
    1414  fclose (f);
    1515  return (status);
    … …  
    1717
    1818/*********************** fits read table ***********************************/
    19 int fits_fread_ftable (FILE *f, FTable *table, char *extname) {
     19int gfits_fread_ftable (FILE *f, FTable *table, char *extname) {
    2020
    2121  int j, found, Nbytes;
    … …  
    2929  for (j = -1; !found; j++) {
    3030    /* load data for this header */
    31     if (!fits_load_header (f, header)) return (FALSE);
     31    if (!gfits_load_header (f, header)) return (FALSE);
    3232
    3333    /* check if this is the correct extension or not */
    3434    bzero (tname, 80);
    35     fits_scan (header, "EXTNAME", "%s", 1, tname);
     35    gfits_scan (header, "EXTNAME", "%s", 1, tname);
    3636    if (!strcmp (tname, extname)) {
    37       if (fits_fread_ftable_data (f, table)) return (TRUE);
    38       fits_free_header (header);
     37      if (gfits_fread_ftable_data (f, table)) return (TRUE);
     38      gfits_free_header (header);
    3939      return (FALSE);
    4040    }
    4141
    4242    /* skip to next header */
    43     Nbytes = fits_matrix_size (header);
     43    Nbytes = gfits_matrix_size (header);
    4444    fseek (f, Nbytes, SEEK_CUR);
    45     fits_free_header (header);
     45    gfits_free_header (header);
    4646  }
    4747  return (TRUE);
    … …  
    4949
    5050/*********************** fits read ftable data ***********************************/
    51 int fits_fread_ftable_data (FILE *f, FTable *table) {
     51int gfits_fread_ftable_data (FILE *f, FTable *table) {
    5252
    5353  int Nbytes, Nread;
    5454
    5555  /* find buffer size */
    56   Nbytes = fits_matrix_size (table[0].header);
     56  Nbytes = gfits_matrix_size (table[0].header);
    5757  ALLOCATE (table[0].buffer, char, Nbytes);
    5858
    5959  Nread = fread (table[0].buffer, sizeof (char), Nbytes, f);
    6060  if (Nread != Nbytes) {
    61     fits_free_table  (table);
     61    gfits_free_table  (table);
    6262    return (FALSE);
    6363  }
    … …  
    6767
    6868/*********************** fits read ftable data ***********************************/
    69 int fits_fread_vtable_range (FILE *f, VTable *table, int start, int Nrows) {
     69int gfits_fread_vtable_range (FILE *f, VTable *table, int start, int Nrows) {
    7070
    7171  int i, Nbytes, Nread, Nskip, Nx, Ny;
    … …  
    7575  Nx = table[0].header[0].Naxis[0];
    7676  Ny = table[0].header[0].Naxis[1];
    77   table[0].size = fits_matrix_size (table[0].header);
     77  table[0].size = gfits_matrix_size (table[0].header);
    7878  table[0].pad = table[0].size - Nx*Ny;
    7979
    … …  
    105105
    106106/*********************** fits read virtual table ***********************************/
    107 int fits_fread_vtable (FILE *f, VTable *table, char *extname, int Nrow, int *row) {
     107int gfits_fread_vtable (FILE *f, VTable *table, char *extname, int Nrow, int *row) {
    108108
    109109  int i, j, Nbytes, Nread;
    … …  
    117117  for (j = -1; TRUE; j++) {
    118118    /* load data for this header */
    119     if (!fits_load_header (f, header)) return (FALSE);
     119    if (!gfits_load_header (f, header)) return (FALSE);
    120120
    121121    /* find buffer size */
    122     Nbytes = fits_matrix_size (header);
     122    Nbytes = gfits_matrix_size (header);
    123123
    124124    /* check if this is the correct extension or not */
    125125    bzero (tname, 80);
    126     fits_scan (header, "EXTNAME", "%s", 1, tname);
     126    gfits_scan (header, "EXTNAME", "%s", 1, tname);
    127127    if (strcmp (tname, extname)) {
    128128      /* skip to next header */
    129129      fseek (f, Nbytes, SEEK_CUR);
    130       fits_free_header (header);
     130      gfits_free_header (header);
    131131      continue;
    132132    }
    … …  
    134134    /* file pointer is at beginning of desired table data */
    135135    start = ftell (f);
    136     fits_scan (header, "NAXIS1", "%d", 1, &Nx);
    137     fits_scan (header, "NAXIS2", "%d", 1, &Ny);
     136    gfits_scan (header, "NAXIS1", "%d", 1, &Nx);
     137    gfits_scan (header, "NAXIS2", "%d", 1, &Ny);
    138138    for (i = 0; i < Nrow; i++) {
    139139      if (row[i] > Ny) { return (FALSE); }
    … …  
    152152    ALLOCATE (table[0].row, int, MAX (1, Nrow));   
    153153    for (i = 0; i < Nrow; i++) table[0].row[i] = row[i];
    154     table[0].size   = fits_matrix_size (table[0].header);
     154    table[0].size   = gfits_matrix_size (table[0].header);
    155155    table[0].pad    = table[0].size - Nx*Ny;
    156156    return (TRUE);
    … …  
    162162/* a poor function: read only the first table, regardless of name */
    163163/* keep until we are sure ohana code does not use it  */
    164 int fits_read_table (char *filename, FTable *table) {
     164int gfits_read_table (char *filename, FTable *table) {
    165165
    166166  FILE *f;
    … …  
    168168  int status, i, nbytes, Nbytes, matrix_size;
    169169
    170   if (!fits_read_header (filename, &header)) {
     170  if (!gfits_read_header (filename, &header)) {
    171171    fprintf (stderr, "error reading header of FITS file %s\n", filename);
    172172    return (FALSE);
    173173  }
    174174
    175   if (!fits_read_Theader (filename, &table[0].header)) {
     175  if (!gfits_read_Theader (filename, &table[0].header)) {
    176176    fprintf (stderr, "error reading Theader of FITS file %s\n", filename);
    177177    return (FALSE);
    … …  
    179179
    180180  if ((f = fopen (filename, "r")) == NULL) {
    181     fits_free_header (&header);
    182     fits_free_header (&table[0].header);
    183     return (FALSE);
    184   }
    185 
    186   matrix_size = fits_matrix_size (&header);
     181    gfits_free_header (&header);
     182    gfits_free_header (&table[0].header);
     183    return (FALSE);
     184  }
     185
     186  matrix_size = gfits_matrix_size (&header);
    187187
    188188  Nbytes = header.size + matrix_size + table[0].header.size;
    189189  fseek (f, Nbytes, 0);
    190   fits_free_header (&header);
     190  gfits_free_header (&header);
    191191 
    192   fits_scan (&table[0].header, "TFIELDS", "%d", 1, &table[0].Nfields);
     192  gfits_scan (&table[0].header, "TFIELDS", "%d", 1, &table[0].Nfields);
    193193
    194194  /* place the field sizes in the structure? */
    … …  
    198198    table[0].Naxis[i] = table[0].header.Naxis[i];
    199199
    200   Nbytes = fits_matrix_size (&table[0].header);
     200  Nbytes = gfits_matrix_size (&table[0].header);
    201201  ALLOCATE (table[0].buffer, char, Nbytes);
    202202  nbytes = fread (table[0].buffer, sizeof(char), Nbytes, f);
    … …  
    213213
    214214
    215 int fits_fread_header_extname (FILE *f, Header *header, char *extname) {
     215int gfits_fread_header_extname (FILE *f, Header *header, char *extname) {
    216216
    217217  int Nbytes;
    … …  
    219219
    220220  fseek (f, 0, SEEK_SET);
    221   fits_fread_header (f, header);
     221  gfits_fread_header (f, header);
    222222
    223223  if (!strcasecmp (extname, "PHU")) return (TRUE);
    224224
    225   Nbytes = fits_matrix_size (header);
     225  Nbytes = gfits_matrix_size (header);
    226226  fseek (f, Nbytes, SEEK_CUR);
    227227
    228   while (fits_fread_header (f, header)) {
    229     fits_scan (header, "EXTNAME", "%s", 1, current);
     228  while (gfits_fread_header (f, header)) {
     229    gfits_scan (header, "EXTNAME", "%s", 1, current);
    230230    if (!strcmp (current, extname)) return (TRUE);
    231     Nbytes = fits_matrix_size (header);
     231    Nbytes = gfits_matrix_size (header);
    232232    fseek (f, Nbytes, SEEK_CUR);
    233233  }
Note: See TracChangeset for help on using the changeset viewer.