Method save_elements

src/IO/MshSaver.cpp:93–153

}

void MshSaver::save_elements(
        const VectorI& elements, MshSaver::ElementType type) {
    size_t nodes_per_element = 0;
    switch (type) {
        case TRI:
            nodes_per_element = 3;
            break;
        case QUAD:
            nodes_per_element = 4;
            break;
        case TET:
            nodes_per_element = 4;
            break;
        case HEX:
            nodes_per_element = 8;
            break;
        default:
            {
                std::stringstream err_msg;
                err_msg << "Unsupported element type " << type;
                throw NotImplementedError(err_msg.str());
            }
    }
    m_num_elements = elements.size() / nodes_per_element;

    // Save elements.
    fout << "$Elements" << std::endl;
    fout << m_num_elements << std::endl;

    if (m_num_elements > 0) {
        int elem_type = type;
        int num_elems = m_num_elements;
        int tags = 0;
        if (!m_binary) {
            for (size_t i=0; i<elements.size(); i+=nodes_per_element) {
                int elem_num = i/nodes_per_element + 1;
                VectorI elem = elements.segment(i, nodes_per_element) +
                    VectorI::Ones(nodes_per_element);

                fout << elem_num << " " << elem_type << " " << tags << " ";
                for (size_t j=0; j<nodes_per_element; j++) {
                    fout << elem[j] << " ";
                }
                fout << std::endl;
            }
        } else {
            fout.write((char*)&elem_type, sizeof(int));
            fout.write((char*)&num_elems, sizeof(int));
            fout.write((char*)&tags, sizeof(int));
            for (size_t i=0; i<elements.size(); i+=nodes_per_element) {
                int elem_num = i/nodes_per_element + 1;
                VectorI elem = elements.segment(i, nodes_per_element) +
                    VectorI::Ones(nodes_per_element);
                fout.write((char*)&elem_num, sizeof(int));
                fout.write((char*)elem.data(), sizeof(int)*nodes_per_element);
            }
        }
    }