Source
349
349
template <class T>
350
350
void AWProjectWBFT::ftWeightImage(Lattice<T>& wtImage,
351
351
const Matrix<Float>& sumWt,
352
352
const Bool& doFFTNorm)
353
353
{
354
354
LogIO log_l(LogOrigin("AWProjectWBFT2", "ftWeightImage[R&D]"));
355
355
if (wtImageFTDone_p) return;
356
356
357
357
//cerr << "From ftWeightImage" << endl;
358
358
359
-
Bool doSumWtNorm=true;
360
-
if (sumWt.shape().nelements()==0) doSumWtNorm=false;
359
+
// Bool doSumWtNorm=true;
360
+
// if (sumWt.shape().nelements()==0) doSumWtNorm=false;
361
361
if ((sumWt.shape().nelements() < 2) ||
362
362
(sumWt.shape()(0) != wtImage.shape()(2)) ||
363
363
(sumWt.shape()(1) != wtImage.shape()(3)))
364
364
log_l << "Sum of weights per poln and chan is required" << LogIO::EXCEPTION;
365
365
//
366
366
// Use only the amplitude of the gridded weights. LatticeExpr
367
367
// classes while useful, appear to be too strict about types (the
368
368
// code below will not compile if the abs(wtImage) is not
369
369
// converted back to a complex type).
370
370