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

Properly initialize b vector in SplineFitting

Browse Source 
InterpolateWithDerivative does not initialize the be vector correctly. This issue is discussed In stackoverflow question 48382939.
This commit is contained in:
Matthew Powelson 2020-01-07 16:20:43 -06:00 committed by Rasmus Munk Larsen
parent 9254974115
commit 2ea5a715cf
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
unsupported/Eigen/src/Splines/SplineFitting.h
View File

@ -385,7 +385,7 @@ namespace Eigen
{ {
const DenseIndex span = SplineType::Span(parameters[i], degree, knots); const DenseIndex span = SplineType::Span(parameters[i], degree, knots);
if (derivativeIndices[derivativeIndex] == i) if (derivativeIndex < derivativeIndices.size() && derivativeIndices[derivativeIndex] == i)
{ {
A.block(row, span - degree, 2, degree + 1) A.block(row, span - degree, 2, degree + 1)
= SplineType::BasisFunctionDerivatives(parameters[i], 1, degree, knots); = SplineType::BasisFunctionDerivatives(parameters[i], 1, degree, knots);
@ -395,8 +395,9 @@ namespace Eigen
} }
else else
{ {
A.row(row++).segment(span - degree, degree + 1) A.row(row).segment(span - degree, degree + 1)
= SplineType::BasisFunctions(parameters[i], degree, knots); = SplineType::BasisFunctions(parameters[i], degree, knots);
b.col(row++) = points.col(i);
} }
} }
b.col(0) = points.col(0); b.col(0) = points.col(0);