LCPQ/dft_tools
3
0
Fork 0
mirror of https://github.com/triqs/dft_tools synced 2024-12-23 04:43:42 +01:00
Code Issues Releases Wiki Activity

fix for bug in wannier90 converter

Browse Source 
* the Fermi level given in the inp file was ignored
* added a new line reporting which energy is subtracted from mu
This commit is contained in:
Alexander Hampel 2020-08-25 08:22:45 -04:00
parent 6a104622d1
commit d61ce6eaae
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
python/triqs_dft_tools/converters/wannier90.py
View File

@ -150,7 +150,7 @@ class Wannier90Converter(ConverterTools):
corr_shells = [{name: int(val) for name, val in zip( corr_shells = [{name: int(val) for name, val in zip(
corr_shell_entries, R)} for icrsh in range(n_corr_shells)] corr_shell_entries, R)} for icrsh in range(n_corr_shells)]
try: try:
self.fermi_energy = R.next() self.fermi_energy = float(next(R))
except: except:
self.fermi_energy = 0. self.fermi_energy = 0.
except StopIteration: # a more explicit error if the file is corrupted. except StopIteration: # a more explicit error if the file is corrupted.
@ -190,6 +190,7 @@ class Wannier90Converter(ConverterTools):
shells_map = [inequiv_to_corr[corr_to_inequiv[ish]] shells_map = [inequiv_to_corr[corr_to_inequiv[ish]]
for ish in range(n_corr_shells)] for ish in range(n_corr_shells)]
mpi.report("Mapping: " + format(shells_map)) mpi.report("Mapping: " + format(shells_map))
mpi.report("Subtracting %f eV from the Fermi level." % self.fermi_energy)
# build the k-point mesh, if its size was given on input (kmesh_mode >= 0), # build the k-point mesh, if its size was given on input (kmesh_mode >= 0),
# otherwise it is built according to the data in the hr file (see # otherwise it is built according to the data in the hr file (see