lsst::jointcal::TwoTransfoMapping Class Reference

The mapping with two transfos in a row. More...

#include <TwoTransfoMapping.h>

Inheritance diagram for lsst::jointcal::TwoTransfoMapping:

Public Member Functions
	TwoTransfoMapping (std::shared_ptr< SimpleGtransfoMapping > chipMapping, std::shared_ptr< SimpleGtransfoMapping > visitMapping)
	TwoTransfoMapping (TwoTransfoMapping const &)=delete
	No copy or move: there is only ever one instance of a given model (i.e.. per ccd+visit) More...
	TwoTransfoMapping (TwoTransfoMapping &&)=delete
TwoTransfoMapping &	operator= (TwoTransfoMapping const &)=delete
TwoTransfoMapping &	operator= (TwoTransfoMapping &&)=delete
unsigned	getNpar () const
	Number of parameters in total. More...
void	getMappingIndices (std::vector< unsigned > &indices) const
	Sets how this set of parameters (of length Npar()) map into the "grand" fit Expects that indices has enough space reserved. More...
void	computeTransformAndDerivatives (FatPoint const &where, FatPoint &outPoint, Eigen::MatrixX2d &H) const
	Actually applies the AstrometryMapping and evaluates the derivatives w.r.t the fitted parameters. More...
void	transformPosAndErrors (FatPoint const &where, FatPoint &outPoint) const
	The same as above but without the parameter derivatives (used to evaluate chi^2) More...
void	offsetParams (Eigen::VectorXd const &delta)
	Remember the error scale and freeze it. More...
Gtransfo const &	getTransfo1 () const
	access to transfos More...
Gtransfo const &	getTransfo2 () const
	access to transfos More...
void	positionDerivative (Point const &where, Eigen::Matrix2d &derivative, double epsilon) const
	Currently not implemented. More...
void	freezeErrorTransform ()
	Currently not implemented. More...

Friends
class	ConstrainedAstrometryModel

Detailed Description

The mapping with two transfos in a row.

Definition at line 38 of file TwoTransfoMapping.h.

Constructor & Destructor Documentation

TwoTransfoMapping() [1/3]

lsst::jointcal::TwoTransfoMapping::TwoTransfoMapping	(	std::shared_ptr< SimpleGtransfoMapping >	chipMapping,
		std::shared_ptr< SimpleGtransfoMapping >	visitMapping
	)

Definition at line 33 of file TwoTransfoMapping.cc.

        : _m1(chipMapping), _m2(visitMapping) {
    /* Allocate the record of temporary variables, so that they are not
       allocated at every call. This is hidden behind a pointer in order
       to be allowed to alter them in a const routine. */
    tmp = std::unique_ptr<tmpVars>(new tmpVars);
    setWhatToFit(true, true);
}

TwoTransfoMapping() [2/3]

lsst::jointcal::TwoTransfoMapping::TwoTransfoMapping ( TwoTransfoMapping const &  )

delete

No copy or move: there is only ever one instance of a given model (i.e.. per ccd+visit)

TwoTransfoMapping() [3/3]

lsst::jointcal::TwoTransfoMapping::TwoTransfoMapping ( TwoTransfoMapping &&  )

delete

Member Function Documentation

computeTransformAndDerivatives()

void lsst::jointcal::TwoTransfoMapping::computeTransformAndDerivatives ( FatPoint const &  where, FatPoint &  outPoint, Eigen::MatrixX2d &  H  ) const

virtual

Actually applies the AstrometryMapping and evaluates the derivatives w.r.t the fitted parameters.

This is grouped into a single call because for most models, evaluating the derivatives w.r.T parameters is not much longer than just transforming

Implements lsst::jointcal::AstrometryMapping.

Definition at line 62 of file TwoTransfoMapping.cc.

                                                                                {
    // not true in general. Will crash if H is too small.
    //  assert(H.cols()==Npar());

    FatPoint pMid;
    // don't need errors there but no Mapping::Transform() routine.

    if (_nPar1) {
        _m1->computeTransformAndDerivatives(where, pMid, tmp->h1);
        // the last argument is epsilon and is not used for polynomials
        _m2->positionDerivative(pMid, tmp->dt2dx, 1e-4);
        H.block(0, 0, _nPar1, 2) = tmp->h1 * tmp->dt2dx;
    } else
        _m1->transformPosAndErrors(where, pMid);
    if (_nPar2) {
        _m2->computeTransformAndDerivatives(pMid, outPoint, tmp->h2);
        H.block(_nPar1, 0, _nPar2, 2) = tmp->h2;
    } else
        _m2->transformPosAndErrors(pMid, outPoint);
}

freezeErrorTransform()

void lsst::jointcal::TwoTransfoMapping::freezeErrorTransform

