File: CompareTwoReferences.cpp

package info (click to toggle)
bornagain 23.0-4
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 103,936 kB
  • sloc: cpp: 423,131; python: 40,997; javascript: 11,167; awk: 630; sh: 318; ruby: 173; xml: 130; makefile: 51; ansic: 24
file content (49 lines) | stat: -rw-r--r-- 1,759 bytes parent folder | download | duplicates (2)
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
//  ************************************************************************************************
//
//  BornAgain: simulate and fit reflection and scattering
//
//! @file      Tests/Functional/Consistence/CompareTwoReferences.cpp
//! @brief     Implements function compareTwoReferences for use in core consistence test.
//!
//! @homepage  http://www.bornagainproject.org
//! @license   GNU General Public License v3 or higher (see COPYING)
//! @copyright Forschungszentrum Jülich GmbH 2018
//! @authors   Scientific Computing Group at MLZ (see CITATION, AUTHORS)
//
//  ************************************************************************************************

#include "BATesting.h"
#include "Base/Util/Assert.h"
#include "Base/Util/PathUtil.h"
#include "Device/Data/Datafield.h"
#include "Device/IO/DiffUtil.h"
#include "Device/IO/IOFactory.h"
#include <cstdlib>
#include <iostream>
#include <memory>

namespace {

Datafield load(const std::string& name)
{
    ASSERT(!name.empty());
    const std::string path = Base::Path::jointPath(BATesting::ReferenceDir_Suite, name + ".int");
    try {
        return IO::readData2D(path, IO::Filetype2D::bornagain2D);
    } catch (const std::exception&) {
        std::cout << "Data file " << path << " could not be read.\n"
                  << "Run the pertinent Core standard test, copy the fresh data file"
                  << " to the reference directory,\n"
                  << "then rerun this test." << std::endl;
        exit(-1);
    }
}

} // namespace

int compareTwoReferences(const std::string& name0, const std::string& name1, const double limit)
{
    Datafield data0 = load(name0);
    Datafield data1 = load(name1);
    return DiffUtil::checkRelativeDifference(data0, data1, limit);
}