/external/eigen/Eigen/src/SparseLU/ |
H A D | SparseLU_column_dfs.h | 43 column_dfs_traits(Index jcol, Index& jsuper, typename SparseLUImpl<Scalar, Index>::GlobalLU_t& glu, SparseLUImpl<Scalar, Index>& luImpl) argument 44 : m_jcol(jcol), m_jsuper_ref(jsuper), m_glu(glu), m_luImpl(luImpl) 66 * \brief Performs a symbolic factorization on column jcol and decide the supernode boundary 76 * \param jcol Current column 93 Index SparseLUImpl<Scalar,Index>::column_dfs(const Index m, const Index jcol, IndexVector& perm_r, Index maxsuper, Index& nseg, BlockIndexVector lsub_col, IndexVector& segrep, BlockIndexVector repfnz, IndexVector& xprune, IndexVector& marker, IndexVector& parent, IndexVector& xplore, GlobalLU_t& glu) argument 96 Index jsuper = glu.supno(jcol); 97 Index nextl = glu.xlsub(jcol); 101 column_dfs_traits<IndexVector, ScalarVector> traits(jcol, jsuper, glu, *this); 103 // For each nonzero in A(*,jcol) do dfs 111 if (kmark == jcol) continu [all...] |
H A D | SparseLU_pivotL.h | 50 * \param jcol The current column of L 60 Index SparseLUImpl<Scalar,Index>::pivotL(const Index jcol, const RealScalar& diagpivotthresh, IndexVector& perm_r, IndexVector& iperm_c, Index& pivrow, GlobalLU_t& glu) argument 63 Index fsupc = (glu.xsup)((glu.supno)(jcol)); // First column in the supernode containing the column jcol 64 Index nsupc = jcol - fsupc; // Number of columns in the supernode portion, excluding jcol; nsupc >=0 69 Scalar* lu_col_ptr = &(glu.lusup.data()[glu.xlusup(jcol)]); // Start of jcol in the supernode 73 Index diagind = iperm_c(jcol); // diagonal index 91 perm_r(pivrow) = jcol; [all...] |
H A D | SparseLU_copy_to_ucol.h | 38 * \param jcol current column to update 50 Index SparseLUImpl<Scalar,Index>::copy_to_ucol(const Index jcol, const Index nseg, IndexVector& segrep, BlockIndexVector repfnz ,IndexVector& perm_r, BlockScalarVector dense, GlobalLU_t& glu) argument 54 Index jsupno = glu.supno(jcol); 58 Index nextu = glu.xusub(jcol); 99 glu.xusub(jcol + 1) = nextu; // close U(*,jcol)
|
H A D | SparseLU_panel_bmod.h | 45 * \param jcol Starting column of the panel 56 void SparseLUImpl<Scalar,Index>::panel_bmod(const Index m, const Index w, const Index jcol, argument 88 for (jj = jcol; jj < jcol + w; jj++) 90 nextl_col = (jj-jcol) * m; 109 for (jj = jcol; jj < jcol + w; jj++) 111 nextl_col = (jj-jcol) * m; 156 for (jj = jcol; jj < jcol [all...] |
H A D | SparseLUImpl.h | 39 Index snode_dfs(const Index jcol, const Index kcol,const MatrixType& mat, IndexVector& xprune, IndexVector& marker, GlobalLU_t& glu); 40 Index snode_bmod (const Index jcol, const Index fsupc, ScalarVector& dense, GlobalLU_t& glu); 41 Index pivotL(const Index jcol, const RealScalar& diagpivotthresh, IndexVector& perm_r, IndexVector& iperm_c, Index& pivrow, GlobalLU_t& glu); 47 void panel_dfs(const Index m, const Index w, const Index jcol, MatrixType& A, IndexVector& perm_r, Index& nseg, ScalarVector& dense, IndexVector& panel_lsub, IndexVector& segrep, IndexVector& repfnz, IndexVector& xprune, IndexVector& marker, IndexVector& parent, IndexVector& xplore, GlobalLU_t& glu); 49 void panel_bmod(const Index m, const Index w, const Index jcol, const Index nseg, ScalarVector& dense, ScalarVector& tempv, IndexVector& segrep, IndexVector& repfnz, GlobalLU_t& glu); 50 Index column_dfs(const Index m, const Index jcol, IndexVector& perm_r, Index maxsuper, Index& nseg, BlockIndexVector lsub_col, IndexVector& segrep, BlockIndexVector repfnz, IndexVector& xprune, IndexVector& marker, IndexVector& parent, IndexVector& xplore, GlobalLU_t& glu); 51 Index column_bmod(const Index jcol, const Index nseg, BlockScalarVector dense, ScalarVector& tempv, BlockIndexVector segrep, BlockIndexVector repfnz, Index fpanelc, GlobalLU_t& glu); 52 Index copy_to_ucol(const Index jcol, const Index nseg, IndexVector& segrep, BlockIndexVector repfnz ,IndexVector& perm_r, BlockScalarVector dense, GlobalLU_t& glu); 53 void pruneL(const Index jcol, const IndexVector& perm_r, const Index pivrow, const Index nseg, const IndexVector& segrep, BlockIndexVector repfnz, IndexVector& xprune, GlobalLU_t& glu);
|
H A D | SparseLU_column_bmod.h | 40 * \param jcol current column to update 53 Index SparseLUImpl<Scalar,Index>::column_bmod(const Index jcol, const Index nseg, BlockScalarVector dense, ScalarVector& tempv, BlockIndexVector segrep, BlockIndexVector repfnz, Index fpanelc, GlobalLU_t& glu) argument 67 jsupno = glu.supno(jcol); 112 nextlu = glu.xlusup(jcol); 140 glu.xlusup(jcol + 1) = nextlu; // close L\U(*,jcol); 150 if (fst_col < jcol) 159 nsupc = jcol - fst_col; // excluding jcol 162 // points to the beginning of jcol i [all...] |
H A D | SparseLU_panel_dfs.h | 41 panel_dfs_traits(Index jcol, Index* marker) argument 42 : m_jcol(jcol), m_marker(marker) 161 //if (marker1(krep) < jcol ) 183 * \brief Performs a symbolic factorization on a panel of columns [jcol, jcol+w) 198 * marker1[i] >= jcol, if i was visited by earlier columns in this panel; 202 * \param[in] jcol Starting column of the panel 219 void SparseLUImpl<Scalar,Index>::panel_dfs(const Index m, const Index w, const Index jcol, MatrixType& A, IndexVector& perm_r, Index& nseg, ScalarVector& dense, IndexVector& panel_lsub, IndexVector& segrep, IndexVector& repfnz, IndexVector& xprune, IndexVector& marker, IndexVector& parent, IndexVector& xplore, GlobalLU_t& glu) argument 227 panel_dfs_traits<IndexVector> traits(jcol, marker1.data()); 230 for (Index jj = jcol; j [all...] |
H A D | SparseLU_pruneL.h | 42 * \param jcol The current column of L 53 void SparseLUImpl<Scalar,Index>::pruneL(const Index jcol, const IndexVector& perm_r, const Index pivrow, const Index nseg, const IndexVector& segrep, BlockIndexVector repfnz, IndexVector& xprune, GlobalLU_t& glu) argument 56 Index jsupno = glu.supno(jcol);
|
H A D | SparseLU.h | 544 Index jcol; local 547 Index nseg1; // Number of segments in U-column above panel row jcol 551 for (jcol = 0; jcol < n; ) 555 for (k = jcol + 1; k < (std::min)(jcol+panel_size, n); k++) 559 panel_size = k - jcol; 564 panel_size = n - jcol; 567 Base::panel_dfs(m, panel_size, jcol, m_mat, m_perm_r.indices(), nseg1, dense, panel_lsub, segrep, repfnz, xprune, marker, parent, xplore, m_glu); 570 Base::panel_bmod(m, panel_size, jcol, nseg [all...] |
/external/cblas/testing/ |
H A D | c_dblas2.c | 192 int i,irow,j,jcol,LDA; local 202 jcol=(*ku)-i; 203 for( j=jcol; j<*n; j++ ) 204 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ]; 212 jcol=i-(*ku); 213 for( j=jcol; j<(*n+*kl); j++ ) 214 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ]; 228 int irow, jcol, i, j, LDA; local 243 jcol=(*k)-i; 244 for( j=jcol; 274 int irow, jcol, i, j, LDA; local 321 int i,j,irow,jcol,LDA; local [all...] |
H A D | c_sblas2.c | 192 int i,irow,j,jcol,LDA; local 202 jcol=(*ku)-i; 203 for( j=jcol; j<*n; j++ ) 204 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ]; 212 jcol=i-(*ku); 213 for( j=jcol; j<(*n+*kl); j++ ) 214 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ]; 228 int irow, jcol, i, j, LDA; local 243 jcol=(*k)-i; 244 for( j=jcol; 274 int irow, jcol, i, j, LDA; local 321 int i,j,irow,jcol,LDA; local [all...] |
H A D | c_cblas2.c | 47 int i,j,irow,jcol,LDA; local 56 jcol=(*ku)-i; 57 for( j=jcol; j<*n; j++ ){ 58 A[ LDA*(j-jcol)+irow ].real=a[ (*lda)*j+i ].real; 59 A[ LDA*(j-jcol)+irow ].imag=a[ (*lda)*j+i ].imag; 70 jcol=i-(*ku); 71 for( j=jcol; j<(*n+*kl); j++ ){ 72 A[ LDA*j+irow ].real=a[ (*lda)*(j-jcol)+i ].real; 73 A[ LDA*j+irow ].imag=a[ (*lda)*(j-jcol)+i ].imag; 181 int i,irow,j,jcol,LD local 299 int irow, jcol, i, j, LDA; local 363 int irow, jcol, i, j, LDA; local [all...] |
H A D | c_zblas2.c | 47 int i,j,irow,jcol,LDA; local 56 jcol=(*ku)-i; 57 for( j=jcol; j<*n; j++ ){ 58 A[ LDA*(j-jcol)+irow ].real=a[ (*lda)*j+i ].real; 59 A[ LDA*(j-jcol)+irow ].imag=a[ (*lda)*j+i ].imag; 70 jcol=i-(*ku); 71 for( j=jcol; j<(*n+*kl); j++ ){ 72 A[ LDA*j+irow ].real=a[ (*lda)*(j-jcol)+i ].real; 73 A[ LDA*j+irow ].imag=a[ (*lda)*(j-jcol)+i ].imag; 181 int i,irow,j,jcol,LD local 299 int irow, jcol, i, j, LDA; local 363 int irow, jcol, i, j, LDA; local [all...] |