Inheritance diagram for lsst::meas::algorithms::shapelet::CrudeSolver:

Public Member Functions
	CrudeSolver (const PixelList &pix, double sigma, double I1, DVector &xinit)

	~CrudeSolver ()

void	calculateF (const DVector &x, DVector &f) const

void	calculateJ (const DVector &x, const DVector &f, DMatrix &df) const

Public Member Functions inherited from lsst::meas::algorithms::shapelet::NLSolver
	NLSolver ()

virtual	~NLSolver ()

virtual bool	solve (DVector &x, DVector &f) const

virtual void	getCovariance (DMatrix &cov) const

virtual void	getInverseCovariance (DMatrix &invcov) const

virtual bool	testJ (const DVector &, DVector &, std::ostream *os=0, double relerr=0.) const

virtual void	useHybrid ()

virtual void	useDogleg ()

virtual void	setFTol (double fTol)

virtual void	setGTol (double gTol)

virtual void	setTol (double fTol, double gTol)

virtual void	setMinStep (double minStep)

virtual void	setMaxIter (int maxIter)

virtual void	setTau (double tau)

virtual void	setDelta0 (double delta0)

virtual double	getFTol ()

virtual double	getGTol ()

virtual double	getMinStep ()

virtual int	getMaxIter ()

virtual double	getTau ()

virtual double	getDelta0 ()

virtual void	setOutput (std::ostream &os)

virtual void	useVerboseOutput ()

virtual void	useExtraVerboseOutput ()

virtual void	noUseVerboseOutput ()

virtual void	useDirectH ()

virtual void	useSVD ()

virtual void	useCholesky ()

virtual void	noUseDirectH ()

virtual void	noUseSVD ()

virtual void	noUseCholesky ()

Private Attributes
double	_sigma

DVector	_I

CDVector	_Z

DDiagMatrix	_W

CDVector	_Z1

DVector	_E

DVector	_Rsq

DVector	_f1

double	_I1

DVector &	_xInit

Detailed Description

Definition at line 36 of file CrudeMeasure.cc.

Constructor & Destructor Documentation

lsst::meas::algorithms::shapelet::CrudeSolver::CrudeSolver	(	const PixelList &	pix,
		double	sigma,
		double	I1,
		DVector &	xinit
	)

Definition at line 62 of file CrudeMeasure.cc.

                                                  :
         _sigma(sigma), _I(pix.size()), _Z(pix.size()), _W(pix.size()), 
         _Z1(pix.size()), _E(pix.size()), _Rsq(pix.size()), _f1(pix.size()), 
         _I1(I1), _xInit(xInit)
     {
         const int nPix = pix.size();
         for(int i=0;i<nPix;++i) {
             _Z(i) = pix[i].getPos();
             _I(i) = pix[i].getFlux();
             double mask = exp(-std::norm(pix[i].getPos()/_sigma)/2.);
             _W(i) = pix[i].getInverseSigma()*mask;
         }
     }

lsst::meas::algorithms::shapelet::CrudeSolver::~CrudeSolver ( )

inline

Definition at line 42 of file CrudeMeasure.cc.

42 {}

Member Function Documentation

void lsst::meas::algorithms::shapelet::CrudeSolver::calculateF	(	const DVector &	x,
		DVector &	f
	)		const

virtual

Implements lsst::meas::algorithms::shapelet::NLSolver.

Definition at line 78 of file CrudeMeasure.cc.

     {
         Assert(x.size() == 4);
         Assert(f.size() == _Z.size());
 
         if ((x-_xInit).TMV_subVector(0,3).TMV_normInf() > 2.) {
             f = 2.e10*_f1;
             f.TMV_addToAll(1.);
             return; 
         }
         if (x(3) < 0.) { 
             f = 2.e10*_f1;
             f.TMV_addToAll(1.);
             return; 
         }
 
         std::complex<double> zc(x[0],x[1]);
         double mu = x[2];
         double I0 = _I1 * x[3];
 
         double m0 = exp(-mu)/_sigma;
         // z' = m*(z-zc)
         //    = m*z - m*zc
         // x' = m x - m xc
         // y' = m y - m yc
         _Z1 = m0*_Z;
         _Z1.TMV_addToAll(-m0*zc);
 
         const int nPix = _Z.size();
         for(int i=0;i<nPix;++i) {
             double rsq = std::norm(_Z1[i]);
             _Rsq[i] = rsq;
             _E[i] = exp(-rsq/2.);
         }
         _f1 = _I - I0*_E;
 #ifdef USE_TMV
         _f1 = _W * _f1;
 #else
         _f1.array() *= _W.array();
 #endif
         f = _f1;
     }

