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
|
// Copyright 2023 Global Phasing Ltd.
#include <cstdio> // for printf, fprintf
#include <algorithm> // for find
#include "gemmi/read_cif.hpp" // for read_cif_gz
#define GEMMI_PROG cifdiff
#include "options.h"
namespace cif = gemmi::cif;
namespace {
enum OptionIndex {
OnlyCategories=4, NoComparison,
};
const option::Descriptor Usage[] = {
{ NoOp, 0, "", "", Arg::None,
"Usage:\n"
" " EXE_NAME " [options] FILE1.cif FILE2.cif\n"
" " EXE_NAME " [options] -n FILE.cif\n\n"
"Compares (or just prints) categories and tags in CIF files.\n"
"First block only." },
CommonUsage[Help],
CommonUsage[Version],
CommonUsage[Verbose],
{ OnlyCategories, 0, "q", "", Arg::None,
" -q \tPrint only categories." },
{ NoComparison, 0, "n", "", Arg::None,
" -n \tNo comparison, just list categories and tags." },
{ 0, 0, 0, 0, 0, 0 }
};
struct DiffItem {
char change;
std::string str;
};
using Diff = std::vector<DiffItem>;
template<typename T>
Diff make_diff(const T& a, const T& b) {
Diff diff;
diff.reserve(a.size() + b.size());
for (const std::string& s : b)
diff.push_back(DiffItem{'+', s});
size_t idx = 0;
for (const std::string& s : a) {
auto pos = std::find_if(diff.begin(), diff.end(),
[&](const DiffItem& d) { return d.str == s; });
if (pos != diff.end()) {
pos->change = ' ';
idx = pos - diff.begin() + 1;
} else {
diff.insert(diff.begin() + idx, DiffItem{'-', s});
idx++;
}
}
return diff;
}
} // anonymous namespace
int GEMMI_MAIN(int argc, char **argv) {
OptParser p(EXE_NAME);
p.simple_parse(argc, argv, Usage);
bool one_file = p.options[NoComparison];
p.require_positional_args(one_file ? 1 : 2);
const char* path1 = p.nonOption(0);
const char* path2 = one_file ? nullptr : p.nonOption(1);
//int verbose = p.options[Verbose].count();
try {
// Starting like an unified diff (with "--- ") enables colordiff.
printf("%sReading %s\n", one_file ? "" : "--- ", path1);
cif::Document doc1 = gemmi::read_cif_or_mmjson_gz(path1);
cif::Block* b1 = &doc1.blocks.at(0);
// NoComparison mode is implemented as comparing Block with itself
// (inefficient, but simple).
cif::Document doc2;
cif::Block* b2 = b1;
if (!one_file) {
printf("+++ Reading %s\n", path2);
doc2 = gemmi::read_cif_or_mmjson_gz(path2);
b2 = &doc2.blocks.at(0);
}
if (b1->name == b2->name) {
printf(" block name: %s\n", b1->name.c_str());
} else {
printf("- block name: %s\n", b1->name.c_str());
printf("+ block name: %s\n", b2->name.c_str());
}
Diff category_diff = make_diff(b1->get_mmcif_category_names(),
b2->get_mmcif_category_names());
for (DiffItem& cat : category_diff) {
cif::Table t1 = b1->find_mmcif_category(cat.str);
cif::Table t2 = b2->find_mmcif_category(cat.str);
size_t len1 = t1.length();
size_t len2 = t2.length();
printf("%c %-37s rows: %5zu", cat.change, cat.str.c_str(), len1);
if (len2 != len1)
printf(" -> %5zu", len2);
putchar('\n');
if (p.options[OnlyCategories])
continue;
size_t prefix_size = cat.str.size();
Diff tag_diff = make_diff(t1.tags(), t2.tags());
for (DiffItem& di : tag_diff)
printf("%c %s\n", di.change, di.str.c_str() + prefix_size);
}
} catch (std::exception& e) {
std::fprintf(stderr, "ERROR: %s\n", e.what());
return 1;
}
return 0;
}
|
