IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Mar 17, 2010, 2:52:38 PM (16 years ago)
Author:
eugene
Message:

more largefile support mods -- code now builds

Location:
branches/eam_branches/largefiles.20100314/Ohana/src/relastro
Files:
10 edited

Legend:

Unmodified
Added
Removed

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/include/relastro.h

    r24308 r27312  
    140140char         *GetPhotnamebyCode   PROTO((PhotCodeData *photcodes, int code));
    141141void          InterpolateGrid     PROTO((float *buffer, int Nx, int Ny, Coords *ccd, Coords *gcoords));
    142 int          *SelectRefMosaic     PROTO((Mosaic **refmosaic, int *Nimage));
     142off_t          *SelectRefMosaic     PROTO((Mosaic **refmosaic, off_t *Nimage));
    143143int           args                PROTO((int argc, char **argv));
    144144int           bcatalog            PROTO((Catalog *subcatalog, Catalog *catalog));
    … …  
    153153void          findImages          PROTO((Catalog *catalog, int Ncatalog));
    154154int           findMosaics         PROTO((Catalog *catalog, int Ncatalog));
    155 Image        *find_images         PROTO((FITS_DB *db, GSCRegion *region, int Nregion, int *Nimage, int **LineNum));
     155Image        *find_images         PROTO((FITS_DB *db, GSCRegion *region, off_t Nregion, off_t *Nimage, off_t **LineNum));
    156156void set_db (FITS_DB *in);
    157157int Shutdown (char *format, ...);
    … …  
    160160int SetSignals (void);
    161161
    162 GSCRegion    *find_regions        PROTO((Image *image, int Nimage, int *Nregions, GSCRegion *fullregion));
     162GSCRegion    *find_regions        PROTO((Image *image, off_t Nimage, int *Nregions, GSCRegion *fullregion));
    163163void          freeGridBins        PROTO((int Ncatalog));
    164164void          freeImageBins       PROTO((int Ncatalog));
    … …  
    166166void          free_catalogs       PROTO((Catalog *catalog, int Ncatalog));
    167167int           gcatalog            PROTO((Catalog *catalog, int FINAL));
    168 Coords       *getCoords           PROTO((int meas, int cat));
    169 float         getMcal             PROTO((int meas, int cat));
    170 float         getMgrid            PROTO((int meas, int cat));
    171 float         getMmos             PROTO((int meas, int cat));
    172 float         getMrel             PROTO((Catalog *catalog, int meas, int cat));
    173 GSCRegion    *get_regions         PROTO((double minRa, double maxRa, double minDec, double maxDec, int *Nregions));
    174 void          getfullregion       PROTO((Image *image, int Nimage, GSCRegion *fullregion));
    175 Image        *getimage            PROTO((int N));
    176 Image        *getimages           PROTO((int *N));
     168Coords       *getCoords           PROTO((off_t meas, int cat));
     169float         getMcal             PROTO((off_t meas, int cat));
     170float         getMgrid            PROTO((off_t meas, int cat));
     171float         getMmos             PROTO((off_t meas, int cat));
     172float         getMrel             PROTO((Catalog *catalog, off_t meas, int cat));
     173GSCRegion    *get_regions         PROTO((double minRa, double maxRa, double minDec, double maxDec, off_t *Nregions));
     174void          getfullregion       PROTO((Image *image, off_t Nimage, GSCRegion *fullregion));
     175Image        *getimage            PROTO((off_t N));
     176Image        *getimages           PROTO((off_t *N));
    177177void          global_stats        PROTO((Catalog *catalog, int Ncatalog));
    178178void          initGrid            PROTO((int dX, int dY));
    179179void          initGridBins        PROTO((Catalog *catalog, int Ncatalog));
    180180void          initImageBins       PROTO((Catalog *catalog, int Ncatalog));
    181 void          initImages          PROTO((Image *input, int N));
     181void          initImages          PROTO((Image *input, off_t N));
    182182void          initMosaicBins      PROTO((Catalog *catalog, int Ncatalog));
    183 void          initMosaicGrid      PROTO((Image *image, int Nimage));
    184 void          initMosaics         PROTO((Image *image, int Nimage));
     183void          initMosaicGrid      PROTO((Image *image, off_t Nimage));
     184void          initMosaics         PROTO((Image *image, off_t Nimage));
    185185void          initMrel            PROTO((Catalog *catalog, int Ncatalog));
    186186void          initialize          PROTO((int argc, char **argv));
    … …  
    189189Catalog      *load_catalogs       PROTO((SkyList *skylist, int *Ncatalog, int subselect));
    190190SkyList      *load_images         PROTO((FITS_DB *db, SkyRegion *region));
    191 Image        *select_images       PROTO((SkyList *skylist, Image *timage, int Ntimage, int **LineNumber, int *Nimage));
     191Image        *select_images       PROTO((SkyList *skylist, Image *timage, off_t Ntimage, off_t **LineNumber, off_t *Nimage));
    192192
    193193void check_permissions (char *basefile);
    … …  
    198198
    199199int           main                PROTO((int argc, char **argv));
    200 void          mark_images         PROTO((Image *image, int Nimage, Image *timage, int Ntimage));
    201 void          matchImage          PROTO((Catalog *catalog, int meas, int cat));
    202 void          matchMosaics        PROTO((Catalog *catalog, int meas, int cat));
    203 GSCRegion    *name_region         PROTO((char *name, int *Nregions));
     200void          mark_images         PROTO((Image *image, off_t Nimage, Image *timage, off_t Ntimage));
     201void          matchImage          PROTO((Catalog *catalog, off_t meas, int cat));
     202void          matchMosaics        PROTO((Catalog *catalog, off_t meas, int cat));
     203GSCRegion    *name_region         PROTO((char *name, off_t *Nregions));
    204204double        opening_angle       PROTO((double x1, double y1, double x2, double y2, double x3, double y3));
    205205void          plot_chisq          PROTO((Catalog *catalog, int Ncatalog));
    … …  
    224224int           setMrelOutput       PROTO((Catalog *catalog, int Ncatalog, int mark));
    225225void          set_ZP              PROTO((double ZERO));
    226 int           setrefcode          PROTO((Image *image, int Nimage));
     226int           setrefcode          PROTO((Image *image, off_t Nimage));
    227227void          skip_measurements   PROTO((Catalog *catalog, int pass));
    228228void          sortA               PROTO((double *X, int N));
    … …  
    268268int UpdateMosaic (Catalog *catalog, int Ncatalog);
    269269int UpdateMeasures (Catalog *catalog, int Ncatalog);
    270 void fixImageRaw (Catalog *catalog, int Ncatalog, int im);
     270void fixImageRaw (Catalog *catalog, int Ncatalog, off_t im);
    271271void FlagOutliers(Catalog *catalog);
    272272int MeasFilterTest(Measure *measure);
    … …  
    278278int FitPMandPar (PMFit *fit, double *X, double *dX, double *Y, double *dY, double *T, double *pR, double *pD, int Npts);
    279279
    280 Mosaic *getMosaicForImage (int N);
    281 
    282 StarData *getImageRef (Catalog *catalog, int Ncatalog, int im, int *Nstars, CoordMode mode);
    283 StarData *getImageRaw (Catalog *catalog, int Ncatalog, int im, int *Nstars, CoordMode mode);
    284 
    285 Mosaic *getmosaics (int *N);
    286 void initMosaics (Image *image, int Nimage);
    287 StarData *getMosaicRaw (Catalog *catalog, int Ncatalog, int mos, int *Nstars);
    288 StarData *getMosaicRef (Catalog *catalog, int Ncatalog, int mos, int *Nstars);
    289 Mosaic *getMosaicForImage (int im);
     280Mosaic *getMosaicForImage (off_t N);
     281
     282StarData *getImageRef (Catalog *catalog, int Ncatalog, off_t im, off_t *Nstars, CoordMode mode);
     283StarData *getImageRaw (Catalog *catalog, int Ncatalog, off_t im, off_t *Nstars, CoordMode mode);
     284
     285Mosaic *getmosaics (off_t *N);
     286void initMosaics (Image *image, off_t Nimage);
     287StarData *getMosaicRaw (Catalog *catalog, int Ncatalog, off_t mos, off_t *Nstars);
     288StarData *getMosaicRef (Catalog *catalog, int Ncatalog, off_t mos, off_t *Nstars);
     289Mosaic *getMosaicForImage (off_t im);
    290290
    291291double getMeanR (Measure *measure, Average *average, SecFilt *secfilt);
    … …  
    297297int relastro_objects (void);
    298298int UpdateObjectOffsets (SkyList *skylist);
     299
     300int relastroVisualPlotRawRef(StarData *raw, StarData *ref, double dRmax, int numObj);
     301int relastroVisualPlotScatter(double values[], double thresh, int npts);
     302int relastroVisualPlotOutliers(Catalog *catalog, int offset, int Nmeasure,
     303                               StatType statsR, StatType statsD, double thresh);
     304

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/src/ImageOps.c

    r24308 r27312  
    44static unsigned int *start;
    55static unsigned int *stop;
    6 static int         **bin;
    7 
    8 static int         **clist;
    9 static int         **mlist;
    10 static int          *Nlist;
    11 static int          *NLIST;
    12 
    13 static Image        *image;
    14 static int          Nimage;
    15 
    16 Image *getimages (int *N) {
     6
     7static off_t       **bin;     // link from catalog,measure to image
     8static int         **clist;   // catalog which supplied measurement on image
     9static off_t       **mlist;   // measure reference for measurement on image
     10static off_t        *Nlist;   // number of measurements on image
     11static off_t        *NLIST;   // allocated number of measurements on image
     12
     13static Image        *image;   // list of available images
     14static off_t        Nimage;   // number of available images
     15
     16Image *getimages (off_t *N) {
    1717  *N = Nimage;
    1818  return (image);
    1919}
    2020
    21 Image *getimage (int N) {
     21Image *getimage (off_t N) {
    2222  return (&image[N]);
    2323}
    2424
    25 void initImages (Image *input, int N) {
    26 
    27   int i;
     25void initImages (Image *input, off_t N) {
     26
     27  off_t i;
    2828
    2929  image = input;
    … …  
    4141void initImageBins (Catalog *catalog, int Ncatalog) {
    4242
    43   int i, j;
    44 
    45   ALLOCATE (bin, int *, Ncatalog);
     43  off_t i, j;
     44
     45  ALLOCATE (bin, off_t *, Ncatalog);
    4646  for (i = 0; i < Ncatalog; i++) {
    47     ALLOCATE (bin[i], int, MAX (catalog[i].Nmeasure, 1));
     47    ALLOCATE (bin[i], off_t, MAX (catalog[i].Nmeasure, 1));
    4848    for (j = 0; j < catalog[i].Nmeasure; j++) bin[i][j] = -1;
    4949  }
    5050
    51   ALLOCATE (Nlist, int, Nimage);
    52   ALLOCATE (NLIST, int, Nimage);
    53   ALLOCATE (clist, int *, Nimage);
    54   ALLOCATE (mlist, int *, Nimage);
     51  ALLOCATE (Nlist, off_t,  Nimage);
     52  ALLOCATE (NLIST, off_t,  Nimage);
     53  ALLOCATE (clist, int *,   Nimage);
     54  ALLOCATE (mlist, off_t *, Nimage);
    5555
    5656  for (i = 0; i < Nimage; i++) {
    … …  
    5858    NLIST[i] = 100;
    5959    ALLOCATE (clist[i], int, NLIST[i]);
    60     ALLOCATE (mlist[i], int, NLIST[i]);
     60    ALLOCATE (mlist[i], off_t, NLIST[i]);
    6161  }
    6262}
    … …  
    6464void freeImageBins (int Ncatalog) {
    6565
    66   int i;
     66  off_t i;
    6767
    6868  for (i = 0; i < Ncatalog; i++) {
    … …  
    8181void findImages (Catalog *catalog, int Ncatalog) {
    8282
    83   int i, j;
     83  off_t i, j;
    8484  char *name;
    8585
    … …  
    9292  for (i = 0; VERBOSE && (i < Nimage); i++) {
    9393    name = GetPhotcodeNamebyCode (image[i].photcode);
    94     fprintf (stderr, "image %d has %d measures (%s, %s)\n", i, Nlist[i],
     94    fprintf (stderr, "image %lld has %lld measures (%s, %s)\n", (long long) i, (long long) Nlist[i],
    9595             ohana_sec_to_date(image[i].tzero), name);
    9696  }
    … …  
    9898
    9999/* modify this function to use the measure->imageID field */
    100 void matchImage (Catalog *catalog, int meas, int cat) {
    101 
    102   int i;
     100void matchImage (Catalog *catalog, off_t meas, int cat) {
     101
     102  off_t i;
    103103  Measure *measure;
    104104
    … …  
    124124    if (Nlist[i] == NLIST[i]) {
    125125      NLIST[i] += 100;
    126       REALLOCATE (clist[i], int, NLIST[i]);
    127       REALLOCATE (mlist[i], int, NLIST[i]);
     126      REALLOCATE (clist[i], int,   NLIST[i]);
     127      REALLOCATE (mlist[i], off_t, NLIST[i]);
    128128    }
    129129    return;
    … …  
    131131}
    132132
    133 Coords *getCoords (int meas, int cat) {
    134 
    135   int i;
     133Coords *getCoords (off_t meas, int cat) {
     134
     135  off_t i;
    136136
    137137  i = bin[cat][meas];
    … …  
    142142void plot_images () {
    143143
    144   int i, bin;
     144  off_t i, bin;
    145145  double *xlist, *Mlist, *dlist;
    146146  Graphdata graphdata;
    … …  
    188188// return StarData values for detections in the specified image, converting coordinates from the
    189189// chip positions: X,Y -> L,M -> P,Q -> R,D
    190 void fixImageRaw (Catalog *catalog, int Ncatalog, int im) {
    191 
    192   int i, m, c, n;
     190void fixImageRaw (Catalog *catalog, int Ncatalog, off_t im) {
     191
     192  off_t i, m, c, n;
    193193  double X, Y, L, M, P, Q, R, D, dR, dD;
    194194
    … …  
    257257// chip positions: X,Y -> L,M -> P,Q -> R,D.  This function is used by the image fitting steps, for
    258258// which the detections have already been filtered when they were loaded (bcatalog)
    259 StarData *getImageRaw (Catalog *catalog, int Ncatalog, int im, int *Nstars, CoordMode mode) {
    260 
    261   int i, m, c, n;
     259StarData *getImageRaw (Catalog *catalog, int Ncatalog, off_t im, off_t *Nstars, CoordMode mode) {
     260
     261  off_t i, m, c, n;
    262262
    263263  Mosaic *mosaic;
    … …  
    332332// the sky positions: R,D -> P,Q -> L,M -> X,Y
    333333
    334 StarData *getImageRef (Catalog *catalog, int Ncatalog, int im, int *Nstars, CoordMode mode) {
    335 
    336   int i, m, c, n;
     334StarData *getImageRef (Catalog *catalog, int Ncatalog, off_t im, off_t *Nstars, CoordMode mode) {
     335
     336  off_t i, m, c, n;
    337337
    338338  Mosaic *mosaic;
    … …  
    399399  return;
    400400
    401   int i, j, k, m, N, Ndel, Nave, Nmax, TOOFEW, Nsecfilt;
     401  int Ndel, Nave;
     402  off_t i, j, k, m, N, Nmax, TOOFEW, Nsecfilt;
    402403  double Ns, theta, x, y;
    403404  double *R, *D, *dR, *dD;
    … …  
    490491void FlagOutliers2D (Catalog *catalog) {
    491492
    492   int i, j, k, m, N, Ndel, Nave, Nmax, TOOFEW, Nsecfilt;
     493  int Ndel, Nave;
     494  off_t i, j, k, m, N, Nmax, TOOFEW, Nsecfilt;
    493495  double *index;
    494496  double Ns, theta, x, y;
    … …  
    569571    // recalculate image center, sigma based on closest 50% of points
    570572    for(k = 0;  k < N; k++) {
    571       int ind = (int) index[k];
     573      off_t ind = (off_t) index[k];
    572574      R[k] = catalog[0].measure[ind].dR;
    573575      D[k] = catalog[0].measure[ind].dD;

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/src/MosaicOps.c

    r21508 r27312  
    22
    33// array of mosaic definition structures
    4 static int    Nmosaic;
     4static off_t    Nmosaic;
    55static Mosaic *mosaic;
    66
    77// list of all images associated with a mosaic
    8 static int   *Nmosaic_own_images; // number of images for this mosaic
    9 static int   *Amosaic_own_images; // size of allocated array
    10 static int   **mosaic_own_images; // array of arrays: mosaic -> images
     8static off_t   *Nmosaic_own_images; // number of images for this mosaic
     9static off_t   *Amosaic_own_images; // size of allocated array
     10static off_t   **mosaic_own_images; // array of arrays: mosaic -> images
    1111
    1212// list of mosaic associated with each image 
    13 static int    Nmosaic_for_images; // number of images (for internal checks)
    14 static int    *mosaic_for_images; // array of: image -> mosaic
     13static off_t    Nmosaic_for_images; // number of images (for off_ternal checks)
     14static off_t    *mosaic_for_images; // array of: image -> mosaic
    1515
    16 Mosaic *getmosaics (int *N) {
     16Mosaic *getmosaics (off_t *N) {
    1717  *N = Nmosaic;
    1818  return (mosaic);
    … …  
    2020
    2121// find mosaic frames (unique time periods & photcode name matches mosaic)
    22 void initMosaics (Image *image, int Nimage) {
     22void initMosaics (Image *image, off_t Nimage) {
    2323
    24   int i, j, found, NMOSAIC;
     24  off_t i, j, found, NMOSAIC;
    2525  unsigned int start, stop;
    2626
    … …  
    2929  ALLOCATE (mosaic, Mosaic, NMOSAIC);
    3030
    31   ALLOCATE (Nmosaic_own_images, int, NMOSAIC);
    32   ALLOCATE (Amosaic_own_images, int, NMOSAIC);
    33   ALLOCATE (mosaic_own_images, int *, NMOSAIC);
     31  ALLOCATE (Nmosaic_own_images, off_t, NMOSAIC);
     32  ALLOCATE (Amosaic_own_images, off_t, NMOSAIC);
     33  ALLOCATE (mosaic_own_images, off_t *, NMOSAIC);
    3434
    3535  /* find the mosaic images (coords.ctype = DIS); generate list of unique mosaics */
    … …  
    5454    Nmosaic_own_images[Nmosaic] = 0;
    5555    Amosaic_own_images[Nmosaic] = 10;
    56     ALLOCATE (mosaic_own_images[Nmosaic], int, Amosaic_own_images[Nmosaic]);
     56    ALLOCATE (mosaic_own_images[Nmosaic], off_t, Amosaic_own_images[Nmosaic]);
    5757
    5858    Nmosaic ++;
    … …  
    6060      NMOSAIC += 10;
    6161      REALLOCATE (mosaic, Mosaic, NMOSAIC);
    62       REALLOCATE (mosaic_own_images, int *, NMOSAIC);
    63       REALLOCATE (Nmosaic_own_images, int, NMOSAIC);
    64       REALLOCATE (Amosaic_own_images, int, NMOSAIC);
     62      REALLOCATE (mosaic_own_images, off_t *, NMOSAIC);
     63      REALLOCATE (Nmosaic_own_images, off_t, NMOSAIC);
     64      REALLOCATE (Amosaic_own_images, off_t, NMOSAIC);
    6565    }
    6666  }
    … …  
    6868  // array to store image->mosaic index
    6969  Nmosaic_for_images = Nimage;
    70   ALLOCATE (mosaic_for_images, int, Nmosaic_for_images);
     70  ALLOCATE (mosaic_for_images, off_t, Nmosaic_for_images);
    7171
    7272  /* now assign the WRP images to these mosaics */
    … …  
    9999    if (Nmosaic_own_images[j] == Amosaic_own_images[j]) {
    100100      Amosaic_own_images[j] += 10;
    101       REALLOCATE (mosaic_own_images[j], int, Amosaic_own_images[j]);
     101      REALLOCATE (mosaic_own_images[j], off_t, Amosaic_own_images[j]);
    102102    }
    103103    assert (Nmosaic_own_images[j] < Amosaic_own_images[j]);
    … …  
    109109// return StarData values for detections in the specified image, converting coordinates from the
    110110// chip positions: X,Y -> L,M -> P,Q -> R,D
    111 StarData *getMosaicRaw (Catalog *catalog, int Ncatalog, int mos, int *Nstars) {
     111StarData *getMosaicRaw (Catalog *catalog, int Ncatalog, off_t mos, off_t *Nstars) {
    112112
    113   int i, j, im, Nraw, Nnew;
     113  off_t i, j, im, Nraw, Nnew;
    114114  StarData *raw, *new;
    115115
    … …  
    142142// return StarData values for averages positions in the specified image, converting coordinates from
    143143// the sky positions: R,D -> P,Q -> L,M -> X,Y
    144 StarData *getMosaicRef (Catalog *catalog, int Ncatalog, int mos, int *Nstars) {
     144StarData *getMosaicRef (Catalog *catalog, int Ncatalog, off_t mos, off_t *Nstars) {
    145145
    146   int i, j, im, Nref, Nnew;
     146  off_t i, j, im, Nref, Nnew;
    147147  StarData *ref, *new;
    148148 
    … …  
    172172}
    173173
    174 Mosaic *getMosaicForImage (int im) {
     174Mosaic *getMosaicForImage (off_t im) {
    175175
    176   int mos;
     176  off_t mos;
    177177
    178178  if (im < 0) abort();

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/src/UpdateChips.c

    r17151 r27312  
    44
    55  /* we can measure new image parameters for each non-mosaic chip independently */
    6   int i, Nimage, Nraw, Nref;
     6  off_t i, Nimage, Nraw, Nref;
    77  Image *image;
    88  StarData *raw, *ref;
    … …  
    2525
    2626    // FitChip does iterative, clipped fitting
    27     fprintf (stderr, "image %d : Nstars: %d\n", i, Nraw);
     27    fprintf (stderr, "image %lld : Nstars: %lld\n", (long long) i, (long long) Nraw);
    2828    FitChip (raw, ref, Nraw, &image[i].coords);
    2929

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/src/UpdateMeasures.c

    r16810 r27312  
    33int UpdateMeasures (Catalog *catalog, int Ncatalog) {
    44
    5   int i, Nimage;
     5  off_t i, Nimage;
    66  Image *image;
    77

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/src/UpdateMosaic.c

    r15600 r27312  
    44
    55  /* we can measure new image parameters for each mosaic independently */
    6   int i, Nmosaic, Nstars;
     6  off_t i, Nmosaic, Nstars;
    77  Mosaic *mosaic;
    88  StarData *raw, *ref;

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/src/UpdateObjects.c

    r24308 r27312  
    11# include "relastro.h"
    22
    3 static int Nmax;
     3static off_t  Nmax;
    44static double *X, *dX;
    55static double *Y, *dY;
    … …  
    1313void initObjectData (Catalog *catalog, int Ncatalog) {
    1414
    15   int i, j;
     15  off_t i, j;
    1616 
    1717  Nmax = 0;
    … …  
    4040int UpdateObjects (Catalog *catalog, int Ncatalog) {
    4141
    42   int i, j, k, m, N, Nsecfilt, found, kp;
     42  int i, j, k, m, N, Nsecfilt;
    4343  StatType statsR, statsD;
    4444  Coords coords;
    … …  
    4747  int mode, Nave, Npm, Npar, Nskip;
    4848  double Tmin, Tmax;
    49   float mag;
    50   int mask;
    51   PhotCode *code;
    5249
    5350  initObjectData (catalog, Ncatalog);
    … …  
    7471  for (i = 0; i < Ncatalog; i++) {
    7572
    76     if (VERBOSE) fprintf (stderr, "astrometrize catalog %d : %d ave, %d meas\n", i, catalog[i].Naverage, catalog[i].Nmeasure);
     73    if (VERBOSE) fprintf (stderr, "astrometrize catalog %d : %lld ave, %lld meas\n", i, (long long) catalog[i].Naverage, (long long) catalog[i].Nmeasure);
    7774
    7875    Nskip = 0;

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/src/bcatalog.c

    r24308 r27312  
    33int bcatalog (Catalog *subcatalog, Catalog *catalog) {
    44
    5   int i, j, k, offset, found;
    6   int NAVERAGE, NMEASURE, Naverage, Nmeasure, Nm, Nsecfilt;
    7   float mag;
    8   int mask;
    9   PhotCode *code;
     5  off_t i, j, offset;
     6  off_t NAVERAGE, NMEASURE, Naverage, Nmeasure, Nm;
     7  int Nsecfilt;
    108
    119  // XXX in the future, use catalog[0].Nsecfilt only?  allow catalogs to have variable Nsecfilt?
    … …  
    106104
    107105  if (VERBOSE) {
    108     fprintf (stderr, "%d: using %d stars (%d measures) for catalog\n", i,
    109              subcatalog[0].Naverage, subcatalog[0].Nmeasure);
     106    fprintf (stderr, "%lld: using %lld stars (%lld measures) for catalog\n", (long long) i, (long long) subcatalog[0].Naverage, (long long) subcatalog[0].Nmeasure);
    110107   }
    111108  return (TRUE);

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/src/load_images.c

    r25757 r27312  
    44
    55  Image     *image, *subset;
    6   int        Nimage, Nsubset;
    7   int       *LineNumber;
     6  off_t      Nimage, Nsubset;
     7  off_t     *LineNumber;
    88
    99  SkyTable *sky = NULL;

  • branches/eam_branches/largefiles.20100314/Ohana/src/relastro/src/select_images.c

    r21508 r27312  
    1313} SkyRegionCoords;
    1414
    15 Image *select_images (SkyList *skylist, Image *timage, int Ntimage, int **LineNumber, int *Nimage) {
     15Image *select_images (SkyList *skylist, Image *timage, off_t Ntimage, off_t **LineNumber, off_t *Nimage) {
    1616 
    1717  Image *image;
    18   int i, j, k, m, found, nimage, NIMAGE;
    19   int InRange;
     18  off_t i, j, k, m, nimage, NIMAGE;
     19  int InRange, found;
    2020  double Ri[5], Di[5], Xi[5], Yi[5], dx, dy;
    21   int *line_number;
     21  off_t *line_number;
    2222  Coords tcoords;
    2323  SkyRegionCoords *skycoords;
    … …  
    7171  NIMAGE = 100;
    7272  ALLOCATE (image, Image, NIMAGE);
    73   ALLOCATE (line_number, int, NIMAGE);
     73  ALLOCATE (line_number, off_t, NIMAGE);
    7474 
    7575  // go through the complete list of images, selecting ones which overlap any region
    … …  
    161161        NIMAGE += 100;
    162162        REALLOCATE (image, Image, NIMAGE);
    163         REALLOCATE (line_number, int, NIMAGE);
     163        REALLOCATE (line_number, off_t, NIMAGE);
    164164      }
    165165    }
    166166  }
    167167     
    168   if (VERBOSE) fprintf (stderr, "found %d images\n", nimage);
     168  if (VERBOSE) fprintf (stderr, "found %lld images\n", (long long) nimage);
    169169
    170170  REALLOCATE (image, Image, MAX (nimage, 1));
    171   REALLOCATE (line_number, int, MAX (nimage, 1));
     171  REALLOCATE (line_number, off_t, MAX (nimage, 1));
    172172  free (skycoords);
    173173
Note: See TracChangeset for help on using the changeset viewer.