Esempio n. 1
0
        internal Tuple <double[], double[]> FunctionApproximationControlPointsAndAuxiliaryKnots(int degree, double[] pointsX, double[] pointsY, double[] knots, int[] knotsIndexes)
        {
            Tuple <double[][], double[]> matrixAndAuxKnots = FunctionApproximationMatrix(pointsX, knots, knotsIndexes, degree);

            double[][]        matrix          = matrixAndAuxKnots.Item1;
            double[][]        transposeMatrix = MatrixOperations.TransposeMatrix(matrix);
            double[][]        transposeAndMatrixMultiplicationInversion = MatrixOperations.MatrixInvert(MatrixOperations.MultiplyMatrices(transposeMatrix, matrix));
            double[][]        YMatrix = MatrixOperations.ArrayToMatrix(pointsY);
            double[][]        finalMatrix;
            double[]          finalMatrixAsArray;
            List <double[][]> matrixList = new List <double[][]>();

            matrixList.Add(transposeAndMatrixMultiplicationInversion);
            matrixList.Add(transposeMatrix);
            matrixList.Add(YMatrix);
            finalMatrix = MatrixOperations.MultiplyListOfMatrices(matrixList);
            try
            {
                finalMatrixAsArray = MatrixOperations.OneColumnMatrixToArray(finalMatrix);
            }
            catch (NumberOfMatrixColumnIsNotOneException)
            {
                System.Diagnostics.Debug.WriteLine("Matrix cannot be converted to an array because number of columns is not one");
                return(null);
            }

            return(new Tuple <double[], double[]>(finalMatrixAsArray, matrixAndAuxKnots.Item2));
        }
Esempio n. 2
0
        internal Tuple <double[], double[]> FunctionControlPointsAndAuxiliaryKnots(int degree, double[] knots, double[] knotsFunctionValues, bool areKnotsService)
        {
            var matrixAndAuxKnots         = FunctionBSplineMatrix(degree, knots, areKnotsService);
            var invertedMatrix            = MatrixOperations.MatrixInvert(matrixAndAuxKnots.Item1);
            var knotsFunctionValuesMatrix = MatrixOperations.ArrayToMatrix(knotsFunctionValues);
            var controlPointsMatrix       = MatrixOperations.MultiplyMatrices(invertedMatrix, knotsFunctionValuesMatrix);

            double[] controlPoints = null;
            try
            {
                controlPoints = MatrixOperations.OneColumnMatrixToArray(controlPointsMatrix);
            }
            catch (NumberOfMatrixColumnIsNotOneException)
            {
                System.Diagnostics.Debug.WriteLine("Matrix cannot be converted to an array because number of columns is not one");
                return(null);
            }
            return(new Tuple <double[], double[]>(controlPoints, matrixAndAuxKnots.Item2));
        }