void lsst::meas::algorithms::shapelet::CrudeSolver::calculateJ	(	const DVector &	x,
		const DVector &	f,
		DMatrix &	df
	)		const

virtual

Reimplemented from lsst::meas::algorithms::shapelet::NLSolver.

Definition at line 121 of file CrudeMeasure.cc.

     {
         //xdbg<<"Start J\n";
         Assert(x.size() == 4);
         Assert(f.size() == _Z.size());
         Assert(df.TMV_colsize() == _Z.size());
         Assert(df.TMV_rowsize() == 4);
 
         double mu = x[2];
         double I0 = _I1 * x[3];
         double m0 = exp(-mu)/_sigma;
 
         // fi = Wi * (Ii - I0 Ei)
         // dfi/dI0 = -Wi Ei
         // dfi/dmu = Wi I0 Ei (x1 dx1/dxc + y1 dy1/dxc) 
         // likewise for the other 4
         //
         // dx1/dxc = -m         dy1/dxc = 0
         // dx1/dyc = 0          dy1/dyc = -m
         // dx1/dmu = -x1        dy1/dmu = -y1
         //
 
 #ifdef USE_TMV
         df.col(0) = -m0 * _Z1.realPart();
         df.col(1) = -m0 * _Z1.imagPart();
         df.col(2) = -_Rsq;
         df.colRange(0,3) = I0 * DiagMatrixViewOf(_E) * df.colRange(0,3);
         df.col(3) = -_I1 * _E;
         df = _W * df;
 #else
         df.col(0).array() = _W.array() * _Z1.real().array();
         df.col(0).array() = _E.array() * df.col(0).array();
         df.col(0) *= -m0 * I0;
         df.col(1).array() = _W.array() * _Z1.imag().array();
         df.col(1).array() = _E.array() * df.col(1).array();
         df.col(1) *= -m0 * I0;
         df.col(2).array() = _W.array() * _Rsq.array();
         df.col(2).array() = _E.array() * df.col(2).array();
         df.col(2) *= -I0;
         df.col(3).array() = _W.array() * _E.array();
         df.col(3) *= -_I1;
 #endif
     }

Member Data Documentation

DVector lsst::meas::algorithms::shapelet::CrudeSolver::_E

mutableprivate

Definition at line 55 of file CrudeMeasure.cc.

DVector lsst::meas::algorithms::shapelet::CrudeSolver::_f1

mutableprivate

Definition at line 57 of file CrudeMeasure.cc.

DVector lsst::meas::algorithms::shapelet::CrudeSolver::_I

private

Definition at line 51 of file CrudeMeasure.cc.

double lsst::meas::algorithms::shapelet::CrudeSolver::_I1

private

Definition at line 58 of file CrudeMeasure.cc.

DVector lsst::meas::algorithms::shapelet::CrudeSolver::_Rsq

mutableprivate

Definition at line 56 of file CrudeMeasure.cc.

double lsst::meas::algorithms::shapelet::CrudeSolver::_sigma

private

Definition at line 50 of file CrudeMeasure.cc.

DDiagMatrix lsst::meas::algorithms::shapelet::CrudeSolver::_W

private

Definition at line 53 of file CrudeMeasure.cc.

DVector& lsst::meas::algorithms::shapelet::CrudeSolver::_xInit

private

Definition at line 59 of file CrudeMeasure.cc.

CDVector lsst::meas::algorithms::shapelet::CrudeSolver::_Z

private

Definition at line 52 of file CrudeMeasure.cc.

CDVector lsst::meas::algorithms::shapelet::CrudeSolver::_Z1

mutableprivate

Definition at line 54 of file CrudeMeasure.cc.

The documentation for this class was generated from the following file:

/home/lsstsw/stack/Linux64/meas_algorithms/11.0-2-gb8b8ce7/src/shapelet/CrudeMeasure.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation