| /external/eigen/Eigen/src/SparseLU/ |
| SparseLU_column_dfs.h | 43 column_dfs_traits(Index jcol, Index& jsuper, typename SparseLUImpl<Scalar, Index>::GlobalLU_t& glu, SparseLUImpl<Scalar, Index>& luImpl)
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)
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) continue
[all...] |
| 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)
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
92 perm_r(pivrow) = jcol;
[all...] |
| 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);
|
| 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)
54 Index jsupno = glu.supno(jcol);
58 Index nextu = glu.xusub(jcol);
99 glu.xusub(jcol + 1) = nextu; // close U(*,jcol)
|
| 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,
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 + w; jj++
[all...] |
| 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)
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 in snode L\U(jsupno)
[all...] |
| SparseLU_panel_dfs.h | 41 panel_dfs_traits(Index jcol, Index* marker)
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)
227 panel_dfs_traits<IndexVector> traits(jcol, marker1.data());
230 for (Index jj = jcol; jj < jcol + w; jj++)
[all...] |
| 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)
56 Index jsupno = glu.supno(jcol);
|
| SparseLU.h | 589 Index jcol; local
592 Index nseg1; // Number of segments in U-column above panel row jcol
596 for (jcol = 0; jcol < n; )
600 for (k = jcol + 1; k < (std::min)(jcol+panel_size, n); k++)
604 panel_size = k - jcol;
609 panel_size = n - jcol;
612 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);
615 Base::panel_bmod(m, panel_size, jcol, nseg1, dense, tempv, segrep, repfnz, m_glu);
[all...] |
| /external/cblas/testing/ |
| 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; j<*n; j++
274 int irow, jcol, i, j, LDA; local
321 int i,j,irow,jcol,LDA; local
[all...] |
| 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; j<*n; j++
274 int irow, jcol, i, j, LDA; local
321 int i,j,irow,jcol,LDA; local
[all...] |
| 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,LDA local
299 int irow, jcol, i, j, LDA; local
363 int irow, jcol, i, j, LDA; local
[all...] |
| 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,LDA local
299 int irow, jcol, i, j, LDA; local
363 int irow, jcol, i, j, LDA; local
[all...] |
