24: 25,
25: 24
}
surface_data = ([{
'l': 1.0,
'd': 0.5
}] * 6 + [{
'l': 1.5,
'd': 1.0 / 4.0
}] * 12 + [{
'l': 1.0,
'd': 1.0 / 3.0
}] * 8)
surface_symmetries = get_surface_symmetries(symmetries, surface_names,
surface_numbers)
def surface_fitting(surfaces):
for i, n1 in enumerate(np.array(surface_names)):
d1 = surface_data[i]['d']
a1 = surfaces[i] * d1
for j, n2 in enumerate(np.array(surface_names)):
d2 = surface_data[j]['d']
a2 = surfaces[j] * d2
nd = np.dot(n1, n2) / (np.linalg.norm(n1) * np.linalg.norm(n2))
if a2 * nd > a1:
surfaces[j] = int(round(a1 / (nd * d2)))
return surfaces
for i, s in enumerate(surface_names):
surface_numbers[s] = i
surface_count = len(surface_names)
surface_mapping = {0: 1, 1: 0, 2: 3, 3: 2, 4: 5, 5: 4,
6: 7, 7: 6, 8: 9, 9: 8, 10: 11, 11: 10,
12: 13, 13: 12, 14: 15, 15: 14, 16: 17, 17: 16,
18: 19, 19: 18, 20: 21, 21: 20, 22: 23, 23: 22,
24: 25, 25: 24}
surface_data = ([{'l': 1.0, 'd': 0.5}] * 6 +
[{'l': 1.5, 'd': 1.0 / 4.0}] * 12 +
[{'l': 1.0, 'd': 1.0 / 3.0}] * 8)
surface_symmetries = get_surface_symmetries(symmetries, surface_names,
surface_numbers)
def surface_fitting(surfaces):
for i, n1 in enumerate(np.array(surface_names)):
d1 = surface_data[i]['d']
a1 = surfaces[i] * d1
for j, n2 in enumerate(np.array(surface_names)):
d2 = surface_data[j]['d']
a2 = surfaces[j] * d2
nd = np.dot(n1, n2) / (np.linalg.norm(n1) * np.linalg.norm(n2))
if a2 * nd > a1:
surfaces[j] = int(round(a1 / (nd * d2)))
return surfaces