IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Apr 9, 2015, 8:43:14 AM (11 years ago)
Author:
eugene
Message:

add lensing ps1-v5-r3, lensobj ps1-v5-r1 (adds R7 apertures)

Location:
branches/eam_branches/ipp-20150405/Ohana/src
Files:
2 added
6 edited
1 moved

Legend:

Unmodified
Added
Removed

  • branches/eam_branches/ipp-20150405/Ohana/src/libautocode/Makefile.Targets

    r37995 r38100  
    4848$(ASRC)/lensing-ps1-v5-r1.$(ARCH).o \
    4949$(ASRC)/lensing-ps1-v5-r2.$(ARCH).o \
     50$(ASRC)/lensing-ps1-v5-r3.$(ARCH).o \
    5051$(ASRC)/lensobj.$(ARCH).o \
    51 $(ASRC)/lensobj-ps1-v5.$(ARCH).o \
     52$(ASRC)/lensobj-ps1-v5-r0.$(ARCH).o \
     53$(ASRC)/lensobj-ps1-v5-r1.$(ARCH).o \
    5254$(ASRC)/starpar.$(ARCH).o \
    5355$(ASRC)/starpar-ps1-v5.$(ARCH).o \
    … …  
    156158$(AINC)/lensing-ps1-v5-r1.h \
    157159$(AINC)/lensing-ps1-v5-r2.h \
     160$(AINC)/lensing-ps1-v5-r3.h \
    158161$(AINC)/lensobj.h \
    159 $(AINC)/lensobj-ps1-v5.h \
     162$(AINC)/lensobj-ps1-v5-r0.h \
     163$(AINC)/lensobj-ps1-v5-r1.h \
    160164$(AINC)/starpar.h \
    161165$(AINC)/starpar-ps1-v5.h \

  • branches/eam_branches/ipp-20150405/Ohana/src/libautocode/def/lensing.d

    r37992 r38100  
    22EXTNAME      DVO_LENSING
    33TYPE         BINTABLE
    4 SIZE         144
     4SIZE         160
    55DESCRIPTION  DVO Lensing Table
    66
    … …  
    3232FIELD  E2_psf,         E2_PSF,       float,          polarization psf stars E2
    3333
    34 FIELD F_ApR5,        FLUX_AP_R5,    float,          Flux inside r = 5
     34FIELD  F_ApR5,           FLUX_AP_R5, float,          Flux inside r = 5
    3535FIELD dF_ApR5,       FLUX_ERR_AP_R5, float,          error on Flux inside r = 5
    3636FIELD sF_ApR5,       FLUX_STD_AP_R5, float,          stdev of Flux inside r = 5
    3737FIELD fF_ApR5,       FLUX_FIL_AP_R5, float,          fill factor for Flux inside r = 5
    3838
    39 FIELD F_ApR6,        FLUX_AP_R6,     float,          Flux inside r = 5
    40 FIELD dF_ApR6,       FLUX_ERR_AP_R6, float,          error on Flux inside r = 5
    41 FIELD sF_ApR6,       FLUX_STD_AP_R6, float,          stdev of Flux inside r = 5
    42 FIELD fF_ApR6,       FLUX_FIL_AP_R6, float,          fill factor for Flux inside r = 5
     39FIELD  F_ApR6,           FLUX_AP_R6, float,          Flux inside r = 6
     40FIELD dF_ApR6,       FLUX_ERR_AP_R6, float,          error on Flux inside r = 6
     41FIELD sF_ApR6,       FLUX_STD_AP_R6, float,          stdev of Flux inside r = 6
     42FIELD fF_ApR6,       FLUX_FIL_AP_R6, float,          fill factor for Flux inside r = 6
     43
     44FIELD  F_ApR7,           FLUX_AP_R7, float,          Flux inside r = 7
     45FIELD dF_ApR7,       FLUX_ERR_AP_R7, float,          error on Flux inside r = 7
     46FIELD sF_ApR7,       FLUX_STD_AP_R7, float,          stdev of Flux inside r = 7
     47FIELD fF_ApR7,       FLUX_FIL_AP_R7, float,          fill factor for Flux inside r = 7
    4348
    4449FIELD   detID,          DET_ID,      int,            detection ID
    … …  
    5055FIELD padding,         PADDING,      int,            padding
    5156
    52 # 28 x float
    53 # 5 x int + padding
    54 # = 136 bytes / detection
     57# 160 bytes / detection
    5558# for 2.5G objects + 60 overlaps, expect 15TB for forced warp

  • branches/eam_branches/ipp-20150405/Ohana/src/libautocode/def/lensobj-ps1-v5-r0.d

    r38099 r38100  
    1 STRUCT       Lensobj_PS1_V5
     1STRUCT       Lensobj_PS1_V5_R0
    22EXTNAME      DVO_LENSOBJ_PS1_V5
    33TYPE         BINTABLE
    … …  
    2929FIELD  E2_sh_psf,      E2_SH_PSF,    float,          mean lensing shear psf stars E2
    3030
    31 FIELD F_ApR5,        FLUX_AP_R5,    float,          Flux inside r = 5
     31FIELD  F_ApR5,           FLUX_AP_R5, float,          Flux inside r = 5
    3232FIELD dF_ApR5,       FLUX_ERR_AP_R5, float,          error on Flux inside r = 5
    3333FIELD sF_ApR5,       FLUX_STD_AP_R5, float,          stdev of Flux inside r = 5
    3434FIELD fF_ApR5,       FLUX_FIL_AP_R5, float,          fill factor for Flux inside r = 5
    3535
    36 FIELD F_ApR6,        FLUX_AP_R6,     float,          Flux inside r = 5
    37 FIELD dF_ApR6,       FLUX_ERR_AP_R6, float,          error on Flux inside r = 5
    38 FIELD sF_ApR6,       FLUX_STD_AP_R6, float,          stdev of Flux inside r = 5
    39 FIELD fF_ApR6,       FLUX_FIL_AP_R6, float,          fill factor for Flux inside r = 5
     36FIELD  F_ApR6,           FLUX_AP_R6, float,          Flux inside r = 6
     37FIELD dF_ApR6,       FLUX_ERR_AP_R6, float,          error on Flux inside r = 6
     38FIELD sF_ApR6,       FLUX_STD_AP_R6, float,          stdev of Flux inside r = 6
     39FIELD fF_ApR6,       FLUX_FIL_AP_R6, float,          fill factor for Flux inside r = 6
    4040
    4141FIELD  gamma,            GAMMA,      float

  • branches/eam_branches/ipp-20150405/Ohana/src/libautocode/def/lensobj.d

    r37807 r38100  
    22EXTNAME      DVO_LENSOBJ
    33TYPE         BINTABLE
    4 SIZE         136
     4SIZE         152
    55DESCRIPTION  DVO Lensobj Table
    66
    … …  
    2929FIELD  E2_sh_psf,      E2_SH_PSF,    float,          mean lensing shear psf stars E2
    3030
    31 FIELD F_ApR5,        FLUX_AP_R5,    float,          Flux inside r = 5
     31FIELD  F_ApR5,           FLUX_AP_R5, float,          Flux inside r = 5
    3232FIELD dF_ApR5,       FLUX_ERR_AP_R5, float,          error on Flux inside r = 5
    3333FIELD sF_ApR5,       FLUX_STD_AP_R5, float,          stdev of Flux inside r = 5
    3434FIELD fF_ApR5,       FLUX_FIL_AP_R5, float,          fill factor for Flux inside r = 5
    3535
    36 FIELD F_ApR6,        FLUX_AP_R6,     float,          Flux inside r = 5
    37 FIELD dF_ApR6,       FLUX_ERR_AP_R6, float,          error on Flux inside r = 5
    38 FIELD sF_ApR6,       FLUX_STD_AP_R6, float,          stdev of Flux inside r = 5
    39 FIELD fF_ApR6,       FLUX_FIL_AP_R6, float,          fill factor for Flux inside r = 5
     36FIELD  F_ApR6,           FLUX_AP_R6, float,          Flux inside r = 6
     37FIELD dF_ApR6,       FLUX_ERR_AP_R6, float,          error on Flux inside r = 6
     38FIELD sF_ApR6,       FLUX_STD_AP_R6, float,          stdev of Flux inside r = 6
     39FIELD fF_ApR6,       FLUX_FIL_AP_R6, float,          fill factor for Flux inside r = 6
     40
     41FIELD  F_ApR7,           FLUX_AP_R7, float,          Flux inside r = 7
     42FIELD dF_ApR7,       FLUX_ERR_AP_R7, float,          error on Flux inside r = 7
     43FIELD sF_ApR7,       FLUX_STD_AP_R7, float,          stdev of Flux inside r = 7
     44FIELD fF_ApR7,       FLUX_FIL_AP_R7, float,          fill factor for Flux inside r = 7
    4045
    4146FIELD      gamma,        GAMMA,      float

  • branches/eam_branches/ipp-20150405/Ohana/src/libdvo/include/ps1_v5_defs.h

    r37995 r38100  
    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);

  • branches/eam_branches/ipp-20150405/Ohana/src/libdvo/src/dvo_convert.c

    r38098 r38100  
    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);
    499511  SKIPPING_FORMAT ("DVO_LENSING_PS1_SIM",         PS1_SIM,         PS1_SIM);
    500512# undef CONVERT_FORMAT
    … …  
    535547//    FORMAT_CASE (PS1_V3,          PS1_V3);
    536548//    FORMAT_CASE (PS1_V4,          PS1_V4);
    537       FORMAT_CASE (PS1_V5,          PS1_V5_R2);
     549      FORMAT_CASE (PS1_V5,          PS1_V5_R3);
    538550# undef FORMAT_CASE
    539551
    … …  
    557569    return (FALSE);
    558570  }
     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); }
    559583
    560584# define SKIPPING_FORMAT(NAME, FORMAT, TYPE)            \
    … …  
    599623  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_V3",          PS1_V3,          PS1_V3);
    600624  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_V4",          PS1_V4,          PS1_V4);
    601   CONVERT_FORMAT  ("DVO_LENSOBJ_PS1_V5",          PS1_V5,          PS1_V5);
     625  CONVERT_FORMAT  ("DVO_LENSOBJ_PS1_V5",          PS1_V5,          PS1_V5_R1);
    602626  SKIPPING_FORMAT ("DVO_LENSOBJ_PS1_SIM",         PS1_SIM,         PS1_SIM);
    603627# undef CONVERT_FORMAT
    … …  
    638662//    FORMAT_CASE (PS1_V3,          PS1_V3);
    639663//    FORMAT_CASE (PS1_V4,          PS1_V4);
    640       FORMAT_CASE (PS1_V5,          PS1_V5);
     664      FORMAT_CASE (PS1_V5,          PS1_V5_R1);
    641665# undef FORMAT_CASE
    642666

  • branches/eam_branches/ipp-20150405/Ohana/src/libdvo/src/dvo_convert_PS1_V5.c

    r38062 r38100  
    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
Note: See TracChangeset for help on using the changeset viewer.