00001 //# SDAlgorithmMSClean.h: Definition for SDAlgorithmMSClean 00002 //# Copyright (C) 1996,1997,1998,1999,2000,2002 00003 //# Associated Universities, Inc. Washington DC, USA. 00004 //# 00005 //# This library is free software; you can redistribute it and/or modify it 00006 //# under the terms of the GNU Library General Public License as published by 00007 //# the Free Software Foundation; either version 2 of the License, or (at your 00008 //# option) any later version. 00009 //# 00010 //# This library is distributed in the hope that it will be useful, but WITHOUT 00011 //# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 00012 //# FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public 00013 //# License for more details. 00014 //# 00015 //# You should have received a copy of the GNU Library General Public License 00016 //# along with this library; if not, write to the Free Software Foundation, 00017 //# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA. 00018 //# 00019 //# Correspondence concerning AIPS++ should be adressed as follows: 00020 //# Internet email: aips2-request@nrao.edu. 00021 //# Postal address: AIPS++ Project Office 00022 //# National Radio Astronomy Observatory 00023 //# 520 Edgemont Road 00024 //# Charlottesville, VA 22903-2475 USA 00025 //# 00026 //# 00027 //# $Id$ 00028 00029 #ifndef SYNTHESIS_SDALGORITHMMSCLEAN_H 00030 #define SYNTHESIS_SDALGORITHMMSCLEAN_H 00031 00032 #include <ms/MeasurementSets/MeasurementSet.h> 00033 #include <synthesis/MeasurementComponents/SkyModel.h> 00034 #include <casa/Arrays/Matrix.h> 00035 #include <images/Images/ImageInterface.h> 00036 #include <images/Images/PagedImage.h> 00037 #include <images/Images/TempImage.h> 00038 #include <casa/Logging/LogMessage.h> 00039 #include <casa/Logging/LogSink.h> 00040 #include <casa/System/PGPlotter.h> 00041 00042 #include<synthesis/ImagerObjects/SDAlgorithmBase.h> 00043 #include <synthesis/MeasurementEquations/MatrixCleaner.h> 00044 00045 namespace casa { //# NAMESPACE CASA - BEGIN 00046 00047 /* Forware Declaration */ 00048 class SIMinorCycleController; 00049 00050 00051 class SDAlgorithmMSClean : public SDAlgorithmBase 00052 { 00053 public: 00054 00055 // Empty constructor 00056 SDAlgorithmMSClean(Vector<Float> scalesizes, 00057 Float smallscalebias=0.6, 00058 // Int stoplargenegatives=-2, 00059 Int stoppointmode=-1 ); 00060 virtual ~SDAlgorithmMSClean(); 00061 00062 // void restore( CountedPtr<SIImageStore> imagestore ); 00063 00064 protected: 00065 00066 // Local functions to be overloaded by various algorithm deconvolvers. 00067 void takeOneStep( Float loopgain, Int cycleNiter, Float cycleThreshold, 00068 Float &peakresidual, Float &modelflux, Int &iterdone ); 00069 // void initializeDeconvolver( Float &peakresidual, Float &modelflux ); 00070 void initializeDeconvolver(); 00071 void finalizeDeconvolver(); 00072 00073 Array<Float> itsMatPsf, itsMatResidual, itsMatModel; 00074 Array<Float> itsMatMask; // Make an array if we eventually use multi-term masks... 00075 00076 MatrixCleaner itsCleaner; 00077 Vector<Float> itsScaleSizes; 00078 Float itsSmallScaleBias; 00079 //Int itsStopLargeNegatives; 00080 Int itsStopPointMode; 00081 00082 private: 00083 Bool itsMCsetup; 00084 00085 }; 00086 00087 } //# NAMESPACE CASA - END 00088 00089 #endif