void write_cmm(const std::string &cmm_filename,
const std::string &marker_set_name, const AnchorsData &ad) {
Floats radii;
// algebra::get_enclosing_sphere(dpa.get_cluster_vectors(i));
radii.insert(radii.begin(), ad.get_number_of_points(), 5.);
std::ofstream out;
out.open(cmm_filename.c_str(), std::ios::out);
write_cmm_helper(out, marker_set_name, ad.points_, ad.edges_, radii);
out.close();
}
algebra::VectorKD get_embedding(const Subset &s,
const Assignment &a,
ParticleStatesTable *pst){
Floats embed;
for (unsigned int i=0; i< s.size(); ++i) {
algebra::VectorKD cur
= pst->get_particle_states(s[i])->get_embedding(a[i]);
embed.insert(embed.end(), cur.coordinates_begin(),
cur.coordinates_end());
}
return embed;
}