lsst::afw::math::LeastSquares::Impl Class Reference

Public Types
enum	StateFlags {   LOWER_FISHER_MATRIX = 0x001, FULL_FISHER_MATRIX = 0x002, RHS_VECTOR = 0x004, SOLUTION_ARRAY = 0x008,   COVARIANCE_ARRAY = 0x010, DIAGNOSTIC_ARRAY = 0x020, DESIGN_AND_DATA = 0x040 }

Public Member Functions
template<typename D >
void	setRank (Eigen::MatrixBase< D > const &values)
void	ensure (int desired)
virtual void	factor ()=0
virtual void	updateRank ()=0
virtual void	updateSolution ()=0
virtual void	updateCovariance ()=0
virtual void	updateDiagnostic ()=0
	Impl (int dimension_, double threshold_=std::numeric_limits< double >::epsilon())
virtual	~Impl ()

Public Attributes
int	state
int	dimension
int	rank
Factorization	factorization
Factorization	whichDiagnostic
double	threshold
Eigen::MatrixXd	design
Eigen::VectorXd	data
Eigen::MatrixXd	fisher
Eigen::VectorXd	rhs
ndarray::Array< double, 1, 1 >	solution
ndarray::Array< double, 2, 2 >	covariance
ndarray::Array< double, 1, 1 >	diagnostic
pex::logging::Debug	log

Detailed Description

Definition at line 37 of file LeastSquares.cc.

Member Enumeration Documentation

enum lsst::afw::math::LeastSquares::Impl::StateFlags

Enumerator
LOWER_FISHER_MATRIX
FULL_FISHER_MATRIX
RHS_VECTOR
SOLUTION_ARRAY
COVARIANCE_ARRAY
DIAGNOSTIC_ARRAY
DESIGN_AND_DATA

Definition at line 40 of file LeastSquares.cc.

                    { 
        LOWER_FISHER_MATRIX     = 0x001,
        FULL_FISHER_MATRIX      = 0x002,
        RHS_VECTOR              = 0x004,
        SOLUTION_ARRAY          = 0x008,
        COVARIANCE_ARRAY        = 0x010,
        DIAGNOSTIC_ARRAY        = 0x020,
        DESIGN_AND_DATA         = 0x040
    };

Constructor & Destructor Documentation

lsst::afw::math::LeastSquares::Impl::Impl ( int  dimension_, double  threshold_ = std::numeric_limits<double>::epsilon()  )

inline

Definition at line 118 of file LeastSquares.cc.

                                                                               : 
        state(0), dimension(dimension_), rank(dimension_), threshold(threshold_), 
        log("afw.math.LeastSquares")
        {}

virtual lsst::afw::math::LeastSquares::Impl::~Impl ( )

inlinevirtual

Definition at line 123 of file LeastSquares.cc.

{}

Member Function Documentation

void lsst::afw::math::LeastSquares::Impl::ensure ( int  desired )

inline

Definition at line 78 of file LeastSquares.cc.

                             {
        if (state & FULL_FISHER_MATRIX) state |= LOWER_FISHER_MATRIX;
        if (desired & FULL_FISHER_MATRIX) desired |= LOWER_FISHER_MATRIX;
        int toAdd = ~state & desired;
        if (toAdd & LOWER_FISHER_MATRIX) {
            assert(state & DESIGN_AND_DATA);
            fisher = Eigen::MatrixXd::Zero(design.cols(), design.cols());
            fisher.selfadjointView<Eigen::Lower>().rankUpdate(design.adjoint());
        }
        if (toAdd & FULL_FISHER_MATRIX) {
            fisher.triangularView<Eigen::StrictlyUpper>() = fisher.adjoint();
        }
        if (toAdd & RHS_VECTOR) {
            assert(state & DESIGN_AND_DATA);
            rhs = design.adjoint() * data;
        }
        if (toAdd & SOLUTION_ARRAY) {
            if (solution.isEmpty()) solution = ndarray::allocate(dimension);
            updateSolution();
        }
        if (toAdd & COVARIANCE_ARRAY) {
            if (covariance.isEmpty()) covariance = ndarray::allocate(dimension, dimension);
            updateCovariance();
        }
        if (toAdd & DIAGNOSTIC_ARRAY) {
            if (diagnostic.isEmpty()) diagnostic = ndarray::allocate(dimension);
            updateDiagnostic();
        }
        state |= toAdd;
    }

virtual void lsst::afw::math::LeastSquares::Impl::factor ( )

pure virtual

template<typename D >

void lsst::afw::math::LeastSquares::Impl::setRank ( Eigen::MatrixBase< D > const &  values )

inline

Definition at line 69 of file LeastSquares.cc.

                                                    {
        double cond = threshold * values[0];
        if (cond <= 0.0) {
            rank = 0;
        } else {
            for (rank = dimension; (rank > 1) && (values[rank-1] < cond); --rank);
        }
    }

virtual void lsst::afw::math::LeastSquares::Impl::updateCovariance ( )

pure virtual

virtual void lsst::afw::math::LeastSquares::Impl::updateDiagnostic ( )

pure virtual

virtual void lsst::afw::math::LeastSquares::Impl::updateRank ( )

pure virtual

virtual void lsst::afw::math::LeastSquares::Impl::updateSolution ( )

pure virtual

Member Data Documentation

ndarray::Array<double,2,2> lsst::afw::math::LeastSquares::Impl::covariance

Definition at line 63 of file LeastSquares.cc.

Eigen::VectorXd lsst::afw::math::LeastSquares::Impl::data

Definition at line 58 of file LeastSquares.cc.

Eigen::MatrixXd lsst::afw::math::LeastSquares::Impl::design

Definition at line 57 of file LeastSquares.cc.

ndarray::Array<double,1,1> lsst::afw::math::LeastSquares::Impl::diagnostic

Definition at line 64 of file LeastSquares.cc.

int lsst::afw::math::LeastSquares::Impl::dimension

Definition at line 51 of file LeastSquares.cc.

Factorization lsst::afw::math::LeastSquares::Impl::factorization

Definition at line 53 of file LeastSquares.cc.

Eigen::MatrixXd lsst::afw::math::LeastSquares::Impl::fisher

Definition at line 59 of file LeastSquares.cc.

pex::logging::Debug lsst::afw::math::LeastSquares::Impl::log

Definition at line 66 of file LeastSquares.cc.

int lsst::afw::math::LeastSquares::Impl::rank

Definition at line 52 of file LeastSquares.cc.

Eigen::VectorXd lsst::afw::math::LeastSquares::Impl::rhs

Definition at line 60 of file LeastSquares.cc.

ndarray::Array<double,1,1> lsst::afw::math::LeastSquares::Impl::solution

Definition at line 62 of file LeastSquares.cc.

int lsst::afw::math::LeastSquares::Impl::state

Definition at line 50 of file LeastSquares.cc.

double lsst::afw::math::LeastSquares::Impl::threshold

Definition at line 55 of file LeastSquares.cc.

Factorization lsst::afw::math::LeastSquares::Impl::whichDiagnostic

Definition at line 54 of file LeastSquares.cc.

---

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

• /home/lsstsw/stack/Linux64/afw/11.0-2-g04d2804/src/math/LeastSquares.cc