Line data Source code
1 : // -*- C++ -*- 2 : //# AzElAperture.cc: Implementation of the AzElAperture class 3 : //# Copyright (C) 1997,1998,1999,2000,2001,2002,2003 4 : //# Associated Universities, Inc. Washington DC, USA. 5 : //# 6 : //# This library is free software; you can redistribute it and/or modify it 7 : //# under the terms of the GNU Library General Public License as published by 8 : //# the Free Software Foundation; either version 2 of the License, or (at your 9 : //# option) any later version. 10 : //# 11 : //# This library is distributed in the hope that it will be useful, but WITHOUT 12 : //# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 13 : //# FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public 14 : //# License for more details. 15 : //# 16 : //# You should have received a copy of the GNU Library General Public License 17 : //# along with this library; if not, write to the Free Software Foundation, 18 : //# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA. 19 : //# 20 : //# Correspondence concerning AIPS++ should be addressed as follows: 21 : //# Internet email: aips2-request@nrao.edu. 22 : //# Postal address: AIPS++ Project Office 23 : //# National Radio Astronomy Observatory 24 : //# 520 Edgemont Road 25 : //# Charlottesville, VA 22903-2475 USA 26 : //# 27 : //# $Id$ 28 : // 29 : #include <synthesis/TransformMachines2/AzElAperture.h> 30 : #include <ostream> 31 : using namespace casacore; 32 : namespace casa{ 33 : namespace refim{ 34 : // 35 0 : void AzElAperture::rotate(const VisBuffer2& vb, CFCell& cfc, const Double& rotAngleIncr) 36 : { 37 0 : LogIO log_l(LogOrigin("AzElAperture", "rotate")); 38 0 : Double R2D = 57.2957795131; 39 : // If the A-Term is a No-Op, the resulting CF is rotationally 40 : // symmetric. 41 0 : if (isNoOp()) return; 42 : 43 0 : Double actualPA = getPA(vb), currentCFPA = cfc.pa_p.getValue("rad"); 44 0 : cerr << "AzEL: ######### " << cfc.pa_p.getValue("deg") << " " << actualPA*R2D << " " << cfc.getStorage()->shape() << endl; 45 0 : Double dPA = currentCFPA-actualPA; 46 : //Double dPA = cachedCFCPA_p-actualPA; 47 0 : if (fabs(dPA) > fabs(rotAngleIncr)) 48 : // if (fabs(actualPA-currentCFPA) > 0.0) 49 : { 50 : //if (cachedCFCPA_p == 400.0) cachedCFCPA_p=currentCFPA; 51 : 52 0 : if (cachedCFC_p.nelements() == 0) 53 : { 54 0 : log_l << "Making a base copy of the current CFC@ PA=" << currentCFPA*R2D << LogIO::POST; 55 0 : cachedCFC_p.assign(*cfc.getStorage()); 56 0 : cachedCFCPA_p = cfc.pa_p.getValue("rad"); 57 : } 58 0 : dPA = cachedCFCPA_p-actualPA; 59 0 : Array<TT> inData; 60 : //inData.assign(*cfc.getStorage()); 61 0 : inData.assign(cachedCFC_p); 62 : // log_l << actualPA << " " << currentCFPA << LogIO::POST; 63 : // cerr << "dPA = " << actualPA << " " << currentCFPA << " " 64 : // << actualPA - currentCFPA << " " << &(*cfc.getStorage()) << endl; 65 0 : log_l << "Rotating the base CFC from PA=" << cachedCFCPA_p*R2D 66 : << " to " << actualPA*R2D << " " << dPA*R2D << " " 67 : << inData.shape() << " " 68 0 : << "SPW: " << vb.spectralWindows() 69 0 : << " Field: " << vb.fieldId() 70 0 : << LogIO::POST; 71 0 : SynthesisUtils::rotateComplexArray(log_l, inData, cfc.coordSys_p, 72 0 : *cfc.getStorage(), 73 : dPA);//,"LINEAR"); 74 : // currentCFPA-actualPA);//,"LINEAR"); 75 : // Update the PA value in the CF-Cache 76 0 : cfc.pa_p=Quantity(actualPA, "rad"); 77 : 78 : // static Int ii=0; 79 : // { 80 : // ostringstream name; 81 : // name << "cfc." << ii << ".im"; 82 : // cfc.makePersistent(name.str().c_str()); 83 : // ii++; 84 : // } 85 : } 86 : }; 87 : 88 0 : void AzElAperture::rotate2(const VisBuffer2& vb, CFCell& baseCFC, CFCell& cfc, const Double& rotAngleIncr) 89 : { 90 0 : LogIO log_l(LogOrigin("AzElAperture", "rotate")); 91 : 92 : // If the A-Term is a No-Op, the resulting CF is rotationally 93 : // symmetric. 94 0 : if (isNoOp()) return; 95 : 96 : (void)baseCFC; 97 : 98 0 : Double actualPA = getPA(vb), currentCFPA = cfc.pa_p.getValue("rad"); 99 : //Double baseCFCPA=baseCFC.pa_p.getValue("rad"); 100 : 101 0 : Double dPA = currentCFPA-actualPA; 102 : 103 0 : if (fabs(dPA) > fabs(rotAngleIncr)) 104 : { 105 0 : Array<TT> inData; 106 : //inData.assign(*baseCFC.getStorage()); 107 : //dPA = baseCFCPA-actualPA; 108 0 : dPA = currentCFPA-actualPA; 109 0 : inData.assign(*cfc.getStorage()); 110 : 111 0 : SynthesisUtils::rotateComplexArray(log_l, inData, cfc.coordSys_p, 112 0 : *cfc.getStorage(), 113 : dPA);//,"LINEAR"); 114 : // currentCFPA-actualPA);//,"LINEAR"); 115 0 : cfc.pa_p=Quantity(actualPA, "rad"); 116 : } 117 : }; 118 : }; 119 : };