Method testLinearWampler1_custom

tests/analysis/fit/FitTest.cpp:532–591

}

void FitTest::testLinearWampler1_custom() {
	// NIST data for Wampler1 dataset
	QVector<int> xData = {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20};
	QVector<int> yData = {1,	  6,	  63,	  364,	  1365,	  3906,	   9331,	19608,	 37449,	  66430,  111111,
						  177156, 271453, 402234, 579195, 813616, 1118481, 1508598, 2000719, 2613660, 3368421};

	// data source columns
	Column xDataColumn(QStringLiteral("x"), AbstractColumn::ColumnMode::Integer);
	xDataColumn.replaceInteger(0, xData);

	Column yDataColumn(QStringLiteral("y"), AbstractColumn::ColumnMode::Integer);
	yDataColumn.replaceInteger(0, yData);

	XYFitCurve fitCurve(QStringLiteral("fit"));
	fitCurve.setXDataColumn(&xDataColumn);
	fitCurve.setYDataColumn(&yDataColumn);

	// prepare the fit
	XYFitCurve::FitData fitData = fitCurve.fitData();
	fitData.modelCategory = nsl_fit_model_custom;
	XYFitCurve::initFitData(fitData);
	fitData.model = QStringLiteral("B0+B1*x+B2*x^2+B3*x^3+B4*x^4+B5*x^5");
	fitData.paramNames << QStringLiteral("B0") << QStringLiteral("B1") << QStringLiteral("B2") << QStringLiteral("B3") << QStringLiteral("B4")
					   << QStringLiteral("B5");
	// fitData.eps = 1.e-12;
	const int np = fitData.paramNames.size();
	// start values
	fitData.paramStartValues << 1. << 1. << 1. << 1. << 1. << 1.;
	for (int i = 0; i < np; i++) {
		fitData.paramLowerLimits << -std::numeric_limits<double>::max();
		fitData.paramUpperLimits << std::numeric_limits<double>::max();
	}
	fitCurve.setFitData(fitData);

	// perform the fit
	fitCurve.recalculate();
	const XYFitCurve::FitResult& fitResult = fitCurve.fitResult();

	// check the results
	QCOMPARE(fitResult.available, true);
	QCOMPARE(fitResult.valid, true);

	QCOMPARE(np, 6);

	for (int i = 0; i < np; i++) {
		const double paramValue = fitResult.paramValues.at(i);
		const double errorValue = fitResult.errorValues.at(i);
		DEBUG(std::setprecision(15) << paramValue << ' ' << paramValue - 1.0);
		FuzzyCompare(paramValue, 1.0, 2.e-12);
		DEBUG(std::setprecision(15) << errorValue);
		FuzzyCompare(errorValue, 0.0, 5.e-12);
	}

	DEBUG(std::setprecision(15) << fitResult.rsd);
	FuzzyCompare(fitResult.rsd, 0.0, 5.e-12);
	QCOMPARE(fitResult.rsquare, 1.0);
	QCOMPARE(fitResult.sse, 0.0);
	QCOMPARE(fitResult.rms, 0.0);