casa6

//# Jones.h: Definition of Jones

//# Copyright (C) 1996,1997,2000,2001,2002,2003

//# Associated Universities, Inc. Washington DC, USA.

//#

//# This library is free software; you can redistribute it and/or modify it

//# under the terms of the GNU Library General Public License as published by

//# the Free Software Foundation; either version 2 of the License, or (at your

//# option) any later version.

//#

//# This library 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 Library General Public

//# License for more details.

//#

//# You should have received a copy of the GNU Library General Public License

//# along with this library; if not, write to the Free Software Foundation,

//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.

//#

//# Correspondence concerning AIPS++ should be adressed 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

//#

//#

​

​

#include <synthesis/MeasurementComponents/Jones.h>

#include <casa/aips.h>

#include <casa/BasicSL/Complex.h>

#include <casa/iostream.h>

#include <casa/Exceptions/Error.h>

#include <casa/namespace.h>

​

using namespace casacore;

namespace casa { //# NAMESPACE CASA - BEGIN

​

// Constructor

Jones::Jones() : 

  j0_(NULL), 

  ok0_(NULL),

  j_(NULL), 

  ji_(NULL),

  ok_(NULL),

  oki_(NULL),

  cOne_(1.0),

  cZero_(0.0),

  scalardata_(false),

  vtmp_(VisVector::Four,true)

{}

​

// In place invert

void Jones::invert() {

​

  if (!ok_) throw(AipsError("Illegal use of Jones::invert()"));

​

  if (ok_[0]&&ok_[1]&&ok_[2]&&ok_[3]) {

    Complex det,tmp;

    tmp=j_[0];

    det=j_[0]*j_[3]-j_[1]*j_[2];

    if (abs(det)>0.0) {