boost/ublas
2
0
Fork 0
mirror of https://github.com/boostorg/ublas.git synced 2026-02-22 03:42:19 +00:00
Code Issues Projects Releases Wiki Activity

lu.hpp: introduced temporary to avoid strange compiler problem.

Browse Source
This commit is contained in:
Gunter Winkler
2007-09-12 21:44:37 +00:00
parent 78d07714d6
commit e7840acdcd
1 changed files with 4 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
include/boost/numeric/ublas/lu.hpp
View File

@@ -101,7 +101,8 @@ namespace boost { namespace numeric { namespace ublas {
matrix_column<M> mci (column (m, i));
matrix_row<M> mri (row (m, i));
if (m (i, i) != value_type/*zero*/()) {
project (mci, range (i + 1, size1)) *= value_type (1) / m (i, i);
value_type m_inv = value_type (1) / m (i, i);
project (mci, range (i + 1, size1)) *= m_inv;
} else if (singular == 0) {
singular = i + 1;
}
@@ -144,7 +145,8 @@ namespace boost { namespace numeric { namespace ublas {
} else {
BOOST_UBLAS_CHECK (pm (i) == i_norm_inf, external_logic ());
}
project (mci, range (i + 1, size1)) *= value_type (1) / m (i, i);
value_type m_inv = value_type (1) / m (i, i);
project (mci, range (i + 1, size1)) *= m_inv;
} else if (singular == 0) {
singular = i + 1;
}