(

)

Currently not implemented.

Definition at line 119 of file TwoTransfoMapping.cc.

                                             {
    throw LSST_EXCEPT(
            pexExcept::TypeError,
            " The routine  TwoTransfoMapping::freezeErrorTransform() was thought to be useless and is "
            "not implemented (yet)");
}

getMappingIndices()

void lsst::jointcal::TwoTransfoMapping::getMappingIndices ( std::vector< unsigned > &  indices ) const

virtual

Sets how this set of parameters (of length Npar()) map into the "grand" fit Expects that indices has enough space reserved.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 45 of file TwoTransfoMapping.cc.

                                                                            {
    unsigned npar = getNpar();
    if (indices.size() < npar) indices.resize(npar);
    // in case we are only fitting one of the two transfos
    if (_nPar1)
        _m1->getMappingIndices(indices);
    else if (_nPar2) {
        _m2->getMappingIndices(indices);
        return;
    }
    // if we get here we are fitting both
    // there is probably a more elegant way to feed a subpart of a std::vector
    std::vector<unsigned> ind2(_nPar2);
    _m2->getMappingIndices(ind2);
    for (unsigned k = 0; k < _nPar2; ++k) indices.at(k + _nPar1) = ind2.at(k);
}

getNpar()

unsigned lsst::jointcal::TwoTransfoMapping::getNpar ( ) const

virtual

Number of parameters in total.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 43 of file TwoTransfoMapping.cc.

{ return _nPar1 + _nPar2; }

getTransfo1()

Gtransfo const& lsst::jointcal::TwoTransfoMapping::getTransfo1 ( ) const

inline

access to transfos

Definition at line 72 of file TwoTransfoMapping.h.

{ return _m1->getTransfo(); }

getTransfo2()

Gtransfo const& lsst::jointcal::TwoTransfoMapping::getTransfo2 ( ) const

inline

access to transfos

Definition at line 75 of file TwoTransfoMapping.h.

{ return _m2->getTransfo(); }

offsetParams()

void lsst::jointcal::TwoTransfoMapping::offsetParams ( Eigen::VectorXd const &  delta )

inlinevirtual

Remember the error scale and freeze it.

Note

this routine is not used when fitting (the Model manages the mappings separately), but can be useful for debugging

Implements lsst::jointcal::AstrometryMapping.

Definition at line 66 of file TwoTransfoMapping.h.

                                                  {
        _m1->offsetParams(delta.segment(_m1->getIndex(), _m1->getNpar()));
        _m2->offsetParams(delta.segment(_m2->getIndex() + _m1->getNpar(), _m2->getNpar()));
    }

operator=() [1/2]

TwoTransfoMapping& lsst::jointcal::TwoTransfoMapping::operator= ( TwoTransfoMapping const &  )

delete

operator=() [2/2]

TwoTransfoMapping& lsst::jointcal::TwoTransfoMapping::operator= ( TwoTransfoMapping &&  )

delete

positionDerivative()

void lsst::jointcal::TwoTransfoMapping::positionDerivative ( Point const &  where, Eigen::Matrix2d &  derivative, double  epsilon  ) const

virtual

Currently not implemented.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 107 of file TwoTransfoMapping.cc.

                                                                 {
    Eigen::Matrix2d d1, d2;  // seems that it does not trigger dynamic allocation
    _m1->positionDerivative(where, d1, 1e-4);
    FatPoint pMid;
    _m1->transformPosAndErrors(where, pMid);
    _m2->positionDerivative(pMid, d2, 1e-4);
    /* The following line is not a mistake. It is a consequence
       of chosing derivative(0,1) = d(y_out)/d x_in. */
    derivative = d1 * d2;
}

transformPosAndErrors()

void lsst::jointcal::TwoTransfoMapping::transformPosAndErrors ( FatPoint const &  where, FatPoint &  outPoint  ) const

virtual

The same as above but without the parameter derivatives (used to evaluate chi^2)

Implements lsst::jointcal::AstrometryMapping.

Definition at line 101 of file TwoTransfoMapping.cc.

                                                                                             {
    FatPoint pMid;
    _m1->transformPosAndErrors(where, pMid);
    _m2->transformPosAndErrors(pMid, outPoint);
}

Friends And Related Function Documentation

ConstrainedAstrometryModel

friend class ConstrainedAstrometryModel

friend

Definition at line 84 of file TwoTransfoMapping.h.

---

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

• /j/snowflake/release/lsstsw/stack/Linux64/jointcal/16.0-25-g07af9f2+1/include/lsst/jointcal/TwoTransfoMapping.h
• /j/snowflake/release/lsstsw/stack/Linux64/jointcal/16.0-25-g07af9f2+1/src/TwoTransfoMapping.cc