Commits

David Mehringer authored 1af00bb9978
debugging

casatools/src/code/mstransform/TVI/StatWtFloatingWindowDataAggregator.cc

Modified
13 13 //# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 14 //# Lesser General Public License for more details.
15 15 //#
16 16 //# You should have received a copy of the GNU Lesser General Public
17 17 //# License along with this library; if not, write to the Free Software
18 18 //# Foundation, Inc., 59 Temple Place, Suite 330, Boston,
19 19 //# MA 02111-1307 USA
20 20
21 21 #include <mstransform/TVI/StatWtFloatingWindowDataAggregator.h>
22 22
23 +// debug
24 +#include <casacore/casa/Logging/LogIO.h>
25 +
23 26 #include <casacore/scimath/StatsFramework/ClassicalStatistics.h>
24 27
25 28 #include <msvis/MSVis/ViImplementation2.h>
26 29 #include <mstransform/TVI/StatWtTypes.h>
27 30
28 31 #ifdef _OPENMP
29 32 #include <omp.h>
30 33 #endif
31 34
32 35 using namespace casacore;
193 196 slStart[2] = subchunkStartRowNum;
194 197 sl.setStart(slStart);
195 198 slEnd = newShape - 1;
196 199 sl.setEnd(slEnd);
197 200 chunkData(sl) = dataCube;
198 201 chunkFlags(sl) = resultantFlags;
199 202 }
200 203 subChunkToTimeStamp.push_back(subchunkTime);
201 204 subchunkStartRowNum += nrows;
202 205 }
206 + LogIO log(LogOrigin("StatWtFloatingWindowDataAggregator", __func__));
207 + log << LogIO::WARN << "Calling _computeWeights at " << __LINE__ << LogIO::POST;
203 208 _computeWeights(
204 209 chunkData, chunkFlags, Vector<Double>(exposures), rowMap, spw
205 210 );
206 211 }
207 212
208 213 void StatWtFloatingWindowDataAggregator::weightSingleChanBin(
209 214 Matrix<Float>& wtmat, Int nrows
210 215 ) const {
211 216 Vector<rownr_t> rowIDs;
212 217 _vii->getRowIds(rowIDs);
291 296 // slice up for correlations and channel binning
292 297 intraChunkSliceEnd[2] = dataShape[2] - 1;
293 298 for (uInt corr=0; corr<ncorr; ++corr) {
294 299 if (! _combineCorr) {
295 300 intraChunkSliceStart[0] = corr;
296 301 intraChunkSliceEnd[0] = corr;
297 302 }
298 303 auto citer = chanBins.begin();
299 304 auto cend = chanBins.end();
300 305 auto iChanBin = 0;
306 + LogIO log(LogOrigin("StatWtFloatingWindowDataAggregator", __func__));
307 + log << LogIO::WARN << "Calling _computeWeight at " << __LINE__ << LogIO::POST;
308 +
301 309 for (; citer!=cend; ++citer, ++iChanBin) {
302 310 intraChunkSliceStart[1] = citer->start;
303 311 intraChunkSliceEnd[1] = citer->end;
304 312 intraChunkSlice.setStart(intraChunkSliceStart);
305 313 intraChunkSlice.setEnd(intraChunkSliceEnd);
306 314 _weights(corr, iChanBin, iRow)
307 315 = _varianceComputer->computeWeight(
308 316 dataArray(intraChunkSlice), flagArray(intraChunkSlice),
309 317 exposureVector, spw, exposures[iRow]
310 318 );

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

Add shortcut