Method write_elem_file

src/IO/NodeWriter.cpp:118–153

}

void NodeWriter::write_elem_file(const std::string& filename, Mesh& mesh) {
    const size_t num_voxels = mesh.get_num_voxels();
    const size_t vertex_per_voxel = mesh.get_vertex_per_voxel();
    if (num_voxels == 0) return;
    if (vertex_per_voxel != 4) {
        throw IOError("Only tet element is supported in .ele file.");
    }

    std::ofstream fout(filename.c_str());
    fout.precision(16);
    if (!is_anonymous()) {
        fout << "# Generated with PyMesh" << std::endl;
    }
    VectorF region;
    if (m_with_region_attribute) {
        if (!mesh.has_attribute("region")) {
            throw IOError("Attribute \"region\" does not exist.");
        } else {
            region = mesh.get_attribute("region");
            assert(num_voxels == region.size());
        }
    }

    const VectorI& voxels = mesh.get_voxels();
    fout << num_voxels << " 4 " << m_with_region_attribute << std::endl;
    for (size_t i=0; i<num_voxels; i++) {
        fout << i;
        for (size_t j=0; j<4; j++) {
            fout << " " << voxels[i*4+ j];
        }
        if (m_with_region_attribute) {
            fout << " " << region[i];
        }
        fout << std::endl;
    }
}