small modif

parent f569c8f9
......@@ -9,7 +9,7 @@ namespace tinygb {
/* Constant values
*****************************************************************************/
// 18 is such that the proba that two randomly generated string with 18 lower
// 18 is such that the proba that two randomly generated strings with 18 lower
// case characters have a probability < 2^(-80) to collide. This is due to the
// inequality 26^18 > 2^80.
......
......@@ -116,7 +116,7 @@ void RowMajMat::MultiplyByInverse(const SpMat& A) {
nb_blocks++;
else
rem_block = SIZE_BLOCK_SPARSE_INV;
// This way, the last block has always size rem_block.
// the last block has always size rem_block.
std::size_t total_nb_blocks = ((nb_blocks*(nb_blocks+1))/2);
entries_blocks.resize(total_nb_blocks);
......@@ -152,6 +152,7 @@ void RowMajMat::MultiplyByInverse(const SpMat& A) {
(it.j_%SIZE_BLOCK_SPARSE_INV)],
it.v_);
// ftrtri computes the inverse of a triangular matrix
FFPACK::ftrtri(GivaroWrapper::k, FFLAS::FflasUpper, FFLAS::FflasUnit,
(ind == 0)?rem_block:SIZE_BLOCK_SPARSE_INV,
diag_block, SIZE_BLOCK_SPARSE_INV);
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment