CholmodSimplicialLDLT2< _MatrixType, _UpLo > Class Template Reference

#include <Eigenstuff.h>

Inheritance diagram for CholmodSimplicialLDLT2< _MatrixType, _UpLo >:

Public Types
using	MatrixType = _MatrixType
using	Index = typename MatrixType::Index
using	RealScalar = typename MatrixType::RealScalar

Public Member Functions
	CholmodSimplicialLDLT2 ()
	CholmodSimplicialLDLT2 (MatrixType const &matrix)
void	update (SparseMatrixD const &H, bool UpOrDown)

Protected Member Functions
void	init ()

Detailed Description

template<typename _MatrixType, int _UpLo = Eigen::Lower>
class CholmodSimplicialLDLT2< _MatrixType, _UpLo >

Definition at line 50 of file Eigenstuff.h.

Member Typedef Documentation

Index

template<typename _MatrixType, int _UpLo = Eigen::Lower>

using CholmodSimplicialLDLT2< _MatrixType, _UpLo >::Index = typename MatrixType::Index

Definition at line 57 of file Eigenstuff.h.

MatrixType

template<typename _MatrixType, int _UpLo = Eigen::Lower>

using CholmodSimplicialLDLT2< _MatrixType, _UpLo >::MatrixType = _MatrixType

Definition at line 56 of file Eigenstuff.h.

RealScalar

template<typename _MatrixType, int _UpLo = Eigen::Lower>

using CholmodSimplicialLDLT2< _MatrixType, _UpLo >::RealScalar = typename MatrixType::RealScalar

Definition at line 58 of file Eigenstuff.h.

Constructor & Destructor Documentation

CholmodSimplicialLDLT2() [1/2]

template<typename _MatrixType, int _UpLo = Eigen::Lower>

CholmodSimplicialLDLT2< _MatrixType, _UpLo >::CholmodSimplicialLDLT2 ( )

inline

Definition at line 60 of file Eigenstuff.h.

: Base() { init(); }

CholmodSimplicialLDLT2() [2/2]

template<typename _MatrixType, int _UpLo = Eigen::Lower>

CholmodSimplicialLDLT2< _MatrixType, _UpLo >::CholmodSimplicialLDLT2 ( MatrixType const & matrix )

inline

Definition at line 62 of file Eigenstuff.h.

                                                     : Base() {
        init();
        this->compute(matrix);
    }

Member Function Documentation

init()

template<typename _MatrixType, int _UpLo = Eigen::Lower>

void CholmodSimplicialLDLT2< _MatrixType, _UpLo >::init ( )

inlineprotected

Definition at line 90 of file Eigenstuff.h.

                {
        m_cholmod.final_asis = 1;
        m_cholmod.supernodal = CHOLMOD_SIMPLICIAL;
        // In CholmodBase::CholmodBase(), the following statement is missing in
        // SuiteSparse 3.2.0.8. Fixed in 3.2.7
        Base::m_shiftOffset[0] = Base::m_shiftOffset[1] = RealScalar(0.0);
    }

update()

template<typename _MatrixType, int _UpLo = Eigen::Lower>

void CholmodSimplicialLDLT2< _MatrixType, _UpLo >::update ( SparseMatrixD const & H, bool UpOrDown )

inline

Definition at line 68 of file Eigenstuff.h.

                                                       {
        // check size
        Index const size = Base::m_cholmodFactor->n;
        EIGEN_UNUSED_VARIABLE(size);
        eigen_assert(size == H.rows());
 
        cholmod_sparse C_cs = viewAsCholmod(H);
        /* We have to apply the magic permutation to the update matrix,
        read page 117 of Cholmod UserGuide.pdf */
        // Using cholmod_l_* functions instead of cholmod_* because index is Eigen::Index instead of int.
        cholmod_sparse *C_cs_perm =
                cholmod_l_submatrix(&C_cs, (Eigen::Index*)Base::m_cholmodFactor->Perm,
                                  Base::m_cholmodFactor->n, nullptr, -1, true, true, &this->cholmod());
        assert(C_cs_perm);
        int isOk = cholmod_l_updown(UpOrDown, C_cs_perm, Base::m_cholmodFactor, &this->cholmod());
        cholmod_l_free_sparse(&C_cs_perm, &this->cholmod());
        if (!isOk) {
            throw(LSST_EXCEPT(lsst::pex::exceptions::RuntimeError, "cholmod_update failed!"));
        }
    }

---

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

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-12.0.0/Linux/jointcal/g29ae962dfc+1c7d47a24f/include/lsst/jointcal/Eigenstuff.h