- Timestamp:
- Jun 17, 2014, 5:58:38 AM (12 years ago)
- Location:
- branches/eam_branches/ipp-20140610
- Files:
-
- 2 edited
-
. (modified) (1 prop)
-
psModules/src/objects/models/pmModel_DEV.c (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20140610
- Property svn:mergeinfo changed
/trunk (added) merged: 36846-36871,36874,36879,36882
- Property svn:mergeinfo changed
-
branches/eam_branches/ipp-20140610/psModules/src/objects/models/pmModel_DEV.c
r36623 r36884 36 36 #include "pmMoments.h" 37 37 #include "pmModelFuncs.h" 38 #include "pmModelClass.h" 38 39 #include "pmModel.h" 39 40 #include "pmModelUtils.h" 40 #include "pmModelClass.h"41 41 #include "pmSourceMasks.h" 42 42 #include "pmSourceExtendedPars.h" … … 59 59 # define PM_MODEL_LIMITS pmModelLimits_DEV 60 60 # define PM_MODEL_RADIUS pmModelRadius_DEV 61 # define PM_MODEL_SET_FWHM pmModelSetFWHM_DEV 61 62 # define PM_MODEL_FROM_PSF pmModelFromPSF_DEV 62 63 # define PM_MODEL_PARAMS_FROM_PSF pmModelParamsFromPSF_DEV … … 123 124 if (radius <= 1.5) { 124 125 // Nsub ~ 10*index^2 + 1 125 psEllipseAxes axes = pmPSF_ModelToAxes(PAR, pmModelClassGetType ("PS_MODEL_DEV"));126 psEllipseAxes axes = pmPSF_ModelToAxes(PAR, true); // DEV uses Reff 126 127 int Nsub = 2 * ((int)(25 / axes.minor)) + 1; 127 128 Nsub = PS_MIN (Nsub, 121); … … 336 337 } 337 338 339 psF64 PM_MODEL_SET_FWHM (const psVector *params, psF64 sigma) { 340 return (NAN); 341 } 342 338 343 bool PM_MODEL_FROM_PSF (pmModel *modelPSF, pmModel *modelFLT, const pmPSF *psf) 339 344 { … … 357 362 // the 2D PSF model fits polarization terms (E0,E1,E2) 358 363 // convert to shape terms (SXX,SYY,SXY) 359 bool useReff = pmModelUseReff (modelPSF->type);364 bool useReff = modelPSF->class->useReff; 360 365 if (!pmPSF_FitToModel (out, 0.1, useReff)) { 361 366 psTrace("psModules.objects", 5, "Failed to fit object at (r,c) = (%.1f,%.1f)", in[PM_PAR_YPOS], in[PM_PAR_XPOS]); … … 411 416 // convert to shape terms (SXX,SYY,SXY) 412 417 // XXX user-defined value for limit? 413 bool useReff = pmModelUseReff (model->type);418 bool useReff = model->class->useReff; 414 419 if (!pmPSF_FitToModel (PAR, 0.1, useReff)) { 415 420 psTrace ("psModules.objects", 3, "Failed to fit object at (r,c) = (%.1f,%.1f)", Xo, Yo);
Note:
See TracChangeset
for help on using the changeset viewer.
