diff options
| author |  Jeff <complexzeros@gmail.com> | 2014-06-23 19:04:52 -0600 |
|---|---|---|
| committer | Jeff <complexzeros@gmail.com> | 2014-06-23 19:04:52 -0600 |
| commit | b59f045c07a41e79df9e27cd56ce06642a7012e7 (patch) | |
| tree | 99b07c7b9777d1ea0eb040273b3c512db13c1869 /unsupported/Eigen/src/Splines | |
| parent | 957c2c291ba0e2d356584d833432a3314094dfee (diff) | |
Using LU decomposition with complete pivoting for better accuracy.
Diffstat (limited to 'unsupported/Eigen/src/Splines')
| -rw-r--r-- | unsupported/Eigen/src/Splines/SplineFitting.h | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/unsupported/Eigen/src/Splines/SplineFitting.h b/unsupported/Eigen/src/Splines/SplineFitting.h index c30129983..3b8e70a18 100644 --- a/unsupported/Eigen/src/Splines/SplineFitting.h +++ b/unsupported/Eigen/src/Splines/SplineFitting.h @@ -17,6 +17,7 @@ #include "SplineFwd.h" +#include <Eigen/LU> #include <Eigen/QR> namespace Eigen @@ -404,8 +405,8 @@ namespace Eigen A(n - 1, n - 1) = 1; // Solve - HouseholderQR<MatrixType> qr(A); - ControlPointVectorType controlPoints = qr.solve(MatrixType(b.transpose())).transpose(); + FullPivLU<MatrixType> lu(A); + ControlPointVectorType controlPoints = lu.solve(MatrixType(b.transpose())).transpose(); SplineType spline(knots, controlPoints); |