Uploaded image for project: 'Data Management'
  1. Data Management
  2. DM-11103

Assertion failure in ScaledPolynomialTransformFitter

    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.
      

        Attachments

          Issue Links

            Activity

            Hide
            price Paul Price added a comment - - edited

            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.

            Show
            price Paul Price added a comment - - edited 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 .

              People

              • Assignee:
                jbosch Jim Bosch
                Reporter:
                price Paul Price
                Watchers:
                Paul Price
              • Votes:
                0 Vote for this issue
                Watchers:
                1 Start watching this issue

                Dates

                • Created:
                  Updated: