#ifdef DAGMC

if (settings::dagmc) {

return;

}

#endif

auto geom_group = create_group(file_id, "geometry");

write_attribute(geom_group, "n_cells", cells.size());

write_attribute(geom_group, "n_surfaces", surfaces.size());

write_attribute(geom_group, "n_universes", universes.size());

write_attribute(geom_group, "n_lattices", lattices.size());

auto cells_group = create_group(geom_group, "cells");

for (Cell* c : cells) c->to_hdf5(cells_group);

close_group(cells_group);

auto surfaces_group = create_group(geom_group, "surfaces");

for (Surface* surf : surfaces) surf->to_hdf5(surfaces_group);

close_group(surfaces_group);

auto universes_group = create_group(geom_group, "universes");

for (Universe* u : universes) u->to_hdf5(universes_group);

close_group(universes_group);

auto lattices_group = create_group(geom_group, "lattices");

for (Lattice* lat : lattices) lat->to_hdf5(lattices_group);

close_group(lattices_group);

close_group(geom_group);