lsst::meas::algorithms::shapelet::EllipseSolver3 Class Reference

#include <EllipseSolver.h>

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

Classes
struct	ESImpl3

Public Member Functions
	EllipseSolver3 (const BVec &b0, int order, bool fixcen=false, bool fixgam=false, bool fixmu=false)
	~EllipseSolver3 ()
void	calculateF (const DVector &x, DVector &f) const
void	calculateJ (const DVector &x, const DVector &f, DMatrix &df) const
void	useNumericJ ()
void	dontZeroB11 ()
void	getCovariance (DMatrix &cov) const
void	getInverseCovariance (DMatrix &invcov) const
void	callF (const DVector &x, DVector &f) const
bool	solve (DVector &x, DVector &f) const
bool	testJ (const DVector &x, DVector &f, std::ostream *os=0, double relerr=0.) const
Public Member Functions inherited from lsst::meas::algorithms::shapelet::NLSolver
	NLSolver ()
virtual	~NLSolver ()
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
ESImpl3 *	_pimpl

Detailed Description

Definition at line 14 of file EllipseSolver.h.

Constructor & Destructor Documentation

lsst::meas::algorithms::shapelet::EllipseSolver3::EllipseSolver3	(	const BVec &	b0,
		int	order,
		bool	fixcen = false,
		bool	fixgam = false,
		bool	fixmu = false
	)

Definition at line 97 of file EllipseSolver.cc.

                                              :
        _pimpl(new ESImpl3(b0,order,fixcen,fixgam,fixmu))
    {}

lsst::meas::algorithms::shapelet::EllipseSolver3::~EllipseSolver3

(

)

Definition at line 103 of file EllipseSolver.cc.

    { delete _pimpl; }

Member Function Documentation

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

virtual

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

Definition at line 106 of file EllipseSolver.cc.

    { _pimpl->calculateF(x,f); }

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

virtual

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

Definition at line 109 of file EllipseSolver.cc.

    { 
#ifdef ALWAYS_NUMERIC_J
        NLSolver::calculateJ(x,f,J);
#else
        if (_pimpl->numeric_j) NLSolver::calculateJ(x,f,J);
        else _pimpl->calculateJ(x,f,J); 
#endif
    }

void lsst::meas::algorithms::shapelet::EllipseSolver3::callF	(	const DVector &	x,
		DVector &	f
	)		const

Definition at line 915 of file EllipseSolver.cc.

    {
        Assert(x.size() == 5);
        Assert(f.size() == 5);
        _pimpl->xinit = x;
        if (_pimpl->fixcen) { _pimpl->fixuc = x[0]; _pimpl->fixvc = x[1]; }
        if (_pimpl->fixgam) { _pimpl->fixg1 = x[2]; _pimpl->fixg2 = x[3]; }
        if (_pimpl->fixmu) { _pimpl->fixm = x[4]; }

        _pimpl->x_short = _pimpl->U * x;

        calculateF(_pimpl->x_short,_pimpl->f_short);

        EIGEN_Transpose(f) = EIGEN_Transpose(_pimpl->f_short)*_pimpl->U;
    }

void lsst::meas::algorithms::shapelet::EllipseSolver3::dontZeroB11

(

)

Definition at line 893 of file EllipseSolver.cc.

{ _pimpl->zerob11 = false; this->useSVD(); }

void lsst::meas::algorithms::shapelet::EllipseSolver3::getCovariance ( DMatrix &  cov ) const

virtual

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

Definition at line 895 of file EllipseSolver.cc.

    {
        DMatrix cov1(_pimpl->x_short.size(),_pimpl->x_short.size());
        NLSolver::getCovariance(cov1);
        cov = _pimpl->U.transpose()*cov1*_pimpl->U;
        dbg<<"getCovariance:\n";
        dbg<<"cov1 = "<<cov1<<std::endl;
        dbg<<"full cov = "<<cov<<std::endl;
    }

void lsst::meas::algorithms::shapelet::EllipseSolver3::getInverseCovariance ( DMatrix &  invcov ) const

virtual

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

