Skip to content

Commit

Permalink
Remove parsing of coordinates and atomic numbers
Browse files Browse the repository at this point in the history
  • Loading branch information
yakutovicha committed Jul 9, 2020
1 parent 3b806eb commit f89ecb6
Showing 1 changed file with 2 additions and 10 deletions.
12 changes: 2 additions & 10 deletions aiida_quantumespresso/parsers/pp.py
Original file line number Diff line number Diff line change
Expand Up @@ -213,9 +213,9 @@ def detect_important_message(logs, line):

# Parse useful data from stdout
for line in stdout_lines:
if 'Check:' in line:
if 'Check:' in line: # QE < 6.5
split_line = line.split('=')
if 'negative/imaginary' in line: # QE6.1
if 'negative/imaginary' in line: # QE6.1-6.3
output_dict['negative_core_charge'] = float(split_line[-1].split()[0])
output_dict['imaginary_core_charge'] = float(split_line[-1].split()[-1])
else: # QE6.4
Expand Down Expand Up @@ -354,12 +354,6 @@ def parse_gaussian(self, data_file_str):
voxel_array = np.array([[x_line[1], x_line[2], x_line[3]], [y_line[1], y_line[2], y_line[3]],
[z_line[1], z_line[2], z_line[3]]],
dtype=np.float64)
atomic_numbers = np.empty(natoms, int)
coordinates = np.empty((natoms, 3))
for i in range(natoms):
line = header[6 + i].split()
atomic_numbers[i] = int(line[0])
coordinates[i] = [float(s) for s in line[2:]]

# Get the volumetric data
data_array = np.empty(xdim * ydim * zdim, dtype=float)
Expand All @@ -378,7 +372,5 @@ def parse_gaussian(self, data_file_str):
arraydata.set_array('data', data_array)
arraydata.set_array('data_units', np.array(data_units))
arraydata.set_array('coordinates_units', np.array(coordinates_units))
arraydata.set_array('coordinates', coordinates)
arraydata.set_array('atomic_numbers', atomic_numbers)

return arraydata

0 comments on commit f89ecb6

Please sign in to comment.