IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Oct 13, 2006, 2:56:37 PM (20 years ago)
Author:
magnier
Message:

support for psfTrendMask in pmPSFAlloc

Location:
trunk/psModules/src/objects
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/psModules/src/objects/pmPSFtry.c

    r9528 r9564  
    55 *  @author EAM, IfA
    66 *
    7  *  @version $Revision: 1.18 $ $Name: not supported by cvs2svn $
    8  *  @date $Date: 2006-10-13 02:29:14 $
     7 *  @version $Revision: 1.19 $ $Name: not supported by cvs2svn $
     8 *  @date $Date: 2006-10-14 00:56:37 $
    99 *
    1010 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
     
    5555
    5656// allocate a pmPSFtry based on the desired sources and the model (identified by name)
    57 pmPSFtry *pmPSFtryAlloc (psArray *sources, char *modelName, bool poissonErrors)
     57pmPSFtry *pmPSFtryAlloc (psArray *sources, char *modelName, bool poissonErrors, psPolynomial2D *psfTrendMask)
    5858{
    5959
     
    6464    // XXX probably need to increment ref counter
    6565    type           = pmModelSetType (modelName);
    66     test->psf      = pmPSFAlloc (type, poissonErrors);
     66    test->psf      = pmPSFAlloc (type, poissonErrors, psfTrendMask);
    6767    test->sources  = psMemIncrRefCounter(sources);
    6868    test->modelEXT = psArrayAlloc (sources->n);
     
    9595// mask values indicate the reason the source was rejected:
    9696
    97 pmPSFtry *pmPSFtryModel (psArray *sources, char *modelName, float RADIUS, bool poissonErrors)
     97pmPSFtry *pmPSFtryModel (psArray *sources, char *modelName, float RADIUS, bool poissonErrors, psPolynomial2D *psfTrendMask)
    9898{
    9999    bool status;
     
    105105    int Npsf = 0;
    106106
    107     pmPSFtry *psfTry = pmPSFtryAlloc (sources, modelName, poissonErrors);
     107    pmPSFtry *psfTry = pmPSFtryAlloc (sources, modelName, poissonErrors, psfTrendMask);
    108108
    109109    // stage 1:  fit an independent model (freeModel) to all sources
     
    280280    psLogMsg ("pmPSFtryMetric", 4, "fit stats: %f +/- %f\n", stats->sampleMedian, stats->sampleStdev);
    281281
    282     pmPSF_MaskApTrend (psfTry->psf->ApTrend, PM_PSF_SKYBIAS);
     282    pmPSFMaskApTrend (psfTry->psf->ApTrend, PM_PSF_APTREND_SKYBIAS);
    283283    psfTry->psf->ApTrend->coeff[0][0][0][0] = poly->coeff[0];
    284284    psfTry->psf->ApTrend->coeff[0][0][1][0] = 0;
  • trunk/psModules/src/objects/pmPSFtry.h

    r7311 r9564  
    66 *  @author EAM, IfA
    77 *
    8  *  @version $Revision: 1.6 $ $Name: not supported by cvs2svn $
    9  *  @date $Date: 2006-06-03 01:02:08 $
     8 *  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
     9 *  @date $Date: 2006-10-14 00:56:37 $
    1010 *
    1111 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
     
    7878 *
    7979 */
    80 pmPSFtry *pmPSFtryAlloc (psArray *sources, char *modelName, bool poissonErrors);
     80pmPSFtry *pmPSFtryAlloc (psArray *sources, char *modelName, bool poissonErrors, psPolynomial2D *psfTrendMask);
    8181
    8282/** pmPSFtryModel()
     
    8787 *
    8888 */
    89 pmPSFtry *pmPSFtryModel (psArray *sources, char *modelName, float RADIUS, bool poissonErrors);
     89pmPSFtry *pmPSFtryModel (psArray *sources, char *modelName, float RADIUS, bool poissonErrors, psPolynomial2D *psfTrendMask);
    9090
    9191/** pmPSFtryMetric()
Note: See TracChangeset for help on using the changeset viewer.