Source
AlwaysAssert(near(wblc[0].getValue("arcmin"), (centerx+radius).getValue("arcmin")), AipsError);
//# Copyright (C) 2000,2001,2003
//# Associated Universities, Inc. Washington DC, USA.
//#
//# This program is free software; you can redistribute it and/or modify it
//# under the terms of the GNU General Public License as published by the Free
//# Software Foundation; either version 2 of the License, or (at your option)
//# any later version.
//#
//# This program is distributed in the hope that it will be useful, but WITHOUT
//# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
//# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for
//# more details.
//#
//# You should have received a copy of the GNU General Public License along
//# with this program; if not, write to the Free Software Foundation, Inc.,
//# 675 Massachusetts Ave, Cambridge, MA 02139, USA.
//#
//# Correspondence concerning AIPS++ should be addressed as follows:
//# Internet email: aips2-request@nrao.edu.
//# Postal address: AIPS++ Project Office
//# National Radio Astronomy Observatory
//# 520 Edgemont Road
//# Charlottesville, VA 22903-2475 USA
//#
using namespace casa;
int main () {
try {
LogIO log;
IPosition shape(4, 400, 400, 1, 5500);
CoordinateSystem csys = CoordinateUtil::defaultCoords4D();
Vector<Double> refpix = csys.referencePixel();
refpix[0] = 200;
csys.setReferencePixel(refpix);
AnnRegion::unitInit();
{
log << LogIO::NORMAL
<<"Test mixed world and pixel coordinates throws exception"
<< LogIO::POST;
Bool thrown = true;
Quantity centerx(0.01, "deg");
Quantity centery(0.01, "pix");
Quantity radius(30, "arcsec");
Quantity beginFreq, endFreq;
String dirTypeString = MDirection::showType(
csys.directionCoordinate().directionType(false)
);
String freqRefFrameString = MFrequency::showType(
csys.spectralCoordinate().frequencySystem()
);
String dopplerString = MDoppler::showType(
csys.spectralCoordinate().velocityDoppler()
);
Quantity restfreq(
csys.spectralCoordinate().restFrequency(), "Hz"
);
Vector<Stokes::StokesTypes> stokes(0);
try {
AnnCircle circle(
centerx, centery, radius, dirTypeString,
csys, shape, beginFreq, endFreq, freqRefFrameString,
dopplerString, restfreq, stokes, false
);
thrown = false;
} catch (const AipsError& x) {
log << LogIO::NORMAL
<< "Exception thrown as expected: " << x.getMesg()
<< LogIO::POST;
}
AlwaysAssert(thrown, AipsError);
}
{
log << LogIO::NORMAL
<< "Test that bad quantity for world direction coordinate throws exception"
<< LogIO::POST;
Bool thrown = true;
Quantity centerx(0.01, "deg");