IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 38436


Ignore:
Timestamp:
Jun 11, 2015, 11:45:55 AM (11 years ago)
Author:
eugene
Message:

various gstar tweaks to make the output more readable

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/eam_branches/ipp-20150419/Ohana/src/opihi/dvo/gstar.c

    r38313 r38436  
    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,
    2124  GSTAR_AVE_NCODE,
    2225  GSTAR_AVE_NUSED,
    … …  
    3134  GSTAR_WARP_FLUX_PSF,
    3235  GSTAR_WARP_FLUX_PSF_ERR,
     36  GSTAR_WARP_FLUX_PSF_CHISQ,
    3337  GSTAR_WARP_FLUX_KRON,
    3438  GSTAR_WARP_FLUX_KRON_ERR,
     39  GSTAR_WARP_FLUX_KRON_CHISQ,
    3540  GSTAR_WARP_FLUX_APER,
    3641  GSTAR_WARP_FLUX_APER_ERR,
     42  GSTAR_WARP_FLUX_APER_CHISQ,
    3743  GSTAR_STACK_MAG_PSF,
    3844  GSTAR_STACK_MAG_PSF_ERR,
    … …  
    4147  GSTAR_STACK_MAG_APER,
    4248  GSTAR_STACK_MAG_APER_ERR,
     49  GSTAR_STACK_NUSED,
    4350  GSTAR_WARP_MAG_PSF,
    4451  GSTAR_WARP_MAG_PSF_ERR,
     52  GSTAR_WARP_MAG_PSF_CHISQ,
    4553  GSTAR_WARP_MAG_KRON,
    4654  GSTAR_WARP_MAG_KRON_ERR,
     55  GSTAR_WARP_MAG_KRON_CHISQ,
    4756  GSTAR_WARP_MAG_APER,
    4857  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,
    4964} GSTAR_SECF_CODES;
    5065
    … …  
    348363        gprint (GP_LOG, "\n");
    349364
     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
    350377        /* average mag errors */
    351378        gprint (GP_LOG, "chp_psf_err: ");
    … …  
    374401          gprint (GP_LOG, "\n");
    375402
     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
    376411          /* average mag Nused */
    377412          gprint (GP_LOG, "chp_ap__N:   ");
    … …  
    391426          gprint (GP_LOG, "\n");
    392427
     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
    393432          /* average mag Nused */
    394433          gprint (GP_LOG, "chp_krn_N:   ");
    … …  
    397436        }
    398437
    399         if (EXTRA_OUTPUT) {
    400           /* Mmin */
    401           gprint (GP_LOG, "chp_psf_min: ");
    402           for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_MAG_MIN);
    403           gprint (GP_LOG, "\n");
    404 
    405           /* Mmax */
    406           gprint (GP_LOG, "chp_psf_max: ");
    407           for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_AVE_MAG_MAX);;
    408           gprint (GP_LOG, "\n");
    409 
    410           /* UCDIST */
    411           gprint (GP_LOG, "ubercal_dst: ");
    412           for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_UCDIST);;
    413           gprint (GP_LOG, "\n");
    414 
    415           /* secfilt flags */
    416           gprint (GP_LOG, "filtflags:   ");
    417           for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_SECF_FLAGS);
    418           gprint (GP_LOG, "\n");
    419         }
    420438        if (STACK_OUTPUT && MAG_OUTPUT) {
    421439          /* FluxPSF */
    … …  
    429447          gprint (GP_LOG, "\n");
    430448 
     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
    431455          if (KRON_OUTPUT) {
    432456            /* MagKron */
    … …  
    464488          gprint (GP_LOG, "\n");
    465489 
     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         
    466508          if (KRON_OUTPUT) {
    467509            /* MagKron */
    … …  
    475517            gprint (GP_LOG, "\n");
    476518
     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            }
    477529          }
    478 
     530         
    479531          if (APER_OUTPUT) {
    480532            /* MagAper */
    … …  
    487539            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_MAG_APER_ERR);
    488540            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            }
    489552          }
    490553        }
    … …  
    514577          if (APER_OUTPUT) {
    515578            /* FluxAper */
    516             gprint (GP_LOG, "stk_ap_ave: ");
     579            gprint (GP_LOG, "stk_ap_ave:  ");
    517580            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_STACK_FLUX_APER);
    518581            gprint (GP_LOG, "\n");
    519582 
    520583            /* dFluxAper */
    521             gprint (GP_LOG, "stk_ap_err: ");
     584            gprint (GP_LOG, "stk_ap_err:  ");
    522585            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_STACK_FLUX_APER_ERR);
    523586            gprint (GP_LOG, "\n");
    524587          }
     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");
    525592        }
    526593        if (WARP_OUTPUT && FLUX_OUTPUT) {
    … …  
    531598 
    532599          /* dFluxPSF */
    533           gprint (GP_LOG, "wrp_psf_ave: ");
     600          gprint (GP_LOG, "wrp_psf_err: ");
    534601          for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_FLUX_PSF_ERR);
    535602          gprint (GP_LOG, "\n");
    536603 
     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
    537620          if (KRON_OUTPUT) {
    538621            /* FluxKron */
    … …  
    545628            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_FLUX_KRON_ERR);
    546629            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");
    547638          }
    548639
    … …  
    557648            for (j = 0; j < Nsecfilt; j++) printPhotcodeSequence (&catalog.average[k], &catalog.secfilt[Nsecfilt*k], j, GSTAR_WARP_FLUX_APER_ERR);
    558649            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");
    559658          }
    560659        }
    561       }
     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
    562673      gprint (GP_LOG, "\n");
    563674     
    … …  
    649760              gprint (GP_LOG, "0x%08x ", catalog.measure[Nv].imageID);
    650761              gprint (GP_LOG, "%5.3f ", catalog.measure[Nv].psfQF);
     762              gprint (GP_LOG, "%5.3f ", catalog.measure[Nv].psfQFperf);
    651763              gprint (GP_LOG, "%7.1f ", catalog.measure[Nv].psfChisq);
    652764              // gprint (GP_LOG, "%3.1f ", catalog.measure[Nv].crNsigma);
    … …  
    659771              gprint (GP_LOG, "%10.3e ", catalog.measure[Nv].FluxKron);
    660772              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);
    661775            }
    662776            gprint (GP_LOG, "\n");
    … …  
    711825void print_double (double value) {
    712826  if (isnan(value))
    713     gprint (GP_LOG, "  NaN    ");
     827    gprint (GP_LOG, "       NaN ");
    714828  else
    715     gprint (GP_LOG, "  %6.3f ", value);
     829    gprint (GP_LOG, "   %7.4f ", value);
    716830}
    717831
    718832void print_double_exp (double value) {
    719833  if (isnan(value))
    720     gprint (GP_LOG, "  NaN    ");
     834    gprint (GP_LOG, "       NaN ");
    721835  else
    722836    gprint (GP_LOG, "%10.3e ", value);
    … …  
    725839void print_short (double value, short int ival) {
    726840  if (ival == NAN_S_SHORT)
    727     gprint (GP_LOG, "  NaN    ");
     841    gprint (GP_LOG, "       NaN ");
    728842  else
    729     gprint (GP_LOG, "  %6.3f ", value);
     843    gprint (GP_LOG, "   %7.4f ", value);
    730844}
    731845
    … …  
    786900        print_double (NAN);
    787901      } else {
    788         print_double (secfilt[seq].Mchisq);
     902        print_double_exp (secfilt[seq].Mchisq);
    789903      }
    790904      break;
    … …  
    792906    case GSTAR_FILTER_NAMES: /* filter names */
    793907      code = GetPhotcodebyNsec (seq);
    794       gprint (GP_LOG, "  %-6s ", code[0].name);
     908      gprint (GP_LOG, "%10s ", code[0].name);
    795909      break;
    796910
    … …  
    805919    case GSTAR_UCDIST: /* ubercal distance */
    806920      if (seq == -1) {
    807         gprint (GP_LOG, "%5d ", 0);
    808       } else {
    809         gprint (GP_LOG, "%5d ", secfilt[seq].ubercalDist);
     921        gprint (GP_LOG, "%10d ", 0);
     922      } else {
     923        gprint (GP_LOG, "%10d ", secfilt[seq].ubercalDist);
    810924      }
    811925      break;
    … …  
    813927    case GSTAR_AVE_NCODE: /* Ncode */
    814928      if (seq == -1) {
    815         gprint (GP_LOG, "%3d ", 0);
    816       } else {
    817         gprint (GP_LOG, "%3d ", secfilt[seq].Ncode);
     929        gprint (GP_LOG, "%10d ", 0);
     930      } else {
     931        gprint (GP_LOG, "%10d ", secfilt[seq].Ncode);
    818932      }
    819933      break;
    … …  
    821935    case GSTAR_AVE_NUSED: /* Nused */
    822936      if (seq == -1) {
    823         gprint (GP_LOG, "%3d ", 0);
    824       } else {
    825         gprint (GP_LOG, "%3d ", secfilt[seq].Nused);
     937        gprint (GP_LOG, "%10d ", 0);
     938      } else {
     939        gprint (GP_LOG, "%10d ", secfilt[seq].Nused);
    826940      }
    827941      break;
    … …  
    829943    case GSTAR_AVE_KRON_NUSED: /* Nused */
    830944      if (seq == -1) {
    831         gprint (GP_LOG, "%3d ", 0);
    832       } else {
    833         gprint (GP_LOG, "%3d ", secfilt[seq].NusedKron);
     945        gprint (GP_LOG, "%10d ", 0);
     946      } else {
     947        gprint (GP_LOG, "%10d ", secfilt[seq].NusedKron);
    834948      }
    835949      break;
    … …  
    837951    case GSTAR_AVE_AP_NUSED: /* Nused */
    838952      if (seq == -1) {
    839         gprint (GP_LOG, "%3d ", 0);
    840       } else {
    841         gprint (GP_LOG, "%3d ", secfilt[seq].NusedAp);
     953        gprint (GP_LOG, "%10d ", 0);
     954      } else {
     955        gprint (GP_LOG, "%10d ", secfilt[seq].NusedAp);
    842956      }
    843957      break;
    … …  
    869983      break;
    870984
     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
    8711001    case GSTAR_AVE_KRON_MAG: /* average ap mags */
    8721002      if (seq == -1) {
    … …  
    8851015      break;
    8861016
     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
    8871025      /******************* STACK : FLUX **********************/
    8881026
    … …  
    9351073      break;
    9361074
     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
    9371083      /******************* STACK : MAG **********************/
    9381084
    … …  
    10031149      break;
    10041150
     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
    10051159    case GSTAR_WARP_FLUX_KRON: /* average ap mags */
    10061160      if (seq == -1) {
    … …  
    10191173      break;
    10201174
     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
    10211183    case GSTAR_WARP_FLUX_APER: /* average ap mags */
    10221184      if (seq == -1) {
    … …  
    10351197      break;
    10361198
     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
    10371207      /******************* WARP : MAG **********************/
    10381208
    … …  
    10531223      break;
    10541224
     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
    10551233    case GSTAR_WARP_MAG_KRON: /* average ap mags */
    10561234      if (seq == -1) {
    … …  
    10691247      break;
    10701248
     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
    10711257    case GSTAR_WARP_MAG_APER: /* average ap mags */
    10721258      if (seq == -1) {
    … …  
    10851271      break;
    10861272
     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;
    10871320  }
    10881321}
Note: See TracChangeset for help on using the changeset viewer.