diff options
Diffstat (limited to 'noao/digiphot/daophot/allstar/dpalmemstar.x')
| -rw-r--r-- | noao/digiphot/daophot/allstar/dpalmemstar.x | 182 |
1 files changed, 182 insertions, 0 deletions
diff --git a/noao/digiphot/daophot/allstar/dpalmemstar.x b/noao/digiphot/daophot/allstar/dpalmemstar.x new file mode 100644 index 00000000..fd9ee80e --- /dev/null +++ b/noao/digiphot/daophot/allstar/dpalmemstar.x @@ -0,0 +1,182 @@ +include "../lib/daophotdef.h" +include "../lib/allstardef.h" + +# DP_ALLSTARSETUP -- Procedure to set up the ALLSTAR parameters. + +procedure dp_allstarsetup (dp) + +pointer dp # pointer to daophot structure + +pointer allstar + +begin + # Allocate memory. + call malloc (DP_ALLSTAR(dp), LEN_ALLSTARSTRUCT, TY_STRUCT) + allstar = DP_ALLSTAR(dp) + + DP_DATA(allstar) = NULL + DP_SUBT(allstar) = NULL + DP_WEIGHTS(allstar) = NULL + DP_DBUF(allstar) = NULL + DP_SBUF(allstar) = NULL + DP_WBUF(allstar) = NULL + + DP_ANUMER(allstar) = NULL + DP_ADENOM(allstar) = NULL + DP_ASUMWT(allstar) = NULL + DP_ARPIXSQ(allstar) = NULL + DP_ASKIP(allstar) = NULL + DP_AXCLAMP(allstar) = NULL + DP_AYCLAMP(allstar) = NULL + DP_AXOLD(allstar) = NULL + DP_AYOLD(allstar) = NULL + DP_AX(allstar) = NULL + DP_AV(allstar) = NULL + DP_AC(allstar) = NULL + DP_ALAST(allstar) = NULL + DP_ANPIX(allstar) = NULL + DP_AIER(allstar) = NULL +end + + +# DP_ALMEMSTAR -- Procedure to allocate sufficient memory for ALLSTAR. + +procedure dp_almemstar (dao, max_star, max_group) + +pointer dao # pointer to daophot structure +int max_star # maximum number of stars +int max_group # maximum group size + +pointer allstar + +begin + allstar = DP_ALLSTAR(dao) + + if (DP_ANUMER (allstar) != NULL) + call mfree (DP_ANUMER (allstar), TY_REAL) + call malloc (DP_ANUMER (allstar), max_star, TY_REAL) + + if (DP_ADENOM (allstar) != NULL) + call mfree (DP_ADENOM (allstar), TY_REAL) + call malloc (DP_ADENOM (allstar), max_star, TY_REAL) + + if (DP_ARPIXSQ (allstar) != NULL) + call mfree (DP_ARPIXSQ (allstar), TY_REAL) + call malloc (DP_ARPIXSQ (allstar), max_star, TY_REAL) + + if (DP_ASUMWT (allstar) != NULL) + call mfree (DP_ASUMWT (allstar), TY_REAL) + call malloc (DP_ASUMWT (allstar), max_star, TY_REAL) + + if (DP_ANPIX (allstar) != NULL) + call mfree (DP_ANPIX (allstar), TY_INT) + call malloc (DP_ANPIX (allstar), max_star, TY_INT) + + if (DP_AIER (allstar) != NULL) + call mfree (DP_AIER (allstar), TY_INT) + call malloc (DP_AIER (allstar), max_star, TY_INT) + + if (DP_ASKIP (allstar) != NULL) + call mfree (DP_ASKIP (allstar), TY_INT) + call malloc (DP_ASKIP (allstar), max_star, TY_INT) + + if (DP_ALAST (allstar) != NULL) + call mfree (DP_ALAST (allstar), TY_INT) + call malloc (DP_ALAST (allstar), max_star, TY_INT) + + if (DP_AXCLAMP (allstar) != NULL) + call mfree (DP_AXCLAMP (allstar), TY_REAL) + call malloc (DP_AXCLAMP (allstar), max_star, TY_REAL) + + if (DP_AYCLAMP (allstar) != NULL) + call mfree (DP_AYCLAMP (allstar), TY_REAL) + call malloc (DP_AYCLAMP (allstar), max_star, TY_REAL) + + if (DP_AXOLD (allstar) != NULL) + call mfree (DP_AXOLD (allstar), TY_REAL) + call malloc (DP_AXOLD (allstar), max_star, TY_REAL) + + if (DP_AYOLD (allstar) != NULL) + call mfree (DP_AYOLD (allstar), TY_REAL) + call malloc (DP_AYOLD (allstar), max_star, TY_REAL) + + # Allocate space for the fitting matrices. Note that nine + # times less space is required if recentering is turned + # off. + if (DP_RECENTER(dao) == YES) { + + if (DP_AX (allstar) != NULL) + call mfree (DP_AX (allstar), TY_REAL) + call malloc (DP_AX (allstar), 3 * max_group + 1, TY_REAL) + + if (DP_AV (allstar) != NULL) + call mfree (DP_AV (allstar), TY_REAL) + call malloc (DP_AV (allstar), 3 * max_group + 1, TY_REAL) + + if (DP_AC (allstar) != NULL) + call mfree (DP_AC (allstar), TY_REAL) + call malloc (DP_AC (allstar), (3 * max_group + 1) * + (3 * max_group + 1), TY_REAL) + + } else { + + if (DP_AX (allstar) != NULL) + call mfree (DP_AX (allstar), TY_REAL) + call malloc (DP_AX (allstar), max_group + 1, TY_REAL) + + if (DP_AV (allstar) != NULL) + call mfree (DP_AV (allstar), TY_REAL) + call malloc (DP_AV (allstar), max_group + 1, TY_REAL) + + if (DP_AC (allstar) != NULL) + call mfree (DP_AC (allstar), TY_REAL) + call malloc (DP_AC (allstar), (max_group + 1) * (max_group + 1), + TY_REAL) + } +end + + +# DP_ALCLOSE -- Procedure to close up the ALLSTAR parameters. + +procedure dp_alclose (dp) + +pointer dp # pointer to daophot structure + +pointer allstar + +begin + allstar = DP_ALLSTAR(dp) + + if (DP_ANUMER (allstar) != NULL) + call mfree (DP_ANUMER (allstar), TY_REAL) + if (DP_ADENOM (allstar) != NULL) + call mfree (DP_ADENOM (allstar), TY_REAL) + if (DP_ARPIXSQ (allstar) != NULL) + call mfree (DP_ARPIXSQ (allstar), TY_REAL) + if (DP_ASUMWT (allstar) != NULL) + call mfree (DP_ASUMWT (allstar), TY_REAL) + if (DP_ASKIP (allstar) != NULL) + call mfree (DP_ASKIP (allstar), TY_INT) + if (DP_ALAST (allstar) != NULL) + call mfree (DP_ALAST (allstar), TY_INT) + if (DP_AXCLAMP (allstar) != NULL) + call mfree (DP_AXCLAMP (allstar), TY_REAL) + if (DP_AYCLAMP (allstar) != NULL) + call mfree (DP_AYCLAMP (allstar), TY_REAL) + if (DP_AXOLD (allstar) != NULL) + call mfree (DP_AXOLD (allstar), TY_REAL) + if (DP_AYOLD (allstar) != NULL) + call mfree (DP_AYOLD (allstar), TY_REAL) + if (DP_AX (allstar) != NULL) + call mfree (DP_AX (allstar), TY_REAL) + if (DP_AV (allstar) != NULL) + call mfree (DP_AV (allstar), TY_REAL) + if (DP_AC (allstar) != NULL) + call mfree (DP_AC (allstar), TY_REAL) + if (DP_ANPIX (allstar) != NULL) + call mfree (DP_ANPIX (allstar), TY_INT) + if (DP_AIER (allstar) != NULL) + call mfree (DP_AIER (allstar), TY_INT) + + call mfree (allstar, TY_STRUCT) +end |