test_parsing.py 12.1 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
# Copyright 2018 Markus Scheidgen
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
#   http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an"AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.

15
16
from io import StringIO
import json
17
import numpy as np
18
import pytest
19
import os
20
21

from nomadcore.local_meta_info import loadJsonFile
22
import nomad_meta_info
23

24
from nomad import utils, files
25
from nomad.parsing import JSONStreamWriter, parser_dict, match_parser, BrokenParser
26
from nomad.parsing import LocalBackend, BadContextURI
27

28
parser_examples = [
29
    ('parsers/random', 'test/data/parsers/random_0'),
30
    ('parsers/template', 'tests/data/parsers/template.json'),
31
32
    ('parsers/exciting', 'tests/data/parsers/exciting/Ag/INFO.OUT'),
    ('parsers/exciting', 'tests/data/parsers/exciting/GW/INFO.OUT'),
Markus Scheidgen's avatar
Markus Scheidgen committed
33
    ('parsers/vasp', 'tests/data/parsers/vasp/vasp.xml'),
34
    ('parsers/vasp', 'tests/data/parsers/vasp_compressed/vasp.xml.gz'),
35
    ('parsers/vaspoutcar', 'tests/data/parsers/vasp_outcar/OUTCAR'),
speckhard's avatar
speckhard committed
36
    ('parsers/fhi-aims', 'tests/data/parsers/fhi-aims/aims.out'),
37
    ('parsers/cp2k', 'tests/data/parsers/cp2k/si_bulk8.out'),
38
39
    ('parsers/crystal', 'tests/data/parsers/crystal/si.out'),
    ('parsers/cpmd', 'tests/data/parsers/cpmd/geo_output.out'),
speckhard's avatar
speckhard committed
40
    ('parsers/nwchem', 'tests/data/parsers/nwchem/single_point/output.out'),
41
    ('parsers/bigdft', 'tests/data/parsers/bigdft/n2_output.out'),
42
    ('parsers/wien2k', 'tests/data/parsers/wien2k/AlN/AlN_ZB.scf'),
43
    # ('parsers/band', 'tests/data/parsers/band_adf.out'),  # causes spglib to segfault
44
    ('parsers/gaussian', 'tests/data/parsers/gaussian/aniline.out'),
Daniel Speckhard's avatar
Daniel Speckhard committed
45
    ('parsers/abinit', 'tests/data/parsers/abinit/Fe.out'),
Daniel Speckhard's avatar
Daniel Speckhard committed
46
    ('parsers/quantumespresso', 'tests/data/parsers/quantum-espresso/benchmark.out'),
47
    ('parsers/orca', 'tests/data/parsers/orca/orca3dot2706823.out'),
Daniel Speckhard's avatar
Daniel Speckhard committed
48
    ('parsers/castep', 'tests/data/parsers/castep/BC2N-Pmm2-Raman.castep'),
Daniel Speckhard's avatar
Daniel Speckhard committed
49
    # ('parsers/dl-poly', 'tests/data/parsers/dl-poly/OUTPUT'), # timeout on Matid System Classification
Daniel Speckhard's avatar
Daniel Speckhard committed
50
    ('parsers/lib-atoms', 'tests/data/parsers/lib-atoms/gp.xml'),
Daniel Speckhard's avatar
Daniel Speckhard committed
51
    ('parsers/octopus', 'tests/data/parsers/octopus/stdout.txt'),
Markus Scheidgen's avatar
Markus Scheidgen committed
52
    ('parsers/phonopy', 'tests/data/parsers/phonopy/control.in'),
Daniel Speckhard's avatar
Daniel Speckhard committed
53
54
    ('parsers/gpaw', 'tests/data/parsers/gpaw/Fe2.gpw'),
    ('parsers/gpaw2', 'tests/data/parsers/gpaw2/H2_lcao.gpw2'),
Daniel Speckhard's avatar
Daniel Speckhard committed
55
    ('parsers/atk', 'tests/data/parsers/atk/Si2.nc'),
56
57
    ('parsers/gulp', 'tests/data/parsers/gulp/example6.got'),
    ('parsers/siesta', 'tests/data/parsers/siesta/Fe/out'),
Markus Scheidgen's avatar
Markus Scheidgen committed
58
    ('parsers/elk', 'tests/data/parsers/elk/Al/INFO.OUT'),
59
60
    ('parsers/elastic', 'dependencies/parsers/elastic/test/examples/2nd/INFO_ElaStic'),  # 70Mb file 2big4git
    ('parsers/turbomole', 'tests/data/parsers/turbomole/acrolein.out'),
Daniel Speckhard's avatar
Daniel Speckhard committed
61
    ('parsers/gamess', 'tests/data/parsers/gamess/exam01.out')
Daniel Speckhard's avatar
Daniel Speckhard committed
62
]
63

