IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Jun 2, 2006, 3:02:08 PM (20 years ago)
Author:
Paul Price
Message:

Applying RHL patch. Generally improves the error handling and traceback. pmConcepts.c and pmConceptsRead.c done by PAP (RHL did this also, but I had already done them). Resolved conflicts, except for pmFPAfile.c, which uses psAbort in some instances where RHL's patch had psError; leave this for Gene to decide.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/psModules/src/objects/pmPSF.c

    r6960 r7311  
    66 *  @author EAM, IfA
    77 *
    8  *  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
    9  *  @date $Date: 2006-04-23 01:07:51 $
     8 *  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
     9 *  @date $Date: 2006-06-03 01:02:08 $
    1010 *
    1111 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
    … …  
    9797    // the ApTrend components are (x, y, r2rflux, flux)
    9898    psf->ApTrend = psPolynomial4DAlloc (PS_POLYNOMIAL_ORD, 2, 2, 1, 1);
    99     pmPSF_MaskApTrend (psf, PM_PSF_SKYBIAS);
     99    pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_SKYBIAS);
    100100
    101101    if (psf->poissonErrors) {
    … …  
    239239 * to enable the specific subset of the coefficients
    240240 **********************************************/
    241 bool pmPSF_MaskApTrend (pmPSF *psf, pmPSF_ApTrendOptions option)
     241bool pmPSF_MaskApTrend (psPolynomial4D *trend,
     242                        pmPSF_ApTrendOptions option)
    242243{
    243244
    244245    switch (option) {
    245246    case PM_PSF_NONE:
    246         maskAllTerms (psf->ApTrend);
     247        maskAllTerms (trend);
    247248        return true;
    248249
    249250    case PM_PSF_CONSTANT:
    250         maskAllTerms (psf->ApTrend);
    251         psf->ApTrend->mask[0][0][0][0] = 0;  // unmask constant
     251        maskAllTerms (trend);
     252        trend->mask[0][0][0][0] = 0;  // unmask constant
    252253        return true;
    253254
    254255    case PM_PSF_SKYBIAS:
    255         maskAllTerms (psf->ApTrend);
    256         psf->ApTrend->mask[0][0][0][0] = 0;  // unmask constant
    257         psf->ApTrend->mask[0][0][1][0] = 0;  // unmask skybias
     256        maskAllTerms (trend);
     257        trend->mask[0][0][0][0] = 0;  // unmask constant
     258        trend->mask[0][0][1][0] = 0;  // unmask skybias
    258259        return true;
    259260
    260261    case PM_PSF_SKYSAT:
    261         maskAllTerms (psf->ApTrend);
    262         psf->ApTrend->mask[0][0][0][0] = 0;  // unmask constant
    263         psf->ApTrend->mask[0][0][1][0] = 0;  // unmask skybias
    264         psf->ApTrend->mask[0][0][0][1] = 0;  // unmask skybias
     262        maskAllTerms (trend);
     263        trend->mask[0][0][0][0] = 0;  // unmask constant
     264        trend->mask[0][0][1][0] = 0;  // unmask skybias
     265        trend->mask[0][0][0][1] = 0;  // unmask skysat
    265266        return true;
    266267
    267268    case PM_PSF_XY_LIN:
    268         maskAllTerms (psf->ApTrend);
    269         psf->ApTrend->mask[0][0][0][0] = 0;  // unmask constant
    270         psf->ApTrend->mask[1][0][0][0] = 0;  // unmask x
    271         psf->ApTrend->mask[0][1][0][0] = 0;  // unmask y
     269        maskAllTerms (trend);
     270        trend->mask[0][0][0][0] = 0;  // unmask constant
     271        trend->mask[1][0][0][0] = 0;  // unmask x
     272        trend->mask[0][1][0][0] = 0;  // unmask y
    272273        return true;
    273274
    274275    case PM_PSF_XY_QUAD:
    275         maskAllTerms (psf->ApTrend);
    276         psf->ApTrend->mask[0][0][0][0] = 0;  // unmask constant
    277         psf->ApTrend->mask[1][0][0][0] = 0;  // unmask x
    278         psf->ApTrend->mask[2][0][0][0] = 0;  // unmask x^2
    279         psf->ApTrend->mask[1][1][0][0] = 0;  // unmask x y
    280         psf->ApTrend->mask[0][1][0][0] = 0;  // unmask y
    281         psf->ApTrend->mask[0][2][0][0] = 0;  // unmask y^2
     276        maskAllTerms (trend);
     277        trend->mask[0][0][0][0] = 0;  // unmask constant
     278        trend->mask[1][0][0][0] = 0;  // unmask x
     279        trend->mask[2][0][0][0] = 0;  // unmask x^2
     280        trend->mask[1][1][0][0] = 0;  // unmask x y
     281        trend->mask[0][1][0][0] = 0;  // unmask y
     282        trend->mask[0][2][0][0] = 0;  // unmask y^2
    282283        return true;
    283284
    284285    case PM_PSF_SKY_XY_LIN:
    285         maskAllTerms (psf->ApTrend);
    286         psf->ApTrend->mask[0][0][0][0] = 0;  // unmask constant
    287         psf->ApTrend->mask[1][0][0][0] = 0;  // unmask x
    288         psf->ApTrend->mask[0][1][0][0] = 0;  // unmask y
    289         psf->ApTrend->mask[0][0][1][0] = 0;  // unmask skybias
     286        maskAllTerms (trend);
     287        trend->mask[0][0][0][0] = 0;  // unmask constant
     288        trend->mask[1][0][0][0] = 0;  // unmask x
     289        trend->mask[0][1][0][0] = 0;  // unmask y
     290        trend->mask[0][0][1][0] = 0;  // unmask skybias
    290291        return true;
    291292
    292293    case PM_PSF_SKYSAT_XY_LIN:
    293         maskAllTerms (psf->ApTrend);
    294         psf->ApTrend->mask[0][0][0][0] = 0;  // unmask constant
    295         psf->ApTrend->mask[1][0][0][0] = 0;  // unmask x
    296         psf->ApTrend->mask[0][1][0][0] = 0;  // unmask y
    297         psf->ApTrend->mask[0][0][1][0] = 0;  // unmask skybias
    298         psf->ApTrend->mask[0][0][0][1] = 0;  // unmask skysat
     294        maskAllTerms (trend);
     295        trend->mask[0][0][0][0] = 0;  // unmask constant
     296        trend->mask[1][0][0][0] = 0;  // unmask x
     297        trend->mask[0][1][0][0] = 0;  // unmask y
     298        trend->mask[0][0][1][0] = 0;  // unmask skybias
     299        trend->mask[0][0][0][1] = 0;  // unmask skysat
    299300        return true;
    300301
    301302    case PM_PSF_ALL:
    302303    default:
    303         maskAllTerms (psf->ApTrend);
    304         psf->ApTrend->mask[0][0][0][0] = 0;  // unmask constant
    305         psf->ApTrend->mask[0][0][1][0] = 0;  // unmask skybias
    306         psf->ApTrend->mask[0][0][0][1] = 0;  // unmask skysat
    307 
    308         psf->ApTrend->mask[1][0][0][0] = 0;  // unmask x
    309         psf->ApTrend->mask[2][0][0][0] = 0;  // unmask x^2
    310         psf->ApTrend->mask[1][1][0][0] = 0;  // unmask x y
    311         psf->ApTrend->mask[0][1][0][0] = 0;  // unmask y
    312         psf->ApTrend->mask[0][2][0][0] = 0;  // unmask y^2
     304        maskAllTerms (trend);
     305        trend->mask[0][0][0][0] = 0;  // unmask constant
     306        trend->mask[0][0][1][0] = 0;  // unmask skybias
     307        trend->mask[0][0][0][1] = 0;  // unmask skysat
     308
     309        trend->mask[1][0][0][0] = 0;  // unmask x
     310        trend->mask[2][0][0][0] = 0;  // unmask x^2
     311        trend->mask[1][1][0][0] = 0;  // unmask x y
     312        trend->mask[0][1][0][0] = 0;  // unmask y
     313        trend->mask[0][2][0][0] = 0;  // unmask y^2
    313314        return true;
    314315    }
Note: See TracChangeset for help on using the changeset viewer.