pub fn to_polyhedron_mesh(
tet2vtx: &[u32],
pyrmd2vtx: &[u32],
prism2vtx: &[u32],
elem2idx_offset: &mut [u32],
idx2vtx: &mut [u32],
) {
let num_tet = tet2vtx.len() / 4;
let num_pyrmd = pyrmd2vtx.len() / 5;
let num_prism = prism2vtx.len() / 6;
let num_elem = num_tet + num_pyrmd + num_prism;
assert_eq!(elem2idx_offset.len(), num_elem + 1);
{
let mut idx = 0;
elem2idx_offset[idx] = 0;
for _i_tet in 0..num_tet {
elem2idx_offset[idx + 1] = elem2idx_offset[idx] + 4;
idx += 1;
}
for _i_pyrmd in 0..num_pyrmd {
elem2idx_offset[idx + 1] = elem2idx_offset[idx] + 5;
idx += 1;
}
for _i_prism in 0..num_prism {
elem2idx_offset[idx + 1] = elem2idx_offset[idx] + 6;
idx += 1;
}
assert_eq!(idx, num_elem);
elem2idx_offset
};
assert_eq!(
idx2vtx.len(),
tet2vtx.len() + pyrmd2vtx.len() + prism2vtx.len()
);
idx2vtx[0..num_tet * 4].copy_from_slice(tet2vtx);
idx2vtx[num_tet * 4..(num_tet * 4 + num_pyrmd * 5)].copy_from_slice(pyrmd2vtx);
idx2vtx[(num_tet * 4 + num_pyrmd * 5)..(num_tet * 4 + num_pyrmd * 5 + num_prism * 6)]
.copy_from_slice(prism2vtx);
}