Commits

Takahiro Tsutsumi authored d267152b02d
Added a tolerance for nigma threshold stopping criterion to use (one

part in 100)

code/synthesis/ImagerObjects/SDAlgorithmBase.cc

Modified
185 185 nsigmathresh = (Float)medians(statsindex) + nsigma * (Float)robustrms(statsindex);
186 186 }
187 187 else {
188 188 nsigmathresh = nsigma * (Float)robustrms(statsindex);
189 189 }
190 190 thresholdtouse = max( nsigmathresh, loopcontrols.getCycleThreshold());
191 191 }
192 192 else {
193 193 thresholdtouse = loopcontrols.getCycleThreshold();
194 194 }
195 - os << LogIO::DEBUG1<<"loopcontrols.getCycleThreshold()="<<loopcontrols.getCycleThreshold()<<LogIO::POST;
196 - os << LogIO::DEBUG1<< "thresholdtouse="<<thresholdtouse<<LogIO::POST;
195 + //os << LogIO::DEBUG1<<"loopcontrols.getCycleThreshold()="<<loopcontrols.getCycleThreshold()<<LogIO::POST;
196 + os << LogIO::NORMAL3<<"current CycleThreshold="<<loopcontrols.getCycleThreshold()<<" nsigma threshold="<<nsigmathresh<<LogIO::POST;
197 197 String thresholddesc = (thresholdtouse == loopcontrols.getCycleThreshold() ? "cyclethreshold" : "n-sigma");
198 + os << LogIO::NORMAL3<< "thresholdtouse="<< thresholdtouse << "("<<thresholddesc<<")"<< LogIO::POST;
198 199
199 200 if (thresholddesc=="n-sigma") {
200 - os << LogIO::DEBUG1<< "Set nsigma thresh="<<nsigmathresh<<LogIO::POST;
201 + //os << LogIO::DEBUG1<< "Set nsigma thresh="<<nsigmathresh<<LogIO::POST;
201 202 loopcontrols.setNsigmaThreshold(nsigmathresh);
202 203 }
203 204 loopcontrols.setPeakResidual( peakresidual );
204 205 loopcontrols.resetMinResidual(); // Set it to current initial peakresidual.
205 206
206 207
207 208 stopCode = checkStop( loopcontrols, peakresidual );
208 209
209 210
210 211 // stopCode=0;

Everything looks good. We'll let you know here if there's anything you should know about.

Add shortcut