diff options
Diffstat (limited to 'Eigen/src/SparseLU/SparseLU_column_bmod.h')
| -rw-r--r-- | Eigen/src/SparseLU/SparseLU_column_bmod.h | 9 |
1 files changed, 2 insertions, 7 deletions
diff --git a/Eigen/src/SparseLU/SparseLU_column_bmod.h b/Eigen/src/SparseLU/SparseLU_column_bmod.h index bf25a33fc..b268c4348 100644 --- a/Eigen/src/SparseLU/SparseLU_column_bmod.h +++ b/Eigen/src/SparseLU/SparseLU_column_bmod.h @@ -46,11 +46,9 @@ * > 0 - number of bytes allocated when run out of space * */ -template <typename IndexVector, typename ScalarVector, typename BlockIndexVector, typename BlockScalarVector> -int LU_column_bmod(const int jcol, const int nseg, BlockScalarVector& dense, ScalarVector& tempv, BlockIndexVector& segrep, BlockIndexVector& repfnz, int fpanelc, LU_GlobalLU_t<IndexVector, ScalarVector>& glu) +template <typename Scalar, typename Index> +int SparseLUBase<Scalar,Index>::LU_column_bmod(const int jcol, const int nseg, BlockScalarVector& dense, ScalarVector& tempv, BlockIndexVector& segrep, BlockIndexVector& repfnz, int fpanelc, GlobalLU_t& glu) { - typedef typename IndexVector::Scalar Index; - typedef typename ScalarVector::Scalar Scalar; int jsupno, k, ksub, krep, ksupno; int lptr, nrow, isub, irow, nextlu, new_next, ufirst; int fsupc, nsupc, nsupr, luptr, kfnz, no_zeros; @@ -95,9 +93,6 @@ int LU_column_bmod(const int jcol, const int nseg, BlockScalarVector& dense, Sca nsupr = glu.xlsub(fsupc+1) - glu.xlsub(fsupc); nrow = nsupr - d_fsupc - nsupc; - // NOTE Unlike the original implementation in SuperLU, the only feature - // available here is a sup-col update. - // Perform a triangular solver and block update, // then scatter the result of sup-col update to dense no_zeros = kfnz - fst_col; |