songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Use CGAL::sqrt() rather than std::sqrt()

This commit is contained in:
andrew morton 2015-04-06 22:44:27 -07:00
parent fa281d882f
commit 21295f07aa
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
Principal_component_analysis/include/CGAL/eigen_2.h
View File

@ -86,8 +86,8 @@ namespace internal {
}
else // generic case
{
FT l1 = (FT)(0.5 * ( -1*std::sqrt(p) + c + a));
FT l2 = (FT)(0.5 * ( std::sqrt(p) + c + a));
FT l1 = (FT)(0.5 * ( -1*CGAL::sqrt(p) + c + a));
FT l2 = (FT)(0.5 * ( CGAL::sqrt(p) + c + a));
// all eigen values of a symmetric positive
// definite matrix must be real and positive
@ -101,15 +101,15 @@ namespace internal {
{
eigen_values.first = l1;
eigen_values.second = l2;
eigen_vectors.first = Vector((FT)1.0, (FT)(-(std::sqrt(p)-c+a) / (2*b)));
eigen_vectors.second = Vector((FT)1.0, (FT)( (std::sqrt(p)+c-a) / (2*b)));
eigen_vectors.first = Vector((FT)1.0, (FT)(-(CGAL::sqrt(p)-c+a) / (2*b)));
eigen_vectors.second = Vector((FT)1.0, (FT)( (CGAL::sqrt(p)+c-a) / (2*b)));
}
else
{
eigen_values.first = l2;
eigen_values.second = l1;
eigen_vectors.first = Vector((FT)1.0, (FT)( (std::sqrt(p)+c-a) / (2*b)));
eigen_vectors.second = Vector((FT)1.0, (FT)(-(std::sqrt(p)-c+a) / (2*b)));
eigen_vectors.first = Vector((FT)1.0, (FT)( (CGAL::sqrt(p)+c-a) / (2*b)));
eigen_vectors.second = Vector((FT)1.0, (FT)(-(CGAL::sqrt(p)-c+a) / (2*b)));
}
} // end generic case
} // end non-degenerate case