Details
-
Type:
Story
-
Status: To Do
-
Resolution: Unresolved
-
Fix Version/s: None
-
Component/s: meas_astrom
-
Labels:None
-
Story Points:1
-
Team:Data Release Production
Description
With lsst_distrib tagged w_2017_25 (meas_astrom 13.0-14-g9415442+2):
pprice@tiger-sumire:/tigress/HSC/users/price/bell $ cat fitter.py
|
from lsst.meas.astrom import MatchPessimisticBTask, FitSipDistortionTask
|
config.calibrate.astrometry.matcher.retarget(MatchPessimisticBTask)
|
config.calibrate.astrometry.wcsFitter.retarget(FitSipDistortionTask)
|
pprice@tiger-sumire:/tigress/HSC/users/price/bell $ processCcd.py DATA --calib /tigress/HSC/HSC/CALIB-LSST-20170105 --rerun test --id visit=19948 ccd=41 -C fitter.py --clobber-config
|
[...]
|
python: src/ScaledPolynomialTransformFitter.cc:320: void lsst::meas::astrom::ScaledPolynomialTransformFitter::fit(int): Assertion `(sxx*fxx + sxy*fxy).isApproxToConstant(1.0)' failed.
|
After changing line 302 to:
if (_keys.outputErr.isValid() && _data[i1].get(_keys.outputErr).allFinite()) {I hit another assertion failure:
python: src/ScaledPolynomialTransformFitter.cc:470: std::pair<double, long unsigned int> lsst::meas::astrom::ScaledPolynomialTransformFitter::rejectOutliers(const lsst::meas::astrom::OutlierRejectionControl&): Assertion `static_cast<std::size_t>(nGood + nClip) == _data.size()' failed.Jim Bosch says this is from multiple NAN values, and that the rankings should be a std::multimap instead of std::map.