Skip to content

tmol

parse_tmol_file(filename)

Convert a provided turbomole coord file into an array of atomic numbers and cartesian positions

Parameters:

Name Type Description Default
filename

path to the turbomole file to read

 required

Returns:

Type Description

array of (N) atomic numbers and (N, 3) Cartesian positions read from

the given file
Source code in chmpy/fmt/tmol.py 
49
50
51
52
53
54
55
56
57
58
59
60
61
def parse_tmol_file(filename):
    """Convert a provided turbomole coord file into an array of
    atomic numbers and cartesian positions

    Args:
        filename: path to the turbomole file to read

    Returns:
        array of (N) atomic numbers and (N, 3) Cartesian positions read from
        the given file
    """
    path = Path(filename)
    return parse_tmol_string(path.read_text(), filename=str(path.absolute()))

parse_tmol_string(contents, filename=None)

Convert provided turbomole coord file contents into an array of atomic numbers and cartesian positions

Parameters:

Name Type Description Default
contents

text contents of the turbomole coord file to read

 required

Returns:

Type Description

array of (N) atomic numbers and (N, 3) Cartesian positions read from

the given file
Source code in chmpy/fmt/tmol.py 
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
def parse_tmol_string(contents, filename=None):
    """Convert provided turbomole coord file contents into an array of
    atomic numbers and cartesian positions

    Args:
        contents: text contents of the turbomole coord file to read

    Returns:
        array of (N) atomic numbers and (N, 3) Cartesian positions read from
        the given file
    """
    lines = contents.splitlines()
    angstroms = "angs" in lines[0]
    elements = []
    positions = []
    for line in lines[1:]:
        stripped = line.strip()
        if (not stripped) or "$end" in line:
            break
        if "$" in line:
            continue
        tokens = line.strip().split()
        xyz = tuple(float(x) for x in tokens[:3])
        positions.append(xyz)
        elements.append(Element[tokens[3]])
    LOG.debug(
        "Found %d atoms lines in %s",
        len(elements),
        "in " + filename if filename else "",
    )
    positions = np.array(positions)
    if not angstroms:
        from chmpy.util.unit import units

        positions = units.angstrom(positions)
    return elements, positions