64
65
# We need to remove some cases with external mainfiles, which might not exist
# in all testing environments (e.g. in the nomad docker image)
Markus Scheidgen's avatar
Markus Scheidgen committed
66
fixed_parser_examples = []
67
68
69
for parser, mainfile in parser_examples:
    if os.path.exists(mainfile) or mainfile.startswith('tests'):
        fixed_parser_examples.append((parser, mainfile))
Markus Scheidgen's avatar
Markus Scheidgen committed
70
parser_examples = fixed_parser_examples
71

72

73
correct_num_output_files = 39
74

75
76
77
78
79

class TestLocalBackend(object):

    @pytest.fixture(scope='session')
    def meta_info(self):
80
81
        file_dir = os.path.dirname(os.path.abspath(nomad_meta_info.__file__))
        path = os.path.join(file_dir, 'all.nomadmetainfo.json')
82
83
84
85
86
87
88
        meta_info, _ = loadJsonFile(path)
        return meta_info

    @pytest.fixture(scope='function')
    def backend(self, meta_info):
        return LocalBackend(meta_info, debug=True)

89
    def test_meta_info(self, meta_info, no_warn):
90
91
        assert 'section_topology' in meta_info

92
    def test_section(self, backend, no_warn):
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
        g_index = backend.openSection('section_run')
        assert g_index == 0
        backend.addValue('program_name', 't0')
        backend.closeSection('section_run', 0)

        g_index = backend.openSection('section_run')
        assert g_index == 1

        g_index = backend.openSection('section_run')
        assert g_index == 2

        backend.addValue('program_name', 't1', 1)
        backend.addValue('program_name', 't2', 2)

        backend.closeSection('section_run', 1)
        backend.closeSection('section_run', 2)

        assert backend.get_sections('section_run') == [0, 1, 2]
        for i in range(0, 3):
            assert backend.get_value('program_name', i) == 't%d' % i

114
115
116
117
118
119
120
121
122
123
124
125
126
127
    def test_section_override(self, backend, no_warn):
        """ Test whether we can overwrite values already in the backend."""
        expected_value = ['Cl', 'Zn']
        backend.openSection('section_run')
        backend.openSection('section_system')
        backend.addArrayValues('atom_labels', np.array(['Al', 'Zn']))
        backend.addArrayValues('atom_labels', np.array(expected_value), override=True)
        backend.closeSection('section_system', 0)

        backend.closeSection('section_run', 0)
        output = StringIO()
        backend.write_json(output)
        assert backend.get_value('atom_labels').tolist() == expected_value

128
129
130
131
132
133
    def test_two_sections(self, backend, no_warn):
        g_index = backend.openSection('section_run')
        assert g_index == 0
        backend.addValue('program_name', 't0')
        backend.closeSection('section_run', 0)

134
        g_index = backend.openSection('section_entry_info')
135
136
        assert g_index == 0
        backend.addValue('parser_name', 'p0')
137
        backend.closeSection('section_entry_info', 0)
138
139

        assert backend.get_sections('section_run') == [0]
140
        assert backend.get_sections('section_entry_info') == [0]
141
142
143
144
145

        output = StringIO()
        backend.write_json(output)
        archive = json.loads(output.getvalue())
        assert 'section_run' in archive
146
        assert 'section_entry_info' in archive
147

148
    def test_subsection(self, backend: LocalBackend, no_warn):
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
        backend.openSection('section_run')
        backend.openSection('section_method')
        backend.closeSection('section_method', -1)

        backend.openSection('section_method')
        backend.closeSection('section_method', -1)

        backend.openSection('section_run')
        backend.closeSection('section_run', 0)
        backend.closeSection('section_run', 1)

        backend.openSection('section_method')
        backend.closeSection('section_method', -1)

        runs = backend.data['section_run']
        assert len(runs) == 2
        assert len(runs[0]['section_method']) == 2
        assert len(runs[1]['section_method']) == 1

168
    def test_context(self, backend: LocalBackend, no_warn):
169
170
171
172
173
        backend.openSection('section_run')
        backend.openSection('section_method')
        backend.closeSection('section_method', -1)
        backend.closeSection('section_run', -1)

174
175
        backend.openSection('section_run')
        backend.closeSection('section_run', -1)
176

177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
        backend.openContext('/section_run/0')
        backend.addValue('program_name', 't1')
        backend.closeContext('/section_run/0')

        backend.openContext('/section_run/1')
        backend.addValue('program_name', 't2')
        backend.closeContext('/section_run/1')

        backend.openContext('/section_run/0/section_method/0')
        backend.closeContext('/section_run/0/section_method/0')

        runs = backend.data['section_run']
        assert runs[0]['program_name'] == 't1'
        assert runs[1]['program_name'] == 't2'

192
    def test_multi_context(self, backend: LocalBackend, no_warn):
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
        backend.openSection('section_run')
        backend.closeSection('section_run', -1)

        backend.openContext('/section_run/0')
        backend.openSection('section_method')
        backend.closeSection('section_method', -1)
        backend.closeContext('/section_run/0')

        backend.openContext('/section_run/0')
        backend.openSection('section_method')
        backend.closeSection('section_method', -1)
        backend.closeContext('/section_run/0')

        assert len(backend.data['section_method']) == 1

