src/netgenreader.cpp

Go to the documentation of this file.

/***************************************************************************
                          netgenreader.cpp  -  description
                             -------------------
    begin                : Fri Dec 12 2003
    copyright            : (C) 2003 by Roman Geus
    email                : roman.geus@psi.ch
***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#include <cassert>

#include "ioerror.h"
#include "netgenreader.h"

using namespace std;

NetgenReader::~NetgenReader() {
}

NetgenReader::NetgenReader(istream* istr, MeshBuilder* builder)
    : _istr(istr), _builder(builder) {}

void NetgenReader::read() {
    const string file_name = "<input stream>";
    istream& istr = *_istr;

    int line = 0;

    // read number of nodes
    int nof_node;
    istr >> nof_node;
    line ++;
    if (istr.fail())
        throw FileIOError(file_name, line, "error reading number of nodes");
    _builder->init_coord(nof_node);

    // read node coordinates
    double x, y, z;
    for (int i = 0; i < nof_node; i ++) {
        istr >> x >> y >> z;
        line ++;
        if (istr.fail())
            throw FileIOError(file_name, line, "error reading node coordinates");
        _builder->set_coord(i, x, y, z);
    }
    _builder->finalize_coord();

    // read number of tets
    int nof_tet;
    istr >> nof_tet;
    line ++;
    if (istr.fail())
        throw FileIOError(file_name, line, "error reading number of elements");
    _builder->init_tet(nof_tet);

    // read tetnode
    int node_ids[4], node_id, mat;
    for (int t = 0; t < nof_tet; t ++) {
        istr >> mat;
        for (int i = 0; i < 4; i ++) {
            istr >> node_id;
            node_id --;
            if (node_id < 0 || node_id >= nof_node)
                throw FileIOError(file_name, line, "node number out of range");
            node_ids[i] = node_id;
        }
        line ++;
        if (istr.fail())
            throw FileIOError(file_name, line, "error reading node numbers");
        _builder->set_tet(t, node_ids[0], node_ids[1], node_ids[2], node_ids[3], mat);
    }
    _builder->finalize_tet();

    // read number of boundary faces
    int nof_bc_face;
    istr >> nof_bc_face;
    line ++;
    if (istr.fail())
        throw FileIOError(file_name, line, "error reading number of boundary faces");
    _builder->init_bc(nof_bc_face);

    // read boundary face data
    int nof_bc[4] = {0, 0, 0, 0};
    int face_node[3];
    for (int i = 0; i < nof_bc_face; i ++) {
        int bc_id;
        istr >> bc_id;
        nof_bc[bc_id] += 1;
        for (int j = 0; j < 3; j ++) {
            istr >> node_id;
            node_id --;
            if (node_id < 0 || node_id >= nof_node)
                throw FileIOError(file_name, line, "node number out of range");
            face_node[j] = node_id;
        }
        line ++;
        if (istr.fail())
            throw FileIOError(file_name, line, "error reading boundary data");
        _builder->set_bc(face_node[0], face_node[1], face_node[2], bc_id);
    }

    // symmetry planes
    int nof_sym;
    for (nof_sym = 3; nof_sym >= 0 && nof_bc[nof_sym] == 0; nof_sym --);
  
    assert(nof_sym >= 0);
    for (int i = 1; i < nof_sym; i ++)
        if (nof_bc[i] == 0)
            throw FileIOError(file_name, "empty symmetry plane detected");
    _builder->finalize_bc(nof_sym);
}

---