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
|
import io
from ase.io import read
header = """
___ ___ ___ _ _ _
| | |_ | | | |
| | | | | . | | | |
|__ | _|___|_____| 21.1.0
|___|_|
"""
atoms = """
Reference energy: -26313.685229
Positions:
0 Al 0.000000 0.000000 0.000000 ( 0.0000, 0.0000, 0.0000)
Unit cell:
periodic x y z points spacing
1. axis: yes 4.050000 0.000000 0.000000 21 0.1929
2. axis: yes 0.000000 4.050000 0.000000 21 0.1929
3. axis: yes 0.000000 0.000000 4.050000 21 0.1929
Energy contributions relative to reference atoms: (reference = -26313.685229)
Kinetic: +23.028630
Potential: -8.578488
External: +0.000000
XC: -24.279425
Entropy (-ST): -0.381921
Local: -0.018721
--------------------------
Free energy: -10.229926
Extrapolated: -10.038965
"""
forces = """
Forces in eV/Ang:
0 Al 0.00000 0.00000 -0.00000
"""
# Three configurations. Only 1. and 3. has forces.
text = header + atoms + forces + atoms + atoms + forces
def test_gpaw_output():
"""Regression test for #896.
"ase.io does not read all configurations from gpaw-out file"
"""
fd = io.StringIO(text)
configs = read(fd, index=':', format='gpaw-out')
assert len(configs) == 3
|
