fix PowerToBezierMatrix bug according to csulijunji.

BIMCoderLiang · BIMCoderLiang · commit 2ebf7fb3ca46 · 2025-12-12T21:00:15.000+08:00
diff --git a/src/LNLib/Algorithm/Polynomials.cpp b/src/LNLib/Algorithm/Polynomials.cpp
@@ -641,14 +641,14 @@ std::vector<std::vector<double>> LNLib::Polynomials::PowerToBezierMatrix(int deg
 
 	for (int k = 1; k < k1; k++)
 	{
-		for (int j = k + 1; j < pk; j++)
+		for (int j = k + 1; j <= pk; j++)
 		{
 			double d = 0.0;
 			for (int i = k; i < j; i++)
 			{
 				d = d - matrix[j][i] * inverseMatrix[i][k];
 			}
-			inverseMatrix[j][k] = d / (inverseMatrix[j][j]);
+			inverseMatrix[j][k] = d / (matrix[j][j]);
 			inverseMatrix[pk][degree - j] = inverseMatrix[j][k];
 		}
 		pk = pk - 1;

Original file line number	Diff line number	Diff line change
`@@ -641,14 +641,14 @@ std::vector<std::vector<double>> LNLib::Polynomials::PowerToBezierMatrix(int deg`
`641`	`641`
`642`	`642`	`for (int k = 1; k < k1; k++)`
`643`	`643`	`{`
`644`		`- for (int j = k + 1; j < pk; j++)`
	`644`	`+ for (int j = k + 1; j <= pk; j++)`
`645`	`645`	`{`
`646`	`646`	`double d = 0.0;`
`647`	`647`	`for (int i = k; i < j; i++)`
`648`	`648`	`{`
`649`	`649`	`d = d - matrix[j][i] * inverseMatrix[i][k];`
`650`	`650`	`}`
`651`		`- inverseMatrix[j][k] = d / (inverseMatrix[j][j]);`
	`651`	`+ inverseMatrix[j][k] = d / (matrix[j][j]);`
`652`	`652`	`inverseMatrix[pk][degree - j] = inverseMatrix[j][k];`
`653`	`653`	`}`
`654`	`654`	`pk = pk - 1;`