As homogeneous_linear_solver() modifies the matrix which is passed as reference

we make sure that rank() is called on the modified matrix.
This commit is contained in:
Andreas Fabri 2013-09-26 16:16:32 +02:00
parent d4c3378e44
commit 5f390a0cd1
1 changed files with 4 additions and 2 deletions

View File

@ -202,7 +202,8 @@ int main(int argc, char* argv[])
Matrix SV;
if (LA::homogeneous_linear_solver(E,x)) {
CGAL_TEST(E*x==Vector(mat_dim));
CGAL_TEST(LA::homogeneous_linear_solver(E,SV)==LA::rank(SV));
int r = LA::homogeneous_linear_solver(E,SV);
CGAL_TEST(r==LA::rank(SV));
}
// inverse
@ -399,7 +400,8 @@ int main(int argc, char* argv[])
Matrix SV;
if (LA::homogeneous_linear_solver(E,x)) {
CGAL_TEST(E*x==Vector(mat_dim));
CGAL_TEST(LA::homogeneous_linear_solver(E,SV)==LA::rank(SV));
int r = LA::homogeneous_linear_solver(E,SV);
CGAL_TEST(r==LA::rank(SV));
}
// inverse
if (LA::inverse(E,D,denom,c)) {