BunchKaufman-methods         package:Matrix         R Documentation

_B_u_n_c_h-_K_a_u_f_m_a_n _D_e_c_o_m_p_o_s_i_t_i_o_n _M_e_t_h_o_d_s

_D_e_s_c_r_i_p_t_i_o_n:

     The Bunch-Kaufman Decomposition of a square symmetric matrix A is
     A = P LDL' P' where P is a permutation matrix, L is _unit_-lower
     triangular and D is _block_-diagonal with blocks of dimension 1 x
     1 or 2 x 2.

_U_s_a_g_e:

     BunchKaufman(x, ...)

_A_r_g_u_m_e_n_t_s:

       x: a symmetric square matrix.

     ...: potentially further arguments passed to methods.

_M_e_t_h_o_d_s:

     Currently, only methods for *dense* numeric symmetric matrices are
     implemented.

     '_x = "_d_s_p_M_a_t_r_i_x"' ...

     '_x = "_d_s_y_M_a_t_r_i_x"' ...

     These rely on the Lapack routines 'dsptrf' and 'dsytrf',
     respectively.

_S_e_e _A_l_s_o:

     The resulting class, 'BunchKaufman'. Related decompositions are
     the LU,  'lu', and the Cholesky, 'chol' (and for _sparse_
     matrices, 'Cholesky').

_E_x_a_m_p_l_e_s:

     data(CAex)
     dim(CAex)
     isSymmetric(CAex)# TRUE
     CAs <- as(CAex, "symmetricMatrix")
     if(FALSE) # no method defined yet for *sparse* :
        bk. <- BunchKaufman(CAs)
     ## does apply to *dense* symmetric matrices:
     bkCA <- BunchKaufman(as(CAs, "denseMatrix"))
     bkCA

