IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Apr 18, 2015, 6:33:51 AM (11 years ago)
Author:
eugene
Message:

merge changes from eam branch ipp-20150405: add R7 to lensing, fix R5,R6,R7 definitions and mean photometry; calculate E1,E2 in dvolens; fix precision for dvopsps FW objects; fix histograms at axis; new colormap options; require primary skycell in diff mean photometry

Location:
trunk
Files:
11 edited

Legend:

Unmodified
Added
Removed

  • trunk

  • trunk/Ohana/src/libdvo/include/dvo.h

    r38062 r38153  
    11501150LENSFIELD(E2_sh_psf)
    11511151
     1152LENSFIELD( F_ApR5);
     1153LENSFIELD(dF_ApR5);
     1154LENSFIELD(sF_ApR5);
     1155LENSFIELD(fF_ApR5);
     1156
     1157LENSFIELD( F_ApR6);
     1158LENSFIELD(dF_ApR6);
     1159LENSFIELD(sF_ApR6);
     1160LENSFIELD(fF_ApR6);
     1161
     1162LENSFIELD( F_ApR7);
     1163LENSFIELD(dF_ApR7);
     1164LENSFIELD(sF_ApR7);
     1165LENSFIELD(fF_ApR7);
     1166
     1167LENSFIELD(E1);
     1168LENSFIELD(E2);
     1169
    11521170# undef LENSFIELD
    11531171

  • trunk/Ohana/src/libdvo/include/dvodb.h

    r38062 r38153  
    8989  MAG_OPTION_F_STDEV_AP_R6,
    9090  MAG_OPTION_F_FILL_AP_R6,
     91
     92  MAG_OPTION_F_AP_R7,
     93  MAG_OPTION_F_ERR_AP_R7,
     94  MAG_OPTION_F_STDEV_AP_R7,
     95  MAG_OPTION_F_FILL_AP_R7,
     96
     97  MAG_OPTION_E1,
     98  MAG_OPTION_E2,
    9199} dvoMagOptionType;
    92100
    … …  
    207215      MEAS_REF_COLOR_BLUE,
    208216      MEAS_REF_COLOR_RED,
     217
    209218      MEAS_X11_SM_OBJ,
    210219      MEAS_X12_SM_OBJ,
    … …  
    227236      MEAS_E1_SH_PSF,
    228237      MEAS_E2_SH_PSF,
     238
     239      MEAS_F_AP_R5,
     240      MEAS_F_ERR_AP_R5,
     241      MEAS_F_STDEV_AP_R5,
     242      MEAS_F_FILL_AP_R5,
     243
     244      MEAS_F_AP_R6,
     245      MEAS_F_ERR_AP_R6,
     246      MEAS_F_STDEV_AP_R6,
     247      MEAS_F_FILL_AP_R6,
     248
     249      MEAS_F_AP_R7,
     250      MEAS_F_ERR_AP_R7,
     251      MEAS_F_STDEV_AP_R7,
     252      MEAS_F_FILL_AP_R7,
     253
    229254      MEAS_E_BV,                      // extinction (mags)
    230255      MEAS_E_BV_ERR,

  • trunk/Ohana/src/libdvo/include/ps1_v5_defs.h

    r37995 r38153  
    88SecFilt_PS1_V5          *SecFiltInternalTo_PS1_V5 (SecFilt *in, off_t Nvalues);
    99
    10 typedef Lensing_PS1_V5_R2 Lensing_PS1_V5;
     10Lensobj                 *Lensobj_PS1_V5_R0_ToInternal (Lensobj_PS1_V5_R0 *in, off_t Nvalues);
     11Lensobj                 *Lensobj_PS1_V5_R1_ToInternal (Lensobj_PS1_V5_R1 *in, off_t Nvalues);
    1112
    12 Lensing                 *Lensing_PS1_V5_ToInternal (Lensing_PS1_V5 *in, off_t Nvalues);
    13 Lensing_PS1_V5_R2       *LensingInternalTo_PS1_V5_R2 (Lensing *in, off_t Nvalues);
    14 Lensobj                 *Lensobj_PS1_V5_ToInternal (Lensobj_PS1_V5 *in, off_t Nvalues);
    15 Lensobj_PS1_V5          *LensobjInternalTo_PS1_V5 (Lensobj *in, off_t Nvalues);
     13Lensobj_PS1_V5_R1       *LensobjInternalTo_PS1_V5_R1 (Lensobj *in, off_t Nvalues);
    1614
    1715Lensing                 *Lensing_PS1_V5_R0_ToInternal (Lensing_PS1_V5_R0 *in, off_t Nvalues);
    1816Lensing                 *Lensing_PS1_V5_R1_ToInternal (Lensing_PS1_V5_R1 *in, off_t Nvalues);
    1917Lensing                 *Lensing_PS1_V5_R2_ToInternal (Lensing_PS1_V5_R2 *in, off_t Nvalues);
     18Lensing                 *Lensing_PS1_V5_R3_ToInternal (Lensing_PS1_V5_R3 *in, off_t Nvalues);
     19
     20Lensing_PS1_V5_R3       *LensingInternalTo_PS1_V5_R3 (Lensing *in, off_t Nvalues);
    2021
    2122StarPar                 *StarPar_PS1_V5_ToInternal (StarPar_PS1_V5 *in, off_t Nvalues);

  • trunk/Ohana/src/libdvo/src/dbExtractAverages.c

    r38062 r38153  
    386386        case MAG_OPTION_E1_SH_PSF:  { value.Flt = LensValue_E1_sh_psf  (field->photcode, lensobj); break; }
    387387        case MAG_OPTION_E2_SH_PSF:  { value.Flt = LensValue_E2_sh_psf  (field->photcode, lensobj); break; }
     388
     389        case MAG_OPTION_F_AP_R5:       { value.Flt = LensValue_F_ApR5  (field->photcode, lensobj); break; }
     390        case MAG_OPTION_F_ERR_AP_R5:   { value.Flt = LensValue_dF_ApR5 (field->photcode, lensobj); break; }
     391        case MAG_OPTION_F_STDEV_AP_R5: { value.Flt = LensValue_sF_ApR5 (field->photcode, lensobj); break; }
     392        case MAG_OPTION_F_FILL_AP_R5:  { value.Flt = LensValue_fF_ApR5 (field->photcode, lensobj); break; }
     393
     394        case MAG_OPTION_F_AP_R6:       { value.Flt = LensValue_F_ApR6  (field->photcode, lensobj); break; }
     395        case MAG_OPTION_F_ERR_AP_R6:   { value.Flt = LensValue_dF_ApR6 (field->photcode, lensobj); break; }
     396        case MAG_OPTION_F_STDEV_AP_R6: { value.Flt = LensValue_sF_ApR6 (field->photcode, lensobj); break; }
     397        case MAG_OPTION_F_FILL_AP_R6:  { value.Flt = LensValue_fF_ApR6 (field->photcode, lensobj); break; }
     398
     399        case MAG_OPTION_F_AP_R7:       { value.Flt = LensValue_F_ApR7  (field->photcode, lensobj); break; }
     400        case MAG_OPTION_F_ERR_AP_R7:   { value.Flt = LensValue_dF_ApR7 (field->photcode, lensobj); break; }
     401        case MAG_OPTION_F_STDEV_AP_R7: { value.Flt = LensValue_sF_ApR7 (field->photcode, lensobj); break; }
     402        case MAG_OPTION_F_FILL_AP_R7:  { value.Flt = LensValue_fF_ApR7 (field->photcode, lensobj); break; }
     403
     404        case MAG_OPTION_E1:            { value.Flt = LensValue_E1      (field->photcode, lensobj); break; }
     405        case MAG_OPTION_E2:            { value.Flt = LensValue_E2      (field->photcode, lensobj); break; }
    388406
    389407        case MAG_OPTION_NONE:

  • trunk/Ohana/src/libdvo/src/dbExtractMeasures.c

    r38062 r38153  
    277277        case MAG_OPTION_F_STDEV_AP_R6:
    278278        case MAG_OPTION_F_FILL_AP_R6:
     279        case MAG_OPTION_F_AP_R7:
     280        case MAG_OPTION_F_ERR_AP_R7:
     281        case MAG_OPTION_F_STDEV_AP_R7:
     282        case MAG_OPTION_F_FILL_AP_R7:
     283        case MAG_OPTION_E1:
     284        case MAG_OPTION_E2:
    279285
    280286        case MAG_OPTION_NONE:
    … …  
    735741      }
    736742      break;
     743
     744    case MEAS_X11_SM_OBJ: { value.Flt = lensing->X11_sm_obj; break; }
     745    case MEAS_X12_SM_OBJ: { value.Flt = lensing->X12_sm_obj; break; }
     746    case MEAS_X22_SM_OBJ: { value.Flt = lensing->X22_sm_obj; break; }
     747    case MEAS_E1_SM_OBJ:  { value.Flt = lensing->E1_sm_obj; break; }
     748    case MEAS_E2_SM_OBJ:  { value.Flt = lensing->E2_sm_obj; break; }
     749    case MEAS_X11_SH_OBJ: { value.Flt = lensing->X11_sh_obj; break; }
     750    case MEAS_X12_SH_OBJ: { value.Flt = lensing->X12_sh_obj; break; }
     751    case MEAS_X22_SH_OBJ: { value.Flt = lensing->X22_sh_obj; break; }
     752    case MEAS_E1_SH_OBJ:  { value.Flt = lensing->E1_sh_obj; break; }
     753    case MEAS_E2_SH_OBJ:  { value.Flt = lensing->E2_sh_obj; break; }
     754    case MEAS_X11_SM_PSF: { value.Flt = lensing->X11_sm_psf; break; }
     755    case MEAS_X12_SM_PSF: { value.Flt = lensing->X12_sm_psf; break; }
     756    case MEAS_X22_SM_PSF: { value.Flt = lensing->X22_sm_psf; break; }
     757    case MEAS_E1_SM_PSF:  { value.Flt = lensing->E1_sm_psf; break; }
     758    case MEAS_E2_SM_PSF:  { value.Flt = lensing->E2_sm_psf; break; }
     759    case MEAS_X11_SH_PSF: { value.Flt = lensing->X11_sh_psf; break; }
     760    case MEAS_X12_SH_PSF: { value.Flt = lensing->X12_sh_psf; break; }
     761    case MEAS_X22_SH_PSF: { value.Flt = lensing->X22_sh_psf; break; }
     762    case MEAS_E1_SH_PSF:  { value.Flt = lensing->E1_sh_psf; break; }
     763    case MEAS_E2_SH_PSF:  { value.Flt = lensing->E2_sh_psf; break; }
     764
     765    case MEAS_F_AP_R5:       { value.Flt = lensing-> F_ApR5; break; }
     766    case MEAS_F_ERR_AP_R5:   { value.Flt = lensing->dF_ApR5; break; }
     767    case MEAS_F_STDEV_AP_R5: { value.Flt = lensing->sF_ApR5; break; }
     768    case MEAS_F_FILL_AP_R5:  { value.Flt = lensing->fF_ApR5; break; }
     769    case MEAS_F_AP_R6:       { value.Flt = lensing-> F_ApR6; break; }
     770    case MEAS_F_ERR_AP_R6:   { value.Flt = lensing->dF_ApR6; break; }
     771    case MEAS_F_STDEV_AP_R6: { value.Flt = lensing->sF_ApR6; break; }
     772    case MEAS_F_FILL_AP_R6:  { value.Flt = lensing->fF_ApR6; break; }
     773    case MEAS_F_AP_R7:       { value.Flt = lensing-> F_ApR7; break; }
     774    case MEAS_F_ERR_AP_R7:   { value.Flt = lensing->dF_ApR7; break; }
     775    case MEAS_F_STDEV_AP_R7: { value.Flt = lensing->sF_ApR7; break; }
     776    case MEAS_F_FILL_AP_R7:  { value.Flt = lensing->fF_ApR7; break; }
    737777
    738778    case MEAS_E_BV:             { value.Flt = starpar->Ebv;      break; }

  • trunk/Ohana/src/libdvo/src/dbFields.c

    r38062 r38153  
    6161    if (fields[i].ID == MEAS_E1_SH_PSF)  return TRUE;
    6262    if (fields[i].ID == MEAS_E2_SH_PSF)  return TRUE;
     63
     64    if (fields[i].ID == MEAS_F_AP_R5)       return TRUE;
     65    if (fields[i].ID == MEAS_F_ERR_AP_R5)   return TRUE;
     66    if (fields[i].ID == MEAS_F_STDEV_AP_R5) return TRUE;
     67    if (fields[i].ID == MEAS_F_FILL_AP_R5)  return TRUE;
     68
     69    if (fields[i].ID == MEAS_F_AP_R6)       return TRUE;
     70    if (fields[i].ID == MEAS_F_ERR_AP_R6)   return TRUE;
     71    if (fields[i].ID == MEAS_F_STDEV_AP_R6) return TRUE;
     72    if (fields[i].ID == MEAS_F_FILL_AP_R6)  return TRUE;
     73
     74    if (fields[i].ID == MEAS_F_AP_R7)       return TRUE;
     75    if (fields[i].ID == MEAS_F_ERR_AP_R7)   return TRUE;
     76    if (fields[i].ID == MEAS_F_STDEV_AP_R7) return TRUE;
     77    if (fields[i].ID == MEAS_F_FILL_AP_R7)  return TRUE;
    6378  }
    6479  return FALSE;
    … …  
    90105    if (fields[i].magOption == MAG_OPTION_E1_SH_PSF)  return TRUE;
    91106    if (fields[i].magOption == MAG_OPTION_E2_SH_PSF)  return TRUE;
     107
     108    if (fields[i].magOption == MAG_OPTION_F_AP_R5)       return TRUE;
     109    if (fields[i].magOption == MAG_OPTION_F_ERR_AP_R5)   return TRUE;
     110    if (fields[i].magOption == MAG_OPTION_F_STDEV_AP_R5) return TRUE;
     111    if (fields[i].magOption == MAG_OPTION_F_FILL_AP_R5)  return TRUE;
     112
     113    if (fields[i].magOption == MAG_OPTION_F_AP_R6)       return TRUE;
     114    if (fields[i].magOption == MAG_OPTION_F_ERR_AP_R6)   return TRUE;
     115    if (fields[i].magOption == MAG_OPTION_F_STDEV_AP_R6) return TRUE;
     116    if (fields[i].magOption == MAG_OPTION_F_FILL_AP_R6)  return TRUE;
     117
     118    if (fields[i].magOption == MAG_OPTION_F_AP_R7)       return TRUE;
     119    if (fields[i].magOption == MAG_OPTION_F_ERR_AP_R7)   return TRUE;
     120    if (fields[i].magOption == MAG_OPTION_F_STDEV_AP_R7) return TRUE;
     121    if (fields[i].magOption == MAG_OPTION_F_FILL_AP_R7)  return TRUE;
     122
     123    if (fields[i].magOption == MAG_OPTION_E1)            return TRUE;
     124    if (fields[i].magOption == MAG_OPTION_E2)            return TRUE;
    92125  }
    93126  return FALSE;
    … …  
    199232  if (!strcasecmp (string, "E1_SH_PSF"))     return MAG_OPTION_E1_SH_PSF;
    200233  if (!strcasecmp (string, "E2_SH_PSF"))     return MAG_OPTION_E2_SH_PSF;
    201   // if (!strcasecmp (string, "F_AP_R5"))       return MAG_OPTION_F_AP_R5;
    202   // if (!strcasecmp (string, "F_ERR_AP_R5"))   return MAG_OPTION_F_ERR_AP_R5;
    203   // if (!strcasecmp (string, "F_STDEV_AP_R5")) return MAG_OPTION_F_STDEV_AP_R5;
    204   // if (!strcasecmp (string, "F_FILL_AP_R5"))  return MAG_OPTION_F_FILL_AP_R5;
    205   // if (!strcasecmp (string, "F_AP_R6"))       return MAG_OPTION_F_AP_R6;
    206   // if (!strcasecmp (string, "F_ERR_AP_R6"))   return MAG_OPTION_F_ERR_AP_R6;
    207   // if (!strcasecmp (string, "F_STDEV_AP_R6")) return MAG_OPTION_F_STDEV_AP_R6;
    208   // if (!strcasecmp (string, "F_FILL_AP_R6"))  return MAG_OPTION_F_FILL_AP_R6;
     234  if (!strcasecmp (string, "F_AP_R5"))       return MAG_OPTION_F_AP_R5;
     235  if (!strcasecmp (string, "F_ERR_AP_R5"))   return MAG_OPTION_F_ERR_AP_R5;
     236  if (!strcasecmp (string, "F_STDEV_AP_R5")) return MAG_OPTION_F_STDEV_AP_R5;
     237  if (!strcasecmp (string, "F_FILL_AP_R5"))  return MAG_OPTION_F_FILL_AP_R5;
     238  if (!strcasecmp (string, "F_AP_R6"))       return MAG_OPTION_F_AP_R6;
     239  if (!strcasecmp (string, "F_ERR_AP_R6"))   return MAG_OPTION_F_ERR_AP_R6;
     240  if (!strcasecmp (string, "F_STDEV_AP_R6")) return MAG_OPTION_F_STDEV_AP_R6;
     241  if (!strcasecmp (string, "F_FILL_AP_R6"))  return MAG_OPTION_F_FILL_AP_R6;
     242  if (!strcasecmp (string, "F_AP_R7"))       return MAG_OPTION_F_AP_R7;
     243  if (!strcasecmp (string, "F_ERR_AP_R7"))   return MAG_OPTION_F_ERR_AP_R7;
     244  if (!strcasecmp (string, "F_STDEV_AP_R7")) return MAG_OPTION_F_STDEV_AP_R7;
     245  if (!strcasecmp (string, "F_FILL_AP_R7"))  return MAG_OPTION_F_FILL_AP_R7;
     246  if (!strcasecmp (string, "E1"))            return MAG_OPTION_E1;
     247  if (!strcasecmp (string, "E2"))            return MAG_OPTION_E2;
    209248
    210249  return MAG_OPTION_NONE;
    … …  
    534573  if (!strcasecmp (fieldName, "E2_SH_PSF"))      ESCAPE (MEAS_E2_SH_PSF,  OPIHI_FLT);
    535574
     575  if (!strcasecmp (fieldName, "F_AP_R5"))        ESCAPE (MEAS_F_AP_R5,       OPIHI_FLT);
     576  if (!strcasecmp (fieldName, "F_ERR_AP_R5"))    ESCAPE (MEAS_F_ERR_AP_R5,   OPIHI_FLT);
     577  if (!strcasecmp (fieldName, "F_STDEV_AP_R5"))  ESCAPE (MEAS_F_STDEV_AP_R5, OPIHI_FLT);
     578  if (!strcasecmp (fieldName, "F_FILL_AP_R5"))   ESCAPE (MEAS_F_FILL_AP_R5,  OPIHI_FLT);
     579  if (!strcasecmp (fieldName, "F_AP_R6"))        ESCAPE (MEAS_F_AP_R6,       OPIHI_FLT);
     580  if (!strcasecmp (fieldName, "F_ERR_AP_R6"))    ESCAPE (MEAS_F_ERR_AP_R6,   OPIHI_FLT);
     581  if (!strcasecmp (fieldName, "F_STDEV_AP_R6"))  ESCAPE (MEAS_F_STDEV_AP_R6, OPIHI_FLT);
     582  if (!strcasecmp (fieldName, "F_FILL_AP_R6"))   ESCAPE (MEAS_F_FILL_AP_R6,  OPIHI_FLT);
     583  if (!strcasecmp (fieldName, "F_AP_R7"))        ESCAPE (MEAS_F_AP_R7,       OPIHI_FLT);
     584  if (!strcasecmp (fieldName, "F_ERR_AP_R7"))    ESCAPE (MEAS_F_ERR_AP_R7,   OPIHI_FLT);
     585  if (!strcasecmp (fieldName, "F_STDEV_AP_R7"))  ESCAPE (MEAS_F_STDEV_AP_R7, OPIHI_FLT);
     586  if (!strcasecmp (fieldName, "F_FILL_AP_R7"))   ESCAPE (MEAS_F_FILL_AP_R7,  OPIHI_FLT);
     587
    536588  if (!strcasecmp (fieldName, "E_BV"))             ESCAPE (MEAS_E_BV          ,   OPIHI_FLT);
    537589  if (!strcasecmp (fieldName, "E_BV_ERR"))         ESCAPE (MEAS_E_BV_ERR      ,   OPIHI_FLT);

  • trunk/Ohana/src/libdvo/src/dvo_catalog_split.c

    r37807 r38153  
    992992    assert (catalog[0].Nlensing_disk >= catalog[0].Nlensing_off);
    993993
    994     // convert to external table format (note that the block above does not damage or free catalog.average)
    995     if (!LensingToFtable (&ftable, &catalog[0].lensing[first], Nrows, catalog[0].catformat)) {
    996       fprintf (stderr, "trouble converting format\n");
    997       goto failure;
    998     }
    999 
    1000     // write out Lensing table
    1001     if (!dvo_catalog_save_subcat (catalog[0].lensing_catalog, &ftable, start, Nrows, catalog[0].Nlensing_disk, Nlensing_disk_new)) {
    1002       fprintf (stderr, "trouble writing Lensing table\n");
    1003       goto failure;
    1004     }
    1005     gfits_free_header (&header);
    1006     gfits_free_table (&ftable);
     994    if (catalog[0].Nlensing) {
     995      // convert to external table format (note that the block above does not damage or free catalog.average)
     996      if (!LensingToFtable (&ftable, &catalog[0].lensing[first], Nrows, catalog[0].catformat)) {
     997        fprintf (stderr, "trouble converting format\n");
     998        goto failure;
     999      }
     1000
     1001      // write out Lensing table
     1002      if (!dvo_catalog_save_subcat (catalog[0].lensing_catalog, &ftable, start, Nrows, catalog[0].Nlensing_disk, Nlensing_disk_new)) {
     1003        fprintf (stderr, "trouble writing Lensing table\n");
     1004        goto failure;
     1005      }
     1006      gfits_free_header (&header);
     1007      gfits_free_table (&ftable);
     1008    }
    10071009  }
    10081010
    … …  
    10191021    assert (catalog[0].Nlensobj_disk >= catalog[0].Nlensobj_off);
    10201022
    1021     // convert to external table format (note that the block above does not damage or free catalog.average)
    1022     if (!LensobjToFtable (&ftable, &catalog[0].lensobj[first], Nrows, catalog[0].catformat)) {
    1023       fprintf (stderr, "trouble converting format\n");
    1024       goto failure;
    1025     }
    1026 
    1027     // write out Lensobj table
    1028     if (!dvo_catalog_save_subcat (catalog[0].lensobj_catalog, &ftable, start, Nrows, catalog[0].Nlensobj_disk, Nlensobj_disk_new)) {
    1029       fprintf (stderr, "trouble writing Lensobj table\n");
    1030       goto failure;
    1031     }
    1032     gfits_free_header (&header);
    1033     gfits_free_table (&ftable);
     1023    if (catalog[0].Nlensobj) {
     1024      // convert to external table format (note that the block above does not damage or free catalog.average)
     1025      if (!LensobjToFtable (&ftable, &catalog[0].lensobj[first], Nrows, catalog[0].catformat)) {
     1026        fprintf (stderr, "trouble converting format\n");
     1027        goto failure;
     1028      }
     1029
     1030      // write out Lensobj table
     1031      if (!dvo_catalog_save_subcat (catalog[0].lensobj_catalog, &ftable, start, Nrows, catalog[0].Nlensobj_disk, Nlensobj_disk_new)) {
     1032        fprintf (stderr, "trouble writing Lensobj table\n");
     1033        goto failure;
     1034      }
     1035      gfits_free_header (&header);
     1036      gfits_free_table (&ftable);
     1037    }
    10341038  }
    10351039
    … …  
    10461050    assert (catalog[0].Nstarpar_disk >= catalog[0].Nstarpar_off);
    10471051
    1048     // convert to external table format (note that the block above does not damage or free catalog.average)
    1049     if (!StarParToFtable (&ftable, &catalog[0].starpar[first], Nrows, catalog[0].catformat)) {
    1050       fprintf (stderr, "trouble converting format\n");
    1051       goto failure;
    1052     }
    1053 
    1054     // write out StarPar table
    1055     if (!dvo_catalog_save_subcat (catalog[0].starpar_catalog, &ftable, start, Nrows, catalog[0].Nstarpar_disk, Nstarpar_disk_new)) {
    1056       fprintf (stderr, "trouble writing StarPar table\n");
    1057       goto failure;
    1058     }
    1059     gfits_free_header (&header);
    1060     gfits_free_table (&ftable);
     1052    if (catalog[0].Nstarpar) {
     1053      // convert to external table format (note that the block above does not damage or free catalog.average)
     1054      if (!StarParToFtable (&ftable, &catalog[0].starpar[first], Nrows, catalog[0].catformat)) {
     1055        fprintf (stderr, "trouble converting format\n");
     1056        goto failure;
     1057      }
     1058
     1059      // write out StarPar table
     1060      if (!dvo_catalog_save_subcat (catalog[0].starpar_catalog, &ftable, start, Nrows, catalog[0].Nstarpar_disk, Nstarpar_disk_new)) {
     1061        fprintf (stderr, "trouble writing StarPar table\n");
     1062        goto failure;
     1063      }
     1064      gfits_free_header (&header);
     1065      gfits_free_table (&ftable);
     1066    }
    10611067  }
    10621068

  • trunk/Ohana/src/libdvo/src/dvo_convert.c

    r38062 r38153  
    455455    return (lensing); }
    456456
     457  if (!strcmp (extname, "DVO_LENSING_PS1_V5") && (ftable[0].header[0].Naxis[0] == 144)) {
     458    Lensing_PS1_V5_R2 *tmpLensing;                                             
     459    tmpLensing = gfits_table_get_Lensing_PS1_V5_R2 (ftable, Nlensing, NULL);
     460    if (!tmpLensing) {                                                 
     461      fprintf (stderr, "ERROR: failed to read lensings\n");             
     462      exit (2);                                                         
     463    }                                                                   
     464    lensing = Lensing_PS1_V5_R2_ToInternal (tmpLensing, *Nlensing);
     465    free (tmpLensing);                                                 
     466    *format = DVO_FORMAT_PS1_V5;                                       
     467    return (lensing); }
     468
     469# define SKIPPING_FORMAT(NAME, FORMAT, TYPE)            \
     470  if (!strcmp (extname, NAME)) {                                        \
     471    fprintf (stderr, "ERROR: format %s not defined for lensing, skipping\n", NAME); \
     472    *Nlensing = 0; \
     473    return NULL; \
     474  }
     475
    457476# define CONVERT_FORMAT(NAME, FORMAT, TYPE)             \
    458477  if (!strcmp (extname, NAME)) {                                        \
    … …  
    468487    return (lensing); }
    469488
    470 # define SKIPPING_FORMAT(NAME, FORMAT, TYPE)            \
    471   if (!strcmp (extname, NAME)) {                                        \
    472     fprintf (stderr, "ERROR: format %s not defined for lensing, skipping\n", NAME); \
    473     *Nlensing = 0; \
    474     return NULL; \
    475   }
    476 
    477489  if (!strcmp (extname, "DVO_LENSING")) {
    478490    lensing = gfits_table_get_Lensing (ftable, Nlensing, NULL);
    … …  
    496508  SKIPPING_FORMAT ("DVO_LENSING_PS1_V3",          PS1_V3,          PS1_V3);
    497509  SKIPPING_FORMAT ("DVO_LENSING_PS1_V4",          PS1_V4,          PS1_V4);
    498   CONVERT_FORMAT  ("DVO_LENSING_PS1_V5",          PS1_V5,          PS1_V5_R2);
     510  CONVERT_FORMAT  ("DVO_LENSING_PS1_V5",          PS1_V5,          PS1_V5_R3);
     511  SKIPPING_FORMAT ("DVO_LENSING_PS1_SIM",         PS1_SIM,         PS1_SIM);
    499512# undef CONVERT_FORMAT
     513# undef SKIPPING_FORMAT
    500514
    501515  fprintf (stderr, "table format unknown: %s\n", extname);
    … …  
    533547//    FORMAT_CASE (PS1_V3,          PS1_V3);
    534548//    FORMAT_CASE (PS1_V4,          PS1_V4);
    535       FORMAT_CASE (PS1_V5,          PS1_V5_R2);
     549      FORMAT_CASE (PS1_V5,          PS1_V5_R3);
    536550# undef FORMAT_CASE
    537551
    … …  
    554568    fprintf (stderr, "EXTNAME missing for lensobj table\n");
    555569    return (FALSE);
     570  }
     571
     572  if (!strcmp (extname, "DVO_LENSOBJ_PS1_V5") && (ftable[0].header[0].Naxis[0] == 136)) {
     573    Lensobj_PS1_V5_R0 *tmpLensobj;                                             
     574    tmpLensobj = gfits_table_get_Lensobj_PS1_V5_R0 (ftable, Nlensobj, NULL);
     575    if (!tmpLensobj) {                                                 
     576      fprintf (stderr, "ERROR: failed to read lensobjs\n");             
     577      exit (2);                                                         
     578    }                                                                   
     579    lensobj = Lensobj_PS1_V5_R0_ToInternal (tmpLensobj, *Nlensobj);
     580    free (tmpLensobj);                                                 
     581    *format = DVO_FORMAT_PS1_V5;                                       
     582    return (lensobj); }
     583
     584# define SKIPPING_FORMAT(NAME, FORMAT, TYPE)            \
     585  if (!strcmp (extname, NAME)) {                                        \
     586    fprintf (stderr, "ERROR: format %s not defined for lensobj, skipping\n", NAME); \
     587    *Nlensobj = 0; \
     588    return NULL; \
    556589  }
    557590
    … …  
    579612  }
    580613
    581 // CONVERT_FORMAT ("DVO_LENSOBJ_PS1_REF",         PS1_REF,         PS1_REF);
    582 // CONVERT_FORMAT ("DVO_LENSOBJ_ELIXIR",          ELIXIR,          Elixir);
    583 // CONVERT_FORMAT ("DVO_LENSOBJ_LONEOS",          LONEOS,          Loneos);
    584 // CONVERT_FORMAT ("DVO_LENSOBJ_PANSTARRS_DEV_0", PANSTARRS_DEV_0, Panstarrs_DEV_0);
    585 // CONVERT_FORMAT ("DVO_LENSOBJ_PANSTARRS_DEV_1", PANSTARRS_DEV_1, Panstarrs_DEV_1);
    586 // CONVERT_FORMAT ("DVO_LENSOBJ_PS1_DEV_1",       PS1_DEV_1,       PS1_DEV_1);
    587 // CONVERT_FORMAT ("DVO_LENSOBJ_PS1_DEV_2",       PS1_DEV_2,       PS1_DEV_2);
    588 // CONVERT_FORMAT ("DVO_LENSOBJ_PS1_V1",          PS1_V1,          PS1_V1);
    589 // CONVERT_FORMAT ("DVO_LENSOBJ_PS1_V2",          PS1_V2,          PS1_V2);
    590 // CONVERT_FORMAT ("DVO_LENSOBJ_PS1_V3",          PS1_V3,          PS1_V3);
    591 // CONVERT_FORMAT ("DVO_LENSOBJ_PS1_V4",          PS1_V4,          PS1_V4);
    592   CONVERT_FORMAT ("DVO_LENSOBJ_PS1_V5",          PS1_V5,          PS1_V5);
     614  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_REF",         PS1_REF,         PS1_REF);
     615  SKIPPING_FORMAT ("DVO_LENSOBJ_ELIXIR",          ELIXIR,          Elixir);
     616  SKIPPING_FORMAT ("DVO_LENSOBJ_LONEOS",          LONEOS,          Loneos);
     617  SKIPPING_FORMAT ("DVO_LENSOBJ_PANSTARRS_DEV_0", PANSTARRS_DEV_0, Panstarrs_DEV_0);
     618  SKIPPING_FORMAT ("DVO_LENSOBJ_PANSTARRS_DEV_1", PANSTARRS_DEV_1, Panstarrs_DEV_1);
     619  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_DEV_1",       PS1_DEV_1,       PS1_DEV_1);
     620  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_DEV_2",       PS1_DEV_2,       PS1_DEV_2);
     621  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_V1",          PS1_V1,          PS1_V1);
     622  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_V2",          PS1_V2,          PS1_V2);
     623  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_V3",          PS1_V3,          PS1_V3);
     624  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_V4",          PS1_V4,          PS1_V4);
     625  CONVERT_FORMAT  ("DVO_LENSOBJ_PS1_V5",          PS1_V5,          PS1_V5_R1);
     626  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_SIM",         PS1_SIM,         PS1_SIM);
    593627# undef CONVERT_FORMAT
     628# undef SKIPPING_FORMAT
    594629
    595630  fprintf (stderr, "table format unknown: %s\n", extname);
    … …  
    627662//    FORMAT_CASE (PS1_V3,          PS1_V3);
    628663//    FORMAT_CASE (PS1_V4,          PS1_V4);
    629       FORMAT_CASE (PS1_V5,          PS1_V5);
     664      FORMAT_CASE (PS1_V5,          PS1_V5_R1);
    630665# undef FORMAT_CASE
    631666
    … …  
    648683    fprintf (stderr, "EXTNAME missing for starpar table\n");
    649684    return (FALSE);
     685  }
     686
     687# define SKIPPING_FORMAT(NAME, FORMAT, TYPE)            \
     688  if (!strcmp (extname, NAME)) {                                        \
     689    fprintf (stderr, "ERROR: format %s not defined for starpar, skipping\n", NAME); \
     690    *Nstarpar = 0; \
     691    return NULL; \
    650692  }
    651693
    … …  
    673715  }
    674716
    675    CONVERT_FORMAT ("DVO_STARPAR_PS1_SIM",         PS1_SIM,         PS1_SIM);
    676 // CONVERT_FORMAT ("DVO_STARPAR_PS1_REF",         PS1_REF,         PS1_REF);
    677 // CONVERT_FORMAT ("DVO_STARPAR_ELIXIR",          ELIXIR,          Elixir);
    678 // CONVERT_FORMAT ("DVO_STARPAR_LONEOS",          LONEOS,          Loneos);
    679 // CONVERT_FORMAT ("DVO_STARPAR_PANSTARRS_DEV_0", PANSTARRS_DEV_0, Panstarrs_DEV_0);
    680 // CONVERT_FORMAT ("DVO_STARPAR_PANSTARRS_DEV_1", PANSTARRS_DEV_1, Panstarrs_DEV_1);
    681 // CONVERT_FORMAT ("DVO_STARPAR_PS1_DEV_1",       PS1_DEV_1,       PS1_DEV_1);
    682 // CONVERT_FORMAT ("DVO_STARPAR_PS1_DEV_2",       PS1_DEV_2,       PS1_DEV_2);
    683 // CONVERT_FORMAT ("DVO_STARPAR_PS1_V1",          PS1_V1,          PS1_V1);
    684 // CONVERT_FORMAT ("DVO_STARPAR_PS1_V2",          PS1_V2,          PS1_V2);
    685 // CONVERT_FORMAT ("DVO_STARPAR_PS1_V3",          PS1_V3,          PS1_V3);
    686 // CONVERT_FORMAT ("DVO_STARPAR_PS1_V4",          PS1_V4,          PS1_V4);
    687   CONVERT_FORMAT ("DVO_STARPAR_PS1_V5",          PS1_V5,          PS1_V5);
    688 //  CONVERT_FORMAT ("DVO_STELLAR_PARAMS_PS1_V5",          PS1_V5,          PS1_V5);
    689   // XXX fix in libauto: STELLAR_PARAMS to STARPAR
     717  CONVERT_FORMAT  ("DVO_STARPAR_PS1_SIM",         PS1_SIM,         PS1_SIM);
     718  SKIPPING_FORMAT ("DVO_STARPAR_PS1_REF",         PS1_REF,         PS1_REF);
     719  SKIPPING_FORMAT ("DVO_STARPAR_ELIXIR",          ELIXIR,          Elixir);
     720  SKIPPING_FORMAT ("DVO_STARPAR_LONEOS",          LONEOS,          Loneos);
     721  SKIPPING_FORMAT ("DVO_STARPAR_PANSTARRS_DEV_0", PANSTARRS_DEV_0, Panstarrs_DEV_0);
     722  SKIPPING_FORMAT ("DVO_STARPAR_PANSTARRS_DEV_1", PANSTARRS_DEV_1, Panstarrs_DEV_1);
     723  SKIPPING_FORMAT ("DVO_STARPAR_PS1_DEV_1",       PS1_DEV_1,       PS1_DEV_1);
     724  SKIPPING_FORMAT ("DVO_STARPAR_PS1_DEV_2",       PS1_DEV_2,       PS1_DEV_2);
     725  SKIPPING_FORMAT ("DVO_STARPAR_PS1_V1",          PS1_V1,          PS1_V1);
     726  SKIPPING_FORMAT ("DVO_STARPAR_PS1_V2",          PS1_V2,          PS1_V2);
     727  SKIPPING_FORMAT ("DVO_STARPAR_PS1_V3",          PS1_V3,          PS1_V3);
     728  SKIPPING_FORMAT ("DVO_STARPAR_PS1_V4",          PS1_V4,          PS1_V4);
     729  CONVERT_FORMAT  ("DVO_STARPAR_PS1_V5",          PS1_V5,          PS1_V5);
    690730# undef CONVERT_FORMAT
     731# undef SKIPPING_FORMAT
    691732
    692733  fprintf (stderr, "table format unknown: %s\n", extname);

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_V5.c

    r38062 r38153  
    448448}
    449449
    450 Lensing *Lensing_PS1_V5_R2_ToInternal (Lensing_PS1_V5_R2 *in, off_t Nvalues) {
     450Lensing *Lensing_PS1_V5_R0_ToInternal (Lensing_PS1_V5_R0 *in, off_t Nvalues) {
    451451
    452452  off_t i;
    … …  
    481481    out[i].E1_sh_psf   = in[i].E1_sh_psf;
    482482    out[i].E2_sh_psf   = in[i].E2_sh_psf;
    483 
    484     out[i].E1_psf      = in[i].E1_psf;
    485     out[i].E2_psf      = in[i].E2_psf;
    486483
    487484    out[i].F_ApR5      = in[i].F_ApR5;
    … …  
    499496    out[i].catID       = in[i].catID;
    500497    out[i].averef      = in[i].averef;
    501 
    502     out[i].imageID     = in[i].imageID;
    503   }
    504   return (out);
    505 }
    506 
    507 Lensing_PS1_V5_R2 *LensingInternalTo_PS1_V5_R2 (Lensing *in, off_t Nvalues) {
    508 
    509   off_t i;
    510   Lensing_PS1_V5_R2 *out;
    511 
    512   ALLOCATE_ZERO (out, Lensing_PS1_V5, Nvalues);
    513 
    514   for (i = 0; i < Nvalues; i++) {
     498  }
     499  return (out);
     500}
     501
     502Lensing *Lensing_PS1_V5_R1_ToInternal (Lensing_PS1_V5_R1 *in, off_t Nvalues) {
     503
     504  off_t i;
     505  Lensing *out;
     506
     507  ALLOCATE_ZERO (out, Lensing, Nvalues);
     508
     509  for (i = 0; i < Nvalues; i++) {
     510    dvo_lensing_init (&out[i]);
    515511
    516512    out[i].X11_sm_obj  = in[i].X11_sm_obj;     
    … …  
    537533    out[i].E1_sh_psf   = in[i].E1_sh_psf;
    538534    out[i].E2_sh_psf   = in[i].E2_sh_psf;
    539 
    540     out[i].E1_psf      = in[i].E1_psf;
    541     out[i].E2_psf      = in[i].E2_psf;
    542535
    543536    out[i].F_ApR5      = in[i].F_ApR5;
    … …  
    561554}
    562555
    563 Lensobj *Lensobj_PS1_V5_ToInternal (Lensobj_PS1_V5 *in, off_t Nvalues) {
     556Lensing *Lensing_PS1_V5_R2_ToInternal (Lensing_PS1_V5_R2 *in, off_t Nvalues) {
     557
     558  off_t i;
     559  Lensing *out;
     560
     561  ALLOCATE_ZERO (out, Lensing, Nvalues);
     562
     563  for (i = 0; i < Nvalues; i++) {
     564    dvo_lensing_init (&out[i]);
     565
     566    out[i].X11_sm_obj  = in[i].X11_sm_obj;     
     567    out[i].X12_sm_obj  = in[i].X12_sm_obj;     
     568    out[i].X22_sm_obj  = in[i].X22_sm_obj;     
     569    out[i].E1_sm_obj   = in[i].E1_sm_obj;     
     570    out[i].E2_sm_obj   = in[i].E2_sm_obj;     
     571                                         
     572    out[i].X11_sh_obj  = in[i].X11_sh_obj;     
     573    out[i].X12_sh_obj  = in[i].X12_sh_obj;     
     574    out[i].X22_sh_obj  = in[i].X22_sh_obj;     
     575    out[i].E1_sh_obj   = in[i].E1_sh_obj;     
     576    out[i].E2_sh_obj   = in[i].E2_sh_obj;     
     577                                         
     578    out[i].X11_sm_psf  = in[i].X11_sm_psf;
     579    out[i].X12_sm_psf  = in[i].X12_sm_psf;
     580    out[i].X22_sm_psf  = in[i].X22_sm_psf;
     581    out[i].E1_sm_psf   = in[i].E1_sm_psf;
     582    out[i].E2_sm_psf   = in[i].E2_sm_psf;     
     583                                         
     584    out[i].X11_sh_psf  = in[i].X11_sh_psf;
     585    out[i].X12_sh_psf  = in[i].X12_sh_psf;
     586    out[i].X22_sh_psf  = in[i].X22_sh_psf;
     587    out[i].E1_sh_psf   = in[i].E1_sh_psf;
     588    out[i].E2_sh_psf   = in[i].E2_sh_psf;
     589
     590    out[i].E1_psf      = in[i].E1_psf;
     591    out[i].E2_psf      = in[i].E2_psf;
     592
     593    out[i].F_ApR5      = in[i].F_ApR5;
     594    out[i].dF_ApR5     = in[i].dF_ApR5;
     595    out[i].sF_ApR5     = in[i].sF_ApR5;
     596    out[i].fF_ApR5     = in[i].fF_ApR5;
     597                                   
     598    out[i].F_ApR6      = in[i].F_ApR6;
     599    out[i].dF_ApR6     = in[i].dF_ApR6;
     600    out[i].sF_ApR6     = in[i].sF_ApR6;
     601    out[i].fF_ApR6     = in[i].fF_ApR6;
     602                                   
     603    out[i].detID       = in[i].detID;
     604    out[i].objID       = in[i].objID;
     605    out[i].catID       = in[i].catID;
     606    out[i].averef      = in[i].averef;
     607
     608    out[i].imageID     = in[i].imageID;
     609  }
     610  return (out);
     611}
     612
     613Lensing *Lensing_PS1_V5_R3_ToInternal (Lensing_PS1_V5_R3 *in, off_t Nvalues) {
     614
     615  off_t i;
     616  Lensing *out;
     617
     618  ALLOCATE_ZERO (out, Lensing, Nvalues);
     619
     620  for (i = 0; i < Nvalues; i++) {
     621    dvo_lensing_init (&out[i]);
     622
     623    out[i].X11_sm_obj  = in[i].X11_sm_obj;     
     624    out[i].X12_sm_obj  = in[i].X12_sm_obj;     
     625    out[i].X22_sm_obj  = in[i].X22_sm_obj;     
     626    out[i].E1_sm_obj   = in[i].E1_sm_obj;     
     627    out[i].E2_sm_obj   = in[i].E2_sm_obj;     
     628                                         
     629    out[i].X11_sh_obj  = in[i].X11_sh_obj;     
     630    out[i].X12_sh_obj  = in[i].X12_sh_obj;     
     631    out[i].X22_sh_obj  = in[i].X22_sh_obj;     
     632    out[i].E1_sh_obj   = in[i].E1_sh_obj;     
     633    out[i].E2_sh_obj   = in[i].E2_sh_obj;     
     634                                         
     635    out[i].X11_sm_psf  = in[i].X11_sm_psf;
     636    out[i].X12_sm_psf  = in[i].X12_sm_psf;
     637    out[i].X22_sm_psf  = in[i].X22_sm_psf;
     638    out[i].E1_sm_psf   = in[i].E1_sm_psf;
     639    out[i].E2_sm_psf   = in[i].E2_sm_psf;     
     640                                         
     641    out[i].X11_sh_psf  = in[i].X11_sh_psf;
     642    out[i].X12_sh_psf  = in[i].X12_sh_psf;
     643    out[i].X22_sh_psf  = in[i].X22_sh_psf;
     644    out[i].E1_sh_psf   = in[i].E1_sh_psf;
     645    out[i].E2_sh_psf   = in[i].E2_sh_psf;
     646
     647    out[i].E1_psf      = in[i].E1_psf;
     648    out[i].E2_psf      = in[i].E2_psf;
     649
     650    out[i].F_ApR5      = in[i].F_ApR5;
     651    out[i].dF_ApR5     = in[i].dF_ApR5;
     652    out[i].sF_ApR5     = in[i].sF_ApR5;
     653    out[i].fF_ApR5     = in[i].fF_ApR5;
     654                                   
     655    out[i].F_ApR6      = in[i].F_ApR6;
     656    out[i].dF_ApR6     = in[i].dF_ApR6;
     657    out[i].sF_ApR6     = in[i].sF_ApR6;
     658    out[i].fF_ApR6     = in[i].fF_ApR6;
     659                                   
     660    out[i].F_ApR7      = in[i].F_ApR7;
     661    out[i].dF_ApR7     = in[i].dF_ApR7;
     662    out[i].sF_ApR7     = in[i].sF_ApR7;
     663    out[i].fF_ApR7     = in[i].fF_ApR7;
     664                                   
     665    out[i].detID       = in[i].detID;
     666    out[i].objID       = in[i].objID;
     667    out[i].catID       = in[i].catID;
     668    out[i].averef      = in[i].averef;
     669
     670    out[i].imageID     = in[i].imageID;
     671  }
     672  return (out);
     673}
     674
     675Lensing_PS1_V5_R3 *LensingInternalTo_PS1_V5_R3 (Lensing *in, off_t Nvalues) {
     676
     677  off_t i;
     678  Lensing_PS1_V5_R3 *out;
     679
     680  ALLOCATE_ZERO (out, Lensing_PS1_V5_R3, Nvalues);
     681
     682  for (i = 0; i < Nvalues; i++) {
     683
     684    out[i].X11_sm_obj  = in[i].X11_sm_obj;     
     685    out[i].X12_sm_obj  = in[i].X12_sm_obj;     
     686    out[i].X22_sm_obj  = in[i].X22_sm_obj;     
     687    out[i].E1_sm_obj   = in[i].E1_sm_obj;     
     688    out[i].E2_sm_obj   = in[i].E2_sm_obj;     
     689                                         
     690    out[i].X11_sh_obj  = in[i].X11_sh_obj;     
     691    out[i].X12_sh_obj  = in[i].X12_sh_obj;     
     692    out[i].X22_sh_obj  = in[i].X22_sh_obj;     
     693    out[i].E1_sh_obj   = in[i].E1_sh_obj;     
     694    out[i].E2_sh_obj   = in[i].E2_sh_obj;     
     695                                         
     696    out[i].X11_sm_psf  = in[i].X11_sm_psf;
     697    out[i].X12_sm_psf  = in[i].X12_sm_psf;
     698    out[i].X22_sm_psf  = in[i].X22_sm_psf;
     699    out[i].E1_sm_psf   = in[i].E1_sm_psf;
     700    out[i].E2_sm_psf   = in[i].E2_sm_psf;     
     701                                         
     702    out[i].X11_sh_psf  = in[i].X11_sh_psf;
     703    out[i].X12_sh_psf  = in[i].X12_sh_psf;
     704    out[i].X22_sh_psf  = in[i].X22_sh_psf;
     705    out[i].E1_sh_psf   = in[i].E1_sh_psf;
     706    out[i].E2_sh_psf   = in[i].E2_sh_psf;
     707
     708    out[i].E1_psf      = in[i].E1_psf;
     709    out[i].E2_psf      = in[i].E2_psf;
     710
     711    out[i].F_ApR5      = in[i].F_ApR5;
     712    out[i].dF_ApR5     = in[i].dF_ApR5;
     713    out[i].sF_ApR5     = in[i].sF_ApR5;
     714    out[i].fF_ApR5     = in[i].fF_ApR5;
     715                                   
     716    out[i].F_ApR6      = in[i].F_ApR6;
     717    out[i].dF_ApR6     = in[i].dF_ApR6;
     718    out[i].sF_ApR6     = in[i].sF_ApR6;
     719    out[i].fF_ApR6     = in[i].fF_ApR6;
     720                                   
     721    out[i].F_ApR7      = in[i].F_ApR7;
     722    out[i].dF_ApR7     = in[i].dF_ApR7;
     723    out[i].sF_ApR7     = in[i].sF_ApR7;
     724    out[i].fF_ApR7     = in[i].fF_ApR7;
     725                                   
     726    out[i].detID       = in[i].detID;
     727    out[i].objID       = in[i].objID;
     728    out[i].catID       = in[i].catID;
     729    out[i].averef      = in[i].averef;
     730
     731    out[i].imageID     = in[i].imageID;
     732  }
     733  return (out);
     734}
     735
     736Lensobj *Lensobj_PS1_V5_R0_ToInternal (Lensobj_PS1_V5_R0 *in, off_t Nvalues) {
    564737
    565738  off_t i;
    … …  
    618791}
    619792
    620 Lensobj_PS1_V5 *LensobjInternalTo_PS1_V5 (Lensobj *in, off_t Nvalues) {
    621 
    622   off_t i;
    623   Lensobj_PS1_V5 *out;
    624 
    625   ALLOCATE_ZERO (out, Lensobj_PS1_V5, Nvalues);
    626 
    627   for (i = 0; i < Nvalues; i++) {
     793Lensobj *Lensobj_PS1_V5_R1_ToInternal (Lensobj_PS1_V5_R1 *in, off_t Nvalues) {
     794
     795  off_t i;
     796  Lensobj *out;
     797
     798  ALLOCATE_ZERO (out, Lensobj, Nvalues);
     799
     800  for (i = 0; i < Nvalues; i++) {
     801    dvo_lensobj_init (&out[i], FALSE);
    628802
    629803    out[i].X11_sm_obj  = in[i].X11_sm_obj;     
    … …  
    660834    out[i].sF_ApR6     = in[i].sF_ApR6;
    661835    out[i].fF_ApR6     = in[i].fF_ApR6;
     836                                   
     837    out[i].F_ApR7      = in[i].F_ApR7;
     838    out[i].dF_ApR7     = in[i].dF_ApR7;
     839    out[i].sF_ApR7     = in[i].sF_ApR7;
     840    out[i].fF_ApR7     = in[i].fF_ApR7;
     841                                   
     842    out[i].gamma       = in[i].gamma;
     843    out[i].E1          = in[i].E1;
     844    out[i].E2          = in[i].E2;
     845
     846    out[i].objID       = in[i].objID;
     847    out[i].catID       = in[i].catID;
     848
     849    out[i].photcode    = in[i].photcode;
     850    out[i].Nmeas       = in[i].Nmeas;
     851  }
     852  return (out);
     853}
     854
     855Lensobj_PS1_V5_R1 *LensobjInternalTo_PS1_V5_R1 (Lensobj *in, off_t Nvalues) {
     856
     857  off_t i;
     858  Lensobj_PS1_V5_R1 *out;
     859
     860  ALLOCATE_ZERO (out, Lensobj_PS1_V5_R1, Nvalues);
     861
     862  for (i = 0; i < Nvalues; i++) {
     863
     864    out[i].X11_sm_obj  = in[i].X11_sm_obj;     
     865    out[i].X12_sm_obj  = in[i].X12_sm_obj;     
     866    out[i].X22_sm_obj  = in[i].X22_sm_obj;     
     867    out[i].E1_sm_obj   = in[i].E1_sm_obj;     
     868    out[i].E2_sm_obj   = in[i].E2_sm_obj;     
     869                                         
     870    out[i].X11_sh_obj  = in[i].X11_sh_obj;     
     871    out[i].X12_sh_obj  = in[i].X12_sh_obj;     
     872    out[i].X22_sh_obj  = in[i].X22_sh_obj;     
     873    out[i].E1_sh_obj   = in[i].E1_sh_obj;     
     874    out[i].E2_sh_obj   = in[i].E2_sh_obj;     
     875                                         
     876    out[i].X11_sm_psf  = in[i].X11_sm_psf;
     877    out[i].X12_sm_psf  = in[i].X12_sm_psf;
     878    out[i].X22_sm_psf  = in[i].X22_sm_psf;
     879    out[i].E1_sm_psf   = in[i].E1_sm_psf;
     880    out[i].E2_sm_psf   = in[i].E2_sm_psf;     
     881                                         
     882    out[i].X11_sh_psf  = in[i].X11_sh_psf;
     883    out[i].X12_sh_psf  = in[i].X12_sh_psf;
     884    out[i].X22_sh_psf  = in[i].X22_sh_psf;
     885    out[i].E1_sh_psf   = in[i].E1_sh_psf;
     886    out[i].E2_sh_psf   = in[i].E2_sh_psf;
     887
     888    out[i].F_ApR5      = in[i].F_ApR5;
     889    out[i].dF_ApR5     = in[i].dF_ApR5;
     890    out[i].sF_ApR5     = in[i].sF_ApR5;
     891    out[i].fF_ApR5     = in[i].fF_ApR5;
     892                                   
     893    out[i].F_ApR6      = in[i].F_ApR6;
     894    out[i].dF_ApR6     = in[i].dF_ApR6;
     895    out[i].sF_ApR6     = in[i].sF_ApR6;
     896    out[i].fF_ApR6     = in[i].fF_ApR6;
     897                                   
     898    out[i].F_ApR7      = in[i].F_ApR7;
     899    out[i].dF_ApR7     = in[i].dF_ApR7;
     900    out[i].sF_ApR7     = in[i].sF_ApR7;
     901    out[i].fF_ApR7     = in[i].fF_ApR7;
    662902                                   
    663903    out[i].gamma       = in[i].gamma;
    … …  
    12711511  }
    12721512  return (data);
    1273 }
    1274 
    1275 Lensing *Lensing_PS1_V5_R0_ToInternal (Lensing_PS1_V5_R0 *in, off_t Nvalues) {
    1276 
    1277   off_t i;
    1278   Lensing *out;
    1279 
    1280   ALLOCATE_ZERO (out, Lensing, Nvalues);
    1281 
    1282   for (i = 0; i < Nvalues; i++) {
    1283     dvo_lensing_init (&out[i]);
    1284 
    1285     out[i].X11_sm_obj  = in[i].X11_sm_obj;     
    1286     out[i].X12_sm_obj  = in[i].X12_sm_obj;     
    1287     out[i].X22_sm_obj  = in[i].X22_sm_obj;     
    1288     out[i].E1_sm_obj   = in[i].E1_sm_obj;     
    1289     out[i].E2_sm_obj   = in[i].E2_sm_obj;     
    1290                                          
    1291     out[i].X11_sh_obj  = in[i].X11_sh_obj;     
    1292     out[i].X12_sh_obj  = in[i].X12_sh_obj;     
    1293     out[i].X22_sh_obj  = in[i].X22_sh_obj;     
    1294     out[i].E1_sh_obj   = in[i].E1_sh_obj;     
    1295     out[i].E2_sh_obj   = in[i].E2_sh_obj;     
    1296                                          
    1297     out[i].X11_sm_psf  = in[i].X11_sm_psf;
    1298     out[i].X12_sm_psf  = in[i].X12_sm_psf;
    1299     out[i].X22_sm_psf  = in[i].X22_sm_psf;
    1300     out[i].E1_sm_psf   = in[i].E1_sm_psf;
    1301     out[i].E2_sm_psf   = in[i].E2_sm_psf;     
    1302                                          
    1303     out[i].X11_sh_psf  = in[i].X11_sh_psf;
    1304     out[i].X12_sh_psf  = in[i].X12_sh_psf;
    1305     out[i].X22_sh_psf  = in[i].X22_sh_psf;
    1306     out[i].E1_sh_psf   = in[i].E1_sh_psf;
    1307     out[i].E2_sh_psf   = in[i].E2_sh_psf;
    1308 
    1309     out[i].F_ApR5      = in[i].F_ApR5;
    1310     out[i].dF_ApR5     = in[i].dF_ApR5;
    1311     out[i].sF_ApR5     = in[i].sF_ApR5;
    1312     out[i].fF_ApR5     = in[i].fF_ApR5;
    1313                                    
    1314     out[i].F_ApR6      = in[i].F_ApR6;
    1315     out[i].dF_ApR6     = in[i].dF_ApR6;
    1316     out[i].sF_ApR6     = in[i].sF_ApR6;
    1317     out[i].fF_ApR6     = in[i].fF_ApR6;
    1318                                    
    1319     out[i].detID       = in[i].detID;
    1320     out[i].objID       = in[i].objID;
    1321     out[i].catID       = in[i].catID;
    1322     out[i].averef      = in[i].averef;
    1323   }
    1324   return (out);
    1325 }
    1326 
    1327 Lensing *Lensing_PS1_V5_R1_ToInternal (Lensing_PS1_V5_R1 *in, off_t Nvalues) {
    1328 
    1329   off_t i;
    1330   Lensing *out;
    1331 
    1332   ALLOCATE_ZERO (out, Lensing, Nvalues);
    1333 
    1334   for (i = 0; i < Nvalues; i++) {
    1335     dvo_lensing_init (&out[i]);
    1336 
    1337     out[i].X11_sm_obj  = in[i].X11_sm_obj;     
    1338     out[i].X12_sm_obj  = in[i].X12_sm_obj;     
    1339     out[i].X22_sm_obj  = in[i].X22_sm_obj;     
    1340     out[i].E1_sm_obj   = in[i].E1_sm_obj;     
    1341     out[i].E2_sm_obj   = in[i].E2_sm_obj;     
    1342                                          
    1343     out[i].X11_sh_obj  = in[i].X11_sh_obj;     
    1344     out[i].X12_sh_obj  = in[i].X12_sh_obj;     
    1345     out[i].X22_sh_obj  = in[i].X22_sh_obj;     
    1346     out[i].E1_sh_obj   = in[i].E1_sh_obj;     
    1347     out[i].E2_sh_obj   = in[i].E2_sh_obj;     
    1348                                          
    1349     out[i].X11_sm_psf  = in[i].X11_sm_psf;
    1350     out[i].X12_sm_psf  = in[i].X12_sm_psf;
    1351     out[i].X22_sm_psf  = in[i].X22_sm_psf;
    1352     out[i].E1_sm_psf   = in[i].E1_sm_psf;
    1353     out[i].E2_sm_psf   = in[i].E2_sm_psf;     
    1354                                          
    1355     out[i].X11_sh_psf  = in[i].X11_sh_psf;
    1356     out[i].X12_sh_psf  = in[i].X12_sh_psf;
    1357     out[i].X22_sh_psf  = in[i].X22_sh_psf;
    1358     out[i].E1_sh_psf   = in[i].E1_sh_psf;
    1359     out[i].E2_sh_psf   = in[i].E2_sh_psf;
    1360 
    1361     out[i].F_ApR5      = in[i].F_ApR5;
    1362     out[i].dF_ApR5     = in[i].dF_ApR5;
    1363     out[i].sF_ApR5     = in[i].sF_ApR5;
    1364     out[i].fF_ApR5     = in[i].fF_ApR5;
    1365                                    
    1366     out[i].F_ApR6      = in[i].F_ApR6;
    1367     out[i].dF_ApR6     = in[i].dF_ApR6;
    1368     out[i].sF_ApR6     = in[i].sF_ApR6;
    1369     out[i].fF_ApR6     = in[i].fF_ApR6;
    1370                                    
    1371     out[i].detID       = in[i].detID;
    1372     out[i].objID       = in[i].objID;
    1373     out[i].catID       = in[i].catID;
    1374     out[i].averef      = in[i].averef;
    1375 
    1376     out[i].imageID     = in[i].imageID;
    1377   }
    1378   return (out);
    13791513}
    13801514

  • trunk/Ohana/src/libdvo/src/dvo_photcode_ops.c

    r37807 r38153  
    18221822LENSFIELD(E1_sh_psf);
    18231823LENSFIELD(E2_sh_psf);
     1824
     1825LENSFIELD( F_ApR5);
     1826LENSFIELD(dF_ApR5);
     1827LENSFIELD(sF_ApR5);
     1828LENSFIELD(fF_ApR5);
     1829
     1830LENSFIELD( F_ApR6);
     1831LENSFIELD(dF_ApR6);
     1832LENSFIELD(sF_ApR6);
     1833LENSFIELD(fF_ApR6);
     1834
     1835LENSFIELD( F_ApR7);
     1836LENSFIELD(dF_ApR7);
     1837LENSFIELD(sF_ApR7);
     1838LENSFIELD(fF_ApR7);
     1839
     1840LENSFIELD(E1);
     1841LENSFIELD(E2);
Note: See TracChangeset for help on using the changeset viewer.