Definition at line 905 of file EllipseSolver.cc.

    {
        DMatrix invcov1(_pimpl->x_short.size(),_pimpl->x_short.size());
        NLSolver::getInverseCovariance(invcov1);
        invcov = _pimpl->U.transpose()*invcov1*_pimpl->U;
        dbg<<"getInverseCovariance:\n";
        dbg<<"invcov1 = "<<invcov1<<std::endl;
        dbg<<"full invcov = "<<invcov<<std::endl;
    }

bool lsst::meas::algorithms::shapelet::EllipseSolver3::solve ( DVector &  x, DVector &  f  ) const

virtual

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

Definition at line 931 of file EllipseSolver.cc.

    {
        Assert(x.size() == 5);
        Assert(f.size() == 5);
        _pimpl->xinit = x;
        _pimpl->bx.TMV_setAllTo(1.); // In case we bail out with f = 2*bx/b(0)
        if (_pimpl->fixcen) { _pimpl->fixuc = x[0]; _pimpl->fixvc = x[1]; }
        if (_pimpl->fixgam) { _pimpl->fixg1 = x[2]; _pimpl->fixg2 = x[3]; }
        if (_pimpl->fixmu) { _pimpl->fixm = x[4]; }
        if (_pimpl->fixcen && (x[0] != 0. || x[1] != 0.)) { 
            calculateZTransform(
                std::complex<double>(x[0],x[1]),_pimpl->bxorder,_pimpl->Taug);
        }
        if (_pimpl->fixgam && (x[2] != 0. || x[3] != 0.)) { 
            calculateGTransform(
                std::complex<double>(x[2],x[3]),_pimpl->bxorder,_pimpl->Saug);
        }
        if (_pimpl->fixmu && x[4] != 0.) { 
            calculateMuTransform(x[4],_pimpl->bxorder,_pimpl->Daug);
        }

        _pimpl->x_short = _pimpl->U * x;

        bool ret;
        try {
            ret = NLSolver::solve(_pimpl->x_short,_pimpl->f_short);
        } catch (...) {
            xdbg<<"Caught exception during NLSolver::solve"<<std::endl;
            ret = false;
        }

        EIGEN_Transpose(x) = EIGEN_Transpose(_pimpl->x_short)*_pimpl->U;
        if (_pimpl->fixcen) { x[0] = _pimpl->fixuc; x[1] = _pimpl->fixvc; }
        if (_pimpl->fixgam) { x[2] = _pimpl->fixg1; x[3] = _pimpl->fixg2; }
        if (_pimpl->fixmu) { x[4] = _pimpl->fixm; }
        EIGEN_Transpose(f) = EIGEN_Transpose(_pimpl->f_short)*_pimpl->U;

        return ret;
    }

bool lsst::meas::algorithms::shapelet::EllipseSolver3::testJ ( const DVector &  x, DVector &  f, std::ostream *  os = 0, double  relerr = 0.  ) const

virtual

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

Definition at line 971 of file EllipseSolver.cc.

    {
        Assert(x.size() == 5);
        Assert(f.size() == 5);
        _pimpl->xinit = x;
        if (_pimpl->fixcen) { _pimpl->fixuc = x[0]; _pimpl->fixvc = x[1]; }
        if (_pimpl->fixgam) { _pimpl->fixg1 = x[2]; _pimpl->fixg2 = x[3]; }
        if (_pimpl->fixmu) { _pimpl->fixm = x[4]; }

        _pimpl->x_short = _pimpl->U * x;

        return NLSolver::testJ(_pimpl->x_short,_pimpl->f_short,os,relerr);
    }

void lsst::meas::algorithms::shapelet::EllipseSolver3::useNumericJ

(

)

Definition at line 891 of file EllipseSolver.cc.

{ _pimpl->numeric_j = true; }

Member Data Documentation

ESImpl3* lsst::meas::algorithms::shapelet::EllipseSolver3::_pimpl

private

Definition at line 40 of file EllipseSolver.h.

---

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

• /home/lsstsw/stack/Linux64/meas_algorithms/11.0-2-gb8b8ce7/include/lsst/meas/algorithms/shapelet/EllipseSolver.h
• /home/lsstsw/stack/Linux64/meas_algorithms/11.0-2-gb8b8ce7/src/shapelet/EllipseSolver.cc