def generate_particle_data_initial(output_path, iconv, sixdump, conf=dict()):
num_iconv = int(len(iconv))
num_dumps = int(len(sixdump.particles))
assert num_iconv > 0
assert num_dumps >= num_iconv
assert (num_dumps % num_iconv) == 0
num_part = num_dumps // num_iconv
assert num_part > 0
assert num_iconv > 0
NORM_ADDR = conf.get("cbuffer_norm_base_addr", 4096)
MAKE_DEMOTRACK = conf.get("make_demotrack_data", False)
MAKE_DEMOTRACK &= st.Demotrack_enabled()
init_particle_idx = 0
initial_p_buffer = create_single_particle_cbuffer(1, num_part, conf)
initial_pset_buffer = create_particle_set_cbuffer(1, num_part, conf)
initial_p_pysix = []
path_initial_pset = os.path.join(output_path, "cobj_initial_particles.bin")
for jj in range(num_part):
kk = num_part * init_particle_idx + jj
assert kk < len(sixdump.particles)
in_p = pysix.Particles(**sixdump[kk].get_minimal_beam())
in_p.state = 1
in_p.turn = 0
in_p.partid = jj
in_p.elemid = iconv[init_particle_idx]
p = st.st_SingleParticle.GET(initial_p_buffer, kk)
pysix_particle_to_single_particle(in_p, p, conf=conf)
pset = st.st_Particles.GET(initial_pset_buffer, 0)
pysix_particle_to_pset(in_p, pset, jj, conf=conf)
initial_p_pysix.append(in_p)
assert len(initial_p_pysix) == num_part
for ii, in_p in enumerate(initial_p_pysix):
assert in_p.state == 1
assert in_p.partid == ii
assert in_p.turn == 0
assert in_p.elemid == iconv[0]
path_init_pset = os.path.join(output_path, "cobj_initial_particles.bin")
if 0 == initial_pset_buffer.tofile_normalised(path_init_pset, NORM_ADDR):
print("**** -> Generated initial particle set data at:\r\n"
f"**** {path_init_pset}")
else:
raise RuntimeError("Unable to generate initial particle set data")
path_init_p = os.path.join(output_path,
"cobj_initial_single_particles.bin")
if 0 == initial_p_buffer.tofile_normalised(path_init_p, NORM_ADDR):
print("**** -> Generated initial single particle data at:\r\n" +
f"**** {path_init_p}")
else:
raise RuntimeError(
"Unable to generate initial single particle set data")
path_init_pysix = os.path.join(output_path,
"pysixtrack_initial_particles.pickle")
with open(path_init_pysix, "wb") as f_out:
pickle.dump(initial_p_pysix, f_out)
print("**** -> Generated initial pysixtrack particle data at:\r\n" +
f"**** {path_init_pysix}")
if MAKE_DEMOTRACK:
dt_particles_buffer = st.st_DemotrackParticle.CREATE_ARRAY(
num_part, True)
assert isinstance(dt_particles_buffer, np.ndarray)
assert len(dt_particles_buffer) == num_part
pset = st.st_Particles.GET(initial_pset_buffer, 0)
assert pset.num_particles == num_part
dt_p = st.st_DemotrackParticle()
for ii in range(0, num_part):
dt_p.from_cobjects(pset, ii)
dt_p.to_array(dt_particles_buffer, ii)
dt_p.clear()
path_init_dt = os.path.join(output_path,
"demotrack_initial_particles.pickle")
with open(path_init_dt, "wb") as f_out:
f_out.write(float_to_bytes(num_part))
f_out.write(dt_particles_buffer.tobytes())
print("**** -> Generated initial demotrack particle data at:\r\n" +
f"**** {path_init_dt}")
return
def generate_particle_data_until_turn(output_path,
line,
iconv,
sixdump,
until_turn,
conf=dict()):
num_iconv = int(len(iconv))
num_belem = int(len(line))
num_dumps = int(len(sixdump.particles))
assert num_iconv > 0
assert num_belem > iconv[num_iconv - 1]
assert num_dumps >= num_iconv
assert (num_dumps % num_iconv) == 0
num_particles = num_dumps // num_iconv
assert num_particles > 0
assert num_belem > 0
assert num_iconv > 0
assert until_turn > 0
start_at_element = iconv[0]
NORM_ADDR = conf.get("cbuffer_norm_base_addr", 4096)
MAKE_DEMOTRACK = conf.get("make_demotrack_data", False)
MAKE_DEMOTRACK &= st.Demotrack_enabled()
pset_buffer = create_particle_set_cbuffer(1, num_particles, conf)
assert pset_buffer.num_objects == 1
pset = st.st_Particles.GET(pset_buffer, 0)
assert pset.num_particles == num_particles
path_initial_pysix_particles = os.path.join(
output_path, "pysixtrack_initial_particles.pickle")
initial_p_pysix = None
with open(path_initial_pysix_particles, "rb") as f_in:
initial_p_pysix = pickle.load(f_in)
assert initial_p_pysix is not None
assert len(initial_p_pysix) == num_particles
for ii, in_p in enumerate(initial_p_pysix):
assert ii == in_p.partid
assert 0 == in_p.turn
assert iconv[0] == in_p.elemid
assert 1 == in_p.state
if MAKE_DEMOTRACK:
dt_p = st.st_DemotrackParticle()
dt_pset_buffer = st.st_DemotrackParticle.CREATE_ARRAY(
num_particles, True)
assert isinstance(dt_pset_buffer, np.ndarray)
assert len(dt_pset_buffer) <= num_particles
for ii, in_p in enumerate(initial_p_pysix):
assert isinstance(in_p, pysix.Particles)
print(f"**** Info :: particle {ii:6d}/{num_particles - 1:6d}")
for jj in range(in_p.turn, until_turn):
for kk, elem in enumerate(line.elements):
assert in_p.elemid == kk
assert in_p.partid == ii
assert in_p.turn == jj
assert in_p.state == 1
if in_p.state == 1:
elem.track(in_p)
if isinstance( elem, pysix.elements.Drift ) and \
in_p.state == 1 and (
in_p.x > 1.0 or in_p.x < -1.0 or
in_p.y > 1.0 or in_p.y < -1.0 ):
in_p.state = 0
if in_p.state == 1:
in_p.elemid += 1
else:
print(
f"lost particle {in_p.partid} at pos {in_p.elemid} : {in_p}"
)
print(f"lost particle {in_p.partid} at elem: {elem}")
break
if in_p.state == 1:
in_p.turn += 1
in_p.elemid = start_at_element
else:
break
if in_p.state != 1:
print(
f"lost particle {in_p.partid} at pos {in_p.elem} / turn {in_p.turn}: {in_p}"
)
pysix_particle_to_pset(in_p, pset, ii, conf=conf)
if MAKE_DEMOTRACK:
dt_p.clear()
dt_p.from_cobjects(pset, ii)
dt_p.to_array(dt_pset_buffer, ii)
path_pset_out = os.path.join(
output_path, f"cobj_particles_until_turn_{until_turn}.bin")
if 0 == pset_buffer.tofile_normalised(path_pset_out, NORM_ADDR):
print("**** -> Generated cbuffer of tracked particle data:\r\n" +
f"**** {path_pset_out}")
else:
raise ValueError(
f"Error during tracking particles until turn {until_turn}")
if MAKE_DEMOTRACK:
path_pset_out = os.path.join(
output_path, f"demotrack_particles_until_turn_{until_turn}.bin")
with open(path_pset_out, "wb") as f_out:
f_out.write(float_to_bytes(len(dt_pset_buffer)))
f_out.write(dt_pset_buffer.tobytes())
print(
"**** -> Generated demotrack data of tracked particles:\r\n" +
f"**** {path_pset_out}")
return
def generate_particle_data_elem_by_elem(output_path,
line,
iconv,
sixdump,
conf=dict()):
num_iconv = int(len(iconv))
num_belem = int(len(line))
num_dumps = int(len(sixdump.particles))
assert num_iconv > 0
assert num_belem > iconv[num_iconv - 1]
assert num_dumps >= num_iconv
assert (num_dumps % num_iconv) == 0
num_particles = num_dumps // num_iconv
assert num_particles > 0
assert num_belem > 0
assert num_iconv > 0
NORM_ADDR = conf.get("cbuffer_norm_base_addr", 4096)
MAKE_DEMOTRACK = conf.get("make_demotrack_data", False)
MAKE_DEMOTRACK &= st.Demotrack_enabled()
pset_buffer = create_particle_set_cbuffer(num_belem + 1, num_particles,
conf)
assert pset_buffer.num_objects == num_belem + 1
path_initial_pysix_particles = os.path.join(
output_path, "pysixtrack_initial_particles.pickle")
initial_p_pysix = None
with open(path_initial_pysix_particles, "rb") as f_in:
initial_p_pysix = pickle.load(f_in)
assert initial_p_pysix is not None
assert len(initial_p_pysix) == num_particles
if MAKE_DEMOTRACK:
dt_p = st.st_DemotrackParticle()
dt_pset_buffer = st.st_DemotrackParticle.CREATE_ARRAY(
(num_belem + 1) * num_particles, True)
assert isinstance(dt_pset_buffer, np.ndarray)
assert len(dt_pset_buffer) <= (num_belem + 1) * num_particles
for ii, in_p in enumerate(initial_p_pysix):
assert isinstance(in_p, pysix.Particles)
assert isinstance(in_p, pysix.Particles)
assert in_p.elemid == iconv[0]
assert in_p.partid == ii
assert in_p.turn == 0
assert in_p.state == 1
print(f"**** Info :: particle {ii:6d}/{num_particles - 1:6d}")
for jj, elem in enumerate(line.elements):
assert in_p.elemid == jj
assert in_p.partid == ii
assert in_p.turn == 0
assert in_p.state == 1
pset = st.st_Particles.GET(pset_buffer, jj)
assert pset.num_particles == num_particles
pysix_particle_to_pset(in_p, pset, ii, conf=conf)
if MAKE_DEMOTRACK:
dt_p.clear()
dt_p.from_cobjects(pset, ii)
kk = jj * num_particles + ii
dt_p.to_array(dt_pset_buffer, kk)
if in_p.state == 1:
elem.track(in_p)
if isinstance( elem, pysix.elements.Drift ) and in_p.state == 1 and \
( in_p.x > 1.0 or in_p.x < -1.0 or
in_p.y > 1.0 or in_p.y < -1.0 ):
in_p.state = 0
if in_p.state == 1:
in_p.elemid += 1
else:
print(
f"lost particle {in_p.partid} at pos {in_p.elemid} : {in_p}"
)
print(f"lost particle {in_p.partid} at elem: {elem}")
break
if in_p.state == 1:
in_p.turn += 1
in_p.elemid = iconv[0]
else:
print(f"lost particle {in_p.partid}")
pset = st.st_Particles.GET(pset_buffer, num_belem)
assert pset.num_particles == num_particles
pysix_particle_to_pset(in_p, pset, ii, conf=conf)
if MAKE_DEMOTRACK:
dt_p.clear()
dt_p.from_cobjects(pset, ii)
dt_p.to_array(dt_pset_buffer, num_particles * num_belem + ii)
path_elem_by_elem = os.path.join(
output_path, "cobj_particles_elem_by_elem_pysixtrack.bin")
if 0 == pset_buffer.tofile_normalised(path_elem_by_elem, NORM_ADDR):
print("**** -> Generated cbuffer of particle elem-by-elem data:\r\n" +
f"**** {path_elem_by_elem}")
else:
raise RuntimeError("Unable to generate cobjects elem-by-elem data")
if MAKE_DEMOTRACK:
path_elem_by_elem = os.path.join(
output_path, "demotrack_elem_by_elem_pysixtrack.pickle")
with open(path_elem_by_elem, "wb") as f_out:
f_out.write(float_to_bytes(len(dt_pset_buffer)))
f_out.write(dt_pset_buffer.tobytes())
print(
"**** -> Generated demotrack particle elem-by-elem data:\r\n" +
f"**** {path_elem_by_elem}")
return
def generate_particle_data_initial(input_path, output_path, conf=dict()):
path_in_particles = os.path.join(input_path,
"pysixtrack_initial_particles.pickle")
with open(path_in_particles, "rb") as f_in:
initial_p_pysix = pickle.load(f_in)
print("**** -> Read input data from:\r\n" +
f"**** {path_in_particles}")
num_part = len(initial_p_pysix)
NORM_ADDR = conf.get("cbuffer_norm_base_addr", 4096)
MAKE_DEMOTRACK = conf.get("make_demotrack_data", False)
MAKE_DEMOTRACK &= st.Demotrack_enabled()
init_particle_idx = 0
start_at_element = 0
initial_p_buffer = create_single_particle_cbuffer(1, num_part, conf)
initial_pset_buffer = create_particle_set_cbuffer(1, num_part, conf)
path_initial_pset = os.path.join(output_path, "cobj_initial_particles.bin")
for jj, in_p in enumerate(initial_p_pysix):
assert jj < num_part
assert isinstance(in_p, pysix.Particles)
in_p.elemid = start_at_element
in_p.turn = 0
in_p.state = 1
in_p.partid = jj
p = st.st_SingleParticle.GET(initial_p_buffer, jj)
pysix_particle_to_single_particle(in_p, p)
pset = st.st_Particles.GET(initial_pset_buffer, 0)
pysix_particle_to_pset(in_p, pset, jj)
path_init_pset = os.path.join(output_path, "cobj_initial_particles.bin")
if 0 == initial_pset_buffer.tofile_normalised(path_init_pset, NORM_ADDR):
print("**** -> Generated initial particle set data at:\r\n"
f"**** {path_init_pset}")
else:
raise RuntimeError("Unable to generate initial particle set data")
path_init_p = os.path.join(output_path,
"cobj_initial_single_particles.bin")
if 0 == initial_p_buffer.tofile_normalised(path_init_p, NORM_ADDR):
print("**** -> Generated initial single particle data at:\r\n" +
f"**** {path_init_p}")
else:
raise RuntimeError(
"Unable to generate initial single particle set data")
path_init_pysix = os.path.join(output_path,
"pysixtrack_initial_particles.pickle")
with open(path_init_pysix, "wb") as f_out:
pickle.dump(initial_p_pysix, f_out)
print("**** -> Generated initial pysixtrack particle data at:\r\n" +
f"**** {path_init_pysix}")
if MAKE_DEMOTRACK:
dt_particles_buffer = st.st_DemotrackParticle.CREATE_ARRAY(
num_part, True)
assert isinstance(dt_particles_buffer, np.ndarray)
assert len(dt_particles_buffer) == num_part
pset = st.st_Particles.GET(initial_pset_buffer, 0)
assert pset.num_particles == num_part
dt_p = st.st_DemotrackParticle()
for ii in range(0, num_part):
dt_p.from_cobjects(pset, ii)
dt_p.to_array(dt_particles_buffer, ii)
dt_p.clear()
path_init_dt = os.path.join(output_path,
"demotrack_initial_particles.pickle")
with open(path_init_dt, "wb") as f_out:
f_out.write(float_to_bytes(num_part))
f_out.write(dt_particles_buffer.tobytes())
print("**** -> Generated initial demotrack particle data at:\r\n" +
f"**** {path_init_dt}")
return
def generate_particle_data_until_turn(input_path,
output_path,
until_turn,
conf=dict()):
path_in_particles = os.path.join(output_path,
"pysixtrack_initial_particles.pickle")
with open(path_in_particles, "rb") as f_in:
initial_p_pysix = pickle.load(f_in)
print("**** -> Read input data from:\r\n" +
f"**** {path_in_particles}")
path_in_line = os.path.join(output_path, "pysixtrack_lattice.pickle")
with open(path_in_line, "rb") as f_in:
line = pickle.load(f_in)
num_belem = len(line)
num_part = len(initial_p_pysix)
start_at_element = 0
NORM_ADDR = conf.get("cbuffer_norm_base_addr", 4096)
MAKE_DEMOTRACK = conf.get("make_demotrack_data", False)
MAKE_DEMOTRACK &= st.Demotrack_enabled()
pset_buffer = create_particle_set_cbuffer(1, num_part, conf)
assert pset_buffer.num_objects == 1
pset = st.st_Particles.GET(pset_buffer, 0)
assert pset.num_particles == num_part
if MAKE_DEMOTRACK:
dt_p = st.st_DemotrackParticle()
dt_pset_buffer = st.st_DemotrackParticle.CREATE_ARRAY(num_part, True)
assert isinstance(dt_pset_buffer, np.ndarray)
assert len(dt_pset_buffer) <= num_part
for ii, in_p in enumerate(initial_p_pysix):
assert isinstance(in_p, pysix.Particles)
assert in_p.elemid == start_at_element
assert in_p.partid == ii
assert in_p.state == 1
assert in_p.turn == 0
print(f"**** Info :: particle {ii:6d}/{num_part - 1:6d}")
start_at_turn = in_p.turn
if start_at_turn < until_turn:
for jj in range(start_at_turn, until_turn):
for kk, elem in enumerate(line):
if conf.get('always_use_drift_exact', False):
assert not isinstance( elem, pysix.elements.Drift ) or \
isinstance( elem, pysix.elements.DriftExact )
if in_p.state == 1:
elem.track(in_p)
if isinstance( elem, pysix.elements.Drift ) and \
in_p.state == 1 and \
( in_p.x > 1.0 or in_p.x < -1.0 or
in_p.y > 1.0 or in_p.y < -1.0 ):
in_p.state = 0
if in_p.state == 1:
in_p.elemid += 1
else:
print(f"lost particle {ii} at pos {jj} : {in_p}")
print(f"lost particle at elem: {elem}")
break
if in_p.state == 1:
in_p.turn += 1
in_p.elemid = start_at_element
else:
break
if in_p.state != 1:
print(
f"lost particle {in_p.partid} at pos {in_p.elem} / turn {in_p.turn}: {in_p}"
)
pysix_particle_to_pset(
in_p,
pset,
ii,
)
if MAKE_DEMOTRACK:
dt_p.clear()
dt_p.from_cobjects(pset, ii)
dt_p.to_array(dt_pset_buffer, ii)
path_pset_out = os.path.join(
output_path, f"cobj_particles_until_turn_{until_turn}.bin")
if 0 == pset_buffer.tofile_normalised(path_pset_out, NORM_ADDR):
print("**** -> Generated cbuffer of tracked particle data:\r\n" +
f"**** {path_pset_out}")
else:
raise ValueError(
f"Error during tracking particles until turn {until_turn}")
if MAKE_DEMOTRACK:
path_pset_out = os.path.join(
output_path, f"demotrack_particles_until_turn_{until_turn}.bin")
with open(path_pset_out, "wb") as f_out:
f_out.write(float_to_bytes(len(dt_pset_buffer)))
f_out.write(dt_pset_buffer.tobytes())
print(
"**** -> Generated demotrack data of tracked particles:\r\n" +
f"**** {path_pset_out}")
return
def generate_particle_data_elem_by_elem(input_path, output_path, conf=dict()):
path_in_particles = os.path.join(output_path,
"pysixtrack_initial_particles.pickle")
with open(path_in_particles, "rb") as f_in:
initial_p_pysix = pickle.load(f_in)
print("**** -> Read input data from:\r\n" +
f"**** {path_in_particles}")
path_in_line = os.path.join(output_path, "pysixtrack_lattice.pickle")
with open(path_in_line, "rb") as f_in:
line = pickle.load(f_in)
num_belem = len(line)
num_part = len(initial_p_pysix)
start_at_element = 0
NORM_ADDR = conf.get("cbuffer_norm_base_addr", 4096)
MAKE_DEMOTRACK = conf.get("make_demotrack_data", False)
MAKE_DEMOTRACK &= st.Demotrack_enabled()
pset_buffer = create_particle_set_cbuffer(num_belem + 1, num_part, conf)
assert pset_buffer.num_objects == num_belem + 1
if MAKE_DEMOTRACK:
dt_p = st.st_DemotrackParticle()
dt_pset_buffer = st.st_DemotrackParticle.CREATE_ARRAY(
(num_belem + 1) * num_part, True)
assert isinstance(dt_pset_buffer, np.ndarray)
assert len(dt_pset_buffer) <= (num_belem + 1) * num_part
for ii, in_p in enumerate(initial_p_pysix):
assert isinstance(in_p, pysix.Particles)
assert in_p.elemid == start_at_element
assert in_p.partid == ii
assert in_p.turn == 0
assert in_p.state == 1
print(f"**** Info :: particle {ii:6d}/{num_part - 1:6d}")
for jj, elem in enumerate(line):
pset = st.st_Particles.GET(pset_buffer, jj)
assert pset.num_particles == num_part
assert in_p.elemid == jj
assert in_p.partid == ii
assert in_p.turn == 0
assert in_p.state == 1
if conf.get('always_use_drift_exact', False):
assert not isinstance( elem, pysix.elements.Drift ) or \
isinstance( elem, pysix.elements.DriftExact )
pysix_particle_to_pset(in_p, pset, ii, conf=conf)
if MAKE_DEMOTRACK:
dt_p.clear()
dt_p.from_cobjects(pset, ii)
kk = jj * num_part + ii
assert kk < len(dt_pset_buffer)
dt_p.to_array(dt_pset_buffer, kk)
if in_p.state == 1:
elem.track(in_p)
if isinstance( elem, pysix.elements.Drift ) and \
in_p.state == 1 and \
( in_p.x > 1.0 or in_p.x < -1.0 or
in_p.y > 1.0 or in_p.y < -1.0 ):
in_p.state = 0
if in_p.state == 1:
in_p.elemid += 1
else:
print(
f"lost particle {in_p.partid} at pos {in_p.elemid} : {in_p}"
)
print(f"lost particle {in_p.partid} at elem: {elem}")
break
if in_p.state == 1:
in_p.turn += 1
in_p.elemid = start_at_element
else:
print(f"lost particle {ii} at pos {in_p.elemid}: {in_p}")
pset = st.st_Particles.GET(pset_buffer, num_belem)
assert pset.num_particles == num_part
pysix_particle_to_pset(in_p, pset, ii, conf=conf)
if MAKE_DEMOTRACK:
dt_p.clear()
dt_p.from_cobjects(pset, ii)
kk = num_belem * num_part + ii
assert len(dt_pset_buffer) > kk
dt_p.to_array(dt_pset_buffer, kk)
path_elem_by_elem = os.path.join(
output_path, "cobj_particles_elem_by_elem_pysixtrack.bin")
if 0 == pset_buffer.tofile_normalised(path_elem_by_elem, NORM_ADDR):
print("**** -> Generated cbuffer of particle elem-by-elem data:\r\n" +
f"**** {path_elem_by_elem}")
else:
raise RuntimeError("Unable to generate cobjects elem-by-elem data")
if MAKE_DEMOTRACK:
path_elem_by_elem = os.path.join(
output_path, "demotrack_elem_by_elem_pysixtrack.pickle")
with open(path_elem_by_elem, "wb") as f_out:
f_out.write(float_to_bytes(len(dt_pset_buffer)))
f_out.write(dt_pset_buffer.tobytes())
print(
"**** -> Generated demotrack particle elem-by-elem data:\r\n" +
f"**** {path_elem_by_elem}")
return