Commit 7cc4a659 authored by Thorbjørn Skovhus's avatar Thorbjørn Skovhus

use np.round when converting to integers for increased roboustness

parent 26a4aeb5
......@@ -462,7 +462,7 @@ def get_pw_coordinates(pd):
# Use cell to get coordinates
B_cv = 2.0 * np.pi * pd.gd.icell_cv
return np.dot(G_Gv, np.linalg.inv(B_cv)).astype(int)
return np.round(np.dot(G_Gv, np.linalg.inv(B_cv))).astype(int)
def write_macroscopic_component(omega_w, chiks_w, chi_w, filename, world):
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment