1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
/*
File : SmoothTest.cpp
Project : LabPlot
Description : Tests for data smoothing
--------------------------------------------------------------------
SPDX-FileCopyrightText: 2020 Stefan Gerlach <stefan.gerlach@uni.kn>
SPDX-License-Identifier: GPL-2.0-or-later
*/
#include "SmoothTest.h"
#include "backend/core/column/Column.h"
#include "backend/worksheet/plots/cartesian/XYSmoothCurve.h"
extern "C" {
#include "backend/nsl/nsl_smooth.h"
}
// ##############################################################################
void SmoothTest::testPercentile() {
// data
QVector<int> xData{1, 2, 3, 4, 5, 6, 7, 8, 9, 10};
QVector<double> yData{47.7, 44., 43., 44.96, 45., 43.73, 38., 47.1, 44., 38.3};
// p=0.5, d=5 (check with LibreOffice and Origin 2020a)
QVector<double> result50_5{47.7, 44, 44.96, 44, 43.73, 44.96, 44, 43.73, 44, 38.3};
// data source columns
Column xDataColumn(QStringLiteral("x"), AbstractColumn::ColumnMode::Integer);
xDataColumn.replaceInteger(0, xData);
Column yDataColumn(QStringLiteral("y"), AbstractColumn::ColumnMode::Double);
yDataColumn.replaceValues(0, yData);
XYSmoothCurve smoothCurve(QStringLiteral("smooth"));
smoothCurve.setXDataColumn(&xDataColumn);
smoothCurve.setYDataColumn(&yDataColumn);
// prepare the smooth
XYSmoothCurve::SmoothData smoothData = smoothCurve.smoothData();
smoothData.type = nsl_smooth_type_percentile;
// default
// smoothData.points = 5;
// smoothData.percentile = 0.5;
smoothCurve.setSmoothData(smoothData);
// perform the smooth
smoothCurve.recalculate();
const auto& smoothResult = smoothCurve.result();
// check the results
QCOMPARE(smoothResult.available, true);
QCOMPARE(smoothResult.valid, true);
const auto* resultXDataColumn{smoothCurve.xColumn()};
const auto* resultYDataColumn{smoothCurve.yColumn()};
const int np{resultXDataColumn->rowCount()};
QCOMPARE(np, 10);
for (int i = 0; i < np; i++) {
QCOMPARE(resultXDataColumn->valueAt(i), (double)i + 1);
QCOMPARE(resultYDataColumn->valueAt(i), result50_5.at(i));
}
}
QTEST_MAIN(SmoothTest)
|
