Skip to content

Commit 7146c3d

Browse files
Added atomic dipole to test.py (#2747)
Added atomic dipole to test.py The current code did not cover this case. --------- Signed-off-by: Hananeh Oliaei <[email protected]> Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>
1 parent 5f8c976 commit 7146c3d

File tree

1 file changed

+30
-8
lines changed

1 file changed

+30
-8
lines changed

deepmd/entrypoints/test.py

Lines changed: 30 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -934,18 +934,40 @@ def test_dipole(
934934

935935
if detail_file is not None:
936936
detail_path = Path(detail_file)
937+
if not atomic:
938+
pe = np.concatenate(
939+
(
940+
np.reshape(test_data["dipole"][:numb_test], [-1, 3]),
941+
np.reshape(dipole, [-1, 3]),
942+
),
943+
axis=1,
944+
)
945+
header_text = "data_x data_y data_z pred_x pred_y pred_z"
946+
else:
947+
pe = np.concatenate(
948+
(
949+
np.reshape(
950+
test_data["atomic_dipole"][:numb_test], [-1, 3 * sel_natoms]
951+
),
952+
np.reshape(dipole, [-1, 3 * sel_natoms]),
953+
),
954+
axis=1,
955+
)
956+
header_text = [
957+
f"{letter}{number}"
958+
for number in range(1, sel_natoms + 1)
959+
for letter in ["data_x", "data_y", "data_z"]
960+
] + [
961+
f"{letter}{number}"
962+
for number in range(1, sel_natoms + 1)
963+
for letter in ["pred_x", "pred_y", "pred_z"]
964+
]
965+
header_text = " ".join(header_text)
937966

938-
pe = np.concatenate(
939-
(
940-
np.reshape(test_data["dipole"][:numb_test], [-1, 3]),
941-
np.reshape(dipole, [-1, 3]),
942-
),
943-
axis=1,
944-
)
945967
np.savetxt(
946968
detail_path.with_suffix(".out"),
947969
pe,
948-
header="data_x data_y data_z pred_x pred_y pred_z",
970+
header=header_text,
949971
)
950972
return {"rmse": (rmse_f, dipole.size)}
951973

0 commit comments

Comments
 (0)