neopdf_capi 0.3.0

C language interface to NeoPDF
Documentation 
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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163

#include <neopdf_capi.h>
#include <NeoPDF.hpp>
#include <cassert>
#include <cmath>
#include <iostream>
#include <string>

using namespace NEOLHAPDF;

const double TOLERANCE = 1e-12;

void test_lhapdf_compatibility_oop() {
    std::cout << "=== Test LHAPDF C++ OOP Compatibility Layer ===\n";

    std::string pdfname = "NNPDF40_nnlo_as_01180";

    // Use the compatibility layer
    setVerbosity(0);
    PDF* neo_lha_pdf = mkPDF(pdfname, 0);

    // Use the native NeoPDF C API for comparison
    NeoPDFWrapper* neo_pdf = neopdf_pdf_load(pdfname.c_str(), 0);

    double x = 1e-4;
    double q2 = 100.0;
    int pid = 21;

    double result = neo_lha_pdf->xfxQ2(pid, x, q2);
    double expected = neopdf_pdf_xfxq2(neo_pdf, pid, x, q2);
    double reldif = std::abs(result - expected) / expected;

    std::cout << "pid=" << pid << ", x=" << x << ", Q2=" << q2 << std::endl;
    std::cout << "Native NeoPDF: " << expected << std::endl;
    std::cout << "NeoPDF (LHAPDF compat): " << result << std::endl;
    std::cout << "Relative difference: " << reldif << std::endl;

    double as_result = neo_lha_pdf->alphasQ2(q2);
    double as_expected = neopdf_pdf_alphas_q2(neo_pdf, q2);
    double as_reldif = std::abs(as_result - as_expected) / as_expected;

    std::cout << "Q2=" << q2 << std::endl;
    std::cout << "Native NeoPDF: " << as_expected << std::endl;
    std::cout << "NeoPDF (LHAPDF compat): " << as_result << std::endl;
    std::cout << "Relative difference: " << as_reldif << std::endl;

    assert(reldif < TOLERANCE);
    assert(as_reldif < TOLERANCE);

    delete neo_lha_pdf;
    neopdf_pdf_free(neo_pdf);

    std::cout << "LHAPDF C++ OOP Compatibility test passed." << std::endl;
}

void test_lhapdf_compatibility_c() {
    printf("=== Test LHAPDF C Compatibility Layer ===\n");

    const char* pdfname = "NNPDF40_nnlo_as_01180";

    // Use the compatibility layer
    initpdfsetbyname(pdfname);
    initpdf(0);
    double x = 0.1;
    double q = 10.0;
    double xfxs[14];
    evolvepdf(x, q, xfxs);
    double a_s = alphaspdf(q);

    // Use the native NeoPDF API for comparison
    NeoPDFWrapper* neo_pdf = neopdf_pdf_load(pdfname, 0);
    double expected_g = neopdf_pdf_xfxq2(neo_pdf, -1, x, q*q);
    double as_expected = neopdf_pdf_alphas_q2(neo_pdf, q*q);

    double reldiff = fabs(xfxs[6] - expected_g) / expected_g;
    double as_reldiff = fabs(a_s - as_expected) / as_expected;

    printf("Comparing gluon PDF at x=0.1, Q=10.0\n");
    printf("Native NeoPDF: %e\n", expected_g);
    printf("NeoPDF (LHAPDF compat): %e\n", xfxs[6]);
    printf("Relative difference: %e\n", reldiff);

    printf("Comparing AlphaS at Q=10.0\n");
    printf("Native NeoPDF: %e\n", as_expected);
    printf("NeoPDF (LHAPDF compat): %e\n", a_s);
    printf("Relative difference: %e\n", as_reldiff);

    assert(reldiff < TOLERANCE);
    assert(as_reldiff < TOLERANCE);

    printf("LHAPDF C Compatibility test passed.\n");
    neopdf_pdf_free(neo_pdf);
}

void test_load_by_lhaid_oop() {
    std::cout << "=== Test Loading PDF by LHAID (NeoPDF OOP) ===\n";

    // NNPDF40_nnlo_as_01180 has base LHAID 331100; member 0 = LHAID 331100
    uint32_t lhaid = 331100;
    std::string pdfname = "NNPDF40_nnlo_as_01180";

    auto pdf_by_id = neopdf::NeoPDF::from_lhaid(lhaid);
    neopdf::NeoPDF pdf_by_name(pdfname, 0);

    double x   = 1e-4;
    double q2  = 100.0;
    int    pid = 21;

    double result   = pdf_by_id->xfxQ2(pid, x, q2);
    double expected = pdf_by_name.xfxQ2(pid, x, q2);
    double reldif   = std::abs(result - expected) / std::abs(expected);

    std::cout << "LHAID " << lhaid << " -> xfxQ2(pid=" << pid << ", x=" << x << ", Q2=" << q2 << ")\n";
    std::cout << "  from_name:  " << expected << "\n";
    std::cout << "  from_lhaid: " << result << "\n";
    std::cout << "  Relative diff: " << reldif << "\n";

    assert(reldif < TOLERANCE);

    double as_result   = pdf_by_id->alphasQ2(q2);
    double as_expected = pdf_by_name.alphasQ2(q2);
    double as_reldif   = std::abs(as_result - as_expected) / std::abs(as_expected);

    assert(as_reldif < TOLERANCE);

    std::cout << "NeoPDF OOP LHAID loading test passed." << std::endl;
}

void test_load_by_lhaid_lhapdf_compat() {
    std::cout << "=== Test Loading PDF by LHAID (NEOLHAPDF mkPDF) ===\n";

    // NNPDF40_nnlo_as_01180 has base LHAID 331100; member 0 = LHAID 331100
    PDF* pdf_by_id   = mkPDF((uint32_t)331100);
    PDF* pdf_by_name = mkPDF(std::string("NNPDF40_nnlo_as_01180"), 0);

    double x   = 1e-4;
    double q2  = 100.0;
    int    pid = 21;

    double result   = pdf_by_id->xfxQ2(pid, x, q2);
    double expected = pdf_by_name->xfxQ2(pid, x, q2);
    double reldif   = std::abs(result - expected) / std::abs(expected);

    std::cout << "mkPDF(331100) -> xfxQ2(pid=" << pid << ", x=" << x << ", Q2=" << q2 << ")\n";
    std::cout << "  mkPDF(name):  " << expected << "\n";
    std::cout << "  mkPDF(lhaid): " << result << "\n";
    std::cout << "  Relative diff: " << reldif << "\n";

    assert(reldif < TOLERANCE);

    delete pdf_by_id;
    delete pdf_by_name;

    std::cout << "NEOLHAPDF mkPDF(LHAID) test passed." << std::endl;
}

int main() {
    test_lhapdf_compatibility_oop();
    test_lhapdf_compatibility_c();
    test_load_by_lhaid_oop();
    test_load_by_lhaid_lhapdf_compat();

    return 0;
}