File indexing completed on 2024-05-05 03:48:44

 /*
     File                 : NSLDiffTest.cpp
     Project              : LabPlot
     Description          : NSL Tests for numerical differentiation
     --------------------------------------------------------------------
     SPDX-FileCopyrightText: 2019 Stefan Gerlach <stefan.gerlach@uni.kn>
 
     SPDX-License-Identifier: GPL-2.0-or-later
 */
 
 #include "NSLDiffTest.h"
 #include <QScopedArrayPointer>
 
 extern "C" {
 #include "backend/nsl/nsl_diff.h"
 }
 
 // ##############################################################################
 // #################  first derivative tests
 // ##############################################################################
 
 const int N = 7;
 double dxdata[] = {1, 2, 4, 8, 16, 32, 64};
 
 void NSLDiffTest::testFirst_order2() {
     double ydata[] = {1, 4, 16, 64, 256, 1024, 4096};
 
     int status = nsl_diff_first_deriv(dxdata, ydata, N, 2);
     QCOMPARE(status, 0);
     for (unsigned int i = 0; i < N; i++)
         QCOMPARE(ydata[i], 2 * dxdata[i]);
 }
 
 void NSLDiffTest::testFirst_order4() {
     double ydata[] = {1, 8, 64, 512, 4096, 32768, 262144};
 
     int status = nsl_diff_first_deriv(dxdata, ydata, N, 4);
     QCOMPARE(status, 0);
     for (unsigned int i = 0; i < N; i++)
         QCOMPARE(ydata[i], 3 * dxdata[i] * dxdata[i]);
 }
 
 void NSLDiffTest::testFirst_avg() {
     double ydata[] = {1, 4, 16, 64, 256, 1024, 4096};
     const double result[] = {3, 4.5, 9, 18, 36, 72, 96};
 
     int status = nsl_diff_first_deriv_avg(dxdata, ydata, N);
     QCOMPARE(status, 0);
     for (unsigned int i = 0; i < N; i++)
         QCOMPARE(ydata[i], result[i]);
 }
 
 // ##############################################################################
 // #################  second derivative tests
 // ##############################################################################
 
 void NSLDiffTest::testSecond_order1() {
     double ydata[] = {1, 4, 16, 64, 256, 1024, 4096};
 
     int status = nsl_diff_second_deriv(dxdata, ydata, N, 1);
     QCOMPARE(status, 0);
     for (double d : ydata)
         QCOMPARE(d, 2.);
 }
 
 void NSLDiffTest::testSecond_order2() {
     double ydata[] = {1, 4, 16, 64, 256, 1024, 4096};
 
     int status = nsl_diff_second_deriv(dxdata, ydata, N, 2);
     QCOMPARE(status, 0);
     for (double d : ydata)
         QCOMPARE(d, 2.);
 }
 
 void NSLDiffTest::testSecond_order3() {
     double ydata[] = {1, 8, 64, 512, 4096, 32768, 262144};
 
     int status = nsl_diff_second_deriv(dxdata, ydata, N, 3);
     QCOMPARE(status, 0);
     for (unsigned int i = 0; i < N; i++)
         QCOMPARE(ydata[i], 6. * dxdata[i]);
 }
 
 // ##############################################################################
 // #################  higher derivative tests
 // ##############################################################################
 
 void NSLDiffTest::testThird_order2() {
     double ydata[] = {1, 8, 64, 512, 4096, 32768, 262144};
 
     int status = nsl_diff_third_deriv(dxdata, ydata, N, 2);
     QCOMPARE(status, 0);
     for (double d : ydata)
         QCOMPARE(d, 6.);
 }
 
 void NSLDiffTest::testFourth_order1() {
     double ydata[] = {1, 8, 64, 512, 4096, 32768, 262144};
 
     int status = nsl_diff_fourth_deriv(dxdata, ydata, N, 1);
     QCOMPARE(status, 0);
     for (double d : ydata)
         QCOMPARE(d, 0.);
 }
 
 void NSLDiffTest::testFourth_order3() {
     double ydata[] = {1, 8, 64, 512, 4096, 32768, 262144};
 
     int status = nsl_diff_fourth_deriv(dxdata, ydata, N, 3);
     QCOMPARE(status, 0);
     for (double d : ydata)
         QCOMPARE(d + 1., 1.);
 }
 
 void NSLDiffTest::testFifth_order2() {
     double ydata[] = {1, 8, 64, 512, 4096, 32768, 262144};
 
     int status = nsl_diff_fifth_deriv(dxdata, ydata, N, 2);
     QCOMPARE(status, 0);
     for (double d : ydata)
         QCOMPARE(d + 1., 1.);
 }
 
 void NSLDiffTest::testSixth_order1() {
     double ydata[] = {1, 8, 64, 512, 4096, 32768, 262144};
 
     int status = nsl_diff_sixth_deriv(dxdata, ydata, N, 1);
     QCOMPARE(status, 0);
     for (double d : ydata)
         QCOMPARE(d + 1., 1.);
 }
 
 // ##############################################################################
 // #################  performance
 // ##############################################################################
 
 void NSLDiffTest::testPerformance_first() {
     const int NN = 1e6;
     QScopedArrayPointer<double> xdata(new double[NN]);
     QScopedArrayPointer<double> ydata(new double[NN]);
 
     for (int i = 0; i < NN; i++)
         xdata[i] = ydata[i] = (double)i;
 
     QBENCHMARK {
         int status = nsl_diff_first_deriv(xdata.data(), ydata.data(), NN, 2);
         QCOMPARE(status, 0);
     }
 }
 
 void NSLDiffTest::testPerformance_second() {
     const int NN = 1e6;
     QScopedArrayPointer<double> xdata(new double[NN]);
     QScopedArrayPointer<double> ydata(new double[NN]);
 
     for (int i = 0; i < NN; i++)
         xdata[i] = ydata[i] = (double)i;
 
     QBENCHMARK {
         int status = nsl_diff_second_deriv(xdata.data(), ydata.data(), NN, 2);
         QCOMPARE(status, 0);
     }
 }
 
 void NSLDiffTest::testPerformance_third() {
     const int NN = 1e6;
     QScopedArrayPointer<double> xdata(new double[NN]);
     QScopedArrayPointer<double> ydata(new double[NN]);
 
     for (int i = 0; i < NN; i++)
         xdata[i] = ydata[i] = (double)i;
 
     QBENCHMARK {
         int status = nsl_diff_third_deriv(xdata.data(), ydata.data(), NN, 2);
         QCOMPARE(status, 0);
     }
 }
 
 QTEST_MAIN(NSLDiffTest)