casa6
Public
Commits
R. V. Urvashi authored 7f978811d1d
| 223 223 | } |
| 224 224 | |
| 225 225 | ///////////////// |
| 226 226 | //findMaxAbs( itsMatResiduals[0], itsPeakResidual, itsMaxPos ); |
| 227 227 | //peakresidual = itsPeakResidual; |
| 228 228 | |
| 229 229 | peakresidual = itsMTCleaner.getpeakresidual(); |
| 230 230 | |
| 231 231 | modelflux = sum( itsMatModels[0] ); // Performance hog ? |
| 232 232 | |
| 233 + | /* // Enable in CAS-13872 |
| 233 234 | // Retrieve residual to be saved to the .residual file in finalizeDeconvolver |
| 234 235 | for(uInt tix=0; tix<itsNTerms; tix++) |
| 235 236 | { |
| 236 237 | casacore::Matrix<Float> tmp; |
| 237 238 | itsMTCleaner.getresidual(tix, tmp); // possible room for optimization here -> get residual without extra tmp copy? maybe change getResidual to accept an array? |
| 238 239 | itsMatResiduals[tix] = tmp; |
| 239 240 | } |
| 241 + | */ |
| 240 242 | } |
| 241 243 | |
| 242 244 | void SDAlgorithmMSMFS::finalizeDeconvolver() |
| 243 245 | { |
| 244 246 | // itsResidual.put( itsMatResidual ); |
| 245 247 | // itsModel.put( itsMatModel ); |
| 246 248 | |
| 247 249 | // Why is this needed ? If the matrices are by reference, then why do we need this ? |
| 248 250 | for(uInt tix=0; tix<itsNTerms; tix++) |
| 249 251 | { |