208
    def test_bad_context(self, backend: LocalBackend, no_warn):
209
210
211
212
213
214
215
216
217
218
219
220
        try:
            backend.openContext('section_run/0')
            assert False
        except BadContextURI:
            pass

        try:
            backend.openContext('dsfds')
            assert False
        except BadContextURI:
            pass

221
222
223
224
225
226
227
228
229

def create_reference(data, pretty):
    if (pretty):
        return json.dumps(data, indent=2)
    else:
        return json.dumps(data, separators=(',', ':'))


@pytest.mark.parametrize("pretty", [False, True])
230
def test_stream_generator(pretty, no_warn):
231
232
233
234
235
236
237
238
239
240
241
242
243
    example_data = [
        {
            'key1': 'value',
            'key2': 1
        },
        {
            'key': {
                'key': 'value'
            }
        }
    ]

    out = StringIO()
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
    writer = JSONStreamWriter(out, pretty=pretty)
    writer.open_array()
    writer.open_object()
    writer.key('key1')
    writer.value('value')
    writer.key('key2')
    writer.value(1)
    writer.close_object()
    writer.open_object()
    writer.key('key')
    writer.open_object()
    writer.key('key')
    writer.value('value')
    writer.close_object()
    writer.close_object()
    writer.close_array()
    writer.close()
261
262
263

    assert create_reference(example_data, pretty) == out.getvalue()

264

265
266
def assert_parser_result(backend):
    status, errors = backend.status
267
268
    assert status == 'ParseSuccess'
    assert errors is None or len(errors) == 0
269
270


Daniel Speckhard's avatar
Daniel Speckhard committed
271
272
273
274
275
def assert_parser_dir_unchanged(previous_wd, current_wd):
    """Assert working directory has not been changed from parser."""
    assert previous_wd == current_wd


276
277
def run_parser(parser_name, mainfile):
    parser = parser_dict[parser_name]
278
    result = parser.run(mainfile, logger=utils.get_logger(__name__))
279
    return add_calculation_info(result, parser_name=parser_name)
Markus Scheidgen's avatar
Markus Scheidgen committed
280
281


282
283
284
@pytest.fixture
def parsed_vasp_example() -> LocalBackend:
    return run_parser(
285
        'parsers/vasp', 'dependencies/parsers/vasp/test/examples/xml/perovskite.xml')
Markus Scheidgen's avatar
Markus Scheidgen committed
286
287


288
289
290
291
292
293
@pytest.fixture
def parsed_template_example() -> LocalBackend:
    return run_parser(
        'parsers/template', 'tests/data/parsers/template.json')


294
295
296
# Function used by normalizer tests.
def parse_file(parser_name_and_mainfile) -> LocalBackend:
    parser_name, mainfile = parser_name_and_mainfile
297
298
299
    return run_parser(parser_name, mainfile)


300
301
302
@pytest.fixture(params=parser_examples, ids=lambda spec: '%s-%s' % spec)
def parsed_example(request) -> LocalBackend:
    parser_name, mainfile = request.param
Markus Scheidgen's avatar
Markus Scheidgen committed
303
304
    result = run_parser(parser_name, mainfile)
    return result
305
306


307
def add_calculation_info(backend: LocalBackend, **kwargs) -> LocalBackend:
308
    backend.openNonOverlappingSection('section_entry_info')
309
    backend.addValue('upload_id', 'test_upload_id')
310
311
    backend.addValue('calc_id', 'test_calc_id')
    backend.addValue('calc_hash', 'test_calc_hash')
312
    backend.addValue('mainfile', 'test/mainfile.txt')
313
314
    for key, value in kwargs.items():
        backend.addValue(key, value)
315
    backend.closeNonOverlappingSection('section_entry_info')
316
    return backend
317
318


319
@pytest.mark.parametrize('parser_name, mainfile', parser_examples)
320
def test_parser(parser_name, mainfile):
Daniel Speckhard's avatar
Daniel Speckhard committed
321
    previous_wd = os.getcwd()  # Get Working directory before parsing.
322
    parsed_example = run_parser(parser_name, mainfile)
323
    assert_parser_result(parsed_example)
Daniel Speckhard's avatar
Daniel Speckhard committed
324
325
    # Check that cwd has not changed.
    assert_parser_dir_unchanged(previous_wd, current_wd=os.getcwd())
Markus Scheidgen's avatar
Markus Scheidgen committed
326
327


328
329
330
331
332
def test_match(raw_files, no_warn):
    example_upload_id = 'example_upload_id'
    upload_files = files.StagingUploadFiles(example_upload_id, create=True, is_authorized=lambda: True)
    upload_files.add_rawfiles('tests/data/parsers')

333
    count = 0
334
335
    for mainfile in upload_files.raw_file_manifest():
        parser = match_parser(mainfile, upload_files)
336
        if parser is not None and not isinstance(parser, BrokenParser):
337
            count += 1
338

339
    assert count == correct_num_output_files