public void IndexerProperty_AtNoneReplenishedGridPoint_BenchmarkResult1(LabelMatrix.MissingValueReplenishment missingValueReplenishment, LabelMatrix.OrderOfInput orderOfInput) { /* some sample data, matrix is provided column-by-column: */ var matrixEntries = new double[] { 1, 2, 3, 4, Double.NaN, 6, 7, 8, 9 }; var xLabels = new double[] { 1.1, 2.12, 3.4 }; var yLabels = new double[] { 1.05, 2.07, 3.65 }; int rowCount = 3; int columnCount = 3; var labelMatrix = LabelMatrix.Create(rowCount, columnCount, matrixEntries, xLabels, yLabels, missingValueReplenishment, orderOfInput); var indexCollection = new Tuple <int, int, int>[] { // (rowIndex, columnIndex, Index of 'matrixEntries') Tuple.Create(0, 0, 0), Tuple.Create(1, 0, 1), Tuple.Create(2, 0, 2), Tuple.Create(0, 1, 3), // skip the 'NaN' value in this test Tuple.Create(2, 1, 5), Tuple.Create(0, 2, 6), Tuple.Create(1, 2, 7), Tuple.Create(2, 2, 8) }; foreach (var indices in indexCollection) { int rowIndex = indices.Item1; int columnIndex = indices.Item2; int expectedMatrixEntryIndex = indices.Item3; double actual = labelMatrix[rowIndex, columnIndex]; double expected = matrixEntries[expectedMatrixEntryIndex]; Assert.That(actual, Is.EqualTo(expected).Within(1E-7), String.Format("Row: {0}, Column: {1}, expected: {2}, actual: {3}", rowIndex, columnIndex, expected, actual)); } }
public void GetValue_TestCase1(GridPointCurve.Interpolator horizontalInterpolator, GridPointCurve.Extrapolator horizontalLeftExtrapolator, GridPointCurve.Extrapolator horizontalRightExtrapolator, GridPointCurve.Interpolator verticalInterpolator, GridPointCurve.Extrapolator verticalAboveExtrapolator, GridPointCurve.Extrapolator verticalBelowExtrapolator, GridPointSurface2d.ConstructionOrder constructionOrder, double x, double y, double expected) { int rowCount = 3; int columnCount = 3; var xLabels = new double[] { 1.1, 2.12, 3.4 }; var yLabels = new double[] { 1.05, 2.07, 3.65 }; double[] matrixEntries = new double[] { 1, 2, 3, 4, 5, 6, 7, 8, 9 }; // matrix is provided column-by-column var labelMatrix = LabelMatrix.Create(rowCount, columnCount, matrixEntries, xLabels, yLabels); var surface = GridPointSurface2d.Create(labelMatrix, horizontalInterpolator, horizontalLeftExtrapolator, horizontalRightExtrapolator, verticalInterpolator, verticalAboveExtrapolator, verticalBelowExtrapolator, constructionOrder); double actual = surface.GetValue(x, y); Assert.That(actual, Is.EqualTo(expected).Within(1E-7), String.Format("expected: {0}, actual: {1}", expected, actual)); }
public void IndexerProperty_yLinear_BenchmarkResult3( [Values(LabelMatrix.OrderOfInput.AscendingOrder, LabelMatrix.OrderOfInput.Disordered, LabelMatrix.OrderOfInput.DisorderedHorizontalLabels, LabelMatrix.OrderOfInput.DisorderedVerticalLabels)] LabelMatrix.OrderOfInput orderOfInput) { var matrixEntries = new double[] { 1, 2, 3, 4, Double.NaN, 6, 7, 8, 9, 10 }; var xLabels = new double[] { 1, 2 }; var yLabels = new double[] { 1, 2, 3, 4, 5 }; int rowCount = 5; int columnCount = 2; var missingValueReplenishment = LabelMatrix.MissingValueReplenishment.WeightedNearestGridPoints.yAxis.Linear; var labelMatrix = LabelMatrix.Create(rowCount, columnCount, matrixEntries, xLabels, yLabels, missingValueReplenishment, orderOfInput); double actual = labelMatrix[4, 0]; double expected = 4.0; // in y-direction only one value is available Assert.That(actual, Is.EqualTo(expected).Within(1E-7), String.Format("Row: {0}, Column: {1}, expected: {2}, actual: {3}", 0, 4, expected, actual)); }
public void IndexerProperty_yLinear_BenchmarkResult1( [Values(LabelMatrix.OrderOfInput.AscendingOrder, LabelMatrix.OrderOfInput.Disordered, LabelMatrix.OrderOfInput.DisorderedHorizontalLabels, LabelMatrix.OrderOfInput.DisorderedVerticalLabels)] LabelMatrix.OrderOfInput orderOfInput) { var matrixEntries = new double[] { 1, 2, 3, 4, Double.NaN, 6, 7, 8, 9 }; var xLabels = new double[] { 1.1, 2.12, 3.4 }; var yLabels = new double[] { 1.05, 2.07, 3.65 }; int rowCount = 3; int columnCount = 3; var missingValueReplenishment = LabelMatrix.MissingValueReplenishment.WeightedNearestGridPoints.yAxis.Linear; var labelMatrix = LabelMatrix.Create(rowCount, columnCount, matrixEntries, xLabels, yLabels, missingValueReplenishment, orderOfInput); double actual = labelMatrix[1, 1]; double expected = (yLabels[1] - yLabels[0]) * (6.0 - 4.0) / (yLabels[2] - yLabels[0]) + 4.0; // a linear interpolation Assert.That(actual, Is.EqualTo(expected).Within(1E-7), String.Format("Row: {0}, Column: {1}, expected: {2}, actual: {3}", 1, 1, expected, actual)); }
public void IndexerProperty_xLinear_BenchmarkResult2( [Values(LabelMatrix.OrderOfInput.AscendingOrder, LabelMatrix.OrderOfInput.Disordered, LabelMatrix.OrderOfInput.DisorderedHorizontalLabels, LabelMatrix.OrderOfInput.DisorderedVerticalLabels)] LabelMatrix.OrderOfInput orderOfInput) { /* some sample data, matrix is provided column-by-column: */ var matrixEntries = new double[] { 1, 2, 3, 4, Double.NaN, 6, 7, 8, 9, 10 }; var xLabels = new double[] { 1, 2, 3, 4, 5 }; var yLabels = new double[] { 1, 2 }; int rowCount = 2; int columnCount = 5; var missingValueReplenishment = LabelMatrix.MissingValueReplenishment.WeightedNearestGridPoints.xAxis.Linear; var labelMatrix = LabelMatrix.Create(rowCount, columnCount, matrixEntries, xLabels, yLabels, missingValueReplenishment, orderOfInput); double actual = labelMatrix[0, 2]; double expected = (xLabels[1] - xLabels[0]) * (7.0 - 3.0) / (xLabels[2] - xLabels[0]) + 3.0; // a linear interpolation Assert.That(actual, Is.EqualTo(expected).Within(1E-7), String.Format("Row: {0}, Column: {1}, expected: {2}, actual: {3}", 0, 2, expected, actual)); }
public void GetValue_AtGridPoint_GridPointValue_TestCase2(GridPointCurve.Interpolator horizontalInterpolator, GridPointCurve.Extrapolator horizontalLeftExtrapolator, GridPointCurve.Extrapolator horizontalRightExtrapolator, GridPointCurve.Interpolator verticalInterpolator, GridPointCurve.Extrapolator verticalAboveExtrapolator, GridPointCurve.Extrapolator verticalBelowExtrapolator, GridPointSurface2d.ConstructionOrder constructionOrder) { int rowCount = 2; int columnCount = 5; var xLabels = new double[] { 1.05, 2.6, 3, 4, 5 }; var yLabels = new double[] { 1.3, 2.8 }; var matrixEntries = new double[] { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 }; // matrix is provided column-by-column var labelMatrix = LabelMatrix.Create(rowCount, columnCount, matrixEntries, xLabels, yLabels); var surface = GridPointSurface2d.Create(labelMatrix, horizontalInterpolator, horizontalLeftExtrapolator, horizontalRightExtrapolator, verticalInterpolator, verticalAboveExtrapolator, verticalBelowExtrapolator, constructionOrder); for (int j = 0; j < rowCount; j++) { for (int k = 0; k < columnCount; k++) { double expected = labelMatrix[j, k]; double actual = surface.GetValue(x: labelMatrix.HorizontalDoubleLabels[k], y: labelMatrix.VerticalDoubleLabels[j]); Assert.That(actual, Is.EqualTo(expected).Within(1E-7), String.Format("Row: {0}, Column: {1}, expected: {2}, actual: {3}", j, k, expected, actual)); } } }