forked from aiidateam/aiida-wannier90-workflows
-
Notifications
You must be signed in to change notification settings - Fork 0
/
example_02.py
executable file
·67 lines (53 loc) · 1.98 KB
/
example_02.py
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
#!/usr/bin/env python
"""Run a ``Wannier90BandsWorkChain`` for Wannier90 band structure.
Usage: ./example_02.py
"""
import click
from aiida import cmdline, orm
from aiida_wannier90_workflows.cli.params import RUN
from aiida_wannier90_workflows.utils.code import check_codes, identify_codes
from aiida_wannier90_workflows.utils.structure import read_structure
from aiida_wannier90_workflows.utils.workflows.builder.serializer import print_builder
from aiida_wannier90_workflows.utils.workflows.builder.setter import set_parallelization
from aiida_wannier90_workflows.utils.workflows.builder.submit import (
submit_and_add_group,
)
from aiida_wannier90_workflows.workflows import Wannier90BandsWorkChain
def submit(
codes: dict,
structure: orm.StructureData,
group: orm.Group = None,
run: bool = False,
):
"""Submit a ``Wannier90BandsWorkChain`` to calculate Wannier bands."""
codes = identify_codes(codes)
check_codes(codes)
builder = Wannier90BandsWorkChain.get_builder_from_protocol(
codes,
structure,
protocol="fast",
run_open_grid=True,
)
# You can change parallelization here
parallelization = {
"num_mpiprocs_per_machine": 1,
"npool": 1,
}
set_parallelization(builder, parallelization, process_class=Wannier90BandsWorkChain)
print_builder(builder)
if run:
submit_and_add_group(builder, group)
@click.command()
@cmdline.utils.decorators.with_dbenv()
@cmdline.params.options.CODES()
@cmdline.params.options.GROUP(help="The group to add the submitted workchain.")
@click.argument("filename", type=click.Path(exists=True))
@RUN()
def cli(filename, codes, group, run):
"""Run a ``Wannier90BandsWorkChain`` to calculate Wannier90 band structure.
FILENAME: a crystal structure file, e.g., ``input_files/GaAs.xsf``.
"""
struct = read_structure(filename, store=True)
submit(codes, struct, group, run)
if __name__ == "__main__":
cli() # pylint: disable=no-value-for-parameter