GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-06-04 18:54:00 +08:00
Code Issues Packages Projects Releases Wiki Activity

update to fix compilation

Browse Source
This commit is contained in:
Benoit Jacob 2008-03-16 21:04:33 +00:00
parent 612350e3f8
commit af131fe770
1 changed files with 9 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
doc/echelon.cpp
View File

@ -4,8 +4,8 @@ USING_PART_OF_NAMESPACE_EIGEN
namespace Eigen { namespace Eigen {
template<typename Scalar, typename Derived> template<typename Derived>
void echelon(MatrixBase<Scalar, Derived>& m) void echelon(MatrixBase<Derived>& m)
{ {
const int N = std::min(m.rows(), m.cols()); const int N = std::min(m.rows(), m.cols());
@ -23,22 +23,22 @@ void echelon(MatrixBase<Scalar, Derived>& m)
} }
} }
template<typename Scalar, typename Derived> template<typename Derived>
void doSomeRankPreservingOperations(MatrixBase<Scalar, Derived>& m) void doSomeRankPreservingOperations(MatrixBase<Derived>& m)
{ {
for(int a = 0; a < 3*(m.rows()+m.cols()); a++) for(int a = 0; a < 3*(m.rows()+m.cols()); a++)
{ {
double d = Eigen::random<double>(-1,1); double d = Eigen::ei_random<double>(-1,1);
int i = Eigen::random<int>(0,m.rows()-1); // i is a random row number int i = Eigen::ei_random<int>(0,m.rows()-1); // i is a random row number
int j; int j;
do { do {
j = Eigen::random<int>(0,m.rows()-1); j = Eigen::ei_random<int>(0,m.rows()-1);
} while (i==j); // j is another one (must be different) } while (i==j); // j is another one (must be different)
m.row(i) += d * m.row(j); m.row(i) += d * m.row(j);
i = Eigen::random<int>(0,m.cols()-1); // i is a random column number i = Eigen::ei_random<int>(0,m.cols()-1); // i is a random column number
do { do {
j = Eigen::random<int>(0,m.cols()-1); j = Eigen::ei_random<int>(0,m.cols()-1);
} while (i==j); // j is another one (must be different) } while (i==j); // j is another one (must be different)
m.col(i) += d * m.col(j); m.col(i) += d * m.col(j);
} }