IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 38463


Ignore:
Timestamp:
Jun 12, 2015, 10:26:46 PM (11 years ago)
Author:
eugene
Message:

merge changes from EAM branch ipp-20150419

Location:
trunk/Ohana/src/opihi
Files:
5 edited
2 copied

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/opihi/cmd.astro/cdensify.c

    r38062 r38463  
    7373  CHECKVAL(dY);
    7474
    75   Nx = (Xmax - Xmin) / dX + 1;
    76   Ny = (Ymax - Ymin) / dY + 1;
     75  Nx = abs((Xmax - Xmin) / dX) + 1;
     76  Ny = abs((Ymax - Ymin) / dY) + 1;
    7777 
    7878  Coords newcoords = graphmode.coords;

  • trunk/Ohana/src/opihi/cmd.data/Makefile

    r38234 r38463  
    2525$(SRC)/book.$(ARCH).o           \
    2626$(SRC)/book_commands.$(ARCH).o  \
     27$(SRC)/bisection.$(ARCH).o              \
    2728$(SRC)/cast.$(ARCH).o           \
    2829$(SRC)/center.$(ARCH).o \

  • trunk/Ohana/src/opihi/cmd.data/init.c

    r38234 r38463  
    99int box              PROTO((int, char **));
    1010int book_command     PROTO((int, char **));
     11int bisection        PROTO((int, char **));
    1112int center           PROTO((int, char **));
    12 int parity           PROTO((int, char **));
    1313int cast             PROTO((int, char **));
    1414int circstats        PROTO((int, char **));
    … …  
    9090int plot             PROTO((int, char **));
    9191int dot              PROTO((int, char **));
     92int parity           PROTO((int, char **));
    9293int point            PROTO((int, char **));
    9394int ps               PROTO((int, char **));
    … …  
    176177  {1, "applyfit2d",   applyfit2d,       "apply 2-d fit to new vector"},
    177178  {1, "applyfit3d",   applyfit3d,       "apply 3-d fit to new vector"},
     179  {1, "bisection",    bisection,        "use bisection to find threshold in vector"},
    178180  {1, "book",         book_command,     "commands to manipulate book/page/word data"},
    179181  {1, "box",          box,              "draw a box on the plot"},

  • trunk/Ohana/src/opihi/cmd.data/tvcolors.c

    r31667 r38463  
    1616  FREE (name);
    1717
     18  int SetNAN = FALSE;
     19  int red, green, blue;
     20  if ((N = get_argument (argc, argv, "-nan"))) {
     21    if (N > argc - 4) {
     22      gprint (GP_ERR, "ERROR: -nan (red) (green) (blue) does not have enough arguments\n");
     23      return FALSE;
     24    }
     25    remove_argument (N, &argc, argv);
     26    red = strtol (argv[N], NULL, 0);
     27    remove_argument (N, &argc, argv);
     28    green = strtol (argv[N], NULL, 0);
     29    remove_argument (N, &argc, argv);
     30    blue = strtol (argv[N], NULL, 0);
     31    remove_argument (N, &argc, argv);
     32    SetNAN = TRUE;
     33  }
     34
    1835  // use the currently-set zero,range values if not supplied
    19   if (argc != 2) {
    20     gprint (GP_ERR, "USAGE: tvcolors (colormap)\n");
     36  if (!SetNAN && (argc != 2)) {
     37    gprint (GP_ERR, "USAGE: tvcolors (colormap) [-nan red green blue]\n");
    2138    gprint (GP_ERR, " colormap options : greyscale, -greyscale, rainbow, heat, fullcolor, ruffcolor (also grayscale, -grayscale)\n");
    2239    return (FALSE);
    2340  }
    2441
    25   KiiSetColormap (kapa, argv[1]);
     42  if (SetNAN) {
     43    KiiSetNanColor (kapa, red, green, blue);
     44  }
     45  if (argc == 2) {
     46    KiiSetColormap (kapa, argv[1]);
     47  }
    2648  return (TRUE);
    2749}

  • trunk/Ohana/src/opihi/dvo/gstar.c

    r38441 r38463  
    88
    99typedef enum {
     10  GSTAR_FILTER_NAMES,
     11  GSTAR_SECF_FLAGS,
     12  GSTAR_UCDIST,
    1013  GSTAR_AVE_MAG,
    1114  GSTAR_AVE_MAG_ERR,
    1215  GSTAR_AVE_MAG_CHISQ,
    13   GSTAR_FILTER_NAMES,
    14   GSTAR_SECF_FLAGS,
    15   GSTAR_UCDIST,
    1616  GSTAR_AVE_AP_MAG,
     17  GSTAR_AVE_AP_MAG_ERR,
     18  GSTAR_AVE_AP_MAG_CHISQ,
    1719  GSTAR_AVE_MAG_MIN,
    1820  GSTAR_AVE_MAG_MAX,
    1921  GSTAR_AVE_KRON_MAG,
    2022  GSTAR_AVE_KRON_MAG_ERR,
     23  GSTAR_AVE_KRON_MAG_CHISQ,
     24  GSTAR_AVE_NCODE,
     25  GSTAR_AVE_NUSED,
     26  GSTAR_AVE_AP_NUSED,
     27  GSTAR_AVE_KRON_NUSED,
    2128  GSTAR_STACK_FLUX_PSF,
    2229  GSTAR_STACK_FLUX_PSF_ERR,
    … …  
    2734  GSTAR_WARP_FLUX_PSF,
    2835  GSTAR_WARP_FLUX_PSF_ERR,
     36  GSTAR_WARP_FLUX_PSF_CHISQ,
    2937  GSTAR_WARP_FLUX_KRON,
    3038  GSTAR_WARP_FLUX_KRON_ERR,
     39  GSTAR_WARP_FLUX_KRON_CHISQ,
    3140  GSTAR_WARP_FLUX_APER,
    3241  GSTAR_WARP_FLUX_APER_ERR,
     42  GSTAR_WARP_FLUX_APER_CHISQ,
    3343  GSTAR_STACK_MAG_PSF,
    3444  GSTAR_STACK_MAG_PSF_ERR,
    … …  
    3747  GSTAR_STACK_MAG_APER,
    3848  GSTAR_STACK_MAG_APER_ERR,
     49  GSTAR_STACK_NUSED,
    3950  GSTAR_WARP_MAG_PSF,
    4051  GSTAR_WARP_MAG_PSF_ERR,
     52  GSTAR_WARP_MAG_PSF_CHISQ,
    4153  GSTAR_WARP_MAG_KRON,
    4254  GSTAR_WARP_MAG_KRON_ERR,
     55  GSTAR_WARP_MAG_KRON_CHISQ,
    4356  GSTAR_WARP_MAG_APER,
    4457  GSTAR_WARP_MAG_APER_ERR,
     58  GSTAR_WARP_MAG_APER_CHISQ,
     59  GSTAR_WARP_NWARP,
     60  GSTAR_WARP_NWARP_GOOD,
     61  GSTAR_WARP_NWARP_USED_PSF,
     62  GSTAR_WARP_NWARP_USED_KRON,
     63  GSTAR_WARP_NWARP_USED_AP,
    4564} GSTAR_SECF_CODES;
    4665
    … …  
    310329        gprint (GP_LOG, "%3d   ",  catalog.average[k].Nmeasure);
    311330        gprint (GP_LOG, "%4.1f ",  catalog.average[k].ChiSqAve);
    312         gprint (GP_LOG, "%5x ",    catalog.average[k].flags);
    313         gprint (GP_LOG, "%x ",    catalog.average[k].objID);
    314         gprint (GP_LOG, "%x ",    catalog.average[k].catID);
     331        gprint (GP_LOG, "0x%08x ", catalog.average[k].flags);
     332        gprint (GP_LOG, "0x%08x ", catalog.average[k].objID);
     333        gprint (GP_LOG, "0x%08x ", catalog.average[k].catID);
    315334       
    316335        if (FULL_OUTPUT) {
    … …  
    344363        gprint (GP_LOG, "\n");
    345364
     365        if (EXTRA_OUTPUT) {
     366          /* Mmin */
     367          gprint (GP_LOG, "chp_psf_min: ");
     368          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_MAG_MIN);
     369          gprint (GP_LOG, "\n");
     370
     371          /* Mmax */
     372          gprint (GP_LOG, "chp_psf_max: ");
     373          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_MAG_MAX);;
     374          gprint (GP_LOG, "\n");
     375        }
     376
    346377        /* average mag errors */
    347378        gprint (GP_LOG, "chp_psf_err: ");
    … …  
    354385        gprint (GP_LOG, "\n");
    355386
     387        /* average mag Ncode */
     388        gprint (GP_LOG, "chp_Ncode:   ");
     389        for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_NCODE);
     390        gprint (GP_LOG, "\n");
     391       
     392        /* average mag Nused */
     393        gprint (GP_LOG, "chp_Nused:   ");
     394        for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_NUSED);
     395        gprint (GP_LOG, "\n");
    356396       
    357397        if (APER_OUTPUT) {
    … …  
    360400          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_AP_MAG);
    361401          gprint (GP_LOG, "\n");
     402
     403          gprint (GP_LOG, "chp_ap__err: ");
     404          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_AP_MAG_ERR);
     405          gprint (GP_LOG, "\n");
     406
     407          gprint (GP_LOG, "chp_ap__chi: ");
     408          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_AP_MAG_CHISQ);
     409          gprint (GP_LOG, "\n");
     410
     411          /* average mag Nused */
     412          gprint (GP_LOG, "chp_ap__N:   ");
     413          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_AP_NUSED);
     414          gprint (GP_LOG, "\n");
    362415        }
    363416
    … …  
    372425          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_KRON_MAG_ERR);
    373426          gprint (GP_LOG, "\n");
     427
     428          gprint (GP_LOG, "chp_krn_chi: ");
     429          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_KRON_MAG_CHISQ);
     430          gprint (GP_LOG, "\n");
     431
     432          /* average mag Nused */
     433          gprint (GP_LOG, "chp_krn_N:   ");
     434          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_KRON_NUSED);
     435          gprint (GP_LOG, "\n");
    374436        }
    375437
    376         if (EXTRA_OUTPUT) {
    377           /* Mmin */
    378           gprint (GP_LOG, "chp_psf_min: ");
    379           for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_MAG_MIN);
    380           gprint (GP_LOG, "\n");
    381 
    382           /* Mmax */
    383           gprint (GP_LOG, "chp_psf_max: ");
    384           for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_MAG_MAX);;
    385           gprint (GP_LOG, "\n");
    386 
    387           /* UCDIST */
    388           gprint (GP_LOG, "ubercal_dst: ");
    389           for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_UCDIST);;
    390           gprint (GP_LOG, "\n");
    391 
    392           /* secfilt flags */
    393           gprint (GP_LOG, "filtflags:   ");
    394           for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_SECF_FLAGS);
    395           gprint (GP_LOG, "\n");
    396         }
    397438        if (STACK_OUTPUT && MAG_OUTPUT) {
    398439          /* FluxPSF */
    … …  
    406447          gprint (GP_LOG, "\n");
    407448 
     449          if (!FLUX_OUTPUT) {
     450            gprint (GP_LOG, "stk_Nstack: ");
     451            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_STACK_NUSED);
     452            gprint (GP_LOG, "\n");
     453          }
     454
    408455          if (KRON_OUTPUT) {
    409456            /* MagKron */
    … …  
    441488          gprint (GP_LOG, "\n");
    442489 
     490          gprint (GP_LOG, "wrp_psf_chi: ");
     491          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_MAG_PSF_CHISQ);
     492          gprint (GP_LOG, "\n");
     493 
     494          if (!FLUX_OUTPUT) {
     495            gprint (GP_LOG, "wrp_Nwarp:   ");
     496            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP);
     497            gprint (GP_LOG, "\n");
     498           
     499            gprint (GP_LOG, "wrp_NGood:   ");
     500            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP_GOOD);
     501            gprint (GP_LOG, "\n");
     502
     503            gprint (GP_LOG, "wrp_NUsed:   ");
     504            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP_USED_PSF);
     505            gprint (GP_LOG, "\n");
     506          }
     507         
    443508          if (KRON_OUTPUT) {
    444509            /* MagKron */
    … …  
    452517            gprint (GP_LOG, "\n");
    453518
     519            /* dMagKron */
     520            gprint (GP_LOG, "wrp_krn_chp: ");
     521            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_MAG_KRON_CHISQ);
     522            gprint (GP_LOG, "\n");
     523
     524            if (!FLUX_OUTPUT) {
     525              gprint (GP_LOG, "wrp_krn_Nuse:");
     526              for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP_USED_KRON);
     527              gprint (GP_LOG, "\n");
     528            }
    454529          }
    455 
     530         
    456531          if (APER_OUTPUT) {
    457532            /* MagAper */
    … …  
    464539            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_MAG_APER_ERR);
    465540            gprint (GP_LOG, "\n");
     541
     542            /* dMagAper */
     543            gprint (GP_LOG, "wrp_ap__chi: ");
     544            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_MAG_APER_CHISQ);
     545            gprint (GP_LOG, "\n");
     546
     547            if (!FLUX_OUTPUT) {
     548              gprint (GP_LOG, "wrp_ap_Nuse: ");
     549              for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP_USED_AP);
     550              gprint (GP_LOG, "\n");
     551            }
    466552          }
    467553        }
    … …  
    491577          if (APER_OUTPUT) {
    492578            /* FluxAper */
    493             gprint (GP_LOG, "stk_ap_ave: ");
     579            gprint (GP_LOG, "stk_ap_ave:  ");
    494580            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_STACK_FLUX_APER);
    495581            gprint (GP_LOG, "\n");
    496582 
    497583            /* dFluxAper */
    498             gprint (GP_LOG, "stk_ap_err: ");
     584            gprint (GP_LOG, "stk_ap_err:  ");
    499585            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_STACK_FLUX_APER_ERR);
    500586            gprint (GP_LOG, "\n");
    501587          }
     588
     589          gprint (GP_LOG, "stk_Nstack:  ");
     590          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_STACK_NUSED);
     591          gprint (GP_LOG, "\n");
    502592        }
    503593        if (WARP_OUTPUT && FLUX_OUTPUT) {
    … …  
    508598 
    509599          /* dFluxPSF */
    510           gprint (GP_LOG, "wrp_psf_ave: ");
     600          gprint (GP_LOG, "wrp_psf_err: ");
    511601          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_FLUX_PSF_ERR);
    512602          gprint (GP_LOG, "\n");
    513603 
     604          gprint (GP_LOG, "wrp_psf_chi: ");
     605          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_FLUX_PSF_CHISQ);
     606          gprint (GP_LOG, "\n");
     607 
     608          gprint (GP_LOG, "wrp_Nwarp:   ");
     609          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP);
     610          gprint (GP_LOG, "\n");
     611         
     612          gprint (GP_LOG, "wrp_NGood:   ");
     613          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP_GOOD);
     614          gprint (GP_LOG, "\n");
     615
     616          gprint (GP_LOG, "wrp_NUsed:   ");
     617          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP_USED_PSF);
     618          gprint (GP_LOG, "\n");
     619
    514620          if (KRON_OUTPUT) {
    515621            /* FluxKron */
    … …  
    522628            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_FLUX_KRON_ERR);
    523629            gprint (GP_LOG, "\n");
     630
     631            gprint (GP_LOG, "wrp_krn_chi: ");
     632            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_FLUX_KRON_CHISQ);
     633            gprint (GP_LOG, "\n");
     634
     635            gprint (GP_LOG, "wrp_krn_Nuse:");
     636            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP_USED_KRON);
     637            gprint (GP_LOG, "\n");
    524638          }
    525639
    … …  
    534648            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_FLUX_APER_ERR);
    535649            gprint (GP_LOG, "\n");
     650
     651            gprint (GP_LOG, "wrp_ap__chi: ");
     652            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_FLUX_APER_CHISQ);
     653            gprint (GP_LOG, "\n");
     654
     655            gprint (GP_LOG, "wrp_ap_Nuse: ");
     656            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_NWARP_USED_AP);
     657            gprint (GP_LOG, "\n");
    536658          }
    537659        }
    538       }
     660        if (EXTRA_OUTPUT) {
     661          /* secfilt flags */
     662          gprint (GP_LOG, "filtflags:   ");
     663          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_SECF_FLAGS);
     664          gprint (GP_LOG, "\n");
     665
     666          /* UCDIST */
     667          gprint (GP_LOG, "ubercal_dst: ");
     668          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_UCDIST);;
     669          gprint (GP_LOG, "\n");
     670        }
     671      }
     672
    539673      gprint (GP_LOG, "\n");
    540674     
    … …  
    626760              gprint (GP_LOG, "0x%08x ", catalog.measure[Nv].imageID);
    627761              gprint (GP_LOG, "%5.3f ", catalog.measure[Nv].psfQF);
     762              gprint (GP_LOG, "%5.3f ", catalog.measure[Nv].psfQFperf);
    628763              gprint (GP_LOG, "%7.1f ", catalog.measure[Nv].psfChisq);
    629764              // gprint (GP_LOG, "%3.1f ", catalog.measure[Nv].crNsigma);
    … …  
    636771              gprint (GP_LOG, "%10.3e ", catalog.measure[Nv].FluxKron);
    637772              gprint (GP_LOG, "%10.3e ", catalog.measure[Nv].dFluxKron);
     773              gprint (GP_LOG, "%10.3e ", catalog.measure[Nv].FluxAp);
     774              gprint (GP_LOG, "%10.3e ", catalog.measure[Nv].dFluxAp);
    638775            }
    639776            gprint (GP_LOG, "\n");
    … …  
    688825void print_double (double value) {
    689826  if (isnan(value))
    690     gprint (GP_LOG, "  NaN    ");
     827    gprint (GP_LOG, "       NaN ");
    691828  else
    692     gprint (GP_LOG, "  %6.3f ", value);
     829    gprint (GP_LOG, "   %7.4f ", value);
    693830}
    694831
    695832void print_double_exp (double value) {
    696833  if (isnan(value))
    697     gprint (GP_LOG, "  NaN    ");
     834    gprint (GP_LOG, "       NaN ");
    698835  else
    699836    gprint (GP_LOG, "%10.3e ", value);
    … …  
    702839void print_short (double value, short int ival) {
    703840  if (ival == NAN_S_SHORT)
    704     gprint (GP_LOG, "  NaN    ");
     841    gprint (GP_LOG, "       NaN ");
    705842  else
    706     gprint (GP_LOG, "  %6.3f ", value);
     843    gprint (GP_LOG, "   %7.4f ", value);
    707844}
    708845
    … …  
    763900        print_double (NAN);
    764901      } else {
    765         print_double (secfilt[seq].Mchisq);
     902        print_double_exp (secfilt[seq].Mchisq);
    766903      }
    767904      break;
    … …  
    769906    case GSTAR_FILTER_NAMES: /* filter names */
    770907      code = GetPhotcodebyNsec (seq);
    771       gprint (GP_LOG, "  %-6s ", code[0].name);
     908      gprint (GP_LOG, "%10s ", code[0].name);
    772909      break;
    773910
    774911    case GSTAR_SECF_FLAGS: /* secfilt flags */
    775912      if (seq == -1) {
    776         gprint (GP_LOG, "%08x ", 0);
    777       } else {
    778         gprint (GP_LOG, "%08x ", secfilt[seq].flags);
     913        gprint (GP_LOG, "0x%08x ", 0);
     914      } else {
     915        gprint (GP_LOG, "0x%08x ", secfilt[seq].flags);
    779916      }
    780917      break;
    … …  
    782919    case GSTAR_UCDIST: /* ubercal distance */
    783920      if (seq == -1) {
    784         gprint (GP_LOG, "%5d ", 0);
    785       } else {
    786         gprint (GP_LOG, "%5d ", secfilt[seq].ubercalDist);
     921        gprint (GP_LOG, "%10d ", 0);
     922      } else {
     923        gprint (GP_LOG, "%10d ", secfilt[seq].ubercalDist);
     924      }
     925      break;
     926
     927    case GSTAR_AVE_NCODE: /* Ncode */
     928      if (seq == -1) {
     929        gprint (GP_LOG, "%10d ", 0);
     930      } else {
     931        gprint (GP_LOG, "%10d ", secfilt[seq].Ncode);
     932      }
     933      break;
     934
     935    case GSTAR_AVE_NUSED: /* Nused */
     936      if (seq == -1) {
     937        gprint (GP_LOG, "%10d ", 0);
     938      } else {
     939        gprint (GP_LOG, "%10d ", secfilt[seq].Nused);
     940      }
     941      break;
     942
     943    case GSTAR_AVE_KRON_NUSED: /* Nused */
     944      if (seq == -1) {
     945        gprint (GP_LOG, "%10d ", 0);
     946      } else {
     947        gprint (GP_LOG, "%10d ", secfilt[seq].NusedKron);
     948      }
     949      break;
     950
     951    case GSTAR_AVE_AP_NUSED: /* Nused */
     952      if (seq == -1) {
     953        gprint (GP_LOG, "%10d ", 0);
     954      } else {
     955        gprint (GP_LOG, "%10d ", secfilt[seq].NusedAp);
    787956      }
    788957      break;
    … …  
    814983      break;
    815984
     985    case GSTAR_AVE_AP_MAG_ERR: /* average ap mags */
     986      if (seq == -1) {
     987        print_double (NAN);
     988      } else {
     989        print_double (secfilt[seq].dMap);
     990      }
     991      break;
     992
     993    case GSTAR_AVE_AP_MAG_CHISQ: /* average ap mags */
     994      if (seq == -1) {
     995        print_double (NAN);
     996      } else {
     997        print_double_exp (secfilt[seq].sMap);
     998      }
     999      break;
     1000
    8161001    case GSTAR_AVE_KRON_MAG: /* average ap mags */
    8171002      if (seq == -1) {
    … …  
    8301015      break;
    8311016
     1017    case GSTAR_AVE_KRON_MAG_CHISQ: /* average ap mags */
     1018      if (seq == -1) {
     1019        print_double (NAN);
     1020      } else {
     1021        print_double_exp (secfilt[seq].sMkron);
     1022      }
     1023      break;
     1024
    8321025      /******************* STACK : FLUX **********************/
    8331026
    … …  
    8801073      break;
    8811074
     1075    case GSTAR_STACK_NUSED: /* Nused */
     1076      if (seq == -1) {
     1077        gprint (GP_LOG, "%10d ", 0);
     1078      } else {
     1079        gprint (GP_LOG, "%10d ", secfilt[seq].Nstack);
     1080      }
     1081      break;
     1082
    8821083      /******************* STACK : MAG **********************/
    8831084
    … …  
    9481149      break;
    9491150
     1151    case GSTAR_WARP_FLUX_PSF_CHISQ: /* average ap mags */
     1152      if (seq == -1) {
     1153        print_double (NAN);
     1154      } else {
     1155        print_double_exp (secfilt[seq].sFpsfWrp);
     1156      }
     1157      break;
     1158
    9501159    case GSTAR_WARP_FLUX_KRON: /* average ap mags */
    9511160      if (seq == -1) {
    … …  
    9641173      break;
    9651174
     1175    case GSTAR_WARP_FLUX_KRON_CHISQ: /* average ap mags */
     1176      if (seq == -1) {
     1177        print_double (NAN);
     1178      } else {
     1179        print_double_exp (secfilt[seq].sFkronWrp);
     1180      }
     1181      break;
     1182
    9661183    case GSTAR_WARP_FLUX_APER: /* average ap mags */
    9671184      if (seq == -1) {
    … …  
    9801197      break;
    9811198
     1199    case GSTAR_WARP_FLUX_APER_CHISQ: /* average ap mags */
     1200      if (seq == -1) {
     1201        print_double_exp (NAN);
     1202      } else {
     1203        print_double_exp (secfilt[seq].sFapWrp);
     1204      }
     1205      break;
     1206
    9821207      /******************* WARP : MAG **********************/
    9831208
    … …  
    9981223      break;
    9991224
     1225    case GSTAR_WARP_MAG_PSF_CHISQ: /* average ap mags */
     1226      if (seq == -1) {
     1227        print_double_exp (NAN);
     1228      } else {
     1229        print_double_exp (secfilt[seq].sFpsfWrp);
     1230      }
     1231      break;
     1232
    10001233    case GSTAR_WARP_MAG_KRON: /* average ap mags */
    10011234      if (seq == -1) {
    … …  
    10141247      break;
    10151248
     1249    case GSTAR_WARP_MAG_KRON_CHISQ: /* average ap mags */
     1250      if (seq == -1) {
     1251        print_double_exp (NAN);
     1252      } else {
     1253        print_double_exp (secfilt[seq].sFkronWrp);
     1254      }
     1255      break;
     1256
    10161257    case GSTAR_WARP_MAG_APER: /* average ap mags */
    10171258      if (seq == -1) {
    … …  
    10301271      break;
    10311272
     1273    case GSTAR_WARP_MAG_APER_CHISQ: /* average ap mags */
     1274      if (seq == -1) {
     1275        print_double_exp (NAN);
     1276      } else {
     1277        print_double_exp (secfilt[seq].sFapWrp);
     1278      }
     1279      break;
     1280
     1281    case GSTAR_WARP_NWARP: /* Nused */
     1282      if (seq == -1) {
     1283        gprint (GP_LOG, "%10d ", 0);
     1284      } else {
     1285        gprint (GP_LOG, "%10d ", secfilt[seq].Nwarp);
     1286      }
     1287      break;
     1288
     1289    case GSTAR_WARP_NWARP_GOOD: /* Nused */
     1290      if (seq == -1) {
     1291        gprint (GP_LOG, "%10d ", 0);
     1292      } else {
     1293        gprint (GP_LOG, "%10d ", secfilt[seq].NwarpGood);
     1294      }
     1295      break;
     1296
     1297    case GSTAR_WARP_NWARP_USED_PSF: /* Nused */
     1298      if (seq == -1) {
     1299        gprint (GP_LOG, "%10d ", 0);
     1300      } else {
     1301        gprint (GP_LOG, "%10d ", secfilt[seq].NusedWrp);
     1302      }
     1303      break;
     1304
     1305    case GSTAR_WARP_NWARP_USED_KRON: /* Nused */
     1306      if (seq == -1) {
     1307        gprint (GP_LOG, "%10d ", 0);
     1308      } else {
     1309        gprint (GP_LOG, "%10d ", secfilt[seq].NusedKronWrp);
     1310      }
     1311      break;
     1312
     1313    case GSTAR_WARP_NWARP_USED_AP: /* Nused */
     1314      if (seq == -1) {
     1315        gprint (GP_LOG, "%10d ", 0);
     1316      } else {
     1317        gprint (GP_LOG, "%10d ", secfilt[seq].NusedApWrp);
     1318      }
     1319      break;
    10321320  }
    10331321}
Note: See TracChangeset for help on using the changeset viewer.